マトロイド
マトロイド(英: matroid)は、ある公理を満たす集合とそのべき集合の部分集合の組である。歴史的には、行列の一次独立・従属を一般化した概念であるが、多くの組合せ最適化問題をマトロイドあるいはより緩い独立性システムとコスト関数で定式化でき、特徴付けを行える等応用範囲は広い。特に組合せ最適化において、マトロイド上の最適化問題には単純な貪欲法によって多項式時間のアルゴリズムとは限らないものの最適解が得られることは非常に重要である。
定義
[編集]有限集合 E とその部分集合族 の組 (E, F) が[注 1]
- (A1)
- (A2)
- (A3)
を満たすとき、マトロイドと呼ばれ、(A1) および (A2) のみを満たすとき独立性システム (independence system) と呼ばれる。さらに (A1) および (A3) を満たすときグリードイド (greedoid) と呼ばれる。
以下に本項で使う用語の定義を挙げる。
- 独立集合 (英: independent set) - の要素
- 従属集合 (英: dependent set) - の要素
- サーキット (英: circuit) - 極小な従属集合
- 基 (英: base, basis) - 極大な独立集合
- の基 - とする の部分集合の中で極大な独立集合
ランク、階数関数
[編集]独立性システム (E, F) の階数 (rank) 関数 は、 に対して
と定義される。マトロイドならば、公理(A3) から E の部分集合 X に対して X のどの基も元数は等しいため、階数関数をの基の大きさと定義しても構わない。
例えば、E = {1, 2, 3}, F={{},{1},{2},{1,2}} というマトロイドならば r({1}) = 1, r({1, 2}) = 2, r({1, 2, 3}) = 2, r({1, 3}) = 1 等となる。
独立性システム の階数関数 は任意の と について次の性質を持つ。
- (R1)
- (R2)
- (R3)
さらに、 がマトロイドならば次の性質も持つ。
- (R4)
- (R5)
- (R6)
特に(R4)に示されているように、階数関数が劣モジュラであることは、マトロイドの極めて重要な性質である。
閉包
[編集]独立性システム の閉包 (closure) 関数 は、 に対して
で定義される。 を X の閉包と呼ぶ。
マトロイド の閉包関数 は次の性質を持つ。
- (L1) 任意の に対して、
- (L2) 任意の に対して、
- (L3) 任意の に対して、
- (L4) 任意の と任意の に対して、
ランク商
[編集]下方ランク (lower rank) 関数 を に含まれる基の最小元数とする。つまり、
- かつ に対し
と定義する。すると、ランク商 (rank quotient) は次のように定義される。
マトロイドのとき である。独立性システムのとき であり、 に対して が高々p 個しかサーキットを持たないならば である[1]ことが知られているのでランク商を見積もることが可能である。
例
[編集]一様マトロイド
[編集]E を有限集合とし、ある整数 r 以下の要素数を持つ 2E の部分集合をFとするとき、(E, F) はマトロイドとなる。これを、一様マトロイド (uniform matroid) と呼び、n=|E|としたとき などと書く。 における基は、要素数がrであるようなEの部分集合であり、サーキットは要素数がr+1であるようなEの部分集合である。また、一様マトロイドの直和は分割マトロイド (partition matroid) と呼ばれる[注 2]。具体的には、を互いに素な集合とし、それぞれに を定めたとき、 それぞれにおいて を満たす の集合を F とすれば、マトロイドとなる。
グラフ理論におけるマトロイド
[編集]E を無向グラフGの辺集合、F を森[注 3]の集合とするとき、(E, F) はマトロイドとなり、グラフ的マトロイド (graphic matroid) または 閉路マトロイド (cycle matroid) と呼ぶ。しばしばグラフGがグラフ的マトロイドであることを M(G) と書く。グラフ的マトロイドにおける従属集合は、閉路を持つグラフの集合であるから、サーキットとは単純閉路となる辺集合である。また(X の)基とは(X の部分集合の中で)できうる最大の森のことで、明らかに(X でカバーされる)点の数 − 1 本の辺で作られる森が最大であり、この森の集合を言う。よって、階数関数は r(X) = (Xがカバーする点数) − 1 と書ける。
他にも、グラフにおけるマトロイドはいくつか知られている。
- E を無向グラフ G の辺集合、F を (G, X) がpseudoforest[注 4]となるようなXの集合とするとき (E, F) はマトロイドとなる。これを、bicircularマトロイド (bicircular matroid) と呼ぶ。
- 点集合 U, V、辺集合 X とする2部グラフ G = (U, V, X) において、マッチングの端点となる点集合 を要素とする集合をFとするとき、(U, F) はマトロイドとなる。これを、横断マトロイド (transversal matroid) と呼ぶ。完全2部グラフ の横断マトロイドは、一様マトロイド である。
- 点集合を V とするグラフにおいて、点の部分集合を とする。U への点素な |U| 本の道が存在する V' の部分集合を F の要素とすると、(V', F) はマトロイドとなる。これを、ガンモイド (gammoid) と呼ぶ。特に、(V, F) を狭義ガンモイド (strict gammoid) と呼ぶ。
線形代数におけるマトロイド
[編集]E を体上の行列の列集合、その体上で線形独立である列の集合を F とするとき、(E, F) はマトロイドとなり、ベクトルマトロイド (vector matroid) と呼ぶ。マトロイドが同等の体 K 上のベクトルマトロイドとして記述できるとき、表現可能であると呼ばれる。任意の体上で表現可能なマトロイドを正則マトロイド (regular matroid) と呼び、位数2の有限体上で表現可能なマトロイドを2値マトロイド (binary matroid) と呼ぶ。これらは、
- マトロイド ⊃ 2値マトロイド ⊃ 正則マトロイド ⊃ グラフ的マトロイド
という包含関係が成り立つ。一方で、Fanoマトロイドは、2値マトロイドであるが(実数体上では表現できないため)正則マトロイドではない。また、Vámosマトロイド (Vámos matroid) は、任意の体上で表現できないマトロイドの最も簡単な例である。
その他のマトロイド
[編集]- 2部マトロイド (bipartite matroid) は、サーキットの大きさがすべて偶数であるマトロイドである。2部グラフのグラフ的マトロイドを一般化しており、グラフ的マトロイドが2部マトロイドであることと、グラフが2部グラフであることは同値である。
- シルベスターマトロイド (Sylvester matroid) は、すべてのサーキットの大きさが3であるようなマトロイドである。例えば、ランク2の一様マトロイド はシルベスターマトロイドである。
- pavingマトロイド (paving matroid) は、すべてのサーキットの大きさがランクよりも大きいマトロイドである。例えば、一様マトロイド のランクはrであり、サーキットの大きさは常に r + 1 であるため、pavingマトロイドである。
- オイラーマトロイド(eulerian matroid) は、要素がサーキットによって分割できるようなマトロイドである。例えば、一様マトロイド は r + 1 が n を割り切るとき、かつそのときのみオイラーマトロイドである。
他の公理系
[編集]集合Eとその部分集合の族Fが(A1)から(A3)を満たすときマトロイドと呼ぶことにし、そこから基やランクを定義した。だが、実はこれらの性質を持つ族あるいは関数からマトロイドとなる F を得ることができる。
基族
[編集]有限集合 E と とする。 がマトロイド (E, F) の基族であるための必要十分条件は次の(B1),(B2)が成り立つことである。
- (B1)
- (B2) 任意の について となる が存在する。
基が1つしかない場合は明らかにマトロイドとなる。そうでない場合、例えば とすると は(B1)および(B2)を満たす。このような基の族を持つマトロイドは、一様マトロイド ただ1つに決まることが分かる。また、基族が のときは(B2)を満たさないため、この基族ではマトロイドにならない。事実、 とすると、この例の基族となるが、(A3)を満たさずマトロイドになっていない。
サーキットの族
[編集]有限集合 E と C ⊆ 2E とする。C がマトロイド (E, F) のサーキットの族であるための必要十分条件は次の(C1),(C2),(C3)が成り立つことである。
- (C1)
- (C2) 任意の C1, C2 ∈ C について C1 ⊆ C2 ならば C1 = C2 である。
- (C3) 任意の C1, C2 ∈ C は C1 ≠ C2 で c ∈ C1 ∩ C2 とするとき、 となる C3 ∈ C が存在する。
ランク関数
[編集]マトロイドのランク関数は(R1)から(R6)を満たすが、逆にEと(R1),(R2),(R4)を満たす[注 5]関数 を与えればFは直ちに決まり (E, F) はマトロイドであり、r はランク関数である。
例えば、E = {1, 2, 3} とし、関数 r を
と定義すると、r は(R1),(R2),(R4)を満たすことが確認できる。すると、r(X) = |X| となる X の族を F と定義すると F = {{1}, {3}, {1,3}} となり、(E, F) はマトロイドになることが確認できる。このように、r を決めれば対応する F がただ1つに決まる。
閉包関数
[編集](L1)から(L4)を満たす関数 はマトロイドの閉包関数となる。
マトロイドの構成法
[編集]ここでは、1つ以上のマトロイド(あるいは独立性システム)から新たなマトロイドを構成する方法について説明する。
双対
[編集]独立性システム (E, F) に対して、F* を となる (E, F) の基 B が存在する の集合とする。このとき、(E, F*) を (E,F) の双対 (dual) と定義する。
双対には次のような性質がある。
- (E, F**) = (E,F)
- (E, F) が独立性システムならば、(E, F*) は独立性システム
- (E, F) はマトロイド (E, F*) はマトロイド[2]。
- マトロイド (E, F) とその双対 (E, F*) とし、それぞれのランク関数を r, r* とすると、r* は任意の X ⊆ E に対して である。
例えば、 E = {1, 2, 3}, F = {{1}, {2}, {1,2}} というマトロイドに対して基の族 B = { {1, 2} } だから F* = { {3} } となる。双対のランク関数も、例えば r*({1, 3}) = 2 + r({2}) - r({1, 2, 3}) = 2 + 1 - 2 = 1 となるように、成り立っていることが分かる。
また、平面グラフに対する双対とグラフ的マトロイドの双対の概念は一致する。つまり、任意の平面的グラフ G の閉路マトロイド M(G) の双対は、G の双対平面グラフ G* のマトロイドと(平面埋め込みの方法によらず)同一である。
交差
[編集]2つの独立性システム (E, F1), (E, F2) とするとき、(E, F1 ∩ F2) を2つの独立性システムの交差 (intersection) と定義する。有限個の独立性システムの交差も同様に定義でき、交差もまた独立性システムとなる。任意の独立性システム (E, F) は、有限個のマトロイドの交差で表せる。さらに、p 個のマトロイドの交差ならば、q(E, F) ≧ 1/p。
例えば2部グラフのマッチング問題の場合、Eを辺集合、Fをマッチング集合とすれば、2部グラフであるから点集合は A ∪ B と書ける。F1 を任意の点 a ∈ A に繋がる辺は高々1本であるという条件とするならば (E, F1) はマトロイドであり、同様に F2 を任意の点 b ∈ B に繋がる辺は高々1本であるという条件とするならば (E, F2) もマトロイドである。F = F1 ∩ F2 なので、(E, F) は2つのマトロイド (E, F1), (E, F2) の交差である。
合併
[編集]2つのマトロイド (E, F1),(E, F2) とする。X1 ∈ F1, X2 ∈ F2 となる X の分割 X = X1 ∪ X2 が存在するとき、X ⊆ E は 分割可能 (partitionable) であると呼ぶ。 とするとき、(E, F) を (E, F1), (E, F2) の合併(union) あるいは和(sum) と呼ぶ。有限個のマトロイドの合併も同様に定義できる。
- マトロイドの合併は、マトロイドとなる。
- k個のマトロイド (E, F1), …, (E, Fk) の各ランク関数を r1, …, rk とすると、これらの合併 (E, F) のランク関数は である[3]。
組合せ最適化
[編集]組合せ最適化問題の多くは、独立性システム とコスト関数 に対して、 を最大(あるいは最小)にする を求める最適化問題に定式化できる。
例えば、以下の中で最小全域木問題はマトロイドになるが、他はマトロイドにはならず、独立性システムとなる。
- 巡回セールスマン問題 - Eをグラフの辺、Fはハミルトン閉路の部分集合
- ナップサック問題 - Eを荷物、Fは規定の重さを超えない荷物の組合せ
- 最小全域木問題 - Eはグラフの辺、Fはグラフの森の集合
貪欲法
[編集]マトロイドにおいては貪欲法で最適解が得られることを示す。なお、マトロイドの貪欲法は組合せ最適化の貪欲法の全てを網羅しているわけではない。たとえば、クラスカル法はマトロイド上の貪欲法で説明できるが、プリム法やダイクストラ法は異なる。
最良選択貪欲法
[編集]最良選択貪欲法は、 の大きい順に を暫定解に追加できるならば追加し追加できない場合は次の要素に移ることを繰り返すアルゴリズムである。つまり、
- となるように をソートする。
- For i = 1 to n : ならば
- を解として出力する。
というアルゴリズムである。
最良選択貪欲法で得られる解のコストを 、最適解のコストを とすると、ランク商 を使って
が任意のコスト関数に対して成立する[4][5]。マトロイドのランク商は1なので、マトロイドである最大化問題は最良選択貪欲法によって最適解を得られる。これは逆も言えるので、独立性システム (E, F) がマトロイドであるための必要十分条件は最良選択貪欲法で全ての に対して最大化問題の最適解を求めることができることである[6][7]。これを Edmonds-Rado 定理という。
証明の概要
[編集]の証明の概要を述べる。
G を最良選択貪欲法で得られる解、 を最適解とする。 は、 となるようにソート済みであるとし、 とする。また、
とし、任意の に対して、 と置く。このとき、次の事実が成り立つ。
以上のことから、
となる。
最悪棄却貪欲法
[編集]独立性システム (E, F) とコスト関数 に対する最小化問題を解く。最悪棄却貪欲法は、都合の悪い e ∈ E を優先して解から除外する。つまり、
- となるように をソートする。
- For i = 1 to n : が基を含むならば、
- を解として出力する。
というアルゴリズムである。
最悪棄却貪欲法で得られる解のコストを c(G)、最適解のコストを c(OPT) とすると、ランク商 q(E, F)、双対な独立性システム (E, F*) の下方ランク関数 ρ*、ランク関数 r* を使って
と書ける[8]。マトロイドならば ρ* = r* なので、常に最適解を得られることが分かる。
オラクル
[編集]組合せ最適化問題において F は明示的に与えられることはまずない。F を列挙しようとすることは無謀であるので、現実には E とコスト関数cのみが与えられる。最良選択貪欲法を実行するには、さらに独立性オラクルを必要となる。独立性オラクルとは、X ⊆ E が与えられたとき X ∈ F であるかどうかを判定するオラクルである。これがないと最良選択貪欲法の3番目のステップは実行できない。同様に最悪棄却貪欲法を実行するためには X ⊆ E が与えられたとき X が基[注 10]を含むかを判定する基拡張集合オラクルを必要とする。
では、独立性オラクルか基拡張集合オラクルどちらか一方が与えられたとき、そのオラクルを使ってもう一方を多項式時間で実行可能(多項式等価)だろうか[注 11]。例えば、巡回セールスマン問題に対する独立性システムの独立性オラクルはつまり、与えられた辺集合がハミルトン閉路の部分集合であるかを判定するものであるが、グラフは完全グラフであるので、多項式時間で実行可能である。対して基拡張集合オラクルは与えられた辺集合からいくらか辺を削除することによってハミルトン閉路になるかということを判定しなくてはならない。それはつまりハミルトン閉路問題と等価であり、ハミルトン閉路問題はNP完全である[9]ため難しいと言える。このように、独立性システムにおいて独立性オラクルと基拡張集合オラクルは必ずしも多項式等価ではない。
マトロイドにおいては独立性オラクル、基拡張集合オラクル、ランク関数を返すランクオラクル、閉包関数を返す閉包オラクル全て多項式等価である。しかし、与えられた部分集合が基かどうかを判定する基決定オラクルは独立性オラクルより弱い[注 12]し、与えられた部分集合の最小元数の従属部分集合を返すオラクルは独立性オラクルより強い[10]。
近似
[編集]最適化問題は厳密解を求めることが現実的でないことが多いために、近似の限界についても研究されている。次の問題が効率的に解ける(入力のサイズと 1/ε の多項式時間で解を出力するアルゴリズムが存在する)ことと、誤差が高々 1+ε 倍の解を出力する多項式時間アルゴリズムが存在することは同値である[11]。
独立性システム (E, F)、コスト関数 、部分集合 S, S' ⊆ E, ε > 0 が
であるとき、S ⊆ B となる基 B が存在して S' ⊆ B' となる基 B' 全てに対して (1 + ε)c(B) ≧ c(B') が成立するか?
つまり、部分的なコスト(c(S) や c(S'))が高々 1 + ε 倍違う程度ならば、それらからできうる最適解も 1 + ε 倍程度しか違わないだろうか、という問題である。部分が最適ならば全体も最適であるという場合はε=0であり、よく知られているように動的計画法が存在する。よって、(E, F) をマトロイドに限定するならば多項式時間アルゴリズムは存在する。
ナップサック問題はこのアルゴリズムが知られている珍しい例で、計算時間が [12][13][14]や [15]で、出力される解の評価が最適解の評価の高々 1 + ε 倍であるアルゴリズムがある。
マトロイドに関する問題
[編集]マトロイド交差問題と分割問題
[編集]マトロイド交差問題 (Matroid Intersection Problem) は、2つのマトロイド (E, F1), (E, F2) が与えられたとき、|F| が最大となるような F ∈ F1 ∩ F2 を求める問題である。マトロイド交差問題は多項式時間で解ける。また、3つ以上のマトロイド交差問題も同様に考えることができるが、これはNP困難問題である。
重み付き版についてもアルゴリズムが知られていて[16]、2つの独立性オラクルの計算量の大きい方を α とすると O(|E|3α) で解ける。
マトロイド分割問題 (Matroid Partitioning Problem) は k 個のマトロイド (E, F1), …, (E, Fk) が与えられたとき |X| が最大になるような分割可能な X ⊆ E を求める問題である。マトロイド交差問題とマトロイド分割問題は等価である。
一般化
[編集]グリードイドはマトロイドと反マトロイドを一般化したものである。グリードイドにも貪欲法が定式化できて、特殊な条件下においては最適解を出力する。だが、グリードイド上での最適化問題はNP困難であることが知られている。
また、マトロイドのランク関数が劣モジュラ関数であることは既に述べたが、有限集合Eと劣モジュラ関数 を用いてポリマトロイド (polymatroid) と呼ばれる有界多面体を定義できる。ポリマトロイドとベクトルに対する分離問題は劣モジュラ関数最小化問題に帰着できる。劣モジュラ関数最小化問題は例えばフローネットワークにおける無向グラフの最小カットを求める問題などを一般化している。劣モジュラ関数最小化問題は楕円体法を用いることで多項式時間で解ける[17]ことが知られて以来、Schrijverのアルゴリズム[18]等が知られている。
脚注
[編集]注釈
[編集]- ^ 記号の意味については「冪集合」「空集合」「集合間の関係を表す記号」「濃度 (数学)」「和集合」「差集合」を参照のこと
- ^ マトロイドの直和もマトロイドになる。
- ^ 閉路のない辺集合
- ^ 各連結成分において、高々1つの閉路を持つグラフ
- ^ (R3),(R5),(R6)を満たす関数と(R1),(R2),(R4)を満たす関数は同値
- ^ は、 のとき、かつそのときのみ1となるから、 が得られる。これを でまとめて右辺を得る。
- ^ Gj は Ej の基であり、ρ の定義より得られる。
- ^ ランク商の定義より明らか。
- ^ であることと、階数関数の定義および性質(R1)より得られる。
- ^ X の基でないことに注意
- ^ 概念は「多項式時間変換」に詳しい
- ^ 最良選択貪欲法を使うことによって(本質的には独立性オラクルを複数回使うことによって)基決定オラクルを作れるが、逆に基決定オラクルを多項式回使っても独立性オラクルを作れない
出典
[編集]- ^ D. Hausmann; B. Korte , T. A. Jenkyns (1980). “Worst case analysis of greedy type algorithms for independence systems”. Mathematical Programming Studies 12: 120-131.
- ^ Hassler Whitney (1935-07). “On the abstract properties of linear dependence” (pdf). American Journal of Mathematics 57: 509-533 2011年3月19日閲覧。.
- ^ Crispin Nash-Williams (1967). “An application of matroids to graph theory”. In P.Rosenstiehl, ed.. Theory of Graphs; proceedings of an international symposium in Rome 1966. New York: Gordon and Breach. pp. 263-265
- ^ T.A. Jenkyns (1976). “The Efficacy of the greedy Algorithm”. Proc. of 7th S-E. Conf. on Combinatorics, Graph Theory, and Computing, Utilitas Mathematica, Winnipeg: 341-350.
- ^ Bernhard Korte; Dirk Hausmann (1978). “An Analysis of the greedy heuristic for independence systems”. In B. Alspach, P. Hell, D.J. Miller, eds.. Aogorithmic Aspects of Combinatorics; Annals of Discrete Mathematics. 2. Amsterdam: North-Holland. pp. 65-74
- ^ R. Rado (1957). “Note on Independence Functions”. Proceedings of the London Mathematical Society 7: 300-320.
- ^ Jack Edmonds (1971). “Matroids and the greedy algorithm”. Mathematical Programming 1: 127-136.
- ^ Bernhard Korte; C.L. Monma (1979). “Some remarks on a classification of oracle-type algorithms”. In L. Collatz, G. Meinardus, W. Wetterling, eds.. Numerische Methoden bei graphentheoretischen und kombinatorischen Problemen. 2. Basel: Birkhäuser. pp. 195-215
- ^ Richard M. Karp (1972). “Reducibility Among Combinatorial Problems”. In R. E. Miller and J. W. Thatcher eds. Complexity of Computer Computations. New York: Plenum. pp. 85-103
- ^ D. Hausmann; B. Korte (1981). “Algorithmic versus axiomatic definitions of matroids”. Mathematical Programming Study 14: 98-111.
- ^ B. Korte; R. Schrader (1981). “On the existence of fast approximation schemes”. In O. Mangaserian, R.R. Meyer, S.M. Robinson, eds.. Nonlinear Programming. New York: Academic Press. pp. 415-437
- ^ Oscar H. Ibarra; Chul E. Kim (1975-10). “Fast Approximation Algorithms for the Knapsack and Sum of Subset Problems”. Journal of the ACM 22 (4): 463-468. doi:10.1145/321906.321909. ISSN 0004-5411.
- ^ Sartaj K. Sahni (1976-01). “Algorithms for Scheduling Independent Tasks”. Journal of the ACM 23 (1): 116-127. doi:10.1145/321921.321934. ISSN 0004-5411.
- ^ G.V. Gens (1979). “Computational complexity of approximation algorithms for combinatorial problems”. In J. Becvar, ed.. Mathematical Foundations of Computer Science. Berlin: Springer. pp. 292-300
- ^ Eugene L. Lawler (1979). “Fast Approximation Algorithms for Knapsack Problems”. Mathematics of Operations Research 4: 339-356. doi:10.1287/moor.4.4.339.
- ^ A. Frank (1981). “A weighted matroid intersection algorithm”. Journal of Algorithms 2: 328-336.
- ^ M. Grötschel; L. Lovász, A. Schrijver (1981). “The ellipsoid method and its consequences in combinatorial optimization” (pdf). Combinatorica 1: 169-197 2011年3月26日閲覧。.
- ^ Alexander Schrijver (2000). “A combinatorial algorithm minimizing submodular functions in strongly polynomial time” (pdf). Journal of Combinatorial Theory, Series B 80 (2): 346-355 2023年5月14日閲覧。.
参考文献
[編集]- B.コルテ、J.フィーゲン 著、浅野孝夫、平田富夫、小野孝男、浅野泰仁 訳『組合せ最適化-理論とアルゴリズム』(第2版)シュプリンガー・ジャパン、2012年2月29日。ISBN 978-4621062029。