讲 Kosaraju 算法之前要先知道什么是强连通分量(SCC)
强连通分量
-
对于一个有向图顶点的子集 S ,如果在 S 内取两个顶点 u 和 v ,都能找到一条从 u 到 v 的路径,那么就称 S 是强连通的。
-
如果在强连通的顶点集合 S 中加入其他任意顶点集合后,它都不再是强连通的,那么就称 S 是原图的一个强连通分量(SCC:strongly connected component)。
-
任意有向图都可以分解成若不相干的强连通分量,这就是强连通分量分解。
-
把分解后的强连通分量缩成一个顶点,就得到了一个DAG(有向无环图)。
如图,图中有 3 个强连通分量(SCC),将三个强连通分量缩成一个顶点,就得到了有向无环图(DAG)。
Kosaraju 算法
-
Kosaraju 算法可以计算出一个有向图的强连通分量(SCC)。
-
Kosaraju 算法基于两个原理:
(1)一个有向图 G,把 G 的所有的边反向,建立反图 rG,反图 rG 不会改变原图 G 的强连通性。也就是说,图 G 的 SCC 数量与 rG 的 SCC 数量相同。
(2)对原图 G 和反图 rG 各做一次 DFS,可以确定 SCC 数量。
-
对原图 G 做一次 DFS 是为了确定点的先后顺序,用到了拓扑排序。在DFS过程中,把递归到最底层的那个点标记为最小,然后回退过程中,其他点的标记逐个递增。即优先级小的点先存,优先级大的点后存。
-
对反图 rG 做一次 DFS ,顺序从标记最大的点开始到最小的点。为什么顺序要相反?因为在反图中,方向相对于原图都相反了,优先级大的点由出度多,变成了入度多,那么 DFS 过程中就会被反边堵住,也就是说,此时能搜索到的点都是同个强连通分量的点。按照这个顺序依次搜索下去。
步骤:
(1)先记录原图 G 和反图 rG
(2)对原图所有的点 DFS 一遍,标记点的先后顺序
(3)根据点的逆序 DFS 一遍,能搜到的点就是在同一个强连通分量。
例题: HDU - 1269 迷宫城堡
题意: 一个有向图,有 n 个点(n <= 10000)和 m 条边(m <= 100000)。判断整个图是否强连通,如果是,输出 “Yes”,否则,输出“No”。
思路: Kosaraju 算法模板题,直接求整个图是否为强连通分量。
Code:
#include <iostream>
#include <cstring>
#include <vector>
using namespace std;
const int N = 100010;
vector<int> G[N],rG[N];
vector<int> S; // 存第一次dfs1()的结果,即标记点的先后顺序,优先级小的点先进
int vis[N]; // vis[i]标记第一次dfs1()点i是否访问过
int sccno[N]; // sccno[i]标记点i属于第几个强连通分量,同时记录dfs2()过程中点i是否访问过
int cnt; //cnt表示强连通分量的个数
void dfs1(int u){
if(vis[u]) return;
vis[u] = 1;
for(int i=0; i<G[u].size(); i++)
dfs1(G[u][i]);
S.push_back(u); //记录点的先后顺序,按照拓扑排序,优先级大的放在S的后面
}
void dfs2(int u){
if(sccno[u]) return;
sccno[u] = cnt;
for(int i=0; i<rG[u].size(); i++)
dfs2(rG[u][i]);
}
void Kosaraju(int n) {
cnt = 0;
S.clear();
memset(vis,0,sizeof(vis));
memset(sccno,0,sizeof(sccno));
for(int i=1; i<=n; i++) //搜索所有点
dfs1(i);
for(int i=n-1; i>=0; i--){
if(!sccno[S[i]]){
cnt++;
dfs2(S[i]);
}
}
}
int main() {
int n,m,u,v;
//坑点,这里 n 或者 m 都不能为 0
while(cin >> n >> m && (n || m)) {
for(int i=0; i<n; i++) {
G[i].clear();
rG[i].clear();
}
for(int i=0; i<m; i++) {
cin >> u >> v;
G[u].push_back(v); // 原图
rG[v].push_back(u); // 反图
}
Kosaraju(n);
cnt==1? cout<<"Yes"<<endl:cout<<"No"<<endl;
}
return 0;
}