Tarjan算法分解强连通分量
算法思路:
算法通过dfs遍历整个连通分量,并在遍历过程中给每个点打上两个记号:一个是时间戳,即首次访问到节点i的时刻,另一个是节点u的某一个祖先被访问的最早时刻。
时间戳用DFN数组存储,最早祖先用low数组来存,每次dfs遍历到一个节点u,即让这两个记号等于当前时刻,在后面回溯或者判断的过程中在来更新low,DNF是一定的,因为第一次访问时刻一定。然后遍历u的子节点,也就是跟u相连的点v,依次看子节点的时间戳有没有打上,也就是看他有没有被访问过。\(1\).没有就继续dfs点v,在后来回溯的时候,如果v的low比u的low小,也就是v的某一祖先比u还早被访问,因为u和v是直接相连的,说明u的low也该更小,便更新u的low值。\(2\).如果有,就看这个点有没有在栈里,在的话就看v被访问的时刻是不是比u的low还早,如果是,更新u的low值。
当遍历完当前节点及其子节点后,检查u的low与DFN是否相等,如果相等,说明已经找到了一个强连通区域,就是现在栈中u及后入栈的节点。挨个pop出来即可处理。
当这一个tarjan执行完了后不能立马退出,因为图本身有可能是不连通的。还要在循环里tarjan多次。
以上只是提了下算法的思路,具体的相关概念、讲解、实例和原理都在参考文章里。
下面是两个版本的tarjan代码,一个是基于邻接矩阵,一个是基于链式前向星(具体介绍参见上一篇博客)。
代码:
版本一:邻接矩阵
#include <iostream>
#include <memory.h>
#define max_n 1005
using namespace std;
int DFN[max_n];//记录dfs访问次序
int low[max_n];//记录节点最早可追溯到的祖先
int cnt = 0;//时间戳
int Stack[max_n];//模拟栈
int top = -1;//栈顶
int flag[max_n];//记录节点是否入栈
int number = 0;//强连通分量的个数
int j;//栈弹出节点
int G[max_n][max_n];//图的邻接矩阵
int n;//图的节点数目
void tarjan(int u)
{
DFN[u] = low[u] = cnt++;//访问到这个节点打上时间戳
Stack[++top] = u;//节点入栈
flag[u] = 1;//节点已在栈内
for(int i = 0;i<n;i++)
{
if(G[u][i])//如果u与i相连
{
if(!DFN[i])//如果i节点未被访问过
{
tarjan(u);//继续遍历
low[u] = min(low[u],low[i]);//回溯时更新u点的最早祖先值
}
else if(flag[i]&&low[u]<DFN[i])//如果i被访问过,并且在栈内
{
low[u] = DFN[i];//更新u点的最早祖先值
}
}
if(DFN[u]==low[u])//如果节点之后全访问完了以后,并且DFN值与最早公共祖先值相等,说明已经得到了一个强连通分量
{
number++;//强连通分量个数加一
do
{
j = Stack[top--];//弹出栈中比该点后入栈的所有点
cout << j << " ";
flag[j] = 0;//节点不在栈中
}while(j!=u);//直到把节点u也弹出
cout << endl;
}
}
}
int main()
{
memset(DFN,0,sizeof(DFN));
memset(flag,0,sizeof(flag));
memset(Stack,0,sizeof(Stack));
for(int i = 1;i<=n;i++)//遍历多遍求出所有的强连通分量
{
if(DFN[i])
{
tarjan(i);
}
}
return 0;
}
版本二:链式前向星
#include <iostream>
#include <cstring>
#include <stack>
#include <memory.h>
#define max_n 1005
using namespace std;
int n,m;//n为节点个数,m为边的数目
int idx=0;//时间戳
int Bcnt=0;//强连通分量的个数
int instack[max_n];//记录节点是否在栈内
int dfn[max_n];//dfs访问顺序标号
int low[max_n];//节点的最早公共祖先
int Belong[max_n];//存储节点属于哪一个强连通分量
stack<int> s;//dfs访问时的栈
//链式前向星结构
int head[max_n];
int cnt = 0;
struct
{
int v;
int next;
}e[max_n<<1];
void add(int u,int v)
{
++cnt;
e[cnt].v = v;
e[cnt].next = head[u];
head[u] = cnt;
}
//读入数据
void readdata()
{
int a,b;
memset(head,0,sizeof(head));
cin >> n >> m;
for(int i = 0;i<m;i++)
{
cin >> a >> b;
add(a,b);
}
}
void tarjan(int u)
{
dfn[u] = low[u] = ++idx;
s.push(u);
instack[u] = 1;
int v;
for(int i = head[u];i;i=e[i].next)//遍历u的相连节点
{
v = e[i].v;
if(!dfn[v])//如果未访问过
{
tarjan(v);
low[u] = min(low[u],low[v]);//尝试更新u的low值
}
else if(instack[v])//如果被访问过且在栈中
{
low[u] = min(low[u],dfn[v]);//尝试更新u的low值
}
}
if(low[u]==dfn[u])//找到一个联通分量
{
Bcnt++;
do
{
v = s.top();
s.pop();
instack[v] = 0;
Belong[v] = Bcnt;
}while(u!=v);
}
}
//结果处理
void solve()
{
for(int i = 1;i<=n;i++)
{
if(!dfn[i])
{
tarjan(i);
}
}
for(int i = 1;i<=n;i++)
{
cout << "d " << dfn[i] << " l " << low[i] << endl;
}
cout << Bcnt << "个强连通分量" << endl;
for(int i = 1;i<=Bcnt;i++)
{
cout << "第 " << i << " 个" << endl;
for(int j = 1;j<=n;j++)
{
if(Belong[j]==i)
{
cout << j << " ";
}
}
cout << endl;
}
}
int main()
{
readdata();
solve();
return 0;
}
参考文章:
我大概整理了一下,这些博客会很有帮助
键盘里的青春,全网最!详!细!tarjan算法讲解,https://blog.csdn.net/qq_34374664/article/details/77488976(主要是原理,实例,我之前看了很多博客云里雾里的,这篇看懂了)
玩人,Tarjan算法详解,https://blog.csdn.net/jeryjeryjery/article/details/52829142?locationNum=4(算法流程比较清楚,实现基于邻接矩阵,代码有少许错误)
five20,浅析强连通分量(Tarjan和kosaraju),https://www.cnblogs.com/five20/p/7594239.html(实例再次,代码不错,基于链式前向星)
BYVoid,有向图强连通分量的Tarjan算法,https://www.byvoid.com/zhs/blog/scc-tarjan (我刚开始看的,由于我本身还不太理解,看不太懂,但渐渐明白后它的讲解还是不错的。偶然发现这个是dalao的网站,dalao网站有繁体字版本(好像鸟哥的),大佬貌似还对中国语言有研究,有很多有意思的文章,本来看算法结果看其他的文章去了(逃)
最后推荐一个有趣的网站,画图!https://csacademy.com/app/graph_editor/ 对,就是给节点画图的那种,感觉好方便
觉得不错要不右下角推荐一下?
原文地址:https://www.cnblogs.com/zhanhonhao/p/11296012.html