tarjan算法模版

#include<iostream>
#define num 5000
#include<cstring>
using namespace std;

int n,m;
int f[num];


typedef struct{
        int u,v,next;
        }z;z edg[num];
int edgs=0,head[10000];
void Addedg(int a,int b){
     edgs++;
     edg[edgs].u=a;edg[edgs].v=b;
     edg[edgs].next=head[a];
     head[a]=edgs;
     }

void Init(){
     cin>>n>>m;
     for(int i=1;i<=n;++i)
     f[i]=i;
     
     memset(head,-1,sizeof(head));
     for(int i=1;i<=m;++i)
     {
       int a,b;
       cin>>a>>b;
       Addedg(i,b);
             }
     }

bool mark[10000],vis[10000];
int cordnode[100000],top=-1,times=0;
int mintime[100000];
void Dfs(int k){
     int now;
     now=mintime[k]=times++;
     vis[k]=1;
     cordnode[++top]=k;
     
     for(int i=head[k];i!=-1;i=edg[i].next)
     {
       int v=edg[i].v;
       if(mark[v]) continue;
       if(!vis[v]) Dfs(v);
       if(mintime[k]>mintime[v]) 
       mintime[k]=mintime[v]; 
             }
     
     if(mintime[k]==now)
     {
       while(cordnode[top]!=k)
       {
         f[cordnode[top]]=k;
         mark[cordnode[top]]=1;
         top--;               
                              } 
       f[cordnode[top]]=k;
       top--; 
                        }
     
     
     }

int main()
{
    Init();
    
    for(int i=1;i<=n;++i)
    if(!mark[i]) Dfs(i);
    
    int p[10000]={0};
    
    for(int i=1;i<=n;++i)
    p[f[i]]=1;
    int ans=0;
    for(int i=1;i<=n;++i)
    if(p[i]) ans++;
    cout<<ans<<endl;
    system("pause");
    
    }