KS求有向图强连通分量模板

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long ll;

 int n,m;
 const int maxn=1e5+;
 const int maxm=*maxn;
 struct nodeg
 {
     int to;
     int nxt;
 }eg[maxm];
 int headg[maxn];
 int headgr[maxn];
 struct nodegr
 {
     int to;
     int nxt;
 }egr[maxn];
 int totg;
 int totgr;
 int id=;
 int num[maxn];
 bool vis[maxn];
 int root;
 int cnt;
 void init()
 {
     memset(headg,-,sizeof(headg));
     memset(headgr,-,sizeof(headgr));
     totg=;
     totgr=;
     cnt=;
 }
 void addg(int u,int v)
 {
     eg[totg].to=v;
     eg[totg].nxt=headg[u];
     headg[u]=totg++;
 }
 void addgr(int u,int v)
 {
     egr[totgr].to=v;
     egr[totgr].nxt=headgr[u];
     headgr[u]=totgr++;
 }
 vector<int> scc[maxn];
 void DFS(int u)
 {
     vis[u]=true;
     for(int i=headg[u];i!=-;i=eg[i].nxt)
     {
         int v=eg[i].to;
         if(!vis[v])
         {
             DFS(v);
         }
     }
     num[++id]=u;
 }
 void RDFS(int u,int k)
 {
     vis[u]=true;
     scc[k].push_back(u);
     for(int i=headgr[u];i!=-;i=egr[i].nxt)
     {
         int v=egr[i].to;
         if(!vis[v])
         {
             RDFS(v,k);
         }
     }
 }
 void SCC()
 {
         memset(vis,false,sizeof(vis));
         for(int i=;i<=n;i++)
         {
             if(!vis[i])
             {
                 DFS(i);
             }
         }
         memset(vis,false,sizeof(vis));
         for(int i=;i<=n;i++)
         {
             scc[i].clear();
         }
         cnt=;
         for(int i=id;i>=;i--)
         {
             if(!vis[num[i]]) RDFS(num[i],++cnt);
         }
 }
 void printSC()
 {
     for(int i=;i<=cnt;i++)
     {
         int sz=scc[i].size();
         for(int j=;j<sz;j++)
             printf("%d ",scc[i][j]);
         puts("");
     }
 }

 int main()
 {
     while(~scanf("%d%d",&n,&m))
     {
         init();
         int u,v;
         for(int i=;i<=m;i++)
         {
             scanf("%d%d",&u,&v);
             addg(u,v);
             addgr(v,u);
         }
         SCC();
         printSC();
     }

     return ;
 } 

Kosaraju-Sharir算法

http://www.cnblogs.com/llhthinker/p/4954082.html