HAOI2006受欢迎的牛——缩点

#include<iostream>
#include<cstdio>
using namespace std;
int n,m,a,b,dfn[10005],low[10005],ct,cnt[10005],hp[10005],p,s[10005],ans,tim;
int col[10005],cr,rd[10005],cd[10005],siz[10005];
bool vis[10005],lin[10005][10005],fl[10005];
struct N{
	int hd;int to;int next;
}edge[50005];
void add(int x,int y)
{
	ct++;
	edge[ct].hd=x;
	edge[ct].to=y;
	edge[ct].next=cnt[x];
	cnt[x]=ct;
}
void tarjan(int x)
{
	vis[x]=1;
	hp[++p]=x;
	dfn[x]=++tim;low[x]=dfn[x];
	for(int i=cnt[x];i;i=edge[i].next)
	{
		int v=edge[i].to;
		if(!dfn[v])
		{
			tarjan(v);
			low[x]=min(low[x],low[v]);
		}
		else if(vis[v])
		{
			low[x]=min(low[x],low[v]);
		}
	}
	if(dfn[x]==low[x])
	{
		cr++;
		vis[x]=0;
		int k=0;
		while(hp[p]!=x)
		{
			vis[hp[p]]=0;
			col[hp[p]]=cr;
			p--;
			k++;
		}
		col[x]=cr;
		p--;k++;
		s[cr]=k;
		siz[cr]=k;
	}
}
void ser(int nw,int jl)
{
	s[nw]+=jl;
	if(!cd[nw])return;
	int k=0;
	for(int i=1;i<=n;i++)
	{
		if(k==cd[nw])break;
		if(lin[nw][i])
		{
			if(!fl[nw])ser(i,s[nw]);//不是jl+s[nw]!
			else ser(i,jl);
			k++;
		}
	}
	fl[nw]=1;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&a,&b);
		add(a,b);
	}
	for(int i=1;i<=n;i++)
		if(!dfn[i])tarjan(i);
	for(int i=1;i<=ct;i++)
	{
		int u=edge[i].hd;
		int v=edge[i].to;
		if(col[u]!=col[v])
			lin[col[u]][col[v]]=1,rd[col[v]]++,cd[col[u]]++;
	}
	int kt=0;int bh;
	for(int i=1;i<=cr;i++)
		if(!cd[i])kt++,bh=i;
	if(kt==1)printf("%d",siz[bh]);
	else printf("0");
//	for(int i=1;i<=cr;i++)
//		if(!rd[i])ser(i,0);
//	for(int i=1;i<=n;i++)
//		if(s[col[i]]>=n)ans++;//>=
//	printf("%d",ans);
	return 0;
}