有向图的强联通分量
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;
}