有向图强连通分量Kosaraju算法

最关键通用部分：强连通分量一定是图的深搜树的一个子树。

反图定义：如果有向图存在边u->v,那么v->u属于改图的反图，记为GT。

基本思路：
这个算法可以说是最容易理解，最通用的算法，其比较关键的部分是同时应用了原图G和反图GT。(步骤1)先用对原图G进行深搜形成森林(树)，（步骤2）然后任选一棵树对其进行深搜(注意这次深搜节点A能往子节点B走的要求是EAB存在于反图GT)，能遍历到的顶点就是一个强连通分量。余下部分和原来的森林一起组成一个新的森林，继续步骤2直到 没有顶点为止。
改进思路：
当然，基本思路实现起来是比较麻烦的(因为步骤2每次对一棵树进行深搜时，可能深搜到其他树上去，这是不允许的，强连通分量只能存在单棵树中(由开篇第一句话可知))，我们当然不这么做，我们可以巧妙的选择第二深搜选择的树的顺序，使其不可能深搜到其他树上去。想象一下，如果步骤2是从森林里选择树，那么哪个树是不连通(对于GT来说)到其他树上的呢？就是最后遍历出来的树，它的根节点在步骤1的遍历中离开时间最晚，而且可知它也是该树中离开时间最晚的那个节点。那么为什么最后遍历出来的树的反图不存在到其他树的边呢？因为他是最后的树，如果它的反图存在到其他树的边那么就是在原图中，存在其他树到这棵树的边，那么这棵树必然在其他树的遍历中别包含进去成为别人的子树，而不是最后遍历出来的树。所以在反图中只存在其他树到这个树的边，因此在对这棵树的反图的深搜中不会到别的树去。这给我们提供了很好的选择，在第一次深搜遍历时，记录时间i离开的顶点j，即numb[i]=j。那么，我们每次只需找到没有找过的顶点中具有最晚离开时间的顶点直接深搜(对于GT来说)就可以了。每次深搜都得到一个强连通分量。

基本步骤： 

1.深搜后序遍历有向图，记录每个点的访问顺序(order);
2.按照1中的顺序(order)的逆顺序，再次深搜遍历原先有向图的转置图(原先的边反向);
3.步骤2中，能在一次dfs中遍历到的点，在同一个强连通分量中;

vector<int> adj[ MAXN ] ; //正向邻接表 vector<int> radj[ MAXN ] ; //反向邻接表 vector<int> ord; //后序访问顺序 void dfs1 ( int v ) { int i,u; for ( i=0 , vis[v]=true; i<adj[v].size() ; i++ ) if ( !vis[u=adj[v][i]] ) dfs1(u); ord.push_back(v); } void dfs2 ( int v ) { int i,u; for ( i=0,map[v]=cnt,vis[v]=true; i<radj[v].size() ; i++ ) if ( !vis[u=radj[v][i]] ) dfs2(u); } void kosaraju () { int i; memset(vis,false,sizeof(vis)); ord.clear(); for ( i=0 ; i<tot ; i++ ) if ( !vis[i] ) dfs1(i); memset(vis,false,sizeof(vis)); for ( cnt=0,i=ord.size()-1 ; i>=0 ; i-- ) if ( !vis[ord[i]] ) { dfs2(ord[i]); cnt++; } }

隐藏性质:

    分析到这里，我们已经知道怎么求强连通分量了。但是，大家有没有注意到我们在第二次深搜选择树的顺序有一个特点呢？它就是：如果我们把求出来的每个强连通分量收缩成一个点，并且用求出每个强连通分量的顺序来标记收缩后的节点，那么这个顺序其实就是强连通分量收缩成点后形成的有向无环图的拓扑序列。为什么呢？首先，应该明确搜索后的图一定是有向无环图呢？如果还有环，那么环上的顶点对应的所有原来图上的顶点构成一个强连通分量，而不是构成环上那么多点对应的独自的强连通分量了。然后就是为什么是拓扑序列，我们在改进分析的时候，不是先选的树不会连通到其他树上（对于反图GT来说），也就是后选的树没有连通到先选的树，也即先出现的强连通分量收缩的点只能指向后出现的强连通分量收缩的点。那么拓扑序列不是理所当然的吗？这就是Kosaraju算法的一个隐藏性质。