POJ2186 Popular Cows（SCC缩点）

题目链接

对于一个有向图，牛(A,B)视为牛A崇拜牛B，具有传递性。定义最受欢迎的牛为其他牛都有路径可以到达它。现在求M组关系的情况下，有多少牛是最受欢迎的牛。

问题分析：
考虑到两头牛都是红牛，那么A,B必然是相互崇拜的。同理若K头牛相互崇拜，那么这K头牛都能相互到达，即K头牛所处同一个联通分量。现在对给定的有向图进行联通分量的分解，缩环成点，构造一个新的图（DAG），统计每个节点的入度出度。可知出度为0的点（可以想象是拓扑排序最后面的点）只有一个时才能满足有解。且解是该节点所包含的节点数。

#include<cstring>
#include<cstdio>
#include<queue>
#define Re register int
const int N=1e4+3,M=1e5+3;
int n,m,x,y,Q_o,ans,A[N],ru[N],gs[N],ip[N],dp[N],cu[N];
inline void in(Re &x){
    int f=0;x=0;char c=getchar();
    while(c<'0'||c>'9')f|=c=='-',c=getchar();
    while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=getchar();
    x=f?-x:x;
}
inline int min(Re a,Re b){return a<b?a:b;}
inline int max(Re a,Re b){return a>b?a:b;}
struct Tarjan{
    int o,t,dfn_o,Q[N],low[N],dfn[N],pan[N],head[N];
    struct QAQ{int x,to,next;}a[M];
    inline void add(Re x,Re y){a[++o].x=x,a[o].to=y,a[o].next=head[x],head[x]=o;}
    inline void tarjan(Re x){
        dfn[x]=low[x]=++dfn_o,Q[++t]=x,pan[x]=1;
        for(Re i=head[x],to;i;i=a[i].next)
            if(!dfn[to=a[i].to])tarjan(to),low[x]=min(low[x],low[to]);
            else if(pan[to])low[x]=min(low[x],dfn[to]);
        if(low[x]==dfn[x]){
            ++Q_o;
            while(1){
                ip[Q[t]]=Q_o,gs[Q_o]+=A[Q[t]],pan[Q[t]]=0;
                if(x==Q[t--])break;
            }
        }
    }
    inline void SuoPoint(){for(Re i=1;i<=n;++i)if(!dfn[i])tarjan(i);}
}T1;
struct Tuopu{
    int o,ru[N],pan[N],head[N];std::queue<int>Q;
    struct QAQ{int to,next;}a[M<<1];
    inline void add(Re x,Re y){a[++o].to=y,a[o].next=head[x],head[x]=o;}
    inline void creat(){
        for(Re i=1;i<=T1.o;++i)
            if((x=ip[T1.a[i].x])!=(y=ip[T1.a[i].to]))
                ++ru[y],add(x,y),++cu[x];
    }
    inline void tuopu(){
        for(Re i=1;i<=Q_o;++i)if(!ru[i])Q.push(i);
        while(!Q.empty()){
            Re x=Q.front();Q.pop();
            for(Re i=head[x],to;i;i=a[i].next){
                to=a[i].to;
                dp[to]=max(dp[to],dp[x]+gs[to]);
                if(!(--ru[to]))Q.push(to);
            }
        }
    }
}T2;
int main(){
    in(n),in(m);
    for(Re i=1;i<=n;++i)A[i]=1;
    while(m--)in(x),in(y),T1.add(x,y);
    T1.SuoPoint();
    T2.creat();
    int zero=0;
    for(Re i=1;i<=Q_o;i++){
    if(!cu[i])zero++,ans=i;
    }
     if(zero!=1){
        printf("0\n");

     }
     else printf("%d\n",gs[ans]);

}