【poj】1904 King's Quest【强连通】

题意：国王有n个王子,王国内有n个妹子,每个王子喜欢一个（些）妹子（们）,给出了一个王子与妹子之间的完美匹配.问一个王子可以匹配的妹子有哪些,要求改王子匹配改妹子之后,其他王子仍有妹子可以匹配,仍能够形成完美匹配

题解：题目真是非常的具有误导性,十分容易考虑到二分匹配上去,实际上确是强连通.

我们考虑转换原问题,,我们从王子向他喜欢的妹子(们)连一条边,对给出的完美匹配,我们对每一对匹配的王子与妹子,从妹子向王子连一条边,在进行强连通,这样每一个强连通分量里面一定会有等数量的王子和妹子

(证明:该图是一个二分图,1如果王子<妹子,那么有妹子无法到达任何一个王子,这与题设相悖2,如果王子>妹子,那么妹子们无法到达某些王子,这与题设相悖,故王子=妹子)，

这样答案就是每个王子所在强连通分量的所有妹子,因为王子选择联通分量中的妹子1.同一个连通分量中的王子可以到达任何一个改连通分量中的妹子2.没有影响不同连通分量的

（证明：王子为什么不能选择不同强连通分量的妹子。反证法：如果强连通分量1中的王子选择了强连通分量2中的妹子,那么势必强连通分量2中的一个王子无法在自己的强连通分量中找到妹子,那么他就会去别的强连通分量找妹子,这样一直循环下去，我们知道最终一定是经过了强连通分量1,2,x1,x2,xn,……,1，王子们才能都找到自己的妹子，这样这些强连通分量1,2,x1,x2,xn,……,1会构成一个强连通分量,与题设在不同强连通分量中找妹子不符)

推荐:hdu4685 该题的加强版就是初始的完美匹配没给出 需要自己来算

#include<queue>
#include<vector>
#include<cstdio>//先用强连通缩点来化简图，然后在图上做拓扑排序，如果排序过程中，出现1个以上的点入度同时为0时，那么就不满足条件。
#include<cstring>
#include<algorithm>
#define PB push_back
#define MP make_pair
using namespace std;
const int N=4010;
const int inf=0x3f3f3f3f;
int dfn[N],low[N],instack[N],belong[N],S[N];
int index,scc,shead;
vector<int>G[N];
vector<pair<int,int> >E[N];
void dfs(int u)
{
    dfn[u]=low[u]=++index;
    S[shead++]=u;
    instack[u]=true;
    for(int i=0;i<G[u].size();i++){
        int v=G[u][i];
        if(!dfn[v]){
            dfs(v);
            low[u]=min(low[v],low[u]);
        }else if(instack[v])
            low[u]=min(dfn[v],low[u]);
    }
    if(low[u]==dfn[u]){
        scc++;
        while(1){
            int tmp=S[--shead];
            belong[tmp]=scc;
            instack[tmp]=false;
            if(tmp==u)break;
        }
    }
}
void tarjan(int n)
{
    index=scc=shead=0;
    memset(dfn,0,sizeof(dfn));
    memset(instack,false,sizeof(instack));
    for(int i=1;i<=n;i++)if(!dfn[i])dfs(i);
}
int ans[N];
int main()
{
    int n;
    while(~scanf("%d",&n)){
        for(int i=1;i<=2*n;i++)G[i].clear();
        for(int i=1;i<=n;i++){
            int x,y;
            scanf("%d",&x);
            while(x--){
                scanf("%d",&y);
                G[i].PB(n+y);
            }
        }
        for(int i=1;i<=n;i++){
            int x;
            scanf("%d",&x);
            G[x+n].PB(i);
        }
        tarjan(2*n);
        for(int i=1;i<=n;i++){
            int cnt=0;
            for(int j=0;j<G[i].size();j++){
                int v=G[i][j];
                if(belong[i]==belong[v]){
                    ans[cnt++]=v-n;
                }
            }
            sort(ans,ans+cnt);
            printf("%d",cnt);
            for(int j=0;j<cnt;j++){
                printf(" %d",ans[j]);
            }
            printf("\n");
        }
    }
    return 0;
}