先上代码:
#include <iostream>
#include <cstring>
#include <vector>
#include <stack>
using namespace std;
int G[][];
int pre[],sccno[],lowlink[]; //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=;i<=G[u][];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=;
memset(sccno,,sizeof(sccno));
memset(pre,,sizeof(pre));
for (int i=;i<=n;i++)
if (!pre[i]) dfs(i); //pre[i]==0说明还没搜索过i点
} int main()
{
cin>>n>>m;
memset(G,,sizeof(G)); for (int i=;i<=m;i++)
{
cin>>x>>y;
G[x][]++; //使用邻接表记录图
G[x][G[x][]]=y;
} find_scc(n);
cout<<scc_cnt<<endl;
for (int i=;i<=n;i++)
cout<<sccno[i]<<" "<<pre[i]<<endl; return ;
}
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