tarjan强联通分量（模板）

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define N 10000
using namespace std;
int x,y,n,m,t,tot,sum,top,time;
int head[N],col[N],stack[N],dfn[N],low[N],a[N][N];
bool vis[N];
struct Edge
{
    int from,next,to;
 }edge[N];
int add(int x,int y)
{
    tot++;
    edge[tot].to=y;
    edge[tot].next=head[x];
    head[x]=tot;
}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return f*x;
}
int tarjan(int now)
{
//stack[]表示递归过程的栈，即用来判断该点是否已经加入到此次递归的栈中，在递归末尾，通过将vis置为false释放所有递归栈的元素
    t=0;
    dfn[now]=low[now]=++time;//初始每一个点的low值dfn等于它的时间戳 
    stack[++top]=now; vis[now]=true;//将该点入栈，标记为在栈中 
    for(int i=head[now];i;i=edge[i].next)//更新于他相连的点的low值 
    {
        x=edge[i].to;
        if(vis[x]) low[now]=min(dfn[x],low[now]);//如果该点已经在栈中,直接更新来到该点的那个点的low，不需要递归查询 
        else if(!dfn[x])
        {
            tarjan(x);
            low[now]=min(low[x],low[now]);//不在栈中，需要从该点继续递归拓展 
        }
    }
    if(low[now]==dfn[now])//说明以这个点结束强连通分量 
    {
        sum++;// 强连通分量的个数加一 
        col[now]=sum;//将该点放在她所属的强连通分量了 
        for(;stack[top]!=now;top--)
        {
            col[stack[top]]=sum;
            vis[stack[top]]=false;
        }
        vis[now]=false;
        top--;
     } 
}
int main()
{
    n=read(),m=read();
    for(int i=1;i<=m;i++)
     {
         x=read();y=read();
         add(x,y);
     }
    for(int i=1;i<=n;i++)
      if(!dfn[i]) tarjan(i);
    printf("%d",sum);
    return 0;
}