poj-2186 强连通分量

//poj 2186
//题意：奶牛之间有仰慕关系，并且这种关系具有传递性，问被所有奶牛仰慕的奶牛个数有多少？
//思路：构造一个图，这些图中有强连通分量，在这强连通分量中的奶牛肯定是被在这个强连通分量里面奶牛所仰慕的
//接下来我们要判断的就是在有没有一个或几个强连通被外面所有的强连通所仰慕。
//如果一个图里面有几个强连通分量的话那么一定至少存在一个出度为0的点，如果不存在那么刚刚的强连通一定不是最大强连通
//又想如果有两个这样的出度为0的点存在，那么这个图中一定没有满足题意的奶牛，因为那两个出度为0的点不会仰慕任何牛
//所以出度为0的点有且只能有一个！如果有多个那么满足题意的牛只有0个，如果没有的话那么整个图都是满足题意的牛
//又证：那个出度为0的点一定有入度（数据中没有涉及）；如果他有入度那么这个出度为0的强连通里面的牛都是满足题意的牛
//因为已经没有了出度为0的点，即除他之外的强连通点都有仰慕的，仰慕关系又可以传递，所以。。。。。
#include <iostream>
using namespace std;
#define MAXM 50010
#define MAXV 10010
#define min(a,b) (a>b?b:a)
typedef struct{
	int s,t,next,next2;
}Edge;
Edge edge[MAXM];
int n,m,headlist[MAXV],headlist2[MAXV];
int order[MAXV],belong[MAXV];
int num,count;
bool vis[MAXV];

void dfs(int x){
	int i,a;
	vis[x]=1;
	for(i=headlist[x];i!=-1;i=edge[i].next){
		a=edge[i].t;
		if(!vis[a]) dfs(a);
	}
	order[++num]=x;
}

void dfst(int x){
	int i,a;
	belong[x]=count;		//记录结点属于哪个连通分量,count从1开始记录
	vis[x]=1;
	for(i=headlist2[x];i!=-1;i=edge[i].next2){		//要将边反过来遍历一遍
		a=edge[i].s;
		if(!vis[a]) dfst(a);
	}
}

void kosaraju(){
	int i;
	memset(vis,0,sizeof(vis));
	num=count=0;
	for(i=1;i<=n;i++)
		if(!vis[i]) dfs(i);
		memset(vis,0,sizeof(vis));
		
	for(i=n;i>=1;i--)
		if(!vis[order[i]]){
			count++;
			dfst(order[i]);
		}
}
int ingree[MAXV];
void solve()
{
	int i,j,k,z;
	int ans;
	int lp=0;
	int sum;
	for(j=1;j<=count;j++)
	{
		ans=0;
		memset(ingree,0,sizeof(ingree));
		for(i=1;i<=n;i++)
		{
			if(belong[i]==j)
			{
				for(k=headlist[i];k!=-1;k=edge[k].next){
					if(belong[edge[k].t]!=j)
					ans++;
					
				}
				
			}
		}
		if(ans==0){
			lp++;
			sum=j;
		}
	}	
	if(lp==1)
	{
		ans=0;
		for(i=1;i<=n;i++)
		{
			if(belong[i]==sum)
				ans++;
		}
		printf("%d\n",ans);
	}
	else
		printf("%d\n",0);
}
int main()
{
	int i,a,b;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		memset(headlist,-1,sizeof(headlist));
		memset(headlist2,-1,sizeof(headlist2));
		for(i=0;i<m;i++)
		{
			scanf("%d%d",&a,&b);
			edge[i].s=a;
			edge[i].t=b;
			edge[i].next=headlist[a];
			edge[i].next2=headlist2[b];
			headlist[a]=i;
			headlist2[b]=i;
		}
		kosaraju();
		solve();
	}
}