首页 > 代码库 > tarjan算法求强连通分量

tarjan算法求强连通分量

2024-07-14 03:35:17   235人阅读

先上代码：

#include <iostream>#include <cstring>#include <vector>#include <stack>using namespace std;int G[1000][1000];int pre[1000],sccno[1000],lowlink[1000];　　　  　//sccno[u]:u节点所属于的强连通分量序号　　
int dfs_clock,scc_cnt,n,m,x,y;　　　　　　　　 　 　//pre[u]:u节点的时间戳stack<int> S;　　　　　　　　　　　　　　　　　　　　 　//lowlink[u]:由u节点所能追溯到的最早的点的时间戳void dfs(int u){    pre[u]=lowlink[u]=++dfs_clock;    S.push(u);    for (int i=1;i<=G[u][0];i++)    {        int v=G[u][i];        if (!pre[v])        {            dfs(v);            lowlink[u]=min(lowlink[u],lowlink[v]);        }        else if (!sccno[v])        {            lowlink[u]=min(lowlink[u],pre[v]);        }    }    if (lowlink[u]==pre[u])　　　　　　//说明找到了一个新的连通分量    {        scc_cnt++;        for (;;)        {            int x=S.top();            S.pop();            sccno[x]=scc_cnt;            if (x==u) break;        }    }}void find_scc(int n){    dfs_clock=scc_cnt=0;    memset(sccno,0,sizeof(sccno));    memset(pre,0,sizeof(pre));    for (int i=1;i<=n;i++)        if (!pre[i])    dfs(i);　　　　　　//pre[i]==0说明还没搜索过i点}int main(){    cin>>n>>m;    memset(G,0,sizeof(G));    for (int i=1;i<=m;i++)    {        cin>>x>>y;        G[x][0]++;　　　　　　　　//使用邻接表记录图        G[x][G[x][0]]=y;    }    find_scc(n);    cout<<scc_cnt<<endl;    for (int i=1;i<=n;i++)        cout<<sccno[i]<<" "<<pre[i]<<endl;　　　　    return 0;}

 

Reference：

http://blog.csdn.net/xinghongduo/article/details/6195337

http://blog.csdn.net/xinghongduo/article/details/6196292

代码中用到了STL栈容器，可参考这里：

http://blog.csdn.net/wallwind/article/details/6858634

 


联系
我们