uestc 方老师和缘分

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

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

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

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

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

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<queue>
#include<algorithm>
#include<map>
#include<iomanip>
#include<climits>
#include<string.h>
#include<cmath>
#include<stdlib.h>
#include<vector>
#include<stack>
#include<set>
using namespace std;
#define INF 1000000007
#define MAXN 4010
#define Mod 1000007
#define N 100007
#define NN 30
#define sigma_size 3
const int maxn = 6e5 + 10;
using namespace std;
typedef long long LL;
const double pi = acos(-1);


vector<int> G[MAXN];
int dfn[MAXN], low[MAXN], instk[MAXN], sccno[MAXN], ans[MAXN];
stack<int> S;
int Time, cnt ,n, res;
void Tarjan(int u)
{
    dfn[u] = low[u] = ++Time;
    S.push(u);
    instk[u] = 1;
    for (int i = 0; i < G[u].size(); ++i) {
        int v = G[u][i];
        if (!dfn[v]) {
            Tarjan(v);
            low[u] = min(low[u], low[v]);
        }
        else if (instk[v]) {
            low[u] = min(low[u], dfn[v]);
        }
    }
    if (low[u] == dfn[u]) {
        cnt++;
        int v;
        do{
            v = S.top();
            S.pop();
            instk[v] = 0;
            sccno[v] = cnt;
        } while (v != u);
    }
}

void find_scc()
{
    for (int i = 1; i <= n; ++i)
        if (!dfn[i]) Tarjan(i);
}

void init()
{
    for (int i = 0; i <= n; ++i)
        G[i].clear();
    Time = cnt = 0;
    memset(dfn, 0, sizeof(dfn));
    memset(low, 0, sizeof(low));
    memset(sccno,0,sizeof(sccno));
    memset(instk, 0, sizeof(instk));
    while (!S.empty()) S.pop();
}

int main()
{
    int k, u, v;
    while (~scanf("%d",&n)) {
        init();
        for (int i = 1; i <= n; ++i) {
            scanf("%d", &k);
            for (int j = 0; j < k; ++j) {
                cin >> v;
                G[i].push_back(v + n);
            }
        }
        for (int i = 1; i <= n; ++i) {
            scanf("%d", &v);
            G[v + n].push_back(i);
        }
        find_scc();
        for (int u = 1; u <= n; ++u) {
            k = 0;
            for (int i = 0; i < G[u].size(); ++i) {
                v = G[u][i];
                if (sccno[u] == sccno[v])
                    ans[k++] = v - n;
            }
            sort(ans, ans + k);
            printf("%d", k);
            for (int i = 0; i < k; ++i)
                printf(" %d",ans[i]);
            puts("");
        }
    }
    return 0;
}