tarjan算法求强连通分量
先上代码:
#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
posted on 2014-07-26 11:23 Pentium.Labs 阅读(197) 评论(0) 编辑 收藏 举报