POJ 2186 Popular Cows

http://poj.org/problem?id=2186

题意：有向图求能被所有其他点到达的点的个数。

题解：强连通缩点，求出度为0的分量，如果只有一个输出分量大小，如果大于一个，则这两个分量不连通，不存在这样的点。

代码：

#include <iostream>
#include<cstdio>
#include<stack>
#include<cstring>
using namespace std;
const int maxn = 10005;
const int maxm = 50005;
int n,m;
struct node{
	int to,next;
}e1[maxm],e2[maxm],e3[maxm];
//bool mp[maxn][maxn];
int head1[maxn],head2[maxn],head3[maxn];
int belong[maxn];
bool vis1[maxn],vis2[maxn];
int tot1,tot2,tot3;
int cnt1,cnt2;
int st[maxn];
int num;
int setnum[maxn];
void add(int u,int v){
	e1[tot1].to=v;
	e1[tot1].next=head1[u];
	head1[u]=tot1++;
	e2[tot2].to=u;
	e2[tot2].next=head2[v];
	head2[v]=tot2++;
}
void add2(int u,int v){
	e3[tot3].to=v;
	e3[tot3].next=head3[u];
	head3[u]=tot3++;
}
void dfs1(int u){
	vis1[u]=true;
	int i;
	for(i=head1[u];i!=-1;i=e1[i].next){
		if(!vis1[e1[i].to]){
			dfs1(e1[i].to);
		}
	}
	st[cnt1++]=u;
}
void dfs2(int u){
	vis2[u]=true;
	num++;
	belong[u]=cnt2;
	int i;
	for(i=head2[u];i!=-1;i=e2[i].next){
		if(!vis2[e2[i].to]){
			dfs2(e2[i].to);
		}
	}
}
void solve(){
	int i;
	cnt1=cnt2=0;
	memset(vis1,false,sizeof(vis1));
	memset(vis2,false,sizeof(vis2));
	for(i=1;i<=n;i++)
		if(!vis1[i])dfs1(i);
	for(i=cnt1-1;i>=0;i--){
		if(!vis2[st[i]]){
			num=0;
			dfs2(st[i]);
			setnum[cnt2++]=num;
		}
	}
}
int out[maxn];
void init(){
	//memset(mp,false,sizeof(mp));
	memset(head1,-1,sizeof(head1));
	memset(head2,-1,sizeof(head2));
	memset(head3,-1,sizeof(head1));
	memset(out,0,sizeof(out));
	tot1=tot2=tot3=0;
}
int main()
{
	int a,b;
	init();
	scanf("%d%d",&n,&m);
	int i,j,k;
	for(i=1;i<=m;i++){
		scanf("%d%d",&a,&b);
		add(a,b);
	}
	solve();
	for(i=1;i<=n;i++){
		for(j=head1[i];j!=-1;j=e1[j].next){
			if(belong[i]==belong[e1[j].to])continue;
			int u=belong[i];
			int v=belong[e1[j].to];
			for(k=head3[u];k!=-1;k=e3[k].next){
				if(e3[k].to==v)break;
			}
			add2(u,v);
			out[u]++;
		}
	}
	int tmp=-1;
	for(i=0;i<cnt2;i++){
		if(out[i]==0){
			if(tmp==-1)tmp=setnum[i];
			else tmp=-2;
		}
	}
	if(tmp>0)printf("%d\n",tmp);
	else printf("0\n");
    return 0;
}