码迷,mamicode.com
首页 > 编程语言 > 详细

Tarjan算法分解强连通分量(附详细参考文章)

时间:2019-08-03 20:02:15      阅读:102      评论:0      收藏:0      [点我收藏+]

标签:https   ext   ret   概念   tarjan   max   htm   浅析   有趣   

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/ 对,就是给节点画图的那种,感觉好方便

觉得不错要不右下角推荐一下?

Tarjan算法分解强连通分量(附详细参考文章)

标签:https   ext   ret   概念   tarjan   max   htm   浅析   有趣   

原文地址:https://www.cnblogs.com/zhanhonhao/p/11296012.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!