目录
Tarjan算法
预备知识
有向图:由有向边的构成的图。需要注意的是这是Tarjan算法的前提和条件。
强连通:如果两个顶点可以相互通达,则称两个顶点 强连通(strongly connected)。如果有向图G的每两个顶点都 强连通,称G是一个强连通图。非 强连通图有向图的极大强连通子图,称为强连通分量(strongly connected components)。
在这个有向图中1、2、3、4四个点可以互相到达,就称这四个点组成的子图为强连通分量。且这四个点两两强连通。然后就可以开始学习神奇的Tarjan算法了!
Tarjan算法是用来求强连通分量的,它是一种基于DFS(深度优先搜索)的算法,每个强连通分量为搜索树中的一棵子树。并且运用了数据结构栈。
在介绍详细原理前,先引入两个非常重要的数组:dfn[ ] 与 low[ ]
dfn[ ]:就是一个时间戳(被搜到的次序),一旦某个点被DFS到后,这个时间戳就不再改变(且每个点只有唯一的时间戳)。所以常根据dfn的值来判断是否需要进行进一步的深搜。
low[ ]:该子树中,且仍在栈中的最小时间戳,像是确立了一个关系,low[ ]相等的点在同一强连通分量中。
注意初始化时 dfn[ ] = low[ ] = ++cnt.
算法思路
首先这个图不一定是一个连通图,所以跑Tarjan时要枚举每个点,若dfn[ ] == 0,进行深搜。
然后对于搜到的点寻找与其有边相连的点,判断这些点是否已经被搜索过,若没有,则进行搜索。若该点已经入栈,说明形成了环,则更新low.
在不断深搜的过程中如果没有路可走了(出边遍历完了),那么就进行回溯,回溯时不断比较low[ ],去最小的low值。如果dfn[x]==low[x]则x可以看作是某一强连通分量子树的根,也说明找到了一个强连通分量,然后对栈进行弹出操作,直到x被弹出。
代码
#include#include #include #include #include #include #define N 75005using namespace std;vector eg[N];//邻接表存边stack sk;//栈int idx = 0;int dfn[N], low[N], ins[N], cnt = 0;void tarjan(int u){ dfn[u] = low[u] = ++idx; sk.push(u); ins[u] = 1; for (int v : eg[u]) { if (!dfn[v]) { tarjan(v); low[u] = min(low[u], low[v]);//回溯时更新low[ ],取最小值 } else if (ins[v]) { low[u] = min(low[u], dfn[v]);//一旦遇到已入栈的点,就将该点作为连通量的根 //这里用dfn[e[i].v]更新的原因是:这个点可能 //已经在另一个强连通分量中了但暂时尚未出栈,所 //以now不一定能到达low[e[i].v]但一定能到达 //dfn[e[i].v]. } } int tem; if (dfn[u] == low[u]) { cnt++;//联通分量个数加 1 do { tem = sk.top(); sk.pop(); ins[tem] = 0; } while (tem != u);//出栈 }}int main(){ int n, m, x, y; while (scanf("%d%d", &n, &m) && n) { idx = 0; cnt = 0; memset(ins, 0, sizeof(ins)); memset(low, 0, sizeof(low)); memset(dfn, 0, sizeof(dfn)); for (int i = 0; i < m; i++) { cin >> x >> y; eg[x].push_back(y); } for(int i=1;i<= n;i++)//防止有根本就不连通的点 if(!dfn[i]) tarjan(i); cout << cnt << endl; for (int i = 1; i <= n; i++) { eg[i].clear(); } }}