Archipelago

Lec 13最小割 Min-Cut 近似算法

讲义完成更新于 2026-09-17

Lec 13 最小割 Min-Cut 近似算法

对应 DGA 第 2.3 节 · 延伸阅读:Ghaffari–Kuhn (DISC 2013)、Ghaffari–Haeupler (SODA 2016)

TL;DR

  • 在 CONGEST 模型中,最小割问题可通过构造连通度为 k′k' 的稀疏证书将网络边数压缩至 O(nk′)O(nk')。
  • 利用人工权重迭代调用分布式 MST 子程序,有效规避了在残差图上直接求森林导致直径失控的难题。
  • 结合稀疏证书与边收缩两情形递归,在 O(k(D+n)log⁡3n)O(k(D+\sqrt n)\log^3 n) 轮内给出了 (2+ε)(2+\varepsilon)-近似最小割,并在平面图上借助低拥塞捷径进一步加速至 O~(D)\tilde O(D)。

1. 本讲要回答的核心问题

最小割(minimum cut) kk 是使图断开所需删除的最少边数。本讲在 CONGEST 模型给出一个 (2+ε)(2+\varepsilon)-近似算法:输出非空顶点子集 SS,使 SS 到 V∖SV\setminus S 的边数 ≤(2+ε)k\le(2+\varepsilon)k。

主结果. 对任意常数 ε>0\varepsilon>0,可在 O ⁣(k(D+n)log⁡3n)O\!\big(k(D+\sqrt n)\log^3 n\big) 轮内算出 (2+ε)(2+\varepsilon)-近似最小割(kk 为最小割大小)。

核心工具是连通性的稀疏证书(sparse certificate),它复用了 Lec 12 的 MST 子程序。

2. 连通性的稀疏证书

直觉:在不改变最小割大小的前提下,把边数压到接近下限 nk/2nk/2 的子图。

连通度

k′k' 的稀疏证书 生成子图 H=(V,EH)H=(V,E_H) 称为 GG 关于连通度 k′k' 的稀疏证书,若:

  1. 对每个非空 S⊂VS\subset V,cutH(S,V∖S)≥min⁡{k′, cutG(S,V∖S)}\mathrm{cut}_H(S,V\setminus S)\ge\min\{k',\,\mathrm{cut}_G(S,V\setminus S)\};
  2. HH 至多 nk′nk' 条边。 :::

连通图 GG 的任意生成树 TT 满足条件 2(n−1n-1 条边)与条件 1(cutT(S)≥1\mathrm{cut}_T(S)\ge1)。更一般地,不连通时任意极大生成森林都是 k′=1k'=1 的稀疏证书。

构造(集中式轮廓). 迭代 k′k' 次:第 ii 次计算当前图的极大生成森林 FiF_i,再令 G←G∖FiG\leftarrow G\setminus F_i。最终 H=⋃i=1k′FiH=\bigcup_{i=1}^{k'}F_i 即所求证书。

引理 2.5. H=⋃iFiH=\bigcup_i F_i 是连通度 k′k' 的稀疏证书。

分布式实现(避免在残图上反复算森林,因残图直径可能很大). 用人工边权:初始全设权 11;第 ii 次算最小权生成树 TiT_i,取其中权为 11 的边作 FiF_i,然后把这些边权改成 ∞\infty(实际取 n2n^2)。

引理 2.6. 在“G∖(∪j<iFj)G\setminus(\cup_{j<i}F_j) 权为 11、其余权为 n2n^2”的加权图上,TiT_i 中权为 11 的边恰构成 G∖(∪j<iFj)G\setminus(\cup_{j<i}F_j) 的极大生成森林。

每个 TiT_i 用 Lec 12 的 MST 算法在 O((D+n)log⁡n)O((D+\sqrt n)\log n) 轮算出,k′k' 次共 O(k′(D+n)log⁡n)O(k'(D+\sqrt n)\log n) 轮。

3. (2+ε)(2+\varepsilon)-近似算法

先假设最小割大小 kk 已知。算法每“层”如下:

  1. 计算连通度 k′=(1+ε/10)kk'=(1+\varepsilon/10)k 的稀疏证书 HH(O(k(D+n)log⁡n)O(k(D+\sqrt n)\log n) 轮)。
  2. 考虑收缩 G∖HG\setminus H 的所有边所得辅助图 H′H'。两条性质:(A) H′H' 的边恰为 HH 的边,故 H′H' 至多 nk′=(1+ε/10)nknk'=(1+\varepsilon/10)nk 条边;(B) H′H' 的最小割 =G=G 的最小割 kk。

按 H′H' 节点数分两情形:

情形 1(度数界与近似割)

若 H′H' 有 ≥n/(1+ε/5)\ge n/(1+\varepsilon/5) 个节点,则平均度 ≤2(1+ε/10)nkn/(1+ε/5)≤2k(1+ε/10)(1+ε/5)≤(2+ε)k,\le\frac{2(1+\varepsilon/10)nk}{n/(1+\varepsilon/5)}\le 2k(1+\varepsilon/10)(1+\varepsilon/5)\le(2+\varepsilon)k, 故存在一个度 ≤(2+ε)k\le(2+\varepsilon)k 的节点。它是 GG 中某些被收缩到一起的顶点子集 SS,于是 cutG(S)≤(2+ε)k\mathrm{cut}_G(S)\le(2+\varepsilon)k——这就是要找的近似割。

情形 2(H′H' 有 <n/(1+ε/5)<n/(1+\varepsilon/5) 个节点): 难以直接找到小割,但节点数已缩小了 (1+ε/5)(1+\varepsilon/5) 倍而最小割不变。令 G←H′G\leftarrow H' 递归。O(log⁡1+ε/5n)O(\log_{1+\varepsilon/5} n) 层后节点数降到 ≤2\le2,必能找到割。

去掉“已知 kk”的假设: 对 11 到 n2n^2 之间所有形如 (1+ε/10)i(1+\varepsilon/10)^i 的估计值各跑一遍,输出所有结果中最小的割。

4. 几点说明与拓展

  • 去掉 kk 因子: 用 Karger 式采样把最小割降到 O(log⁡n)O(\log n),可把轮复杂度里的 kk 换成 O(log⁡n)O(\log n),得 O~(D+n)\tilde O(D+\sqrt n)。
  • 加权图: 把权 ww 边视为 ww 条平行边即可推广。
  • 近最优性: O~(D+n)\tilde O(D+\sqrt n) 已近最优——任何非平凡最小割近似都需 Ω~(D+n)\tilde\Omega(D+\sqrt n) 轮(下界见 Lec 14)。

5. 平面图:低拥塞捷径下的 O~(D)\tilde O(D)

Ghaffari–Haeupler (SODA 2016) 把 Lec 12 的低拥塞捷径框架系统化,证明平面图(及有界亏格、排除稠密 minor 等)总存在 max⁡{α,β}≤O~(D)\max\{\alpha,\beta\}\le\tilde O(D) 的捷径且 O~(D)\tilde O(D) 轮可构造。由此在平面网络上同时得到 O~(D)\tilde O(D)-轮的 MST 与 Min-Cut——突破了一般图 Ω~(n)\tilde\Omega(\sqrt n) 的壁垒。

6. 本讲小结与近似割体系

内容结论
稀疏证书保最小割、边数 ≤nk′\le nk';生成树为 k′=1k'{=}1 证书
分布式构造加权 MST 技巧,O(k′(D+n)log⁡n)O(k'(D+\sqrt n)\log n) 轮
(2+ε)(2+\varepsilon)-近似证书 + 收缩 + 两情形递归,O(k(D+n)log⁡3n)O(k(D+\sqrt n)\log^3 n)
去 kk / 改进采样 ⇒ O~(D+n)\tilde O(D+\sqrt n)
平面图低拥塞捷径 ⇒ O~(D)\tilde O(D) MST 与 Min-Cut

Min-Cut 与 MST 共享同一套机制(MST 子程序、稀疏证书、低拥塞捷径),并同样止步于 O~(D+n)\tilde O(D+\sqrt n)——下一讲证明这是全局问题的硬下界。