强连通分量Tarjan模板

 #include<iostream>
 #include<stdio.h>
 #include<string.h>
 #include<stack>
 #include<vector>
 #include<algorithm>

 #define N 1000

 using namespace std;

 int pre[N],lowlink[N],sccno[N],dfn_clock,scc_cnt;

 stack<int> stk;

 vector<int> G[N];

 void DFN(int u)
 {
     pre[u] = lowlink[u] = ++dfn_clock;  //时间戳
     stk.push(u);
     int i;
     for(i=;i<G[u].size();i++)
     {
         int v = G[u][i];
         if(!pre[v])  //没有访问过
         {
             DFN(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++;  //强连通分量的个数计算
         while()
         {
             int x = stk.top();
             stk.pop();
             sccno[x] = scc_cnt;
             if(x==u)    break;
         }
     }
 }

 void find_scc(int n)
 {
     dfn_clock = scc_cnt = ;
     memset(sccno,,sizeof(sccno));
     memset(pre,,sizeof(pre));
     int i;
     for(i=;i<=n;i++)
         if(!pre[i]) DFN(i);
 }

 int main()
 {
     int n,m;
     scanf("%d%d",&n,&m);
     int i;
     for(i=;i<m;i++)
     {
         int a,b;
         scanf("%d%d",&a,&b);
         G[a].push_back(b);   //构造图
     }
     find_scc(n);  //查找强连通图
     for(i=;i<=n;i++)
         printf("%d ",sccno[i]);
     putchar();
     return ;
 }

 /*
 //测试用图
 6 8
 1 3
 3 5
 5 6
 3 4
 4 6
 4 1
 1 2
 2 4
 */