POJ 2186 Popular Cows强连通分量分解

题目链接http://poj.org/problem?id=2186

 假设两头牛A和B都被其他牛认为是红人，则A认为B是，B认为A是，及存在一个包含A和B的圈，或者A和B属于同一个强连通分量。反之，如果一个牛被其他所有牛认为是红人，那么其所在的强连通分量的所有牛都被其他所有牛认为是红人，显然这样的强连通分量最多只有一个，否则强连通分量还能继续扩大，不满足强连通分量的定义。我们把图进行强连通分量分解后，可以得到各个强连通分量的拓扑排序后的顺序，显然能够成为解的只有可能是拓扑序最后的强连通分量，只需判断这个强连通分量能否从其他任意点到达即可。时间复杂度为O(N+M)。

AC代码

#include <iostream>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn =10000+100;

int V;
vector<int> G[maxn];
vector<int> rG[maxn];//反向后的图的邻接表
vector<int> vs;//后序遍历顺序的顶点列表
bool used[maxn];//访问标记
int cmp[maxn];//所属的强连通分量的拓扑序

void add_edge(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 v=0;v<V;v++){
        if(!used[v]){
            dfs(v);
        }
    }

    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;
}

int N,M;
int A[5*maxn],B[5*maxn];

void solve(){
    V=N;
    for(int i=0;i<M;i++){
        add_edge(A[i]-1,B[i]-1);
    }
    int n=scc();

    int u=0,num=0;
    for(int v=0;v<V;v++){
        if(cmp[v]==n-1){
            u=v;
            num++;
        }
    }

    memset(used,0,sizeof(used));
    rdfs(u,0);
    for(int v=0;v<V;v++){
        if(!used[v]){
            num=0;
            break;
        }
    }
    printf("%d\n",num);
}

int main(){
    while(scanf("%d%d",&N,&M)==2){
        for(int i=0;i<M;i++){
            scanf("%d%d",&A[i],&B[i]);
        }
        solve();
    }
    return 0;
}