poj 1904 强连通图

题目链接：http://poj.org/problem?id=1904

题目意思：有n个女生和n个男生,给定一些关系表示男生喜欢女生(即两个人可以结婚),再给定一个初始匹配,表示这个男生和哪个女生结婚,初始匹配必定是合法的.求每个男生可以和哪几个女生可以结婚且能与所有人不发生冲突。
思路：将男生从1到n编号,女生从(n+1)到2*n编号,输入时如果男生u可以和女生v结婚,那么就做一条从u到v的边(u,v),对于输入的初始匹配(u,v)(表示男生u和女生v结婚),那么从v做一条到u的边(v,u),然后求解图的强连通分量,如果男身i和女生j在同一个强连通分量内,且i到j有直接的边。则他们可以结婚。

给定2*n个顶点的二分图以及一个完美匹配，确定所有的点对(u,v)，使去除u,v之后的二分图仍有完美匹配。

首先，把原图的边改成有向的（从男到女）。然后，对于给定的匹配(u,v)，添加边v->u。如果两个点在同一个强连通分量里，他们就是合法的点对。

对于一个有完美匹配的二分图。
1:任意找出一个完美匹配。
2:在原来的边不变的情况下。对于求出的任一个完美匹配i（左），j（右），则连接反向边j到i。
3：求强联通。
4：对于每一个联通块。
按二分图的左右分成两堆A和B。
所有在A中点i，如果和B中的点j有直接的边i到j（不是j到i）
那么i就可以和j匹配。
5：所有满足条件的就是i的所有的可能的解。

#include<cstdio> #include<cstring> #include<stack> using namespace std; #define N 2005 #define Max 4010 //失误，一直开的都是N，wrong了两次 int dfn[Max],low[Max],belong[Max]; int instack[Max],index,scc,n,tmp[N]; stack<int>S; typedef struct node { int to; node *next; }node; node *g[Max]; int cmp(const void *a,const void *b) { return *(int *)a-*(int *)b; } void tarjan(int x) { int i,v; node *p; dfn[x]=low[x]=index; index++; S.push(x); instack[x]=1; p=g[x]; while(p) { i=p->to; if(dfn[i]==-1) { tarjan(i); low[x]=low[x]>low[i]?low[i]:low[x]; } else if(instack[i]) low[x]=low[x]>dfn[i]?dfn[i]:low[x]; p=p->next; //第一次忘写了，导致死循环 } if(low[x]==dfn[x]) { while(1) { v=S.top(); S.pop(); belong[v]=scc; instack[v]=0; if(v==x) break; } scc++; } } int main() { int i,j,k,t; node *p; while(scanf("%d",&n)!=EOF) { scc=index=0; memset(g,0,sizeof(g)); memset(dfn,-1,sizeof(dfn)); memset(instack,0,sizeof(instack)); for(i=1;i<=n;i++) { scanf("%d",&t); while(t--) { p=new node; scanf("%d",&j); p->to=n+j; p->next=g[i]; g[i]=p; } } for(i=1;i<=n;i++) { p=new node; scanf("%d",&j); p->to=i; p->next=g[n+j]; //注意存储是n+j g[n+j]=p; } for(i=1;i<=n*2;i++) //此处注意，应到n*2 if(dfn[i]==-1) tarjan(i); for(i=1;i<=n;i++) { p=g[i]; k=0; while(p) { j=p->to; if(belong[j]==belong[i]) tmp[k++]=j-n; p=p->next; } qsort(tmp,k,sizeof(tmp[0]),cmp); //结果需排序 printf("%d",k); for(j=0;j<k;j++) printf(" %d",tmp[j]); printf("/n"); } } return 0; }