POJ 2186(Popular Cows)

题意：有n个奶牛，它们之间互相膜拜，并且如果A膜拜B，B膜拜C则，A也膜拜C，求有多少头奶牛被所有奶牛都膜拜；
思路：强连通分量+缩点，刚学SCC，这题是书上的例题，顺手敲下熟悉下，用的是Kosaraju算法，这个算法比较明了；
对于这个题目，先求出强连通分量的个数，如果强连通分量个数只有一个，就输出n，否则对其进行缩点，然后统计各点
缩完点之后的度数，若缩完点后度数为0的点有1个，则答案就是这个强连通分量的奶牛个数，如果大于1个则说明至少有
两群奶牛不相互膜拜，答案为0。
第一题强连通分量题目，贴个代码

#include <set>
#include <map>
#include <queue>
#include <vector>
#include <cstdio>
#include <cstring>
#include <iostream>
#define pub push_back
using namespace std;

const int N = 10010;
//vis:访问数组 st:访问时间数组 scc:各强连通分量编号数组 oe:出入度素组 sum:各强连通分量点的个数 
int vis[N], st[N], scc[N], oe[N], sum[N];
vector<int> v1[N], v2[N];//v1为原图,v2为转置图,也叫逆图 

void dfsOne(int to, int &sig){
	vis[to] = 1;
	for(int i = 0;i < v1[to].size();i++)
		if(!vis[v1[to][i]])
			dfsOne(v1[to][i], sig);
	st[sig++] = to;
}

void dfsTwo(int to, int sig){
	vis[to] = 1;
	scc[to] = sig;
	for(int i = 0;i < v2[to].size();i++)
		if(!vis[v2[to][i]])
			dfsTwo(v2[to][i], sig);
}

int Kosaraju(int n){
	int sig = 1, i;
	memset(scc, 0, sizeof(scc));
	memset(st, 0, sizeof(st));
	memset(vis, 0, sizeof(vis));
	for(i = 1;i <= n;i++)
		if(!vis[i])
			dfsOne(i, sig);
	memset(vis, 0, sizeof(vis));
	sig = 0;
	for(i = n;i >= 1;i--)
		if(!vis[st[i]])
			dfsTwo(st[i], sig++);
	return sig;
}

void solve(int n, int sig){
	memset(oe, 0, sizeof(oe));//出度 
	memset(sum, 0, sizeof(sum));//各连通分量点的数量
	//缩点 
	for(int i = 1;i <= n;i++){
		sum[scc[i] + 1]++;
		for(int j = 0;j < v1[i].size();j++)
			if(scc[v1[i][j]] != scc[i])
				oe[scc[i] + 1]++;
	} 
	int res = 0, cnt = 0;
	for(int i = 1;i <= sig;i++)
		if(oe[i] == 0){
			cnt++;
			res += sum[i];
		}
	printf("%d\n", cnt > 1?0:res);
}

int main(){
	int n, m;
	while(cin >> n >> m){
		int i, a, b;
		for(i = 0;i < m;i++){
			scanf("%d%d", &a, &b);
			v1[a].pub(b);
			v2[b].pub(a);
		}
		int tot = Kosaraju(n);
		if(tot == 1)
			printf("%d\n", n);
		else
			solve(n, tot);
		for(i = 1;i <= n;i++){
			v1[i].clear();
			v2[i].clear();
		}
	}
	return 0;
}

这个是用tarjan实现的，居然比kosaraju要慢。。估计写挫了。。

#include <set>
#include <map>
#include <queue>
#include <vector>
#include <cstdio>
#include <cstring>
#include <iostream>
#define pub push_back
using namespace std;

const int N = 10010;
int oe[N], sum[N];

/********************Tarjan********************/
int dfn[N], low[N], idx, a, scc[N];
int vis[N], stack[N], tot;
vector<int> v[N];

int tarjanDFS(int to){
	vis[to] = 1;
	stack[++idx] = to;
	low[to] = dfn[to] = ++a;
	for(int i = 0;i < v[to].size();i++){
		if(!vis[v[to][i]]){
			tarjanDFS(v[to][i]);
			low[to] = min(low[to], low[v[to][i]]);
		}
		if(vis[v[to][i]] == 1)
			low[to] = min(dfn[v[to][i]], low[to]);
	}
 	if(dfn[to] == low[to]){
 		tot++;
 		int e;
	 	do{
		 	e = stack[idx--];
	 		scc[e] = tot;
			vis[e] = 2;
	 	}while(e != to);
	 }
}

int tarjan(int n){
	idx = 0, tot = 0, a = 0;
	memset(vis, 0, sizeof(vis));
	memset(dfn, 0, sizeof(dfn));
	for(int i = 1;i <= n;i++)
		if(!vis[i])
			tarjanDFS(i);
	return tot;
}
/*********************END**********************/

void solve(int n){
	tarjan(n);
	if(tot == 1){
		printf("%d\n", n);
		return ;
	}
	memset(oe, 0, sizeof(oe));//出度 
	memset(sum, 0, sizeof(sum));//各连通分量点的数量
	//缩点
	for(int i = 1;i <= n;i++){
		sum[scc[i]]++;
		for(int j = 0;j < v[i].size();j++)
			if(scc[v[i][j]] != scc[i])
				oe[scc[i]]++;
	} 
	int res = 0, cnt = 0;
	for(int i = 1;i <= tot;i++)
		if(oe[i] == 0){
			cnt++;
			res += sum[i];
		}
	printf("%d\n", cnt > 1?0:res);
}

int main(){
	int n, m;
	while(cin >> n >> m){
		int i, a, b;
		for(i = 0;i < m;i++){
			scanf("%d%d", &a, &b);
			v[a].pub(b);
		}
		solve(n);
		for(int i = 1;i <= n;i++)
			v[i].clear();
	}
	return 0;
}