浅谈 Tarjan 算法

在了解 Tarjan 算法之前,我们先来了解 dfs 搜索树。

1 dfs 生成树

定义: dfs 遍历整张图,按照 dfs 序构成一棵树。

1.1 有向图的 dfs 生成树

在这里插入图片描述

有向图的 dfs 生成树包括四种边:

  • 树边(tree edge):图中黑色边表示。表示搜索访问到未访问过的结点。
  • 回边(back edge):图中橙色边表示。表示该边指向先前访问过的祖先。
  • 叉边(cross edge):图中蓝色边表示。表示该边指向先前访问过的非祖先结点。
  • 前向边(forward edge):图中红色边表示。表示该边指向先前访问过的孩子结点。

1.2 无向图的 dfs 生成树

在这里插入图片描述
无向图的 dfs 生成树仅包含两种边:树边和回边。

证明没有叉边和前向边:

  • 如果存在叉边 u → v u \to v uv,其中 u u u 为当前搜索访问到的结点, v v v 为之前访问过的非祖先结点,那么在访问 v v v 时,就应该访问了 u u u,与 u u u 为当前搜索访问到的结点矛盾。

  • 如果存在前向边 u → v u \to v uv ,那么在首次访问 v v v 时,该边就会作为一条回边指向 u u u ,不会作为 u → v u \to v uv 一条前向边。


2 Tarjan 算法

2.1 桥

题目

2.1.1 桥的定义

无向图中,若删去一条边会使得这个图的极大连通分量数增加,则该边被称为桥。

也可以理解为无向图的一个连通块中,若删除一条边会使得至少两点之间无法相互到达,该边被称为桥。

2.1.2 实现

样例:
在这里插入图片描述
样例中, 1 → 2 1 \to 2 12 5 → 6 5 \to 6 56 就是该图的桥。

  • 容易发现,只有 树边 才可能作为桥,而 dfs 生成树中红色的回边 4 → 2 4 \to 2 42 标记了 2 → 3 → 5 → 4 2 \to 3 \to 5 \to 4 2354 沿途上的树边都不可能作为桥。因此我们得出结论:一条边是桥,当且仅当它是 不被回边标记的树边

在 Tarjan 算法中,我们对每个点 u u u 按 dfs 序打上时间戳 d f n u dfn_u dfnu,同时记录 l o w u low_u lowu 表示 u u u 不经过其父亲能到达的最小时间戳,即按照 dfs 生成树中有向边能到达的最小时间戳。

那么一条树边 u → v u \to v uv 不被回边标记的充要条件即为: l o w v > d f n u low_v > dfn_u lowv>dfnu。如例中, 1 → 2 1 \to 2 12low[2]=2,dfn[1]=1,可以得到 1 → 2 1 \to 2 12 就是一座桥。

void Tarjan(int u, int fa)
{
	dfn[u] = low[u] = ++ idx;
	for(auto v : G[u])
	{
		if(v == fa) continue;
		
		if(!dfn[v]) // tree edge 
		{
			Tarjan(v, u);
			low[u] = min(low[u], low[v]);
			if(low[v] > dfn[u]) // is bridge
				bridge[u][v] = 1;
		}
		
		else // back edge 
			low[u] = min(low[u], dfn[v]);
	}
}

2.2 割点

题目

2.2.1 割点的定义

无向图中,若删去一个点会使得这个图的极大连通分量数增加,这个点被称为割点。

也可以理解为无向图的一个连通块中,若删除一个点会使得至少两点之间无法相互到达,该点被称为割点。

2.2.2 实现

在这里插入图片描述
首先,由于图不连通,我们大体的决策是将图分成若干连通块,在每个连通块对应的 dfs 生成树上进行搜索。

  • 对于结点 5 5 5,删去 5 5 5 显然不会对儿子结点 3 , 4 3,4 3,4 造成影响,因为它们有回边可以连到结点 5 5 5 的祖先,而对于儿子结点 6 6 6 会有影响,其没有回边可以连到可以连到结点 5 5 5 的祖先。

    我们推广到一般情况:对于结点 u u u ,若存在儿子结点的 l o w v ≥ d f n u low_v \ge dfn_u lowvdfnu ,即 v v v 没有回边可以连到结点 u u u 的祖先,那么 u u u 一定是个割点。

  • 然而 1 1 1 号虽然满足上述条件,然而 1 1 1 并不是割点。

    因此我们还需判断一种情况,即 u u u 作为了 dfs 搜索树的根节点,不论如何其任何儿子 v v v l o w v low_v lowv 都不可能小于 d f n u dfn_u dfnu ,因此需要分开讨论:当 u u u 作为根节点时,若 u u u 有两个及以上的儿子时, u u u 为割点。

int dfn[N], low[N], idx, cnt;
bool cut[N];
void Tarjan(int u, int fa)
{
	dfn[u] = low[u] = ++ idx;
	
	int son = 0;
	
	for(auto v : G[u])
	{	
		if(v == fa) continue;
		
		if(!dfn[v])
		{	
			Tarjan(v, u);
			low[u] = min(low[u], low[v]);
			
			son ++ ;
			if(fa != 0 and low[v] >= dfn[u] and !cut[u]) cnt ++ , cut[u] = 1;
		}
		
		else low[u] = min(low[u], dfn[v]);
	}
	
	if(fa == 0 and son >= 2 and !cut[u]) cnt ++ , cut[u] = 1;
	
}

代码

2.3 强连通分量

题目

2.3.1 强连通分量的定义

注意,强连通分量(SCC,Strongly Connected Components)是在单向图中的。

强联通子图,定义为:在 单向图 中,该子图上的任意两点之间能互相到达。

强联通分量,定义为:极大连通子图

2.3.2 实现

以下图为例:

在这里插入图片描述

该例中强连通分量为 { 1 , 2 , 3 , 4 , 5 , 6 , 7 } \lbrace 1,2,3,4,5,6,7 \rbrace {1,2,3,4,5,6,7}

前向边没用。

  • u u u 是当前 scc 第一个被访问的结点,那么该 scc 的结点一定在 dfs 生成树中以 u u u 为根的子树内。

    证明:若存在该 scc 中的点 v v v 不在 u u u 子树内,则 u u u v v v 的路径中一定存在叉边或回边,根据叉边、回边的定义,则 v v v u u u 之前已经被访问过,与 u u u 是当前 scc 第一个被访问的结点矛盾。

  • 如果我们把单个节点也算作是强连通分量,那么所有 scc 的根应 u u u 当都满足 l o w u = d f n u low_u=dfn_u lowu=dfnu

具体实现

考虑使用一个栈 s s s 存储遍历到的结点,每当访问到一个 l o w u = d f n u low_u=dfn_u lowu=dfnu 的结点,就将栈内元素从栈顶 弹出直至 u u u ,这些结点即为以 u u u 为根的强连通分量。

在遍历 u u u 的子树过程中,确实有可能会遍历到其他强连通分量的点,根据树的结构,一定会先遍历到其他 scc 的根 v v v,而在遍历 v v v 的过程中,发现了 l o w v = d f n v low_v=dfn_v lowv=dfnv ,便会将栈内的元素弹出直到 v v v,因此在遍历完 u u u 的子树后仍然留在栈中的元素一定是以 u u u 为根的 scc 中的元素。

而在搜索过程中,我们会遍历到以下三种结点:

  • 未访问的结点:该边为树边,将 v v v 入栈,对 v v v 进行 dfs,更新 l o w low low

  • 访问过并且不在栈中:该边为叉边或前向边,说明已经是其他 scc 中的结点,对当前 scc 不会产生任何影响,忽略。

  • 访问过并且在栈中:该边为叉边、前向边或回边,更新 l o w low low 值。

int low[N], dfn[N], idx;
bool ins[N];
int tot, belong[N];

stack <int> s;

void Tarjan(int u)
{
	low[u] = dfn[u] = ++ idx;
	
	s.push(u); ins[u] = 1;
	
	for(auto v : G[u])
	{
		if(!dfn[v])
		{
			Tarjan(v);
			low[u] = min(low[u], low[v]);
		}
		if(dfn[v] and ins[v]) 
			low[u] = min(low[u], dfn[v]);
	}
	
	if(low[u] == dfn[u]) // 当前 scc 的根
	{
		++ tot;
		while(s.top() != u) // 出栈,标记结点所属 scc
		{
			int x = s.top(); s.pop();
			belong[x] = tot, ins[x] = 0;
		}
		s.pop(); belong[u] = tot, ins[u] = 0;
	}
	
}

代码

2.3.3 强连通分量缩点

在无向图的一些问题中,若将所有 scc 合并为一个点,那么这张无向图就会变成一张有向无环图(DAG,Directed Acyclic Graph)。

因此我们只需保留连接两个不同 scc 的边来构成缩点后的新图。

	for(int i=1; i<=n; i++)
	{
		if(!dfn[i])
			idx = 0, Tarjan(i);
	}
	
	for(int u=1; u<=n; u++) // 缩点后的图 e
	{
		for(auto v : G[u])
		{
			int x = belong[u], y = belong[v];
			if(x != y)
				e[x].push_back(y);
		}
	}

[模板] 缩点

题目

由于题目存在环且能多次经过同一点,直接在图上进行搜索显然不可行。

由于点权均为非负数,在经过一个环时,显然应当经过还上的所有点,此时我们考虑缩点,将图变成一张 DAG ,然后用记忆化搜索或者拓扑排序求解最长路径即可。

2.4 双联通分量

无向图中,双联通分量包含点 点双联通分量边双联通分量

2.4.1 双联通分量的定义
  • 点双联通:对于点 u u u v v v ,删去图上任意一个点, u u u v v v 始终连通,则称 u u u v v v 点双联通。

    容易发现,若一个子图点双连通,那么这个子图上一定不存在割点。

  • 边双连通:对于点 u u u v v v ,仅删去图上任意一条边, u u u v v v 始终连通,则称 u u u v v v 双联通。

    容易发现,若一个子图边双连通,那么这个子图上一定不存在桥。

  • 点双连通分量:极大点双连通子图。

  • 边双连通分量:极大边双连通子图。

2.4.2 边双的实现

题目

在找出原图上的桥后,将其在原图中去掉,剩下的连通块即为边双连通分量。

int cnt; 
bool vis[N];
vector <int> dcc[N];
void dfs(int u)
{
	vis[u] = 1;
	dcc[cnt].push_back(u);
	for(int i=head[u]; i; i=e[i].nxt)
	{
		if(vis[e[i].v] or bridge[i]) continue;
		dfs(e[i].v);
	}
}

代码

2.4.3 边双缩点

我们将每个边双看作一个新节点,原图上的桥看作连接边双之间的边,这样新图会构成一棵树。

for(int i=2; i<=tot1; i++)
	{
		int u = e1[i].v, v = e1[i ^ 1].v;
		int x = belong[u], y = belong[v];
		if(x != y) add_edge2(x, y);
	}
[NOIP2022] 建造军营

题目

在边双中的任意两点之间都不需要保护道路,考虑边双缩点。由于原图为连通图,边双缩点后变为一棵树,进行树 dp 统计答案。

2.4.4 点双的实现

详见圆方树章节。

  • 6
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
Tarjan算法和Kosaraju算法都是求解有向图强连通分量的算法,它们的时间复杂度都为O(N+M),其中N为图中节点数,M为图中边数。 Tarjan算法的基本思想是通过DFS遍历图中的节点,并在遍历的过程中维护一个栈,用于存储已经遍历过的节点。在遍历的过程中,对于每个节点,记录它被遍历到的时间戳和能够到达的最小时间戳,当一个节点的最小时间戳等于它自身的时间戳时,说明这个节点及其之前遍历到的节点构成了一个强连通分量,将这些节点从栈中弹出即可。 Kosaraju算法的基本思想是先对原图进行一次DFS,得到一个反向图,然后再对反向图进行DFS。在第二次DFS的过程中,每次从未被访问过的节点开始遍历,遍历到的所有节点构成一个强连通分量。 两种算法的具体实现可以参考以下代码: ```python # Tarjan算法 def tarjan(u): dfn[u] = low[u] = timestamp timestamp += 1 stk.append(u) for v in graph[u]: if not dfn[v]: tarjan(v) low[u] = min(low[u], low[v]) elif v in stk: low[u] = min(low[u], dfn[v]) if dfn[u] == low[u]: scc = [] while True: v = stk.pop() scc.append(v) if v == u: break scc_list.append(scc) # Kosaraju算法 def dfs1(u): vis[u] = True for v in graph[u]: if not vis[v]: dfs1(v) stk.append(u) def dfs2(u): vis[u] = True scc.append(u) for v in reverse_graph[u]: if not vis[v]: dfs2(v) # 构建图和反向图 graph = [[] for _ in range(n)] reverse_graph = [[] for _ in range(n)] for u, v in edges: graph[u].append(v) reverse_graph[v].append(u) # Tarjan算法求解强连通分量 dfn = [0] * n low = [0] * n timestamp = 1 stk = [] scc_list = [] for i in range(n): if not dfn[i]: tarjan(i) # Kosaraju算法求解强连通分量 vis = [False] * n stk = [] scc_list = [] for i in range(n): if not vis[i]: dfs1(i) vis = [False] * n while stk: u = stk.pop() if not vis[u]: scc = [] dfs2(u) scc_list.append(scc) ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值