POJ 1904（强连通分量）Tarjan

题意描述：有n个王子和n个美女，每个王子有自己喜欢的美女，保证每个王子都能娶到一个自己喜欢的美女，现在巫师给出了一个王子和美女匹配的方案，但是国王想知道每个王子能够娶到的美女的名单（保证每个王子都能娶到美女）

 

思路：一看就是二分图，求所有可行的完全匹配的边集，但是想想二分图里面没有可以直接或者间接套用的算法。

不知道大家看到数据提供的一组匹配时是什么感想，应该要猜测到这组数据肯定有用。

那么假设现在我们按照数据提供的匹配已经匹配好了，这是对于u王子，如果他和v美女（非原配）匹配的话，到最后，则必然要有另一个王子和u的原配美女匹配

想想，发现这是一个环，想到强连通了么？？

先别看以下的内容，先思考再看！

构图：每个王子向他喜欢的美女连一条有向边，对于数据中的匹配，每个美女向其匹配的王子连一条有向边，然后强连通就行了。

PS：注意在同一个强连通分量中的美女并不是都能和在该连通分量王子匹配（王子不一定喜欢她，为此狂WA，单恋就是不好。。。）

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#define N 5000
#define M 700000 
using namespace std;
int n,t,divg,cnt,p,stack[N],head[N],to[M],next[M],dfn[N],low[N],belong[N],map[2001][2001];
bool fg[N];
inline void add(int u,int v) 
{
    to[cnt]=v; next[cnt]=head[u]; head[u]=cnt++;
}
void read()
{
    memset(fg,0,sizeof fg);
    memset(head,-1,sizeof head); cnt=0;
    memset(dfn,0,sizeof dfn);
    memset(belong,0,sizeof belong);
    memset(map,0,sizeof map);
    t=0; divg=0; p=0;
    for(int i=1,a,b;i<=n;i++)
    {
        scanf("%d",&a);
        while(a--) scanf("%d",&b),add(i,n+b),map[i][b]=true;
    }
    for(int i=1,a;i<=n;i++) scanf("%d",&a),add(a+n,i);
}
void dfs(int u)
{
    t++;
    dfn[u]=low[u]=t;
    stack[++p]=u; fg[u]=true;
    for(int i=head[u];~i;i=next[i])
    {
        if(!dfn[to[i]])
        {
            dfs(to[i]);
            low[u]=min(low[u],low[to[i]]);
        }
        else if(fg[to[i]]) low[u]=min(low[u],dfn[to[i]]);
    }
    if(dfn[u]==low[u])
    {
        divg++;
        int tmp=-1;
        while(tmp!=u)
        {
            tmp=stack[p--];
            belong[tmp]=divg;
            fg[tmp]=false;
        }
    }
}
void go()
{
    for(int i=1;i<=n;i++)
        if(!dfn[i]) dfs(i);
    for(int i=1,num;i<=n;i++)
    {
        num=0;
        for(int j=n+1;j<=n+n;j++)
            if(belong[j]==belong[i]&&map[i][j-n]) num++;
        printf("%d ",num);
        for(int j=n+1;j<=n+n;j++)
            if(belong[j]==belong[i]&&map[i][j-n]) printf("%d ",j-n);
        printf("\n");
    }
}
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        read();
        go();
    }
    return 0;
}