POJ 2186 Popular Cows

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

题目大意：

    每头牛的梦想就是成为牛群中最受欢迎的牛。 在一群N（1 <= N <= 10,000）母牛中，你可以得到M（1 <= M <= 50,000）有序的形式对（A，B），告诉你母牛A认为母牛 B很受欢迎。 由于流行是传递性的，如果A认为B很受欢迎，B认为C受欢迎，那么A也会认为C是受欢迎的，即使这不是输入中有序对明确规定的。 你的任务是计算每头奶牛认为受欢迎的奶牛数量。

解题思路：

    题目中的突破口为这一句话： 如果A认为B很受欢迎，B认为C很受欢迎，那A也会认为C受欢迎。 故需要先利用强连通分量进行缩点，然后分情况讨论：

    <1>若图中只有1个强连通分量（即为强连通图），则所有牛都受欢迎，输出牛的总数。

    <2>若图中有多个强连通分量，则统计缩点后图中 出度为0的点的数量（因为出度为0肯定是公认受欢迎的，假设出度不为0，就肯定还有比他更受欢迎的）：

    1.出度为0的点的数量为0：输出牛的总数

    2.出度为0的点的数量为1：输出该连通分量中点的总数（这个连通分量内的牛是最受欢迎的）

    3.出度为0的点的数量为2及以上：输出0（一山不容二虎）

ps:缩点后，还需要判断一下图是否连通（所有顶点的出度之和是否等于边数），若不等，输出也是0。

pps:刚开始写缩点后的讨论代码时，有点顾虑把图遍历一遍会不会超时，仔细看了看题目，发现最多50000条边，然后就不要大意的上了。

ppps:感觉有向图求强连通分量稍微有那么一点点像并查集~~~

pppps：现在才刚学Tarjan算法，以后有时间其他算法也写一遍。

下面是用Tarjan算法写的代码： 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
const int maxn=10000+5;
vector <int> pic[maxn];
int low[maxn],dfn[maxn],vis[maxn];
int st[maxn],top=0;
int cnt=0,co=0;
int indeg[maxn];//统计出度
int lian[maxn];//标记属于哪个连通分量 
void tarjan(int begin)
{
	int i;
	dfn[begin]=low[begin]=++cnt;//初始化 
	vis[begin]=1;//标记在栈内 
	st[top++]=begin;//入栈 
	for(i=0;i<pic[begin].size();i++)//如果出度为0，则不更新 
	{
		int tp=pic[begin][i];
		if(dfn[tp]==0)//从来没访问过 
		{
			tarjan(tp);
			low[begin]=min(low[begin],low[tp]);//回溯时更新 
		}
		else if(vis[tp])//在栈中 
			low[begin]=min(low[begin],dfn[tp]);
	}
	if(dfn[begin]==low[begin])//关键点 
	{
		co++;//计数器 
		do{
			i=st[--top];
			vis[i]=0;
			lian[i]=co;		
		}while(begin!=i);
	}
}
int main()
{
	int n,m;
	int begin,end;
	while(~scanf("%d%d",&n,&m)&&(n+m))
	{
		memset(dfn,0,sizeof(dfn));
		memset(low,0,sizeof(low));
		memset(vis,0,sizeof(vis));
		memset(indeg,0,sizeof(indeg));
		co=0,cnt=0,top=0;
		
		for(int i=0;i<=n;i++)//清空vector 
			pic[i].clear();
		while(m--)//加边 
		{
			scanf("%d%d",&begin,&end);
			pic[begin].push_back(end);
		}
		for(int i=1;i<=n;i++)//每一个点都搜一次 
		{
			if(dfn[i]==0)
			tarjan(i);
		}
		int sum=0,sum1=0,sign=0;
		if(co==1)//只有一个强连通分量 
			printf("%d\n",n);
		else
		{
			int tol=0;
			for(int i=1;i<=n;i++)
			{
				for(int j=0;j<pic[i].size();j++)
				{
					if(lian[i]!=lian[pic[i][j]])//两个点分属不同连通分量
					{
						tol++;//边数 
						indeg[lian[i]]++;//统计出度
					} 
				}
			}
			for(int i=1;i<=co;i++)
			{
				sum1+=indeg[i];//所有点出度 
				if(indeg[i]==0)//出度为0的点的个数 
				{
					sign=i;
					sum++;
				}
			}
			if(sum1==tol)//所有顶点的出度之和等于边数 
			{
				if(sum==0)//没有出度为0的，就是所有人 
				printf("%d\n",n);
				else if(sum==1)//有一个是该连通分量里的所有人 
				{
					int c=0;
					for(int i=1;i<=n;i++)
					{
						if(lian[i]==sign)
						c++;
					}
					printf("%d\n",c);
				}
				else//2个及以上为0 
				printf("0\n");
			}
			else
				printf("0\n");
		}
	}
	return 0;
}