poj 2186 Popular Cows Kosaraju求强连通分量

题意：有一群牛，每对牛之间可以有崇拜关系，同时，崇拜关系满足传递性。问被除自己以外的所有牛崇拜的牛有多少头？

第一次做强连通分量的题。这道题想想其实很简单。

思路：先按题意建图，显然，每个强联通分量中的牛都是相互崇拜的，所以我们把所有强连通分量找出来，然后把它们都看成一个点。两个强连通分量之间如果存在边，则必存在边（V1,V2）∈E，其中V1∈强连通分量1， V2∈强连通分量2。点有了，边有了，新图就出来了。然后找到这个新图中度数为0的点，如果有1个以上，则不可能存在一头牛被所有牛崇拜。如果为1，则输出度数为0的强连通分量中的点的数量，如果为0，则输出n，即所有牛都满足条件。

代码如下：

#include<stdio.h>
#include<string.h>

const int N = 10005;
const int M = 50005;

int n, m;

struct Edge {
    int v1, v2;
    int next1, next2;
};
Edge edge[M];

int fir1[N], fir2[N];
int tstamp[N], vis[N];
int scc[N];
int outdeg[N];
int cnt, p;

void init() {
    cnt = p = 0;
    memset(fir1, -1, sizeof(fir1));
    memset(fir2, -1, sizeof(fir2));
    memset(outdeg, 0, sizeof(outdeg));
}

void addEdge(int u, int v) {
    edge[cnt].v1 = v;
    edge[cnt].next1 = fir1[u];
    fir1[u] = cnt;

    edge[cnt].v2 = u;
    edge[cnt].next2 = fir2[v];
    fir2[v] = cnt++;
}

void dfs1(int u) {
    vis[u] = 1;

    for(int i = fir1[u]; i != -1; i = edge[i].next1) {
        int v = edge[i].v1;
        if(!vis[v])
            dfs1(v);
    }
    tstamp[p++] = u;
}

void dfs2(int u, int st) {
    vis[u] = 1;
    scc[u] = st;
    for(int i = fir2[u]; i != -1; i = edge[i].next2) {
        int v = edge[i].v2;
        if(!vis[v])
            dfs2(v, st);
    }
}

int kosaraju() {
    memset(vis, 0, sizeof(vis));
    for(int i = 1; i <= n; i++)
        if(!vis[i])
            dfs1(i);
    memset(vis, 0, sizeof(vis));
    int sccNum = 1;
    for(int i = n-1; i >= 0; i--) {
        if(!vis[tstamp[i]])
            dfs2(tstamp[i], sccNum++);
    }
    return sccNum;
}

void solve() {

    int sccNum = kosaraju();
    int ans = 0;
    for(int i = 1; i <= n; i++)
        for(int j = fir1[i]; j != -1; j = edge[j].next1) {
            int v = edge[j].v1;
            if(scc[i] != scc[v]) {
                outdeg[scc[i]]++;
            }
        }
    int tmp = 0;
    int mark;
    for(int i = 1; i < sccNum; i++)
        if(outdeg[i] == 0){
            tmp++;
            mark = i;
        }
    if(tmp == 0)
        printf("%d\n", n);
    else if(tmp > 1)
        printf("0\n");
    else {
        for(int i = 1; i <= n; i++)
            if(scc[i] == mark)
                ans++;
        printf("%d\n", ans);
    }
}

int main () {

    while(~scanf("%d%d", &n, &m)) {
        int u, v;
        init();
        for(int i = 0; i < m; i++) {
            scanf("%d%d", &u, &v);
            addEdge(u, v);
        }
        solve();
    }
    return 0;
}