这是一个漫(jian)长(nan)的过程
请大家耐心读完,相信你一定能学会
首先来介绍一下强连通分量
神奇海螺指引你:
有向图强连通分量:在有向图G中,如果两个顶点vi,vj间(vi>vj)有一条从vi到vj的有向路径,同时还有一条从vj到vi的有向路径,则称两个顶点强连通(strongly connected)。如果有向图G的每两个顶点都强连通,称G是一个强连通图。有向图的极大强连通子图,称为强连通分量(strongly connected components)。——百度百科
Tarjan是基于迪法师(DFS)的一种算法,可以说是一种流(du)批(liu)的操作,本蒟蒻苦哈哈的学了好几天才学会强连通分量和缩点;
在这里推荐给大家几道题:受欢迎的牛,消息扩散;
这些题我会在后期的博客中进行详细的讲解
步入正轨:
图中{6}是第一个被发现的强连通分量,其次是{5},最后被发现的是{3,4,1,2}(顺序无所谓)
Tarjan算法是解决强连通分量和缩点问题的一种比较常见的方法
接下来开始🔛一步一步的
首先给你一个图(求其中的强连通分量)(注:以下图片摘自:~~)
在这里我们假设从点1开始遍历,并且将访问的点压入栈中
在进行操作之前,我们要明确Tarjan中的两个至关重要的数组:DFN和LOW;
DFN[x]:表示这个点x几次被遍历到;
LOW[x]:表示点x或点x的子树能够追溯到的最早的栈中节点的次序号;
即LOW[x]=min{LOW[x],LOW[next]};
其中next为x的子节点;
接下来,我们来一步一步的找;
第一遍从1直接一直遍历到6,另访问到的点的初始low和dfn都为num++;num为计数器;
但当我们遍历到6是发现6没有子节点了,所以low[6]=dfn[6]=4;发现一个强连通分量{6},6出栈;
退回到5,low[5]=min(low[5],low[6]);
low[6]=4;而low[5]=dfn[5]=3;
所以low[5]=3;
low[5]=dfn[5];所以5也是一个强连通分量;5出栈;
返回3发现3有子节点4,搜索4,并且4入栈;发现4也有子节点,6已经为强连通分量且不再栈中,跳过;
遍历1,发现一在栈中且已经被访问过,那么4值得被发现;
所以low[4]=min(low[4],dfn[1]);所以low[4]=1;
4遍历完,又因为dfn[4]!=low[4],所以4不退站;回到3,low[3]=min(low[3],low[4]);
所以low[3]=1;3访问完,又因为dfn[3]!=low[3],所以3不退站;
回到1,发现1还有子节点2,2入栈;发现2的子节点4已经被访问且在栈中,所以low[2]=min(dfn[4],low[2]);
所以low[2]=5;访问完2,回到一,发现low[1]=dfn[1];所以将栈中的所有点出栈,并且这些点为一个强连通分量;
至此,所有点已经都被访问,算法结束。
时间复杂度为O(N+M);
下面插上代码:
#include<iostream>
#include<cstdio>
#include<vector>
#include<string.h>
#include<stack>
#define Size 101
using namespace std;
inline int read()
{
int x=0,f=1;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')
f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=(x<<1)+(x<<3)+(ch^48);
ch=getchar();
}
return x*f;
}
stack<int> atack;
vector<int> map[Size];
int a,b,c,d,num=0,low[Size],dfn[Size];
bool inst[Size];
void tarjan(int now)
{
dfn[now]=low[now]=++num;
atack.push(now);
inst[now]=true;
for(int i=0;i<map[now].size();i++)
{
int next=map[now][i];
if(dfn[next]==0)
{
tarjan(next);
low[now]=min(low[now],low[next]);
}else if(inst[next]==1)
{
low[now]=min(low[now],dfn[next]);
}
}
if(dfn[now]==low[now])
{
int p;
do
{
p=atack.top();
atack.pop();
cout<<p<<" ";
inst[p]=0;
}while(p!=now);
cout<<endl;
}
}
int main()
{
a=read();b=read();
for(int i=1;i<=b;i++)
{
c=read();d=read();
map[c].push_back(d);
}
memset(dfn,0,sizeof(dfn));
memset(inst,false,sizeof(inst));
for(int i=1;i<=a;++i)
{
if(dfn[i]==0)
tarjan(i);
}
return 0;
}
写一篇文章不容易,可否来个赞👍。
来源:CSDN
作者:123-KB
链接:https://blog.csdn.net/dfsdsdasa/article/details/104510023