强连通分量 模板

忘了哪道题了，直接发全部代码，需要时改着用。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<algorithm>
using namespace std;
int m,n,k,cas;
int x,y;
vector<int> ans;
int t;
vector<int> e[11000];///
vector<int> be[11000];
int sco[6000];
int belong[11000];///
int st[11000];///
int dep[6000],low[6000];///
bool vis[6000];///
int cc[11000];
int de,cnt;
int bb;///连通分量个数
int out[11000],in[11000];
int dfs(int a)
{
    dep[a]=low[a]=de++;
    st[cnt++]=a;
    vis[a]=1;
    for(int i=0;i<e[a].size();i++)
    {
        int v=e[a][i];
        if(dep[v]==-1)
        {
            dfs(v);
            low[a]=min(low[a],low[v]);
        }
        else if(vis[v])
        {
            low[a]=min(low[a],dep[v]);
        }
    }
    if(dep[a]==low[a])
    {
	int tt;
        do{
            tt=st[--cnt];
            //cout<<tt<<" "<<bb<<endl;
            belong[tt]=bb;
            cc[bb]++;
            vis[tt]=0;
        }while(tt!=a);
        bb++;
    }
}

/*
void fdfs(int root,int u)
{
    sco[root]+=cc[u];
    vis[u]=1;
    int num=unique(be[u].begin(),be[u].end())-be[u].begin();
    for(int i=0;i<num;i++)
    {
        int v=be[u][i];
        if(!vis[v])
        {
            fdfs(root,v);
        }
    }
    return ;
}
*/
int main()
{
    int t;
    scanf("%d",&t);
    for(int cas=1;cas<=t;cas++)
    {
        scanf("%d%d",&n,&m);
        for(int i=0;i<=n;i++)
            {e[i].clear();be[i].clear();}
        for(int i=0;i<m;i++)
        {
            scanf("%d%d",&x,&y);
            e[x].push_back(y);
        }
        de=1;cnt=0;bb=0;
        memset(dep,-1,sizeof(dep));memset(low,0,sizeof(low));
        memset(vis,0,sizeof(vis));memset(cc,0,sizeof(cc));
        memset(sco,0,sizeof(sco));
        for(int i=0;i<n;i++)///
        {

         if(dep[i]==-1)
        dfs(i);
        }
        memset(out,0,sizeof(out));
        memset(in,0,sizeof(in));
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<e[i].size();j++)
            {
                int v=e[i][j];
                if(belong[i]==belong[v]) continue;
                  out[belong[i]]++;
                  in[belong[v]]++;
                be[belong[v]].push_back(belong[i]);
            }
        }
        for(int i=0;i<bb;i++)
        {
            if(out[i]==0)
            {
                memset(vis,0,sizeof(vis));
                fdfs(i,i);
                sco[i]--;
            }
        }
        int maxn=-1;
        for(int i=0;i<bb;i++)
        {
            maxn=max(sco[i],maxn);
        }
        memset(vis,0,sizeof(vis));
        ans.clear();
        for(int i=0;i<bb;i++)
        {
            if(sco[i]==maxn)
            {
                vis[i]=1;
            }
        }
        for(int i=0;i<n;i++)
        {
            if(vis[belong[i]]==1)
            {
                ans.push_back(i);
            }
        }
        sort(ans.begin(),ans.end());
        printf("Case %d: %d\n",cas,maxn);
        for(int i=0;i<ans.size();i++)
        {
            if(i)printf(" ");
            printf("%d",ans[i]);
        }
        puts("");

    }
}