POJ 3180 The Cow Prom（SCC）

【題目連結】 http://poj.org/problem?id=3180

【題目大意】

　　N頭牛，M條有向繩子，能組成幾個歌舞團？要求順時針逆時針都能帶動舞團内所有牛。

【題解】

　　等價于求點數大于1的SCC數量。

#include <cstdio>
#include <algorithm>
#include <vector>
#include <cstring> 
using namespace std;
const int MAX_V=10000;
int V; //頂點數
vector<int> G[MAX_V]; //圖的鄰接表表示
vector<int> rG[MAX_V]; //反向圖
vector<int> vs; //後序周遊
bool used[MAX_V];
int cmp[MAX_V]; //所屬強連通分量的拓撲序
void add_edge(int from,int to){
    G[from].push_back(to);
    rG[to].push_back(from);
} 
void dfs(int v){
    used[v]=1;
    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]=1;
    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;
}
const int MAX_M=50000;
int N,M;
int A[MAX_M],B[MAX_M];
int cnt[MAX_M];
void solve(){
    V=N;
    for(int i=0;i<M;i++){
        add_edge(A[i]-1,B[i]-1);
    }int n=scc();
    int num=0;
    memset(cnt,0,sizeof(cnt));
    for(int i=0;i<V;i++)++cnt[cmp[i]];
    for(int i=0;i<n;i++)if(cnt[i]>1)num++;
    printf("%d\n",num);
}
int main(){
    while(~scanf("%d%d",&N,&M)){
        for(int i=0;i<M;i++)scanf("%d%d",&A[i],&B[i]);
        solve();
    }return 0;
}