poj-1904

题目来源：http://poj.org/problem?id=1904

大意：N个王子和N个美女配对，每个王子喜欢一个或多个美女，大臣已经配好了一个列表，但是国王想知道每个王子和美女的可行匹配。就是求出每个王子可以选择哪些美女

思路：对于每对男女u,v，连边(u,u+n)，而输入的最后给出一个可行匹配，则对于每个美女v和配对的王子u，连边(v+n,u)，然后求出图的强联通分量，最后计数输出就是了。

这是用vector储存的程序，5000ms+

#include <stdio.h>
#include <string>
#include <vector>
#include <algorithm>
using namespace std;

#define NN 5005
vector <int> a[NN];
vector <int> b[NN];
int real[NN],low[NN],p[NN],color[NN];
bool f[NN];
int n,top,total,timee;

void tarjan(int h)
{
     p[++top]=h;
     f[h]=true;
     real[h]=low[h]=++timee;
     for (int i=0;i<a[h].size();i++)
     {
         if (!real[a[h][i]])
         {
                         tarjan(a[h][i]);
                         low[h]=min(low[h],low[a[h][i]]);
         }
         else if (f[a[h][i]]) low[h]=min(low[h],real[a[h][i]]);
     }
     
     if (real[h]==low[h])
     {
                         total++;
                         while (p[top+1]!=h)
                         {
                               color[p[top]]=total;
                               f[p[top]]=false;
                               top--;
                         }
     }
}

int main()
{
    scanf("%d",&n);
    int m,x;
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&m);
        for (int j=0;j<m;j++)
        {
            scanf("%d",&x);
            a[i].push_back(x+n);
        }
    }
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&x);
        a[x+n].push_back(i);
    }
    
    memset(real,0,sizeof(real));
    memset(low,0,sizeof(low));
    memset(p,0,sizeof(p));
    memset(f,false,sizeof(f));
    top=total=timee=0;
    for (int i=1;i<=n;i++)
        if (!real[i]) 
           tarjan(i);
    
    for (int i=1;i<=n;i++)
        b[i].clear();
    for (int i=1;i<=n;i++)
    {
        for (int j=0;j<a[i].size();j++)
            if (color[i]==color[a[i][j]])
            b[i].push_back(a[i][j]-n);
        sort(b[i].begin(),b[i].end());
    }
    for (int i=1;i<=n;i++)
    {
        printf("%d",b[i].size());
        for (int j=0;j<b[i].size();j++)
            printf(" %d",b[i][j]);
        printf("\n");
    }
    return 0;
}

这是用池子法（前向星法）写的程序，3000ms+

#include <stdio.h>
#include <cstring>
#include <algorithm>
using namespace std;

#define maxn 20005
#define maxm 2000005
struct node {
       int y,next;
}p[maxm];

int head[maxn],a[maxm];
int n,m,cl;
bool instack[maxn];
int real[maxn],low[maxn],stack[maxn],color[maxn];
int depth,top,total;

void add(int x,int y)
{
     p[++cl].y=y;
     p[cl].next=head[x];
     head[x]=cl;
}

void dfs(int x)
{
     int y;
     real[x]=low[x]=++depth;
     instack[x]=true;
     stack[++top]=x;
     for (int i=head[x];i!=-1;i=p[i].next)
     {
         y=p[i].y;
         if (!real[y])
         {
                      dfs(y);
                      low[x]=min(low[x],low[y]);
         }
         else if (instack[y])
              low[x]=min(low[x],real[y]);
     }
     
     if (real[x]==low[x])
     {
                         total++;
                         while (stack[top+1]!=x)
                         {
                               y=stack[top--];
                               instack[y]=false;
                               color[y]=total;
                         }
     }
}
     

void tarjan()
{
     memset(real,0,sizeof(real));
     memset(low,0,sizeof(low));
     memset(instack,false,sizeof(instack));
     depth=top=total=0;
     for (int i=1;i<=n*2;i++)
         if (!real[i])
            dfs(i);
}

void solve()
{
     int y;
     for (int i=1;i<=n;i++)
     {
         top=0;
         for (int j=head[i];j!=-1;j=p[j].next)
         {
             y=p[j].y;
             if (color[y]==color[i])
                a[top++]=y-n;
         }
         printf("%d",top);
         sort(a,a+top);
         for (int j=0;j<top;j++)
             printf(" %d",a[j]);
         printf("\n");
     }
}

int main()
{
    scanf("%d",&n);
    int x,y;
    memset(head,-1,sizeof(head));
    cl=0;
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&y);
        while (y--)
        {
              scanf("%d",&x);
              add(i,x+n);
        }
    }
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&x);
        add(x+n,i);
    }
    tarjan();
    solve();
    return 0;
}