P2341 [HAOI2006] 受欢迎的牛 (【模板】强连通分量 )

题目链接：https://www.luogu.org/problem/P2341

求有向图的强连通分量–Tarjan算法
代码（模板）：

void Tarjan(int u)
{
    dfn[u]=low[u]=++num;/*初始化*/
    st[++top]=u;/*将u入栈*/
    for(int i=fir[u]; i; i=nex[i])
    {
        int v=to[i];
        if(!dfn[v])/*如果结点v未被访问过*/
        {
            Tarjan(v);/*继续向下找*/
            low[u]=min(low[u],low[v]);
        }
        else if(!co[v])/*如果结点v还在栈内即v不属于任何强联通分量*/
            low[u]=min(low[u],dfn[v]);
    }
    if(low[u]==dfn[u])
    {
        /*找到一个关键结点,即找到一个强联通分量,加上下面的while循环,把一个强联通分量内所有的节点标记为相同的col*/
        co[u]=++col;
        while(st[top]!=u)
        {
            co[st[top]]=col;/*将st[top]退栈,为该强联通分量中的一个顶点*/
            --top;
        }
        --top;/*将u退栈*/
    }
}

解析Tarjan算法的博客
https://blog.csdn.net/mengxiang000000/article/details/51672725

分析：
喜欢具有传递性，在有向图上找出强联通子图，可以发现每个子图的奶牛一定是相互喜欢的。
进行缩点
整个图变成了有向无环图
同一头奶牛需要受到其它所有奶牛的喜欢，那么它不能喜欢任意除了自己这个连通子图的牛的其它牛(图中已经不存在环)，统计出度为0的点（强连通子图）的个数
注：存在两个或两个以上的出度为0的点，输出0（子图中的奶牛各自都无法喜欢）
代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int MaxN=1e4+11;
const int MaxM=5e4+10;
int to[MaxM],nex[MaxM],fir[MaxN];
int col,num,dfn[MaxN],low[MaxN],de[MaxN],si[MaxN];
int tot=0,co[MaxN],n,m;
int top,st[MaxN];
void Ins(int x,int y)
{
    to[++tot]=y;
    nex[tot]=fir[x];
    fir[x]=tot;
}
void Tarjan(int u)
{
    dfn[u]=low[u]=++num;
    st[++top]=u;
    for(int i=fir[u]; i; i=nex[i])
    {
        int v=to[i];
        if(!dfn[v])
        {
            Tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else
        {
            if(!co[v])/*结点v还在栈中*/
                low[u]=min(low[u],dfn[v]);
        }

    }
    if(low[u]==dfn[u])
    {
        co[u]=++col;
        ++si[col];/*统计强联通子图中有多少个点*/
        while(st[top]!=u)
        {
            ++si[col];
            co[st[top]]=col;
            --top;
        }
        --top;
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1,x,y; i<=m; i++)
    {
        scanf("%d%d",&x,&y);/*将图反向,变成统计入度*/
        Ins(y,x);
    }
    for(int i=1; i<=n; i++)/*进行强联通子图的缩点*/
    {
        if(!dfn[i])
            Tarjan(i);
    }
    for(int i=1; i<=n; i++)/*统计入度,只用计算不同强联通子图之间存在边*/
        for(int j=fir[i]; j; j=nex[j])
            if(co[i]!=co[to[j]])
                de[co[to[j]]]++;/*这个地方++的是强联通子图的标号,从而实现缩点的功能*/
    int ans=0,u=0;
    for(int i=1; i<=col; i++)
    {
        if(!de[i])
        {
            ans=si[i];
            u++;
        }
    }
    if(u==1)
        printf("%d\n",ans);
    else
        printf("0\n");
    return 0;
}