Tarjan强连通(自学版)

强连通分量

　　首先我们把强连通分量看做一个齿轮或环（他会转啊），不考虑其他的限制则可认为分量结点可以互换（就是转一下）不会影响分量中包含的结点（为什么这么想呢？你理解tarjan时中的low值有帮助）

Tarjan算法是基于对图深度优先搜索的算法，每个强连通分量为搜索树中的一棵子树。搜索时，把当前搜索树中未处理的节点加入一个堆栈，回溯时可以判断栈顶到栈中的节点是否为一个强连通分量。
定义DFN(u)为节点u搜索的次序编号(时间戳)，Low(u)为u或u的子树能够追溯到的最早的栈中节点的次序号。由定义可以得出，
Low(u)=Min

{    DFN(u),    Low(v),(u,v)为树枝边，u为v的父节点    DFN(v),(u,v)为指向栈中节点的后向边(非横叉边)}当DFN(u)=Low(u)时，以u为根的搜索子树上所有节点是一个强连通分量。
伪代码：tarjan(u){        
DFN[u]=Low[u]=++Index                      // 为节点u设定次序编号和Low初值             
Stack.push(u)                              // 将节点u压入栈中      
for each (u, v) in E                       // 枚举每一条边                
if (v is not visted)               // 如果节点v未被访问过                        
tarjan(v)                  // 继续向下找                        
Low[u] = min(Low[u], Low[v])                
else if (v in S)                   // 如果节点v还在栈内                        
Low[u] = min(Low[u], DFN[v])        
if (DFN[u] == Low[u])                      // 如果节点u是强连通分量的根               
repeat                        
v = S.pop                  // 将v退栈，为该强连通分量中一个顶点                       
print v                
until (u== v)}
cpp： void tarjan(int i) { int j; DFN[i]=LOW[i]=++Dindex; instack[i]=true; Stap[++Stop]=i; for (edge *e=V[i];e;e=e->next) { j=e->t; if (!DFN[j]) { tarjan(j); if (LOW[j]<LOW[i]) LOW[i]=LOW[j]; } else if (instack[j] && DFN[j]<LOW[i]) LOW[i]=DFN[j]; } if (DFN[i]==LOW[i]) { Bcnt++; do { j=Stap[Stop--]; instack[j]=false; Belong[j]=Bcnt; } while (j!=i); } } void solve() { int i; Stop=Bcnt=Dindex=0; memset(DFN,0,sizeof(DFN)); for (i=1;i<=N;i++) if (!DFN[i]) tarjan(i); }

pascal：procedure dfs(s:longint); var ne:longint; begin view[s]:=1; //view[i]表示点i的访问状态.未访问,正访问,已访问的点,值分别为0,1,2 inc(top);stack[top]:=s; //当前点入栈 inc(time);rea[s]:=time;low[s]:=time; //记录访问该点的真实时间rea和最早时间low ne:=head[s]; while ne<>0 do begin if view[e[ne]]=0 then dfs(e[ne]); //如果扩展出的点未被访问,继续扩展 if view[e[ne]]<2 then low[s]:=min(low[s],low[e[ne]]); //如果扩展出的不是已访问的点,更新访问源点s的最早时间. //容易理解,如果一个点能到达之前访问过的点,那么路径中存在一个环使它能更早被访问 ne:=next[ne]; end; if rea[s]=low[s] then begin //如果s的最早访问时间等于其实际访问时间,则可把其视作回路的"始点" inc(tot); //连通块编号 while stack[top+1]<>s do begin //将由s直接或间接扩展出的点标记为同一连通块,标记访问后出栈 lab[stack[top]]:=tot; //lab[i]表示点i所属的连通块 view[stack[top]]:=2; dec(top); end; end; end;