poj 2186 Popular Cows (强联通分量+缩点)

给出一些牛自己认为popular的牛，问那些牛是被所有牛认为popular的。

裸的强联通分量+缩点。

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
typedef long long ll;
void cmax(int& a,int b){ if(b>a)a=b; }
void cmin(int& a,int b){ if(b<a)a=b; }
const int oo=0x3f3f3f3f;
const int MOD=1000000007;
const int maxn=10050;
const int maxm=50050;
struct EDGE{
    int v,next;
}E[maxm<<1];
int head[maxn],tol;
int dfn[maxn],low[maxn],instack[maxn],Stack[maxn];
int id[maxn],d[maxn];
int g_cnt,ID,top;
int U[maxm],V[maxm];

void Init(){

    memset(dfn,0,sizeof dfn);
    memset(instack,0,sizeof instack);
    memset(low,0,sizeof low);
    memset(d,0,sizeof d);
    memset(id,0,sizeof id);
    memset(head,-1,sizeof head);
    g_cnt=ID=top=tol=0;
}

void add_edge(int u,int v){
    E[tol].v=v;
    E[tol].next=head[u];
    head[u]=tol++;
}

void Tarjan(int u){
    dfn[u]=low[u]=++g_cnt;
    instack[u]=1;
    Stack[++top]=u;
    int v;
    for(int i=head[u];i!=-1;i=E[i].next){
        v=E[i].v;
        if(!dfn[v]){
            Tarjan(v);
            if(low[v]<low[u])
                low[u]=low[v];
        }else if(instack[v]&&dfn[v]<low[u])
            low[u]=dfn[v];
    }
    if(dfn[u]==low[u]){
        ID++;
        do{
            v=Stack[top--];
            instack[v]=0;
            id[v]=ID;
        }while(u!=v);
    }
}

void gao(int n,int m){
    for(int i=1;i<=m;i++){
        int u=id[U[i]];
        int v=id[V[i]];
        if(u!=v)
            d[u]++;
    }
    int cnt=0,pos=1;
    for(int i=1;i<=ID;i++){
        if(d[i]==0){
            cnt++;
            pos=i;
        }
    }
    if(cnt>1){
        printf("0\n");
        return ;
    }
    cnt=0;
    for(int i=1;i<=n;i++){
        if(id[i]==pos)
            cnt++;
    }
    printf("%d\n",cnt);
}

int main(){
    int n,m,u,v;
    //freopen("G:\\read.txt","r",stdin);
    while(scanf("%d %d",&n,&m)!=EOF){
        Init();
        for(int i=1;i<=m;i++){
            scanf("%d %d",&u,&v);
            U[i]=u;V[i]=v;
            add_edge(u,v);
        }
        for(int i=1;i<=n;i++)if(!dfn[i])Tarjan(i);
        gao(n,m);
    }
    return 0;
}