[模板]tarjan求强连通分量

大约是今年4月学的算法了，后来5月的时候做题还写了一个退化的tarjanQAQ。

时间复杂度：O(n+m)

用途：有向图缩环

 #include<set>
 #include<cmath>
 #include<ctime>
 #include<queue>
 #include<stack>
 #include<cstdio>
 #include<vector>
 #include<cstring>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>
 #define N 10001
 #define M 100001
 using namespace std;
 struct graph{
     int nxt,to;
 }e[M];
 int g[N],f[N]/*father*/,dfn[N],low[N],sta[N]/*stack;stack[0]:the top of the stack*/,m,n,cnt;
 /*dfn[u]:遍历到u点的时间; low[u]:u点可到达的各点中最小的dfn[v],即最高层的点*/
 bool ins[N];//in stack
 inline void addedge(int x,int y){
     e[++cnt].nxt=g[x];g[x]=cnt;e[cnt].to=y;
 }
 inline void tarjan(int u){
     dfn[u]=low[u]=++cnt;
     sta[++sta[]]=u;ins[u]=true;
     for(int i=g[u];i;i=e[i].nxt)
         if(!dfn[e[i].to]){
             tarjan(e[i].to);
             low[u]=min(low[u],low[e[i].to]);
         }
         else if(ins[e[i].to])
             low[u]=min(low[u],low[e[i].to]);
     if(dfn[u]==low[u]){
         while(sta[sta[]+]!=u){
             f[sta[sta[]]]=u;
             ins[sta[sta[]--]]=false;
         }
     }
 }
 inline void init(){
     scanf("%d%d",&n,&m);
     for(int i=,j,k;i<=m;i++){
         scanf("%d%d",&j,&k);
         addedge(j,k);
     }
     cnt=;
     for(int i=;i<=n;i++)
         if(!dfn[i]) tarjan(i);
     for(int i=;i<=n;i++)
         printf("%d's father is %d\n",i,f[i]);
 }
 int main(){
     freopen("tarjan.in","r",stdin);
     freopen("tarjan.out","w",stdout);
     init();
     fclose(stdin);
     fclose(stdout);
     return ;
 }