强连通分量做题笔记

有向图的强联通分量

tarjan算法基础板子

tarjan缩点后就是一个拓扑序，没必要再拓扑排序

//链式前向星存图
int h[N],e[M],ne[M],idx;

/*
timestamp时间戳变量
dfn dfs搜索到当前点的时间戳
low 当前点能到达点的最小时间戳
if dfn==low 说明当前点不能回到前方点

*/
int dfn[N],low[N],timestamp;

//模拟栈
int stk[N],top;
bool in_stk[N];//是否入栈

//id当前点缩点后的点指向（即缩点后的点）
//siz强连通分量的大小
int id[N],scc_cnt,siz[N];



void add(int a, int b)
{
    idx++;
    e[idx]=b;
    ne[idx]=h[a];
    h[a]=idx;
}

void tarjan(int x)
{
    dfn[x]=low[x]=++timestamp;
    stk[++top]=x;in_stk[x]=true;
    for(int i=h[x];i;i=ne[i])
    {
        int j=e[i];
        if(!dfn[j])
        {
            tarjan(j);
            low[x]=min(low[j],low[x]);
        }
        else if(in_stk[j])low[x]=min(low[x],dfn[j]);
    }
    if(dfn[x]==low[x])
    {
        ++scc_cnt;
        int y;
        do
        {
            y=stk[top--];
            in_stk[y]=false;
            id[y]=scc_cnt;
            siz[scc_cnt]++;
        }while(y!=x);
    }
}

题目练习

acwing1174受欢迎的牛

思路：tarjan缩点后若只有一个点的出度为0，则该点的强连通分量大小即为所求，否则为0.

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=10010;
const int M=50010;

int n,m;
int h[N],e[M],ne[M],idx;
int dfn[N],low[N],timestamp;
int stk[N],top;
bool in_stk[N];
int id[N],scc_cnt,siz[N];
int dout[N];

void add(int a, int b)
{
    idx++;
    e[idx]=b;
    ne[idx]=h[a];
    h[a]=idx;
}

void tarjan(int x)
{
    dfn[x]=low[x]=++timestamp;
    stk[++top]=x;in_stk[x]=true;
    for(int i=h[x];i;i=ne[i])
    {
        int j=e[i];
        if(!dfn[j])
        {
            tarjan(j);
            low[x]=min(low[j],low[x]);
        }
        else if(in_stk[j])low[x]=min(low[x],dfn[j]);
    }
    if(dfn[x]==low[x])
    {
        ++scc_cnt;
        int y;
        do
        {
            y=stk[top--];
            in_stk[y]=false;
            id[y]=scc_cnt;
            siz[scc_cnt]++;
        }while(y!=x);
    }
}


int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);cout.tie(0);

    cin>>n>>m;
    for(int i=1;i<=m;i++)
    {
        int a,b;cin>>a>>b;
        add(a,b);
    }
    for(int i=1;i<=n;i++)
    {
        if(!dfn[i])tarjan(i);
    }

    for(int i=1;i<=n;i++)
    {
        for(int j=h[i];j;j=ne[j])
        {
            int a=id[i];
            int b=id[e[j]];
            if(a!=b)dout[a]++;
        }
    }

    int res=0,ans=0;
    for(int i=1;i<=scc_cnt;i++)
    {
        if(!dout[i])
        {
            res++;
            ans+=siz[i];
            if(res>1)
            {
                ans=0;
                break;
            }
        }
    }
    cout<<ans<<"\n";

    return 0;
}