POJ 1904 King's Quest

http://poj.org/problem?id=1904

方法同HDU 4685 2013多校8 1010，只是那道题需要自己找匹配，这道题已经给定了匹配了。

根据初始给定的匹配直接找强连通分量，在同一强连通分量的人可以配对。

这题试了一下std::ios_base::sync_with_stdio(false);，加了这句cin比scanf竟然还快。。。。。。

//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<iostream>
#include<cstring>
#include<string>
#include<cmath>
#include<set>
#include<list>
#include<map>
#include<iterator>
#include<cstdlib>
#include<vector>
#include<queue>
#include<stack>
#include<algorithm>
#include<functional>
using namespace std;
typedef long long LL;
#define ROUND(x) round(x)
#define FLOOR(x) floor(x)
#define CEIL(x) ceil(x)
const int maxn=4040;
const int maxm=400010;
const int inf=0x3f3f3f3f;
const LL inf64=0x3f3f3f3f3f3f3f3fLL;
const double INF=1e30;
const double eps=1e-6;
int N,g[maxn/2][maxn/2];
int Mx[maxn];
struct Edge
{
    int v,next;
}edge[maxm];
int edgeNum,head[maxn];
void addSubEdge(int u,int v)
{
    edge[edgeNum].v=v;
    edge[edgeNum].next=head[u];
    head[u]=edgeNum++;
}
int sid[maxn];
int mark[maxn],low[maxn];
int check[maxn];
int sstack[maxn],top;
int dfn,ssn;
int n,m;
void dfs(int k)
{
    int i,j;
    check[k]=1;
    low[k]=mark[k]=dfn++;
    sstack[top++]=k;
    for(int i=head[k]; i!=-1; i=edge[i].next)
    {
        int j=edge[i].v;
        if(mark[j]==0)
        {
            dfs(j);
            low[k]=min(low[k],low[j]);
        }
        else if(check[j])
            low[k]=min(low[k],mark[j]);
    }
    if(mark[k]==low[k])
    {
        while(sstack[--top]!=k)
        {
            check[sstack[top]]=0;
            sid[sstack[top]]=ssn;
        }
        sid[k]=ssn;
        check[k]=0;
        ++ssn;
    }
    return;
}
void tarjan()
{
    ssn=1;
    dfn=1;
    top=0;
    memset(check,0,sizeof(check));
    memset(mark,0,sizeof(mark));
    for(int i=0; i<n; ++i) if(mark[i]==0) dfs(i);
}
void init()
{
    memset(g,0,sizeof(g));
    memset(head,-1,sizeof(head));
    edgeNum=0;
}
void input()
{
    for(int i=0; i<N; i++)
    {
        int k;
//        scanf("%d",&k);
        cin>>k;
        while(k--)
        {
            int x;
            cin>>x;
//            scanf("%d",&x);
            g[i][x-1]=1;
            addSubEdge(i,x+N-1);
        }
    }
    for(int i=0;i<N;i++)
    {
        int x;
        cin>>x;
//        scanf("%d",&x);
        addSubEdge(x+N-1,i);
    }
}
void solve()
{
    n=N;
    tarjan();
}
void output()
{
    vector<int> ans;
    for(int i=0;i<N;i++)
    {
        ans.clear();
        for(int j=N;j<2*N;j++)
        {
            if(g[i][j-N]&&sid[i]==sid[j])
            {
                ans.push_back(j-N+1);
            }
        }
        cout<<ans.size();
//        printf("%d",ans.size());
        for(int j=0;j<ans.size();j++)
        {
            cout<<" "<<ans[j];
//            printf(" %d",ans[j]);
        }
        cout<<endl;
//        puts("");
    }
}
int main()
{
    std::ios_base::sync_with_stdio(false);
//    freopen("in.cpp","r",stdin);
    while(cin>>N)
    {
        init();
        input();
        solve();
        output();
    }
    return 0;
}