Tarjan算法求有向图强连通分量1.2

详细解释参考:https://blog.csdn.net/lanshan1111/article/details/83279153

一个模板:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<stack>
using namespace std;
const int maxn=1000+10;
int n,m; //点数,边数
int dfs_clock;//时钟
int scc_cnt;//强连通分量总数
vector<int> G[maxn];//G[i]表示i节点指向的所有点,vector很方便;
int pre[maxn]; //时间戳
int low[maxn]; //u以及u的子孙能到达的祖先pre值
int sccno[maxn];//sccno[i]==j表示i节点属于j连通分量
stack<int> S;
 
void dfs(int u)
{
    pre[u]=low[u]=++dfs_clock;//初始化;
    S.push(u);
    for(int i=0;i<G[u].size();i++)
    {
        int v=G[u][i];
        if(!pre[v])
        {
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if(!sccno[v])
        {
            low[u]=min(low[u],pre[v]);
        }
    }
    if(low[u] == pre[u])//u为当前强连通分量的入口
    {
        scc_cnt++;
        while(true)
        {
            int x=S.top(); S.pop();
            sccno[x]=scc_cnt;
            if(x==u) break;
        }
    }
}
 
//求出有向图所有连通分量
void find_scc(int n)
{
    scc_cnt=dfs_clock=0;//时钟强连通分量数初始化为0;
    memset(sccno,0,sizeof(sccno));
    memset(pre,0,sizeof(pre));
    for(int i=0;i<n;i++)
        if(!pre[i]) dfs(i);
}
int main()
{
    while(scanf("%d%d",&n,&m)==2&&n)
    {
        for(int i=0;i<n;i++) G[i].clear();
        while(m--)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            G[u].push_back(v);//记录u可到达的点
        }
        find_scc(n);//n个点;
        for(int i=0;i<n;i++)
            printf("%d号点属于%d分量\n",i,sccno[i]);
    }
}
/*
刘汝佳 训练指南P319测试图
输入:
12 17
0 1
1 2
1 3
1 4
4 1
2 5
5 2
4 5
4 6
5 7
6 7
8 6
6 9
9 8
7 10
10 11
11 9
输出:
0号点属于5分量
1号点属于4分量
2号点属于2分量
3号点属于3分量
4号点属于4分量
5号点属于2分量
6号点属于1分量
7号点属于1分量
8号点属于1分量
9号点属于1分量
10号点属于1分量
11号点属于1分量
*/

猜你喜欢

转载自blog.csdn.net/lanshan1111/article/details/83315536
今日推荐