[模板]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[0]]=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[0]+1]!=u){
            f[sta[sta[0]]]=u;
            ins[sta[sta[0]--]]=false;
        }
    }
}
inline void init(){
    scanf("%d%d",&n,&m);
    for(int i=1,j,k;i<=m;i++){
        scanf("%d%d",&j,&k);
        addedge(j,k);
    }
    cnt=0;
    for(int i=1;i<=n;i++)
        if(!dfn[i]) tarjan(i);
    for(int i=1;i<=n;i++)
        printf("%d's father is %d
",i,f[i]);
}
int main(){
    freopen("tarjan.in","r",stdin);
    freopen("tarjan.out","w",stdout);
    init();
    fclose(stdin);
    fclose(stdout);
    return 0;
}