uestc 方老师和缘分

关于怎么建图，自己还真是想不出来。

我觉得就是找到每个方老师所在的scc里的所有缘分吧。。

转自http://www.cnblogs.com/whatbeg/p/3765621.html

做法：建图，每个方老师和它想要的缘分之间连一条有向边，然后，在给出的初始匹配中反向建边，即如果第i个方老师现在找到的是缘分u，则建边u->i。这样求出所有的强连通分量，每个强连通分量中方老师和缘分的数目一定是相等的，所以每个方老师一定可以找到与他在同一个强连通分量里的缘分，因为强连通分量中每个点都是可达的，某个方老师找到了其强连通分量中的非原配点，则该原配缘分一定可以在强连通分量中找到"新欢"。可以画个图看看。

由于要构造非二分图，缘分的编号从n+1开始，到2n。

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstdlib>
  4 #include<cstring>
  5 #include<string>
  6 #include<queue>
  7 #include<algorithm>
  8 #include<map>
  9 #include<iomanip>
 10 #include<climits>
 11 #include<string.h>
 12 #include<cmath>
 13 #include<stdlib.h>
 14 #include<vector>
 15 #include<stack>
 16 #include<set>
 17 using namespace std;
 18 #define INF 1000000007
 19 #define MAXN 4010
 20 #define Mod 1000007
 21 #define N 100007
 22 #define NN 30
 23 #define sigma_size 3
 24 const int maxn = 6e5 + 10;
 25 using namespace std;
 26 typedef long long LL;
 27 const double pi = acos(-1);
 28 
 29 
 30 vector<int> G[MAXN];
 31 int dfn[MAXN], low[MAXN], instk[MAXN], sccno[MAXN], ans[MAXN];
 32 stack<int> S;
 33 int Time, cnt ,n, res;
 34 void Tarjan(int u)
 35 {
 36     dfn[u] = low[u] = ++Time;
 37     S.push(u);
 38     instk[u] = 1;
 39     for (int i = 0; i < G[u].size(); ++i) {
 40         int v = G[u][i];
 41         if (!dfn[v]) {
 42             Tarjan(v);
 43             low[u] = min(low[u], low[v]);
 44         }
 45         else if (instk[v]) {
 46             low[u] = min(low[u], dfn[v]);
 47         }
 48     }
 49     if (low[u] == dfn[u]) {
 50         cnt++;
 51         int v;
 52         do{
 53             v = S.top();
 54             S.pop();
 55             instk[v] = 0;
 56             sccno[v] = cnt;
 57         } while (v != u);
 58     }
 59 }
 60 
 61 void find_scc()
 62 {
 63     for (int i = 1; i <= n; ++i)
 64         if (!dfn[i]) Tarjan(i);
 65 }
 66 
 67 void init()
 68 {
 69     for (int i = 0; i <= n; ++i)
 70         G[i].clear();
 71     Time = cnt = 0;
 72     memset(dfn, 0, sizeof(dfn));
 73     memset(low, 0, sizeof(low));
 74     memset(sccno,0,sizeof(sccno));
 75     memset(instk, 0, sizeof(instk));
 76     while (!S.empty()) S.pop();
 77 }
 78 
 79 int main()
 80 {
 81     int k, u, v;
 82     while (~scanf("%d",&n)) {
 83         init();
 84         for (int i = 1; i <= n; ++i) {
 85             scanf("%d", &k);
 86             for (int j = 0; j < k; ++j) {
 87                 cin >> v;
 88                 G[i].push_back(v + n);
 89             }
 90         }
 91         for (int i = 1; i <= n; ++i) {
 92             scanf("%d", &v);
 93             G[v + n].push_back(i);
 94         }
 95         find_scc();
 96         for (int u = 1; u <= n; ++u) {
 97             k = 0;
 98             for (int i = 0; i < G[u].size(); ++i) {
 99                 v = G[u][i];
100                 if (sccno[u] == sccno[v])
101                     ans[k++] = v - n;
102             }
103             sort(ans, ans + k);
104             printf("%d", k);
105             for (int i = 0; i < k; ++i)
106                 printf(" %d",ans[i]);
107             puts("");
108         }
109     }
110     return 0;
111 }