HAOI2006受欢迎的牛——缩点

题目：http://poj.org/problem?id=2186

本题是缩点模板题，将强连通分量缩成一个点，从而形成一个有向无环图，当仅有一个出度为0的点时答案即此点的大小，否则无解。

代码如下：

#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;
}

　　其中注释掉的部分以及dfs的部分是不使用上述结论的暴力枚举；

　　但有一出不解：即第96行判断条件若是“==n”则有些有解的情况输出0导致WA，不知何故。