[HDOJ5876]Sparse Graph（补图最短路）

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

题意：求一个图的补图中的单源最短路。

题解说：补图上的 BFS 是非常经典的问题。一般的做法是用链表（或者偷懒用 std::set）维护还没 BFS 过的点。当要扩展点 u 的时候，遍历一次还没访问过的点 v，如果 uv 没边，那么将 v 入队。否则将 v 留在未扩展点中。很明显，后者只会发生 m 次，前者只会发生 n 次，所以复杂度是

总得说，求补图的最短路就是要存原图，并且需要有一个额外的set来维护未扩展的点集合。

在bfs的时候从队列取出一个点u，访问该点u邻接的点v，假如uv之间有边则补图无边，否则v会入队。

最后的时候要把入队的点从未扩展点集合中删掉。

  1 #include <algorithm>
  2 #include <iostream>
  3 #include <iomanip>
  4 #include <cstring>
  5 #include <climits>
  6 #include <complex>
  7 #include <fstream>
  8 #include <cassert>
  9 #include <cstdio>
 10 #include <bitset>
 11 #include <vector>
 12 #include <deque>
 13 #include <queue>
 14 #include <stack>
 15 #include <ctime>
 16 #include <set>
 17 #include <map>
 18 #include <cmath>
 19 using namespace std;
 20 #define fr first
 21 #define sc second
 22 #define cl clear
 23 #define BUG puts("here!!!")
 24 #define W(a) while(a--)
 25 #define pb(a) push_back(a)
 26 #define Rint(a) scanf("%d", &a)
 27 #define Rll(a) scanf("%I64d", &a)
 28 #define Rs(a) scanf("%s", a)
 29 #define Cin(a) cin >> a
 30 #define FRead() freopen("in", "r", stdin)
 31 #define FWrite() freopen("out", "w", stdout)
 32 #define Rep(i, len) for(int i = 0; i < (len); i++)
 33 #define For(i, a, len) for(int i = (a); i < (len); i++)
 34 #define Cls(a) memset((a), 0, sizeof(a))
 35 #define Clr(a, x) memset((a), (x), sizeof(a))
 36 #define Full(a) memset((a), 0x7f7f7f, sizeof(a))
 37 #define lrt rt << 1
 38 #define rrt rt << 1 | 1
 39 #define pi 3.14159265359
 40 #define RT return
 41 #define lowbit(x) x & (-x)
 42 #define onecnt(x) __builtin_popcount(x)
 43 typedef long long LL;
 44 typedef long double LD;
 45 typedef unsigned long long ULL;
 46 typedef pair<int, int> pii;
 47 typedef pair<string, int> psi;
 48 typedef pair<LL, LL> pll;
 49 typedef map<string, int> msi;
 50 typedef vector<int> vi;
 51 typedef vector<LL> vl;
 52 typedef vector<vl> vvl;
 53 typedef vector<bool> vb;
 54 
 55 const int maxn = 200200;
 56 const int inf = 10000010;
 57 set<int> G[maxn];
 58 set<int> x;
 59 int n, m, s;
 60 int d[maxn];
 61 
 62 void bfs(int s) {
 63     Clr(d, -1);
 64     d[s] = 0;
 65     queue<int> q;
 66     q.push(s); x.erase(s);
 67     set<int>::iterator v;
 68     while(!q.empty()) {
 69         int u = q.front(); q.pop();
 70         set<int> y;
 71         for(v = x.begin(); v != x.end(); v++) {
 72             if(G[u].find(*v) == G[u].end()) {
 73                 if(d[*v] == -1) d[*v] = d[u] + 1;
 74                 else d[*v] = min(d[*v], d[u] + 1);
 75                 y.insert(*v);
 76                 q.push(*v);
 77             }
 78         }
 79         for(v = y.begin(); v != y.end(); v++)    {
 80             x.erase(*v);
 81         }
 82     }
 83 }
 84 
 85 int main() {
 86     // FRead();
 87     int T, u, v;
 88     Rint(T);
 89     W(T) {
 90         Rint(n); Rint(m);
 91         x.clear();
 92         For(i, 1, n+1) {
 93             x.insert(i);
 94             G[i].clear();
 95         }
 96         Rep(i, m) {
 97             Rint(u); Rint(v);
 98             G[u].insert(v); G[v].insert(u);
 99         }
100         Rint(s);
101         bfs(s);
102         For(i, 1, n+1) {
103             if(i == s) continue;
104             printf("%d%c", d[i], i == n ? '
' : ' ');
105         }
106     }
107     RT 0;
108 }