form:Christopher Yan
概念
如果两个顶点可以相互通达，则称两个顶点强连通
如果有向图 (G) 的每两个顶点都强连通，称 (G) 是一个强连通图.
非强连通图有向图的极大强连通子图，称为强连通分量
(Tarjan) 算法是用来求强连通分量的，它是一种基于(DFS)（深度优先搜索）的算法，每个强连通分量为搜索树中的一棵子树。并且运用了数据结构栈。
在介绍详细原理前，先引入两个非常重要的数组：(dfn[ ]) 与 (low[ ])
(dfn[ ])：就是一个时间戳（被搜到的次序），一旦某个点被 (DFS) 到后，这个时间戳就不再改变（且每个点只有唯一的时间戳）。所以常根据 (dfn) 的值来判断是否需要进行进一步的深搜
(low[ ])：该子树中，且仍在栈中的最小时间戳，像是确立了一个关系，(low[ ]) 相等的点在同一强连通分量中。
注意初始化时 (dfn[ ] = low[ ] = ++cnt.)
算法思路：
首先这个图不一定是一个连通图，所以跑 (Tarjan) 时要枚举每个点，若 (dfn[ ] == 0)，进行深搜。
然后对于搜到的点寻找与其有边相连的点，判断这些点是否已经被搜索过，若没有，则进行搜索。若该点已经入栈，说明形成了环，则更新 (low).
在不断深搜的过程中如果没有路可走了（出边遍历完了），那么就进行回溯，回溯时不断比较 (low[ ])，去最小的 (low) 值。如果 (dfn[x]==low[x]) 则 (x) 可以看作是某一强连通分量子树的根，也说明找到了一个强连通分量，然后对栈进行弹出操作，直到 (x) 被弹出
手动模拟一下过程：
从1进入 dfn［1］＝ low［1］＝ ＋＋cnt = 1
入栈 1
由1进入2 dfn［2］＝low［2］＝ ＋＋cnt = 2
入栈 1 2
之后由2进入4 dfn［4］＝low［4］＝ ＋＋cnt = 3
入栈 1 2 4
之后由4进入 6 dfn［6］＝low［6］＝＋＋cnt = 4
入栈 1 2 4 6
6无出度，之后判断 (dfn［6］＝＝low［6］)
说明6是个强连通分量的根节点：(6) 及 (6) 以后的点出栈并输出。
回溯到4后发现4找到了一个已经在栈中的点 (1)，更新 low [ 4 ] = min ( low [ 4 ] , dfn [ 1 ] )
于是 low [ 4 ] = 1 .
由4继续回到2 Low[2] = min ( low [ 2 ] , low [ 4 ] ).
low［2］＝1；
由2继续回到1 判断 low[1] = min ( low [ 1 ] , low [ 2 ] ).
low［1］还是 1
然后更新3的过程省略，大家可以自己手动模拟一下。
省略了 (1->3) 的更新过程之后，(1)的所有出边就跑完了
于是判断： (low [ 1 ] == dfn [ 1 ]) 说明以 (1) 为根节点的强连通分量已经找完了。
将栈中 (1) 以及 (1) 之后进栈的所有点，都出栈并输出
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
inline int read(){
int x = 0, f = 1;char ch = getchar();
while (ch < '0' || ch > '9') {if(ch == '-') f = -1; ch = getchar();}
while (ch >= '0' && ch <= '9') {x = x * 10 + ch - '0';ch = getchar();}
return x * f;
}
int n,m,x,y,top=0,cnt=0,t,col;
int ans1=-1,ans2=-1,ans3=-1;
int d[200020];
int a[200020];
int c[200020];
int f[200020];
int dfn[200020];
int low[200020];
int stack[200020];
bool v[200020];
struct edge{
int u;
int v;
int w;
int next;
}e[1000020];
void Add(int u,int v,int w){
e[++top].v=v;
e[top].w=w;
e[top].next=f[u];
f[u]=top;
}
void tarjan(int now)
{
dfn[now]=low[now]=++cnt;
stack[++t]=now;
v[now]=1;
for(int i=f[now];i!=-1;i=e[i].next)
if(!dfn[e[i].v]) {
tarjan(e[i].v);
low[now]=min(low[now],low[e[i].v]);
}
else if(v[e[i].v])
low[now]=min(low[now],dfn[e[i].v]);
int cur;
if(dfn[now]==low[now]){
do
{
cur=stack[t--];
v[cur]=false;
printf("%d ",cur);
}while(now!=cur);
printf("
");
}
}
int main()
{
n=read();
m=read();
memset(f,-1,sizeof f);
for(int i=1;i<=n;++i)
a[i]=read();
for(int i=1;i<=m;++i)
{
x=read();
y=read();
Add(x,y,0);
}
for(int i=1;i<=n;++i)
if(!dfn[i]) tarjan(i);
return 0;
}