メインコンテンツへスキップ

【月刊組合せ論 Natori】マトロイドに入門しよう【2024 年 7 月号】

·
箱星
著者
箱星
のんびり暮らしたい。
目次

月刊組合せ論 Natori は面白そうな組合せ論のトピックを紹介していく企画です。今回は競技プログラミングにも登場するマトロイドに入門していきます。タイトルの「入門しよう」には、筆者と共に入門しようという気持ちが込められています。つまり筆者も素人です。

ベクトル空間
#

ベクトル a=(1,1,0),b=(0,1,1),c=(3,2,1),d=(1,1,1)\vec{a}=(1,-1,0), \vec{b}=(0,1,-1),\vec{c}=(3,-2,-1),\vec{d}=(1,1,1) を考えます。c\vec{c}c=3a+b\vec{c}=3\vec{a}+\vec{b} と表せる一方、d\vec{d}pa+qbp\vec{a}+q\vec{b} の形で表すことができません。このような関係を表す概念として、線形独立と線形従属があります。

ベクトルの集合 X={x1,,xn}X=\{\vec{x}_1,\ldots,\vec{x}_n\}線形独立であるとは、a1x1++anxn=0a_1\vec{x}_1+\cdots+a_n\vec{x}_n=\vec{0} ならば a1==an=0a_1=\cdots=a_n=0 が成り立つことをいいます。そうでないとき XX線形従属であるといいます。例えば {a,b,d}\{\vec{a},\vec{b},\vec{d}\} は線形独立ですが、{a,b,c}\{\vec{a},\vec{b},\vec{c}\}3a+bc=03\vec{a}+\vec{b}-\vec{c}=\vec{0} をみたすので線形従属です。

XX が線形独立ならば、部分集合 YXY\subset X も線形独立です。さらに次の性質が成り立ちます。

命題: X,YX,Y はベクトル空間 VV の線形独立な部分集合で X>Y|X|>|Y| をみたすとする。このときある xXYx\in X\setminus Y であって Y{x}Y\cup\{x\} が線形独立であるものが存在する。

証明: もしすべての xXYx\in X\setminus Y について Y{x}Y\cup\{x\} が線形従属であったとすると、XX の元はすべて YY の元の線形結合として表せる。YY を基底とする部分空間の次元は Y|Y| で、X>Y|X|>|Y| であることから XX は線形従属となる。これは矛盾である。

グラフ
#

グラフは頂点と辺からなるものです。厳密にはグラフ GG は頂点集合 V(G)V(G) と辺集合 E(G)E(G) の組で、E(G)E(G)V(G)V(G) の非順序対 {x,y}\{x,y\} からなります。ここでは xyx\ne y とします。

グラフ GGサイクルとは、頂点の列 (v1,v2,,vn)(v_1,v_2,\ldots,v_n) (ただし n3n\ge 3 とする) であって {v1,v2},{v2,v3},,{vn1,vn},{vn,v1}\{v_1,v_2\},\{v_2,v_3\},\ldots,\{v_{n-1},v_n\},\{v_n,v_1\} が辺であるものをいいます。とは、サイクルをもたない部分グラフのことをいいます。

命題: X,YX,Y はグラフ GG の森で、E(X)>E(Y)|E(X)|>|E(Y)| をみたすとする。このときある辺 eE(X)E(Y)e\in E(X)\setminus E(Y) であって Y{e}Y\cup\{e\} が森であるものが存在する。

証明: すべての eE(X)E(Y)e\in E(X)\setminus E(Y) について Y{e}Y\cup\{e\} が森ではないと仮定する。また V(X)=V(Y)=V(G)V(X)=V(Y)=V(G) とする。r(X)r(X)XX の連結成分の個数とするとき、r(X)=V(X)E(X)r(X)=|V(X)|-|E(X)| より r(X)<r(Y)r(X)<r(Y) が成り立つ。一方、e={u,v}E(X)e=\{u,v\}\in E(X) に対して、頂点 u,vu,vYY の同じ連結成分に属する。このことから r(X)r(Y)r(X)\ge r(Y) となる。これは矛盾である。

マトロイド
#

マトロイドはベクトル空間とグラフに対する上述の性質を抽象化したものです。

EE を有限集合、I\mathcal{I}EE の部分集合からなる集合とします。(E,I)(E,\mathcal{I})マトロイドであるとは

  1. I\emptyset\in \mathcal{I}
  2. XI,YXX\in \mathcal{I}, Y\subset X ならば YIY\in \mathcal{I}
  3. X,YI,X>YX,Y\in \mathcal{I}, |X|>|Y| ならばある xXYx\in X\setminus Y について Y{x}IY\cup\{x\}\in \mathcal{I}

をみたすことをいいます。XIX\in \mathcal{I} であるとき、XX は独立であるといいます。

1, 2 をみたす場合、(E,I)(E,\mathcal{I}) は独立性システムであるといいます。

ランク関数
#

マトロイド (E,I)(E,\mathcal{I}) に対してランク関数 r ⁣:2EZ+r\colon 2^E\to \mathbb{Z}_+r(X)=max{YYX,YI}r(X)=\max\{|Y| \mid Y\subseteq X, Y\in \mathcal{I}\} により定めます。すると

  • 0r(X)X0\le r(X)\le |X|
  • YXY\subseteq X ならば r(Y)r(X)r(Y)\le r(X)
  • r(X)+r(Y)r(XY)+r(XY)r(X)+r(Y)\ge r(X\cup Y)+r(X\cap Y)

が成り立ちます。最初の 2 つは明らかです。三番目の不等式は劣モジュラ性と呼ばれています。

逆を考えることもできます。EE は有限集合で関数 r ⁣:2EZ+r\colon 2^E\to \mathbb{Z}_+ は上の 3 つの条件をみたすとします。

I={XEr(X)=X} \mathcal{I}=\{X\subseteq E\mid r(X)=|X|\}

とおくと、(E,I)(E,\mathcal{I}) はマトロイドになります。証明は省略します。

このように、マトロイドはランク関数によって導入することもできます。他にも様々なマトロイドの特徴づけが知られています。次のセクションではアルゴリズムによってマトロイドを特徴づける話を紹介します。

貪欲法
#

「マトロイドは貪欲法」というフレーズを聞いたことがある人もいるかもしれません。実は、貪欲法が上手く動くこととマトロイドであることが同値になります。これを正確に述べます。

(E,I)(E,\mathcal{I}) を独立性システムとします。各 eEe \in E に対して重み w(e)R0w(e)\in\mathbb{R} _ {\ge 0} が定まっているとします。求めたいものは、極大独立集合 BB であって重みの和 eBw(e)\sum_{e\in B}w(e) が最小となるものです。この問題に対して、次のような貪欲法を考えます。

w(e1)w(en)w(e_1)\le\cdots\le w(e_n) となるように EE の元に e1,,ene_1,\ldots,e_n と名前を付けます。まず B=B=\emptyset とします。i=1,2,,ni=1,2,\ldots,n の順に次の操作を行います。

  • B{ei}IB\cup\{e_i\}\in\mathcal{I} ならば BBeie_i を加える。

最終的な BB が出力です。

このアルゴリズムは本当に最小コストの極大独立集合を出力するでしょうか?

定理: (E,I)(E,\mathcal{I}) を独立性システムとする。任意の関数 w ⁣:ER0w\colon E\to\mathbb{R} _ {\ge 0} に対して貪欲法が最小コストの極大独立集合を出力することと、(E,I)(E,\mathcal{I}) がマトロイドであることは同値である。

なお極大独立集合のことは基ともいいます。基の要素数は等しいことが知られており、以下の証明でも用いられます。

まずはマトロイドならば貪欲法が上手く動くことを示します。貪欲法の出力が B={x1,,xr}B=\{x_1,\ldots,x_r\} で、x1,,xrx_1,\ldots,x_r の順に加えられたとします。もちろん w(x1)w(xr)w(x_1)\le\cdots\le w(x_r) です。B={y1,,yr}IB^{\prime}=\{y_1,\ldots,y_r\}\in\mathcal{I} (w(y1)w(yr)w(y_1)\le\cdots\le w(y_r)) が最適解であるとします。このとき任意の ii に対して w(xi)w(yi)w(x_i)\le w(y_i) が成り立つことを示します。そうでないと仮定すると、w(xk)>w(yk)w(x_k)>w(y_k) となる kk が存在します。このような kk のうち最小のものをとり、S={x1,,xk1},T={y1,,yk}S=\{x_1,\ldots,x_{k-1}\}, T=\{y_1,\ldots,y_k\} を考えます。T=S+1|T|=|S|+1 なので、マトロイドの性質よりある j (1jr)j \ (1\le j\le r) について S{yj}IS\cup\{y_j\}\in\mathcal{I} となります。しかし w(yj)w(yk)<w(xk)w(y_j)\le w(y_k)<w(x_k) なので、貪欲法の動かし方と合いません。

次に貪欲法が上手く動くならばマトロイドであることを示します。マトロイドでないとすると、ある独立集合 A,BA,BA=B+1|A|=|B|+1 をみたすが、任意の aABa\in A\setminus B について B{a}B\cup\{a\} が独立ではありません。重みを

w(x)={w1xABw2xBAw3xABw4otherwise w(x)=\begin{cases} w_1 & x\in A\cap B \\ w_2 & x\in B\setminus A \\ w_3 & x\in A\setminus B \\ w_4 & \text{otherwise} \end{cases}

とします。まず w1<w2<w3<w4w_1<w_2<w_3<w_4 とします。貪欲法は、まず ABA\cap B の元をすべて選び、BAB\setminus A の元をすべて選びます。仮定より ABA\setminus B の元は選びません。BB は極大独立集合でないので、X(AB)X\setminus (A\cup B) から選ばれます。mm を極大独立集合の要素数としたとき、コストは

ABw1+BAw2+(mB)w4 |A\cap B|w_1+|B\setminus A|w_2+(m-|B|)w_4

となります。これが最適でないことを示したいので、A,BA,B の役割を変えて得られるコスト

ABw1+ABw3+(mA)w4 |A\cap B|w_1+|A\setminus B|w_3+(m-|A|)w_4

の方が小さいようにしたいです。この不等式は

w4>ABw3BAw2 w_4>|A\setminus B|w_3-|B\setminus A|w_2

と変形できます。これが成り立つように w1,w2,w3,w4w_1,w_2,w_3,w_4 を具体的に定めることができますが、これは読者の演習問題とします。

また、最小を最大に置き換えても類似の定理が成り立ちます。

競技プログラミングにおける具体例
#

競技プログラミングにおいてマトロイドが現れる例をいくつか紹介します。

代表的な例は最小全域木です。マトロイドはグラフの一般化だったことを思い出すと、サイクルのないグラフのうち辺の数が極大なものが基 (極大独立集合) です。グラフが連結であるとすると、基は全域木です。するとマトロイドに対する貪欲法を用いることで最小全域木が求められます。このアルゴリズムにはクラスカル法という名前がついています。

他にも、AtCoder Beginner Contest 236 F - Spices がマトロイドの考え方で解けます。その他のマトロイドが現れる問題はぜひ詳しい人に聞いてみてください。

競技プログラミングとマトロイドに関する記事は色々あります。いくつかピックアップして紹介します。

おわりに
#

マトロイドの基礎を紹介しました。筆者も素人なので、詳しい人からのコメントをお待ちしております。

高度な話題も、勉強が進めば記事にしたいと思っています。

今後も月刊組合せ論 Natori では組合せ論の様々な話題をお届けする予定なので、応援のほどよろしくお願いいたします。

参考文献
#