Tarjan+缩点【强连通分量】【模板】

 #include<iostream>
 #include<cstring>
 #include<algorithm>
 #include<queue>
 #include<stack>
 #include<vector>
 using namespace std;
 const int maxn=;//边数
 const int maxn1=;//顶点数
 struct edge{
     int from;
     int to;
     int next;
 }EDGE[maxn];
 vector<int>vc[maxn1];
 int head[maxn1],dfn[maxn1],vis[maxn1],low[maxn1],col[maxn1],out[maxn1],in[maxn1],en[maxn1],stk[maxn1];//各个变量的意义可参照上篇博客
 int edge_cnt=,tot1=,tot2=,scc_cnt=,tot0=;
 void add(int x,int y)
 {
     EDGE[edge_cnt].from=x;
     EDGE[edge_cnt].to=y;
     EDGE[edge_cnt].next=head[x];
     head[x]=edge_cnt++;
 }
 void Tarjan(int u)
 {
     low[u]=dfn[u]=++tot1;//注意tot1的初值必须是1【因为dfn必须为正数】，所以这里使用++tot1而不用tot1++;
     vis[u]=;
     stk[++tot2]=u;
     for(int i = head[u]; i != -  ; i = EDGE[i].next)
     {
         if(!dfn[EDGE[i].to]){
             Tarjan(EDGE[i].to);
             low[u]=min(low[u],low[EDGE[i].to]);
         }
         else if(vis[EDGE[i].to]){
             low[u]=min(low[u],low[EDGE[i].to]);
         }
     }
             if(low[u]==dfn[u]){
             int xx;
             scc_cnt++;//注意scc_cnt也是从1开始的，因为要染色，区别于为染色的0
             do{
                 xx=stk[tot2--];
                 vc[scc_cnt].push_back(xx);
                 col[xx]=scc_cnt;
                 vis[xx]=;
             }while(xx!=u);
         }
 }
 void INIT()
 {
     for(int i =  ; i < maxn1 ; i++)
     vc[i].clear();
     edge_cnt=,tot1=,tot2=,scc_cnt=,tot0=;
     memset(head,-,sizeof(head));
     memset(stk,,sizeof(stk));
     memset(in,,sizeof(in));
     memset(out,,sizeof(out));
     memset(dfn,,sizeof(dfn));
     memset(low,,sizeof(low));
     memset(col,,sizeof(col));
  }
  void suodian()//缩点
  {
      for(int i =  ; i < edge_cnt ; i++)
     {
         if(col[EDGE[i].from]!=col[EDGE[i].to])
         {
             in[col[EDGE[i].to]]++;//缩点
                out[col[EDGE[i].from]]++;
         }
     }
   }
 int main()
 {
     int n,m;
     scanf("%d%d",&n,&m);
     INIT();
     while(m--)
     {
         int a,b;
         scanf("%d%d",&a,&b);
         add(a,b);
     }
     for(int i =  ; i <= n; i++)
     {
         if(!dfn[i])Tarjan(i);
     }
     suodian();
     return ;
 }
 /*4 5
 1 3
 2 4
 4 2
 1 4
 2 1*/