hdu 3836 Equivalent Sets 强联通 tarjan缩点

求强连通分量，然后统计每一个分量的入度和初度

这两个量中较大的那一个就是答案了

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stack>
#include<vector>
using namespace std;
const int N = 20010;
int n,index,ans;
vector<int>map[N];
stack<int>S;
int low[N],dfn[N];
int instack[N];
int in[N],out[N],belong[N];
void init()
{
    memset(low,0,sizeof(int)*(n+1));
    memset(dfn,0,sizeof(int)*(n+1));
    memset(in,0,sizeof(int)*(n+1));
    memset(out,0,sizeof(int)*(n+1));
    ans=0;
    index=0;
    for(int i=0;i<=n;i++)
    {
        map[i].clear();
        instack[i]=0;
    }
}
void tarjan(int u)
{
    dfn[u]=low[u]=++index;
    S.push(u);
    instack[u]=1;
    for(int i=0;i<map[u].size();i++)
    {
        int v=map[u][i];
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if(instack[v])
            low[u]=min(low[u],dfn[v]);
    }
    if(dfn[u]==low[u])
    {
        ans++;
        int v;
        do{
            v=S.top();
            S.pop();
            instack[v]=0;
            belong[v]=ans;
        }while(u!=v);
    }
}
void solve()
{
    for(int i=1;i<=n;i++)
    {
        if(!dfn[i])
        tarjan(i);
    }
    for(int u=1;u<=n;u++)
    {
        for(int j=0;j<map[u].size();j++)
        {
            int v=map[u][j];
            if(belong[u]!=belong[v])
            {
                in[belong[v]]++;
                out[belong[u]]++;
            }
        }
    }
    int t1=0,t2=0;

    for(int i=1;i<=ans;i++)
    {
        if(!in[i]) t1++;
        if(!out[i]) t2++;
    }
    if(ans<=1)
        printf("0\n");
    else
    printf("%d\n",max(t1,t2));

}
int main()
{
    int T,m;
        while(scanf("%d%d",&n,&m)!=EOF)
        {

        init();
        for(int i=1;i<=m;i++)
        {
            int a,b;
            scanf("%d%d",&a,&b);
            map[a].push_back(b);
        }
        solve();

    }
    return 0;
}