1 回 グラフの定義:数理 岡本吉央 [email protected]/okamotoy/lect/2019/gn/lect01.pdfグラフとネットワーク 第1 回 グラフの定義:数理
2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ...
Transcript of 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ...
![Page 1: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/1.jpg)
1
一般グラフのマッチングとその周辺
高澤兼二郎
京都大学数理解析研究所
第 11 回組合せ最適化セミナー2014 年 7 月 30 日
![Page 2: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/2.jpg)
今日の目標 2
Goemans 2013
(1) 最大マッチングをひとつ示せ
(2)それが最大マッチングであることを示せ(= それよりも大きいマッチングが存在しないことを示せ)
![Page 3: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/3.jpg)
目次 3
一般グラフのマッチング
一般グラフの重みつきマッチング
2部グラフのマッチング
Tutte-Berge 公式 Edmonds のアルゴリズム Edmonds-Gallai 分解
周辺の話題
Kőnig の定理 増加道アルゴリズム
2部グラフの重みつきマッチング マッチング多面体 完全単模性
マッチング多面体 完全双対整数性
(Cunningham-Marsh公式)
![Page 4: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/4.jpg)
|M| 最大のマッチング M を求めよ
M⊆E がマッチング
Mの枝は端点を共有しない
マッチング問題 4
多項式時間可解 Edmonds 1965 など
無向グラフ G = (V, E)頂点集合 枝集合
6 12
11
5
2
42
1
8
3
枝重み w: ER+
w(M) = 2+8+5 = 15
1
定義
問題
w(M) = ∑e∈Mw(e) 最大のマッチング M を求めよ
![Page 5: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/5.jpg)
交互道 5
Mに関する交互道: M と E-M の枝が交互に現れる道
Mに関する増加道: 交互道で,両端点がMに接続しない定義
P が M-増加道 M’ = M∆P = (M - P)∪(P - M)はマッチングで|M’| = |M|+1
M P MΔP
![Page 6: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/6.jpg)
増加道アルゴリズム (?) 6
M が最大マッチング M-増加道が存在しない定理
【演習】証明せよ
M: 任意のマッチング
M-増加道 P を探す
Pが見つかった
Pが存在しない
M := M∆P
Return M
アルゴリズム(?)
![Page 7: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/7.jpg)
存在しないことの証明 7
「○○が存在しない」ことを説明するには・・・??
(^_^;) 「探したけど増加道ありません」(゚⊿゚) 「は? ちゃんと探せや」
(^O^) 「全探索したけどありません!(u v の全頂点列ドバー)」
(゚⊿゚) 「こんなにたくさんチェックできん, やり直し」
u
v
M が最大マッチング M-増加道が存在しない定理
![Page 8: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/8.jpg)
2部グラフの場合
E – M の枝: UV の向き付け M の枝: VU へ向き付け
M-交互道 有向道 M-増加道 U°-V° 有向道
G = (U∪V, E)
U°
V°
M-増加道が存在 U°からV°へ到達可能 R
R := U°から到達可能な頂点の集合R∩V° = Ø M-増加道は存在しない
M は最大マッチング
![Page 9: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/9.jpg)
2部グラフのアルゴリズム
E – M の枝: UV の向き付け M の枝: VU へ向き付け
G = (U∪V, E)
(*)
U°
V°
RM: 任意のマッチング
U°-V° 有向道 P を探す
Pが見つかった
Pが存在しない
M := M∆P
Return M
E の枝に, (*) の向き付け
アルゴリズム
![Page 10: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/10.jpg)
マッチング と 点被覆 10
C⊆U∪V が点被覆
E のすべての枝は C に端点をもつ
定義
M はマッチング, C は点被覆 |M| ≤ |C|
命題
|M| = |C| が成立すればM は最大マッチングC は最小点被覆
点被覆 C |C| はマッチングサイズの上限
![Page 11: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/11.jpg)
C は点被覆 |M| = |C|
最大マッチング と 最小点被覆11
Rアルゴリズム終了時: U°-V° 有向道なし
R⊆U∪V : U° から有向道で到達可能C := (U - R)∪(V∩R)
命題
max{|M| : M はマッチング} = min{|C| : C は点被覆}
定理 [Kőnig 1931]
![Page 12: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/12.jpg)
最小点被覆全体の構造: DM分解 12
点被覆 C1 点被覆 C2
C1, C2 はともに最小点被覆
C⊆U∪V が最小点被覆 ∈C・・・
Dulmage-Mendelsohn 分解
浅野: 情報の構造 上. 日本評論社, 1994 伊理, 藤重, 大山: グラフ・ネットワーク・マトロイド, 産業図書, 1986 K. Murota: Matrices and Matroids for Systems Analysis,
Springer, 2010
![Page 13: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/13.jpg)
完全マッチングの存在判定 13
Γ(X) := { v∈V : ∃u∈X, uv∈E}
G = (U∪V, E), X⊆U
G が完全マッチングをもつ |Γ(X)| ≥ |X| ∀X⊆U
定理 [Hall 1935]
X Γ(X)
【演習】証明せよ
![Page 14: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/14.jpg)
ここまでのまとめ 14
最大マッチング 増加道が存在しない
存在しないことの証明は (一般には) 難しい
2部グラフ: 向き付け + パス探索で存在判定可能
最大最小定理による最適性の証明
max{|M| : M はマッチング} = min{|C| : C は点被覆}
max ≤ min を示す |M| = |C| をみたす M, C 求めるアルゴリズムを設計
max = min の証明 & アルゴリズムの正しさの保証
![Page 15: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/15.jpg)
目次 15
一般グラフのマッチング
一般グラフの重みつきマッチング
2部グラフのマッチング
Tutte-Berge 公式 Edmonds のアルゴリズム Edmonds-Gallai 分解
周辺の話題
Kőnig の定理 増加道アルゴリズム
2部グラフの重みつきマッチング マッチング多面体 完全単模性
マッチング多面体 完全双対整数性
(Cunningham-Marsh公式)
![Page 16: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/16.jpg)
一般グラフのマッチング 16
M が最大マッチング M-増加道が存在しない定理
一般グラフでも成立
増加道の探索・存在判定が難しい
1 2 3
4 5
6
7 8
• 向き付けができない• 各枝はどちら向きにも通れる
u
v
探索: 1234567 … 増加道なし??
![Page 17: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/17.jpg)
最大最小定理 17
Tutte 1947, Berge 1958
X⊆V
12
odd(V-X)=4
X
max {|M| : M はマッチング}
= min {|V| + |X| - odd(V - X)}
定理 [Tutte-Berge 公式]
頂点数奇数の成分数
G が完全マッチングをもつodd(V - X) ≤ |X| ∀X⊆V
M, X は最適解
定理 [Tutte の定理]
max ≤ min
等号でみたす M, X
![Page 18: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/18.jpg)
冒頭の問題 の 答え 18
X = {2, 15, 16} odd(V-X) = 5
18
1 17
4
6 9 7
10 11 12
13 14
2 15 16
3
5 8
X
∀M, |M| ≤ (|V| + |X| - odd(V-X))/2 = (18 + 3 – 5)/2 = 8
|M|=8 のマッチングMが見つかっていれば最適解(Xも最適解)
![Page 19: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/19.jpg)
1 2 3
4 5
6
87
ブロッサム・アルゴリズム: 縮約・増加 19
Edmonds 1965
1 2 3
4 5
6
7 8
1 2 83~7
1 2 83~7
∀奇閉路 C は factor-critical命題
∀v, C – v は完全マッチングをもつ
![Page 20: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/20.jpg)
ブロッサム・アルゴリズム: 終了時 20
X = { から奇数長の交互道でたどり着く頂点}
|M| = 1
2(|V| + |X| - odd(V-X))
Edmonds 1965
M, X はそれぞれ最適解Tutte-Berge 公式の構成的証明
![Page 21: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/21.jpg)
Edmonds-Gallai 分解 21
D = {v∈V : v を被覆しない最大マッチングが存在}A = {v∈V-D : vはDの頂点に隣接}C = V – D - A
(1) G[D] の連結成分は factor-critical
(2) G[C] は完全マッチングをもつ
(3) C を削除, E[A] を削除, G[D] の連結成分を縮約して得られる2部グラフにおいて, Γ(X) - |X| ≥ 1 ∀X⊆A
(4) M: G の最大マッチング
• M[D]: 各連結成分の near-perfect matching• M[C]: G[C] の完全マッチング• M[A]: A の頂点は D の異なる連結成分とマッチ
(5) max{|M|} = (|V| - #component(D) + |A|)/2
A
D
C
D
A
定理 [Edmonds-Gallai 分解]
![Page 22: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/22.jpg)
ここまでのまとめ 22
最大マッチング 増加道が存在しない
一般グラフ: 増加道の存在判定は難しい
最大最小定理による最適性の証明
max{|M| : M はマッチング} = 1
2min{|V| + |X| - odd (V-X)}
max ≤ min を示す 左辺 = 右辺をみたす M, X 求めるアルゴリズムを設計
max = min の証明 & アルゴリズムの正しさの保証
![Page 23: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/23.jpg)
目次 23
一般グラフのマッチング
一般グラフの重みつきマッチング
2部グラフのマッチング
Tutte-Berge 公式 Edmonds のアルゴリズム Edmonds-Gallai 分解
周辺の話題
Kőnig の定理 増加道アルゴリズム
2部グラフの重みつきマッチング マッチング多面体 完全単模性
マッチング多面体 完全双対整数性
(Cunningham-Marsh公式)
![Page 24: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/24.jpg)
24
G = (U∪V, E) w : E R≥0
w(M) 最大のマッチング M w(M) 最小の完全マッチング M
問題
多項式時間可解 (Hungarian method など)
3
1
2
1
2
1
12
w(M) = ∑e∈Mw(e)
2部グラフの重みつきマッチング
![Page 25: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/25.jpg)
最小重み完全マッチングの整数計画表現 25
x(e) = 1 e∈Mx(e) = 0 e∈E - M
3
122
1
2変数 x∈{0,1}E : マッチングの特性ベクトル
x = (1 0 1 0 0 1)
minimize ∑e∈E w(e) x(e)
subject to ∑e∈δ(v) x(e) = 1 ∀v∈U∪Vx(e) ∈ {0,1} ∀e∈E
IPmin
w(M) = ∑e∈Mw(e)x(e)
vδ(v)
![Page 26: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/26.jpg)
線形計画緩和 と 整数最適解 26
minimize ∑e∈E w(e) x(e)
subject to ∑e∈δ(v) x(e) = 1 ∀v∈U∪Vx(e) ≥ 0 ∀e∈E
LPmin
OPT(LPmin) ≤ OPT(IPmin)
G が 2部グラフのとき,(LPmin) は整数最適解 x* をもつ
定理 [Birkhoff 1946 など]
(LPmin) を解いて整数最適解 x* を求める x* は (IPmin) の最適解
(最小重み完全マッチングの特性ベクトル)
![Page 27: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/27.jpg)
整数性の証明 27
G が 2部グラフのとき,(LPmin) の制約式は完全マッチング多面体を定める
定理 [Birkhoff 1946 など]
【証明】
• x : 制約式が定める多面体 Q の端点
∑e∈δ(v) x(e) = 1 ∀v∈U∪Vx(e) ≥ 0 ∀e∈E
• F = {e∈E | x(e) > 0}
• F が閉路 C を含むとする
C の長さは偶数
C はマッチング M, N の disjoint union
x + ε(χM - χN), x - ε(χM - χN) ∈ Q
x が Q の端点であることに矛盾
【証明終】
• F は閉路を含まない x(e) = 1 ∀e∈F
![Page 28: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/28.jpg)
最大重みマッチングの線形計画表現 28
maximize ∑e∈E w(e) x(e)
subject to ∑e∈δ(v) x(e) ≤ 1 ∀v∈U∪Vx(e) ≥ 0 ∀e∈E
LP
G が 2部グラフのとき,(LP) は整数最適解 x* をもつ, i.e., (LP) の制約式はマッチング多面体を定める
定理
【演習】証明せよHint: (LPmin) の整数性を使ってよい
![Page 29: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/29.jpg)
LP の最適性保証 : 双対定理 29
maximize ∑e∈E w(e) x(e)subject to ∑e∈δ(v) x(e) ≤ 1 ∀v∈U∪V
x(e) ≥ 0 ∀e∈E
LP
minimize ∑v∈U∪V y(v)
subject to y(u) + y(v) ≥ w(e) ∀e = uv∈Ey(v) ≥ 0 ∀v∈U∪V
Dual-LP y∈RU∪V
OPT(LP) = OPT(Dual-LP)
(LP) の許容解 x が最適解
相補性条件をみたす (Dual-LP) の許容解 y が存在: • x(e) > 0 y(u) + y(v) = w(e)• y(v) > 0 ∑e∈δ(v) x(e) = 1
![Page 30: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/30.jpg)
双対最適解 = “最小点被覆” 30
minimize ∑v∈U∪V y(v)
subject to y(u) + y(v) ≥ 1 ∀e = uv∈Ey(v) ≥ 0 ∀v∈U∪V
Dual-LP
y∈{0,1}U∪V y は点被覆の特性ベクトル
最小点被覆問題の線形緩和
y( ) = 1y( ) = 0
w(e) = 1 ∀e の場合
(Dual-LP) は整数最適解をもつ命題
相補性条件
x: 最大マッチングy: 最小点被覆• y(v) = 1 ∑e∈δ(v) x(e) = 1
• x(uv) = 1 y(u) + y(v) = 1
【略証】• OPT(IP) = OPT(LP) = OPT(Dual-LP) ≤ OPT(Dual-IP)• OPT(IP) = OPT(Dual-IP)
![Page 31: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/31.jpg)
双対最適解 = “最小点被覆” 31
minimize ∑v∈U∪V y(v)
subject to y(u) + y(v) ≥ w(e) ∀e = uv∈Ey(v) ≥ 0 ∀v∈U∪V
Dual-LP
最小点被覆問題の一般化
w(e)∈Z ∀e の場合
w が整数ならば, (Dual-LP) は整数最適解をもつ
定理[Egerváry 1931]
相補性条件
x: 最大重みマッチングy: “最小点被覆”• y(v)∈Z>0 ∑e∈δ(v) x(e) = 1
• x(uv) = 1 y(u) + y(v) = w(e)
4
121
1
2
0 1
1 0
22
![Page 32: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/32.jpg)
係数行列 の 完全単模性 32Total unimodularityTotally unimodular
A : n×m 行列
A が 完全単模 (完全単模性をもつ)
A の任意の正方小行列の行列式が 0, +1 or -1
定義
A が完全単模 n×m 行列, b∈Zn
P = {x | Ax ≤ b} は整数多面体定理
整数格子点の凸包
A が完全単模 n×m 行列, b∈Zn, c∈Zm
max {cx | Ax ≤ b, x ≥ 0} = min {yb | yA ≥ c, y ≥ 0} 両辺は 整数最適解 をもつ(※両辺は実行可能とする)
定理
![Page 33: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/33.jpg)
接続行列の完全単模性 33
A(G) が完全単模 G は2部グラフ
定理
G = (V,E)A(G): G の接続行列 e1
e2
e5e3
e4v1
v2 v3
v4
e1 e2 e3 e4 e5
v1
v2
v3
v4
1100
0110
0011
1001
1010
maximize ∑e∈E w(e) x(e)
subject to ∑e∈δ(v) x(e) ≤ 1 ∀v∈U∪Vx(e) ≥ 0 ∀e∈E
LP
A(G) x ≤ 1
【演習】証明せよ
![Page 34: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/34.jpg)
ここまでのまとめ 34
最大重みマッチング
min. 1ys.t. yA ≥ w,
y ≥ 0
max. wxs.t. Ax ≤ 1,
x∈{0,1}E
(IP)
max. wxs.t. Ax ≤ 1,
x ≥ 0
(LP) (LP-Dual)
G: 2部グラフ
min. 1ys.t. yA ≥ w,
y∈ZU∪V
(IP-Dual) “最小点被覆“
線形緩和
=(双対定理)
≦≦= =
A の完全単模性
=
![Page 35: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/35.jpg)
目次 35
一般グラフのマッチング
一般グラフの重みつきマッチング
2部グラフのマッチング
Tutte-Berge 公式 Edmonds のアルゴリズム Edmonds-Gallai 分解
周辺の話題
Kőnig の定理 増加道アルゴリズム
2部グラフの重みつきマッチング マッチング多面体 完全単模性
マッチング多面体 完全双対整数性
(Cunningham-Marsh公式)
![Page 36: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/36.jpg)
36
min. ∑ w(e) x(e)
sub. to ∑e∈δ(v) x(e) = 1 ∀v∈Vx(e) ∈ {0,1} ∀e∈E
OPT(IPmin) > OPT(×LPmin)
e1 e2
e3
(IP): 許容解なし
(LP): x = (½ ½ ½)
(IPmin)
min. ∑ w(e) x(e)
sub. to ∑e∈δ(v) x(e) = 1 ∀v∈Vx(e) ≥ 0 ∀e∈E
(×LPmin)
最小重み完全マッチング
![Page 37: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/37.jpg)
min. ∑ w(e) x(e)
sub. to ∑e∈δ(v) x(e) = 1 ∀v∈V
∑e∈δ(U) x(e) ≥ 1 ∀U⊆V, |U| is odd
x(e) ≥ 0 ∀e∈E
U
最小重み完全マッチングの線形計画表現 37
(LPmin)
(LPmin) は整数最適解 x* をもつi.e., (LPmin) の制約式は整数多面体を定める
定理 [Edmonds 1965]
U
![Page 38: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/38.jpg)
整数性の証明 38
(LPmin) の制約式は完全マッチング多面体を定める
定理 [Edmonds 1965]
∑e∈δ(v) x(e) = 1 ∀v∈V∑e∈δ(U) x(e) ≥ 1 ∀U⊆V, |U| is oddx(e) ≥ 0 ∀e∈E
【証明】P=完全マッチング多面体, Q=制約式が定める多面体
P⊆Q : 簡単 Q⊆P を示す.x ∈ P をみたす Q の頂点 x が存在したとする
• |V|+|E| 最小の例をとる
0 < x(e) < 1 ∀e∈E deg(v) ≥ 2 ∀v∈V |E| ≥ |V|
|E| = |V| の場合閉路の disjoint union
α
1-α
完全マッチングの凸結合 二つ目の制約に矛盾
![Page 39: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/39.jpg)
整数性の証明 39∑e∈δ(v) x(e) = 1 ∀v∈V∑e∈δ(U) x(e) ≥ 1 ∀U⊆V, |U| is oddx(e) ≥ 0 ∀e∈E
|E| > |V| の場合
• x は |E| 本の制約を等号でみたす
∃U⊆V, x(δU) = 1 and 3 ≤ |U| ≤ |V| - 3
• x’, x’’ : x の G/Ū, G/U への射影
x’∈Q(G/Ū) = P(G/Ū)x’’∈Q(G/U) = P(G/U)
Q(G)
x’ = χM’(i), x’’ = χM’’(i)1
𝑘
𝑖=1
𝑘1
𝑘
𝑖=1
𝑘
∀i∈[k], M’(i)∩δU = M’’(i)∩δU としてよい
∀i∈[k], M(i) = M’(i)∪M’’(i) は G の完全マッチングで
1
𝑘
𝑖=1
𝑘
x = χM(i) となって矛盾 【証明終】
U Ū
ŪU
![Page 40: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/40.jpg)
max. ∑ w(e) x(e)
sub. to ∑e∈δ(v) x(e) ≤ 1 ∀v∈V
∑e∈E[U] x(e) ≤𝑈 −1
2∀U⊆V, |U| is odd
x(e) ≥ 0 ∀e∈E
U
最大重みマッチングの線形計画表現 40
(LP)
(LP) は整数最適解 x* をもつi.e., (LP) の制約式は整数多面体を定める
定理 [Edmonds 1965]U
【演習】証明せよ
![Page 41: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/41.jpg)
双対問題 41
min. ∑v∈V y(v) + ∑U
𝑈 −1
2z(U)
sub. to y(u) + y(v) + ∑U: u,v∈U z(U) ≥ w(e) ∀e=uv∈Ey(v) ≥ 0 ∀v∈V z(U) ≥ 0 ∀U⊆V,|U| is odd
(Dual-LP)
U2
U1
U3
![Page 42: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/42.jpg)
完全双対整数性 42y(u) + y(v) + ∑U: u,v∈U z(U) ≥ w(e) y(v) ≥ 0z(U) ≥ 0
w が整数のとき, (Dual-LP) は整数最適解 (y, z) をもつさらに, F = {U⊆V | z(U) > 0} がラミナー集合族をなす整数最適解が存在する
定理 [Cunningham & Marsh 1978]
U1, U2∈F U1⊆U2 or U1⊇U2 or U1∩U2=Ø
(LP) の制約式は完全双対整数性をもつ
定理
Rational system Ax ≤ b が完全双対整数性をもつmin{ yb | y ≥ 0, yA = c} が任意の整数ベクトル c に対して整数最適解をもつ
![Page 43: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/43.jpg)
ここまでのまとめ 43
max. wxs.t. Ax ≤ 1,
x(E[U]) ≤ 𝑈 −1
2
x ≥ 0
(LP)
min. ∑vy(v) + ∑U𝑈 −1
2z(U)
s.t. y(u) + y(v) + ∑U∋u,vz(U) ≥ w, y,z ≥ 0
(Dual-LP)
G: 一般グラフ
(双対定理)
=
マッチング多面体を定める 完全双対整数性をもつ
(LP), (Dual-LP) ともに整数最適解をもつ
![Page 44: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/44.jpg)
目次 44
一般グラフのマッチング
一般グラフの重みつきマッチング
2部グラフのマッチング
Tutte-Berge 公式 Edmonds のアルゴリズム Edmonds-Gallai 分解
周辺の話題
Kőnig の定理 増加道アルゴリズム
2部グラフの重みつきマッチング マッチング多面体 完全単模性
マッチング多面体 完全双対整数性
(Cunningham-Marsh公式)
![Page 45: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/45.jpg)
参考文献 45
高澤: マッチング. 薩摩順吉, 大石進一, 杉原正顕 (編),応用数理ハンドブック, 朝倉書店, 2013, pp. 288-289.
A. Schrijver: Combinatorial Optimization---Polyhedra and Efficiency, Springer, 2003.
L. Lovász, M.D. Plummer: Matching Theory, AMS Chelsea Publishing, 2009.
G.L. Nemhauser, L.A. Wolsey: Integer and Combinatorial Optimization, Wiley, 1988.
浅野: 情報の構造 上. 日本評論社, 1994.
K. Murota: Matrices and Matroids for Systems Analysis, Springer, 2010.
伊理, 藤重, 大山: グラフ・ネットワーク・マトロイド. 産業図書, 1986.
![Page 46: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/46.jpg)
46
![Page 47: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/47.jpg)
Finding 2-Factors Closer to TSP Tours in Cubic Graphs
11th COSS
July 30, 201447
Sylvia Boyd (U. Ottawa)
Satoru Iwata (U. Tokyo)
Kenjiro Takazawa (Kyoto U.)
Related Topic
![Page 48: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/48.jpg)
Petersen’s Theorem 48
Every bridgeless cubic graph has a perfect matching
Every bridgeless cubic graph has a 2-factor
[1891]
G=(V,E): Bridgeless Cubic Graph
Thm.
= 2-edge-connected deg(v) = 3 for every v in V
![Page 49: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/49.jpg)
Schönberger’s Theorem 49
G has a perfect matching including e*
G has a 2-factor excluding e*e*
[1935]
G=(V,E): Bridgeless Cubic Graphe* in E
O(n log4 n) algorithm[Biedl, Bose, Demaine, Lubiw 2001] n = |V|
Thm.
![Page 50: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/50.jpg)
Kaiser & Škrekovski’s Theorem 50
[2008]G=(V,E): Bridgeless Cubic e* in E
G has a 2-factor excluding e* and covering all 3- and 4-edge cuts
3-edge cut 4-edge cut Not a 4-edge cut
Thm.
![Page 51: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/51.jpg)
2-factors and TSP Tours51
TSP tour = 2-factor of one cycle of length n
2-factor without cycles of length k or less :C≤k-free 2-factor (in simple graphs)
✓ C≤3-free ✓C≤4-free
k = n/2 TSP tour
Relax
![Page 52: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/52.jpg)
Complexity of C≤k-free 2-factors52
Unweighted Weighted
k ≥ 5 NP-hard [Papadimitriou ’80] NP-hard
k = 4 (a) OPEN (b) NP-hard[Vornberger ’80]
k = 3 (c) P [Hartvigsen ’84] (d) OPEN
k = 2 P P
Bipartite graphs(a) : P [Hartvigsen ’06, Pap ’07]
(b) : NP-hard for general weight [Király 00]
P if the weight hass a special property [Makai ’07, T. ’09]
Subcubic graphs(a) : P [Bérczi & Végh ’10]
(c) : P [Bérczi & Végh ’10, Hartvigsen & Li ’11]
(d) : P [Vornberger ’80, Kobayashi ’10, Hartvigsen & Li ’13]
![Page 53: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/53.jpg)
2-factors Covering Cuts53
TSP tour = 2-factor covering all edge cuts
2-factor covering all 3-edge cuts C≤3-free
2-factor covering all 3,4-edge cuts C≤4-free
G: 3-edge-connected cubic
G: Cubic
2-factor covering prescribed edge cutsRelax
![Page 54: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/54.jpg)
Our Results54
(2) An O(n3)-algorithm for finding a 2-factor covering all 3-, 4-edge cuts in bridgeless cubic graphs
(3) A 6/5-approx. algorithm for the minimum 2-edge-connected subgraph problemin 3-edge-connected cubic graphs
Constructive proof for [Kaiser, Škrekovski 2008]
Start with the 2-factor found by Algorithm (2) Previous ratio: 5/4 for 3-edge-connected cubic graphs [Huh 2004]
(1) - An O(n3)-algorithm for finding a min.-weight 2-factor covering all 3-edge cuts in bridgeless cubic graphs
- Polyhedral description
Application
![Page 55: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/55.jpg)
Contents 55
Introduction
Summary
(1) - An O(n3)-algorithm for finding a min.-weight 2-factor covering all 3-edge cuts in bridgeless cubic graphs
- Polyhedral description
(1) - 2辺連結3正則グラフにおける, すべての3辺カットと丁度 1 本の辺で交わる完全マッチングで最小重みのものを求めるアルゴリズム
- 多面体的表現3-カット完全マッチング
![Page 56: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/56.jpg)
多面体的表現 56
∑e∈δ(v) x(e) = 1 ∀v∈V∑e∈δ(U) x(e) ≥ 1 ∀U⊆V, |U| is odd
x(e) ≥ 0 ∀e∈E∑e∈δ(S) x(e) = 1 ∀S⊆V, δ(S) is a proper 3-edge cut
S
A 3-edge cut δ(S) is proper2 ≤|S|≤ n - 2
V-S
上記の制約式は 3-カット完全マッチング多面体を定める
定理
【証明】
P⊆Q : 簡単
P=3-カット完全マッチング多面体, Q=制約式が定める多面体
![Page 57: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/57.jpg)
多面体的表現 57
Q⊆P を示す: x*∈Q を3-カット完全マッチングの凸結合で表す
∑e∈δ(v) x(e) = 1 ∀v∈V∑e∈δ(U) x(e) ≥ 1 ∀U⊆V, |U| is odd
x(e) ≥ 0 ∀e∈E∑e∈δ(S) x(e) = 1 ∀δ(S): proper 3-edge cut
Q
x*∈Q⊆[完全マッチング多面体]
x* = ∑i λi χMi : 完全マッチング {Mi} の凸結合
∑e∈δ(S) x*(e) = 1
∑e∈δ(S) χMi(e) ≥ 1 ∑e∈δ(S) χMi(e) = 1 ∀i
i.e., Mi は 3-カット完全マッチング ∀i 【証明終】
![Page 58: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/58.jpg)
アルゴリズム 58
(1) Find a proper 3-edge cut δ(S)
S
(2) Contract V – S, S
(3) For i = 1,2,3, find a min.-weight perfect matching Mi including ei
(4) For i = 1,2,3, let w(ei) := w(Mi)
Smallerbridgeless cubic graphs
G2
S e3
G1
e2
e1
(5) Recurse in G2 N return M = Mi∪N
M intersects all 3-cuts with 1 edge : Proven by submodularity
![Page 59: 2部グラフにおける 制約付き最小重み t-因子の組合せ的アルゴリ …ds.ws.hosei.ac.jp/slides/140730COSS.pdf · 1 一般グラフのマッチングと その周辺](https://reader035.fdocument.pub/reader035/viewer/2022070719/5edf6000ad6a402d666ab978/html5/thumbnails/59.jpg)
For bridgeless cubic graphs: A 2-factor covering all 3- and 4-edge cuts: Algorithm A min-weight 2-factor covering all 3-edge cuts:
Algorithm Polyhedral description
Summary59
For 3-edge-connected cubic graphs 6/5-approx. algorithm for
the min. 2-edge-connected subgraph problem
Min-weight 2-factor covering all 3- and 4-edge cuts in bridgeless cubic graphs
6/5–approx. algorithm for the min. 2-edge-connected subgraph problem in bridgeless cubic graphs
Open Problems