Tarjan算法
一种求解有向图强连通分量的线性时间的算法。
也是一个基于DFS深度优先搜索的一个有向图
强连h通：如果两个结点能够相互到达，则称这两个结点强连通
强连接图：如果有向图G中任意两个结点强连通，则称这个有向图G为强连接图
强连通分量：对有向图的一个子图，这个子图任意两个点满足强连通，则称这个子图为有向图G的强连通分量
栗子：
对于这个图，1，2两个点是强连通关系，因为1可以直接到2，2可以通过3到1；
对1，2，3，这三个结点是整个图的一个子图，任意两个点可以相互到达，是强连通分量；
一： 要搞懂tarjan算法就需要理解三个知识点 （dfn[]数组，len[]数组，堆栈）
1：dfn[]数组：存储的是当前结点是第几个被搜索的，作为这个点搜索的时间戳，（每个点的时间戳都不一样，且不能修改）
2：len[]数组：存储的是该子图中，且仍在栈中的最小的时间戳（当所有点都更新完后，同意强连通分量中每个len[i]都相等，且是最小的）
3：为了存储整个强连通分量，需要用到堆栈为容器；每到一个新的结点，就入栈，遍历该节点后是否有后继节点（这里可以用链式前向星的方法，百度上有详解），并不断更新low[]数组的值，直到dfn[u]==low[u]，就说明这个结点是这个强连通分量的根结点，找到根结点后，就将栈中所有结点出栈，组成一个新的强连通分量。
// 可以用一个数组sta[]来实现堆栈
//scc[]数组存储每个强连通分量的编号
void tarjan(int u)
{
dfn[u]=low[u]=++index; //先对u的dfn[]和low[]赋初值，这里的u为新的结点，可在主函数里判断一下
sta[top++]=u; //把u存入栈中
for(int i=head[u];i!=-1;i=edge[i].next) //链式前向星
{
int t=edge[i].v;
if(!dfn[t]) // 如果t这个结点还没有去过
{
tarjan(t); //继续往下找
low[u]=min(low[u],low[t]); //更新low[u] 数组
}
else if(!scc[t]) //scc[]数组存储的是该结点在第几个强连通分量中
{
low[u]=min(low[u],low[t]); //同样也需要更新
}
}
if(dfn[u]==low[u]) // 如果u是这个强连通分量的根节点
{
sig++; //sig是每个强连通分量的编号
while(1)
{
int j=sta[--top]; //一个个的出栈
scc[j]=sig; //每个j所在的强连通分量的编号都是sig
if(j==u) //知道根节点u出栈就借书
break;
}
}
二： 算法演示
借鉴：http://www.cnblogs.com/shadowland/p/5872257.html
以1为Tarjan 算法的起始点，如图
顺次DFS搜到节点6
回溯时发现LOW[ 5 ]==DFN[ 5 ] , LOW[ 6 ]==DFN[ 6 ] ,则{ 5 } , { 6 } 为两个强连通分量。回溯至3节点，拓展节点4.
拓展节点1 ， 发现1再栈中更新LOW[ 4 ]，LOW[ 3 ] 的值为1
回溯节点1，拓展节点2
自此，Tarjan 算法就结束了，{1 , 2 , 3 , 4 } , { 5 } , { 6 } 为图中的三个强连通分量。
不难发现，Tarjan算法的时间复杂度为O(E+V).
例题：A - 迷宫城堡 https://vjudge.net/contest/284357#problem
为了训练小希的方向感，Gardon建立了一座大城堡，里面有N个房间(N<=10000)和M条通道(M<=100000)，每个通道都是单向的，就是说若称某通道连通了A房间和B房间，只说明可以通过这个通道由A房间到达B房间，但并不说明通过它可以由B房间到达A房间。Gardon需要请你写个程序确认一下是否任意两个房间都是相互连通的，即：对于任意的i和j，至少存在一条路径可以从房间i到房间j，也存在一条路径可以从房间j到房间i。
Input
输入包含多组数据，输入的第一行有两个数：N和M，接下来的M行每行有两个数a和b，表示了一条通道可以从A房间来到B房间。文件最后以两个0结束。
Output
对于输入的每组数据，如果任意两个房间都是相互连接的，输出"Yes"，否则输出"No"。
Sample Input
3 3 1 2 2 3 3 1 3 3 1 2 2 3 3 2 0 0
Sample Output
Yes No
已AC代码
#include<stdio.h>
#include<string.h>
#include<algorithm>
typedef long long ll;
const int xmax=1e6+100;
const int INF=1e9+7;
using namespace std;
int cnt,index;
int dfn[xmax];
int scc[xmax],sig;
int low[xmax];
int sta[xmax],top;
struct node
{
int v;
int next;
}edge[xmax];
int head[xmax];
void init()
{
memset(head,-1,sizeof(head));
memset(dfn,0,sizeof(dfn));
memset(sta,0,sizeof(sta));
memset(scc,0,sizeof(scc));
cnt=0;index=0;top=0;sig=0;
}
void addedge(int u,int v)
{
edge[cnt].v=v;
edge[cnt].next=head[u];
head[u]=cnt++;
}
void tarjan(int u)
{
dfn[u]=low[u]=++index;
sta[top++]=u;
for(int i=head[u];i!=-1;i=edge[i].next)
{
int t=edge[i].v;
if(!dfn[t])
{
tarjan(t);
low[u]=min(low[u],low[t]);
}
else if(!scc[t])
{
low[u]=min(low[u],low[t]);
}
}
if(dfn[u]==low[u])
{
sig++;
while(1)
{
int j=sta[--top];
scc[j]=sig;
if(j==u)
break;
}
}
}
int main()
{
int n,m;
while(scanf("%d%d",&n,&m)!=EOF&&(n||m))
{
init();
int a,b;
while(m--)
{
scanf("%d%d",&a,&b);
addedge(a,b);
}
for(int i=1;i<=n;i++)
{
if(!dfn[i]) //判断i是否是新结点
tarjan(i);
}
if(sig!=1)
printf("No\n");
else
printf("Yes\n");
}
}