POJ 2186 Popular Cows 强连通分量 Kosaraju算法

题意：有N头牛,他们之间有M个关系(A,B)，即A认为B是欢迎的，同时这种欢迎关系是传递的。现在要求出有多少牛被其他的所有牛认为是欢迎的。

思路：首先这是一个有向图。可以想到，对于一个强连通分量里的牛，每个牛都是被其他牛认为是欢迎的。这样，通过强连通缩点后，就变成了一个DAG。这样，能被其他牛认为是欢迎的，必然是拓扑排序的最后一个点。这样的话，我们就要用强连通分量的Kosaraju算法算法了。

注意：需要注意的一点是，最后我们还要判断其他所有的点能否到达最后一个点，否则就没有满足条件的牛。

代码如下：

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>

using namespace std;

template<class T>
inline bool read(T &n){
	T x = 0, tmp = 1; char c = getchar();
	while ((c < '0' || c > '9') && c != '-' && c != EOF) c = getchar();
	if (c == EOF) return false;
	if (c == '-') c = getchar(), tmp = -1;
	while (c >= '0' && c <= '9') x *= 10, x += (c - '0'), c = getchar();
	n = x*tmp;
	return true;
}

template <class T>
inline void write(T n) {
	if (n < 0) {
		putchar('-');
		n = -n;
	}
	int len = 0, data[20];
	while (n) {
		data[len++] = n % 10;
		n /= 10;
	}
	if (!len) data[len++] = 0;
	while (len--) putchar(data[len] + 48);
}

#define MAX_N 10001
#define MAX_M 50001

int v;
vector <int> g[MAX_N];
vector <int> rg[MAX_N];
vector <int> vs;
bool used[MAX_M];
int cmp[MAX_M];
int n,m;
int a[MAX_M];
int b[MAX_M];
void addedge(int from,int to)
{
    g[from].push_back(to);
    rg[to].push_back(from);
}
void dfs(int v)
{
    used[v] = true;
    for(int i = 0 ; i < g[v].size(); ++i)
        if(!used[g[v][i]])
            dfs(g[v][i]);
    vs.push_back(v);
}
void rdfs(int v, int k)
{
    used[v] = true;
    cmp[v] = k;
    for(int i = 0; i < rg[v].size(); ++i)
        if(!used[rg[v][i]])
            rdfs(rg[v][i],k);
}
int scc()
{
    memset(used,0,sizeof(used));
    vs.clear();
    for(int i = 0 ; i < v ; ++i)
        if(!used[i])
            dfs(i);
    memset(used,0,sizeof(used));
    int  k = 0;
    for(int i  = vs.size() -1; i >= 0; i--)
        if(!used[vs[i]])
            rdfs(vs[i],k++);
    return k;
}
void solve()
{
    v= n;
    for(int i = 0 ; i <m; ++i)
        addedge(a[i]-1,b[i]-1);
    int nn = scc();
    int u = 0 , num= 0;
    for(int i = 0 ; i < v ; ++i)
        if(cmp[i] == nn-1)
    {
        u = i;
        num++;
    }
    memset(used,0,sizeof(used));
    rdfs(u,0);
    for(int  i = 0 ; i < v; ++i)
        if(!used[i])
    {
        num = 0;
        break;
    }
    printf("%d\n",num);
}
int main(void)
{
    read(n),read(m);
    for(int i = 0; i <m ; ++i)
        read(a[i]),read(b[i]);
    solve();
    return 0;
}