【BZOJ 1051】[HAOI2006]受欢迎的牛

【链接】 我是链接,点我呀:)
【题意】

在这里输入题意

【题解】

Tarjan算法强连通缩点 。 最后出度为0的点。 如果只有一个。 那么这个“大点”所包含的点的个数就是答案了。

【代码】

/*
    n个点,m条有向边.
    把有向图G的环进行缩点;
    缩完之后的图存在vector <int> g[N]里面;
    n变为缩完点之后的图的节点的个数了。
*/
#include <bits/stdc++.h>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define LL long long
#define rep1(i,a,b) for (int i = a;i <= b;i++)
#define rep2(i,a,b) for (int i = a;i >= b;i--)
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define ms(x,y) memset(x,y,sizeof x)
#define ri(x) scanf("%d",&x)
#define rl(x) scanf("%lld",&x)
#define rs(x) scanf("%s",x)
#define oi(x) printf("%d",x)
#define ol(x) printf("%lld",x)
#define oc putchar(' ')
#define os(x) printf(x)
#define all(x) x.begin(),x.end()
#define Open() freopen("F:\\rush.txt","r",stdin)
#define Close() ios::sync_with_stdio(0)


typedef pair<int,int> pii;
typedef pair<LL,LL> pll;


const int dx[9] = {0,1,-1,0,0,-1,-1,1,1};
const int dy[9] = {0,0,0,-1,1,-1,1,-1,1};
const double pi = acos(-1.0);
const int N = 1e4;//节点个数


vector <int> G[N+10],g[N+10];
int n,m,tot = 0,top = 0,dfn[N+10],low[N+10],z[N+10],totn,in[N+10];
int bh[N+10];


void dfs(int x){
    dfn[x] = low[x] = ++ tot;
    z[++top] = x;
    in[x] = 1;
    int len = G[x].size();
    rep1(i,0,len-1){
        int y = G[x][i];
        if (!dfn[y]){
            dfs(y);
            low[x] = min(low[x],low[y]);
        }else
        if (in[y] && dfn[y]<low[x]){
            low[x] = dfn[y];
        }
    }
    if (low[x]==dfn[x]){
        int v = 0;
        totn++;
        while (v!=x){
            v = z[top];
            in[v] = 0;
            bh[v] = totn;
            top--;
        }
    }
}

bool bo[N+10];

int main(){
	#ifdef LOCAL_DEFINE
	    freopen("rush_in.txt", "r", stdin);
	#endif
    ms(dfn,0);
    ms(in,0);
    tot = 0,totn = 0;
    ri(n),ri(m);
    rep1(i,1,n) G[i].clear(),g[i].clear();
    rep1(i,1,m){
        int x,y;
        ri(x),ri(y);
        G[x].pb(y);
    }


    rep1(i,1,n)
        if (dfn[i]==0)
            dfs(i);


    rep1(i,1,n){
        int len = G[i].size();
        int xx = bh[i];
        rep1(j,0,len-1){
            int y = G[i][j];
            int yy = bh[y];
            if (xx!=yy)
                g[xx].pb(yy);
        }
    }

    int cnt = 0,idx = -1;
    for (int i = 1;i <= totn;i++)
        if (g[i].empty()){
            cnt++;
            idx = i;
        }
    if (cnt!=1){
        oi(0);
        puts("");
    }else{
        int ans = 0;
        for (int i=1;i <= n;i++)
        if (bh[i]==idx){
            ans++;
        }
        oi(ans);
        puts("");
    }
    return 0;
}