【模板】tarjan算法

#include<bits/stdc++.h> 
#define MAX 10005 
using namespace std;  
int dfn[MAX],low[MAX],inst[MAX],st[MAX],tot,sum,top;
vector<int>G[MAX];
void tarjan(int u)
{
    int v;
    dfn[u]=low[u]=++tot;
    st[++top]=u;
    inst[u]=1;
    for(int i=0;i<G[u].size();i++)
    {
        v=G[u][i];
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]); //找出父子结点对应关系 
        }
        else if(inst[v])               //如果访问过，并且还在栈里 
            low[u]=min(low[u],dfn[v]); //连接父子结点对应关系 
    }
    if(dfn[u]==low[u])                 //找到强连通分量子树里的最小根 
    {
        sum++;                         //强连通分量个数sum 
        do{
            v=st[top--];
            inst[v]=0;
        }while(u!=v);
    }
}
void init()
{
    for(int i=0;i<MAX;i++)
        G[i].clear();
    top=-1;
    sum=tot=0;
    memset(inst,0,sizeof(inst)); 
    memset(dfn,0,sizeof(dfn)); 
    memset(low,0,sizeof(low));     
}
int main()  
{  
    int n,m,x,y,i;  
    while(cin>>n>>m) 
    {
        init();
        for(i=1;i<=m;i++)  
        {  
            scanf("%d%d",&x,&y);  
            G[x].push_back(y);
        }  
        for(i=1;i<=n;i++) 
            if(!dfn[i])
                tarjan(i);
    }
    return 0;  
}