强连通算法模板

#include<stdio.h>
#include<string.h>
#include<vector>
#include<algorithm>
using namespace std;
vector<int>vec;
const int N=1e5+10;
int head[N],book[N],cnt,rcnt,rhead[N],cmp[N];
struct l
{
	int next,to;
}dis[N],rdis[N];
void add(int from,int to)
{
	dis[++cnt].to=to;
	dis[cnt].next=head[from];
	head[from]=cnt;
}
void radd(int from,int to)
{
	rdis[++rcnt].to=to;
	rdis[rcnt].next=rhead[from];
	rhead[from]=rcnt;
}
void dfs(int u)
{
	int i;
	book[u]=1;
	for(i=head[u];i!=-1;i=dis[i].next)
	{
		if(book[dis[i].to]==0) dfs(dis[i].to);
	}
	vec.push_back(u);
}
void rdfs(int u,int k)
{
	int i;
	book[u]=1;
	cmp[u]=k;
	for(i=rhead[u];i!=-1;i=rdis[i].next)
	{
		if(book[rdis[i].to]==0) rdfs(rdis[i].to,k);
	}
}
int main()
{
	int n,m,i,j;
	cnt=-1,rcnt=-1;
	memset(head,-1,sizeof(head));
	memset(rhead,-1,sizeof(rhead));
	scanf("%d %d",&n,&m);
	for(i=0;i<m;i++)
	{
		int a,b;
		scanf("%d %d",&a,&b);
		add(a,b);
		radd(b,a);
	}
	for(i=1;i<=n;i++)
	{
		if(!book[i]) dfs(i);
	}
	memset(book,0,sizeof(book));
	int k=0;
	for(i=vec.size()-1;i>=0;i--)
	{
		if(!book[vec[i]]) rdfs(vec[i],k++);
	}
	printf("%d\n",k);
}