HDU 5876 Sparse Graph

题目：Sparse Graph

链接：http://acm.split.hdu.edu.cn/showproblem.php?pid=5876

题意：给出一个图（V<=20万，E<=2万），要求先化为补图（每两个点，原本有边删去边，原本没边添加边），然后问指定点S到其他每个点的最短距离。

思路：

　　普通的广搜应该解决不了。。。O（n*m）太大，不会很难，比赛时没做出来有点可惜，当知道过了也进不了时，就安慰了许多。

　　变化一下思路，从起点出发，因为原本边<=2万，那么大部分的点都已经可以到达了，也就是len=1（len最短距离），现在用集合st 来存放还没到达的点，用ad[i]=true来表示i 已经到达，用co表示已经到达的点的数量。每一次循环遍历st，对于某一个还没到达的点x ，假定和x 相邻的点有m 个，如果m<co，那说明x 肯定能到达。因为原本y 和x 没有边，补图里肯定就有，也就是说最坏情况是m个点对应的都是已经遍历过的，那么co-m这些点肯定能到达x。如果m>=co，并不能说明x 一定到不了，现在就要遍历m个点，如果某个点是没有遍历过的，m可以-1，如果最终m<co，说明x 还是可以遍历到的。每一次循环，就是找出若干个x ，x能到达，如果没有这样的x 就可以退出了，如果有，集合删去这些x，记录len，记录ad，记录co ，进入下次循环。

　　时间复杂度计算：集合最初最多2万个点，每次删掉就不恢复了，这里的时间为E，m>=co要找点时也就是遍历所有的边，每个边最多遍历1次，时间E。也就是说循环里看似套了很多层，其实很快，外面大概就是V，时间复杂度V+E。

AC代码：

  1 #include<stdio.h>
  2 #include<string.h>
  3 #include<stdlib.h>
  4 #include<math.h>
  5 #include<set>
  6 #include<map>
  7 #include<list>
  8 #include<stack>
  9 #include<queue>
 10 #include<vector>
 11 #include<string>
 12 #include<iostream>
 13 #include<algorithm>
 14 using namespace std;
 15 #define lson rt<<1
 16 #define rson rt<<1|1
 17 #define N 200010
 18 #define M 100010
 19 #define Mod 1000000007
 20 #define LL long long
 21 #define INF 0x7fffffff
 22 #define FOR(i,f_start,f_end) for(int i=f_start;i<=f_end;i++)
 23 #define For(i,f_start,f_end) for(int i=f_start;i<f_end;i++)
 24 #define REP(i,f_end,f_start) for(int i=f_end;i>=f_start;i--)
 25 #define Rep(i,f_end,f_start) for(int i=f_end;i>f_start;i--)
 26 #define MT(x,i) memset(x,i,sizeof(x))
 27 #define gcd(x,y) __gcd(x,y)
 28 const double PI = acos(-1);
 29
 30 bool ad[N];
 31 int len[N];
 32 vector<int> v[N];
 33 vector<int> zan;
 34 set<int> st;
 35 int main()
 36 {
 37   int t,n,m,x,y;
 38   scanf("%d",&t);
 39   while(t--)
 40   {
 41     scanf("%d%d",&n,&m);
 42     while(m--)
 43     {
 44       scanf("%d%d",&x,&y);
 45       v[x].push_back(y);
 46       v[y].push_back(x);
 47     }
 48     int s,co=1;
 49     scanf("%d",&s);
 50     st.clear();
 51     MT(ad,0);
 52     ad[s]=true;
 53     For(i,0,v[s].size()){
 54       ad[v[s][i]]=true;
 55     }
 56
 57     MT(len,-1);
 58     FOR(i,1,n){
 59       if(i==s) continue;
 60       if(ad[i]==false){
 61         len[i]=1;
 62         co++;
 63         ad[i]=true;
 64       }
 65       else{
 66         st.insert(i);
 67         ad[i]=false;
 68       }
 69     }
 70     int ll=2;
 71     while(1)
 72     {
 73       zan.clear();
 74       bool ff=false;
 75       for(set<int>::iterator it=st.begin();it!=st.end();it++)
 76       {
 77         int x=*it;
 78         bool flag=false;
 79         if(v[x].size()<co)
 80         {
 81           flag=true;
 82         }
 83         else
 84         {
 85           int k=v[x].size();
 86           For(i,0,v[x].size())
 87           {
 88             int son=v[x][i];
 89             if(ad[son]==false) k--;
 90             if(k<co) break;
 91           }
 92           if(k<co)
 93           {
 94             flag=true;
 95           }
 96         }
 97         if(flag==true)
 98         {
 99           ff=true;
100           zan.push_back(x);
101         }
102       }
103       if(ff==false) break;
104       For(i,0,zan.size())
105       {
106         ad[zan[i]]=true;
107         len[zan[i]]=ll;
108         co++;
109       }
110       ll++;
111       for(set<int>::iterator it=st.begin();it!=st.end();)
112       {
113         int x=*it;
114         if(ad[x]==true) st.erase(it++);
115         else it++;
116       }
117     }
118     bool ff=false;
119     FOR(i,1,n){
120       if(i!=s){
121         if(ff!=false) printf(" ");
122         else ff=true;
123         printf("%d",len[i]);
124       }
125     }
126     printf("\n");
127     FOR(i,1,n){
128       v[i].clear();
129     }
130   }
131   return 0;
132 }