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