DFS性质的应用——利用Tarjan算法求割顶、BCC、SCC
整理自《算法竞赛入门经典——训练指南》以及网络
DFS (depth first search)深度优先搜索算法
dfs森林:按照dfs的执行顺序,将图的所有边重新梳理,分为四个类别:前向边、反向边、交叉边和树边。在无向图中不存在交叉边,前向边与后向边等价。
关键变量:
pre[u]:记录u点被访问到的次序。 pre[u] = ++dfs_clock;
low[u]:在DFS过程中,u及其后代能连回的最早祖先的pre值。 low[u] = min{low[v] | u -> v}
$1 计算割顶和桥
- 割顶:对于连通图,删除后使图不再连通的点。
- 桥:对于连通图,删除后使图不再连通的边。
计算割顶的方法:在DFS过程中,如果一个点u存在一个子节点v,使得v及其后代都没有反向边连回u的祖先(不包括u),即lowv >= pre[u],则u是割顶。
计算桥的方法:如果v的后代只能连回v自己(即low(v) > pre(u))则u-v是桥。
注意:
- 对于已访问点,只处理反向边(条件pre[v] < pre[u])。前向边的pre已被传递过,不需要处理。
- 根节点需要特判:当DFS树根只有一个孩子时不是割顶,需手动取消割顶标记。
$2 计算BCC和SCC
- 无向图的(点)双连通分量(BCC):内部无割顶,即任意两条边都在一个简单环中。
- 无向图的边-双连通分量:所有边都不是桥,即每条边都至少在一个简单环中。
计算BCC的方法:将边入栈,在找到割顶后将这条边及其之后入栈的边(子树的边)出栈,对点进行标记。
计算edge-BCC的方法:找到桥后再做一次不经过桥的DFS即可。
- 有向图的强连通分量(SCC):分量内的点相互可达
计算SCC的方法:将点入栈,DFS过程中记录lowlink,当lowlink(u) == pre(u)时将其中的点出栈标记。
注意:
- 当遇到已访问点时,要忽略已经确定编号的点(条件!sccno[v])
- lowlinku的最小值取任意一个都可以(都比树根的小):
lowlinku = min(lowlinku, pre[v] OR lowlinkv);
- BCC算法会将单独连接的两个点也标记为BCC(两个点虽然不能称为环,但他们都不是割顶,符合BCC的定义。需要根据实际情况特判)
- SCC算法会将每一个点都划分到一个SCC,即单独的点也会被scc_count编号,这使得缩点的过程更为方便。
贴上代码
计算BCC
#include <iostream>
#include <stack>
#include <vector>
#include <algorithm>
using namespace std;
const int maxn = 2e4 + 1;
int dfs_clock, pre[maxn], iscut[maxn], bccno[maxn], bcc_count;
struct Edge {
int u, v;
};
vector <int> bcc[maxn], G[maxn];
stack <Edge> S;
void AddE(int u, int v) {
G[u].push_back(v);
}
int tarjan(int u, int fa) {
int lowu = pre[u] = ++dfs_clock, child = 0;
for (auto v: G[u]) {
if (!pre[v]) {
child++;
S.push((Edge) {u, v});
int lowv = tarjan(v, u);
lowu = min(lowu, lowv);
if (lowv >= pre[u]) {
iscut[u] = 1;
++bcc_count;
bcc[bcc_count].clear();
//cout<<dfs_clock<<' '<<u<<'>'<<v<<','<<lowv<<'|'<<bcc_count<<endl;
while (1) {
Edge e = S.top(); S.pop();
if (bccno[e.u] != bcc_count) {
bcc[bcc_count].push_back(e.u);
bccno[e.u] = bcc_count;
}
if (bccno[e.v] != bcc_count) {
bcc[bcc_count].push_back(e.v);
bccno[e.v] = bcc_count;
}
if (e.u == u && e.v == v) break;
}
}
} else if (pre[v] < pre[u] && v != fa) {
S.push((Edge) {u, v});
lowu = min(lowu, pre[v]); //Err
}
}
if (fa == -1 && child == 1) iscut[u] = 0; //Err
return lowu;
}
signed main() {
int n, m, u, v, i;
cin>>n>>m;
while (m--) {
cin>>u>>v;
AddE(u, v);
AddE(v, u);
}
for (i = 1; i <= n; ++i)
if (!pre[i])
tarjan(i, -1);
cout<<count(iscut + 1, iscut + n + 1, 1)<<endl;
for (i = 1; i <= n; ++i)
if (iscut[i])
cout<<i<<' ';
return 0;
}
Tarjan求SCC+缩点+拓扑排序
#include <iostream>
#include <vector>
#include <stack>
using namespace std;
const int maxn = 1e4 + 1;
int n;
int pre[maxn], dfs_clock, sccno[maxn], scc_count;
int c[maxn], topo[maxn];
vector <int> G[maxn], GA[maxn], scc[maxn];
stack <int> S;
void AddE(int u, int v) {
G[u].push_back(v);
}
void AddED(int u, int v) {
//cout<<"Add"<<u<<' '<<v<<endl;
GA[u].push_back(v);
}
int tarjan(int u) {
int lowlinku = pre[u] = ++dfs_clock;
S.push(u);
for (auto x: G[u])
if (!pre[x]) {
int lowlinkv = tarjan(x);
lowlinku = min(lowlinku, lowlinkv);
} else if (!sccno[x])
lowlinku = min(lowlinku, pre[x]);
if (lowlinku == pre[u]) {
++scc_count;
scc[scc_count].clear();
//cout<<endl<<scc_count<<":";
while (1) {
int v = S.top(); S.pop();
scc[scc_count].push_back(v);
sccno[v] = scc_count;
//cout<<v<<' ';
if (v == u) break;
}
}
return lowlinku;
}
int dfs(int u) {
static int cnt = scc_count;
if (c[u] == -1)
return false;
c[u] = -1;
for (auto v: GA[u])
if (!c[v] && !dfs(v)) return false;
c[u] = 1;
topo[cnt--] = u;
return true;
}
int toposort(int n) {
for (int i = 1; i <= n; ++i)
if (!c[i])
if (!dfs(i))
return false;
return true;
}
signed main() {
int m, i, u, v;
cin>>n>>m;
while (m--) {
cin>>u>>v;
AddE(u, v);
}
for (i = 1; i <= n; ++i) {
if (!pre[i])
tarjan(i);
}
for (i = 1; i <= n; ++i) //Err Impppp
for (int v: G[i]) if (sccno[i] != sccno[v]) //Err2.
AddED(sccno[i], sccno[v]);
if (!toposort(scc_count))
cout<<"Shit!"<<endl;
return 0;
}
来源:oschina
链接:https://my.oschina.net/u/4302647/blog/3400535