计蒜客习题：受欢迎的蒜头

---

问题描述

蒜厂除了蒜头君还有很多小蒜头。
每只蒜头的梦想是成为最受欢迎的蒜头。有 N 只蒜头，有 M 对二元关系 (A,B)，告诉你蒜头 A 认为蒜头 B 是受欢迎的。如果 A 认为 B 是受欢迎的，B 认为 C 是受欢迎的，则 A 也认为 C 是受欢迎的。你的任务是计算被其余蒜头都认为是受欢迎的蒜头数量。
输入格式
第一行两个正整数 N 和 M，分别表示一共有 N 只小蒜头和 M 对二元关系(1≤N≤10000,1≤M≤50000)。
接下来 M 行，每行两个整数 A 和 B，表示 A 认为 B 是受欢迎的 (1≤A,B≤N)。
输出格式
一行一个整数，表示答案。
样例输入
4 6
1 2
1 3
1 4
2 3
2 4
3 4
样例输出
1

---

AC代码

#include <cstdio>  
#include <iostream>  
#include <cstring>  
#include <queue>  
#include <stack>  
#include <vector>  


using namespace std;    
const int maxn=2e5+10;    
vector<int>v[maxn];//邻接表,存原图    
vector<int>scc[maxn];//缩点后的图    
stack<int>s;    
int dfn[maxn],low[maxn],tot,ins[maxn];    
//dfn是否为0可以判断点是否访问过，ins数组用来判断点是否在栈中    
//dfn数组表示顶点dfs的时间戳，low[]为u能够追溯到的最早的栈中顶点的次序号    
int scc_cnt;//强联通分量的个数    
int sccnum[maxn];//缩点数组，表示某个点对应的缩点值,即缩完点之后的下标    
int in[maxn],out[maxn];//出度入度    

int n,m;    
void readin()//因为初始化问题，莫名其妙的wa，真的醉    
{    
    int x,y;    
    tot=0;    
    scc_cnt=0;    
    fill(ins,ins+maxn,0);    
    fill(sccnum,sccnum+maxn,0); fill(dfn,dfn+maxn,0);fill(low,low+maxn,0);   
    fill(in,in+maxn,0);    
    fill(out,out+maxn,0);    
    scanf("%d%d",&n,&m);    
    for(int i=1;i<=n;i++) v[i].clear(),scc[i].clear();    while(!s.empty()) s.pop();    

    for(int i=1;i<=m;i++){    
        scanf("%d%d",&x,&y);    
        v[x].push_back(y);    
    }    
}    
void tarjan(int x)    
{    
    low[x]=dfn[x]=++tot;//    
    s.push(x);ins[x]=1;    
    for(int i=0;i<v[x].size();i++){    
        int p=v[x][i];    
        if(!dfn[p]){//判断点是否被访问过    
            tarjan(p);    
            low[x]=min(low[x],low[p]);    
        }    
        else if(ins[p]) low[x]=min(low[x],dfn[p]);//判断点是否在栈中    
    }    

    if(low[x]==dfn[x]){//dfs到叶子节点，开始判断，如果==就到达了某个强联通分量的根节点,dfs回溯到了这个点    
        scc_cnt++;    
        scc[scc_cnt].clear();    
        while(1){//得出某个强联通分量的集合    
            int now=s.top();    
            s.pop();    
            ins[x]=0;//出栈ins置0    
            if(sccnum[now]!=scc_cnt) scc[scc_cnt].push_back(x);    
            sccnum[now]=scc_cnt;//将该连通分量中的每个点都赋予相同的值    
            if(now==x) break;//通过栈找到那个点    
        }    
    }    
}    

int main()    
{    
        readin();    
        for(int i=1;i<=n;i++){    
            if(!dfn[i]){    
                tarjan(i);    
            }    
        }    

        //出度入度的统计方法    
        for(int u=1;u<=n;u++){    
            for(int i=0;i<v[u].size();i++){    
                int uu=v[u][i];    
                if(sccnum[u]!=sccnum[uu]){    
                    in[sccnum[uu]]++;    
                    out[sccnum[u]]++;    
                }    
            }    
        }    


        int ansp=-1;    
        int flag=-1;    
        for(int i=1;i<=scc_cnt;i++){    
            if(out[i]==0){    
                flag++;    
                ansp=i;    
            }    
        }    
        if(flag==0) cout<<scc[ansp].size()<<endl;    
        else cout<<"0"<<endl;    
}