poj2186 Popular Cows(强连通分量tarjan)

题中给出一张有向图，不保证图连通，如果对于某个点，其他点都有至少一条路径能到达该点，该点就是符合要求的（ Popular Cows），问图中有多少这样的点。

解法：求强连通分量，对强连通分量缩点，缩点后是一棵树，如果出度为0的缩点个数大于1，则符合要求的点为0个，否则出度为0的缩点（强连通分量）包含的点的个数就是答案。

求强连通分量时和双连通分量时，tarjan算法还是不同的！无向图中tarjan回边不包括反向边，而有向图，反两边就是另外的一条边，可以作为回边优化low值。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<queue>
#include<stack>
#include<cstring>
#include<map>
#include<string>
#define LL long long
#define clean(a) memset(a,0,sizeof(a))
#define maxn 11000
#define eps 1e-8

using namespace std;
int dfn[maxn],low[maxn],_dfn;
bool vis[maxn];
typedef vector<int>ve;
vector<ve>head(maxn);
stack<int>sta;
int area_num;
int area[maxn];
int deg[maxn];
void tarjan(int rt,int father)
{
	sta.push(rt);
	dfn[rt]=low[rt]=_dfn++;
	vis[rt]=true;
	for(int i=0,len=head[rt].size();i<len;++i)
	{
		int v=head[rt][i];
		if(!vis[v])
		{
			tarjan(v,rt);
			low[rt]=min(low[rt],low[v]);
		}
		else//有向图不用判father
			low[rt]=min(low[rt],dfn[v]);
	}
	if(low[rt]==dfn[rt])
	{
		while(sta.top()!=rt)
		{
			area[sta.top()]=area_num;
			sta.pop();
		}
		sta.pop();
		area[rt]=area_num++;
	}
}
int main()
{
	int n,m;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		for(int i=0;i<=n;++i)
			head[i].clear();
		for(int i=0;i<m;++i)
		{
			int a,b;
			scanf("%d%d",&a,&b);
			head[a].push_back(b);
		}
		while(!sta.empty()) sta.pop();
		area_num=1;
		_dfn=0;
		memset(area,0,sizeof(area));
		memset(deg,0, sizeof(deg));
		memset(dfn,0,sizeof(dfn));
		memset(low,0,sizeof(low));
		memset(vis,false,sizeof(vis));
		for(int i=1;i<=n;++i)
		{
			if(!vis[i])
			{
				tarjan(i,1);
			}
		}
		for(int i=1;i <=n;++i)
		{
			vector<int>&temp=head[i];
			for(int j=0,len=temp.size();j<len;++j)
			{
				int v=temp[j];
				if(area[i]!=area[v])
				{
					deg[area[i]]++;
				}
			}
		}
		int num=0;
		int area_name=0;
		for(int i=1;i<area_num;++i)
		{
			if(deg[i]==0)
			{
				area_name=i;
				num++;
			}
		}
		if(num>1 )
			printf("0\n");
		else
		{
			int ans=0;
			for(int i=1;i<=n;++i)
			{
				if(area[i]==area_name)
					ans++;
			}
			printf("%d\n",ans);
		}
	}
}