次数制約付き最小全域木問題

次数制約付き最小全域木問題

このページでは、次数制約付きの最小全域木問題というNP困難な問題をイジングモデルにより表現する方法について述べます。

問題の定式化 #

無向グラフ G=(V,E)G = (V, E) が与えられます。 また各辺 uvEuv \in E について、コスト cuvc_{uv} が定まっています。 さらに、正の自然数の定数 Δ\Delta が与えられます。

このとき、 GG の全ての頂点を含む木 TGT \subseteq G であって、 TT のコスト c(T)=uvETcuv c(T) = \sum_{uv \in E_T} c_{uv} が最小になるようなもので、 TT の各頂点の次数が Δ\Delta 以下となるような物を求める問題を考えましょう。

問題の具体例 #

格子点上に7個の頂点が以下の図のような位置関係で存在しています。 このときの最小全域木は図の左側です。 この最小全域木は最大次数が3であるので、 次数制限 Δ3\Delta\geq3 のとき、次数制約付きの最小全域木は図の左側です。 次数制限 Δ=2\Delta=2 のとき、次数制約付きの最小全域木は図の右側です。 ( Δ1\Delta\leq1 のときは解がありません。) 問題の具体例

イジングモデルへの変形 #

この問題を解くためのイジングモデルのエネルギー関数は以下のようになります。

まず、必要な変数をいくつか導入します。まず、 N=VN = |V| としましょう。 各辺 eEe \in E について、それが TT の辺に含まれるなら1を、そうでないなら0を取るようなバイナリ変数 yey_e を導入します。 また、各頂点 vVv \in V と自然数 i{0,1,,N}i \in \{ 0, 1, \ldots, N \} に対して、頂点 vv が木 TT において根から深さ ii のところにあるときに1を取り、それ以外の時に0を取るようなバイナリ変数 xv,ix_{v,i} を導入します。 また、各辺 uvEuv \in E と自然数 i{0,1,,N}i \in \{ 0, 1, \ldots, N \} に対して、木 TT において辺 uvuv が根から深さ ii の所に含まれて、かつ uu の方が vv よりも根に近い時に1を、そうでない時に0を取るバイナリ変数 xuv,ix_{uv, i} を導入します。 xuv,ix_{uv,i} xvu,ix_{vu,i} は別のスピンであることに注意してください。 さらに、各頂点 vVv \in V と自然数 j{1,,Δ}j \in \{ 1, \ldots, \Delta \} に対して、木 TT において頂点 vv の次数が kk の時に j=1Δzv,j=k\sum_{j=1}^\Delta z_{v,j}=k となるようなバイナリ変数 zv,jz_{v,j} を導入します。 AとBを正の定数とします。

これらを元にして、エネルギー関数は以下のようになります。 HA1=A(1vVxv,0)2 H_{A_1} = A \left( 1 - \sum_{v \in V} x_{v,0} \right)^2 HA2=AvV(1i=0N/2xv,i)2 H_{A_2} = A \sum_{v \in V} \left( 1 - \sum_{i = 0}^{N/2} x_{v,i} \right)^2 HA3=AuvE(yuvi=1N/2(xuv,i+xvu,i))2 H_{A_3} = A \sum_{uv \in E} \left( y_{uv} - \sum_{i=1}^{N/2} (x_{uv,i} + x_{vu,i}) \right)^2 HA4=AvVi=1N(xv,iu:uvExuv,i)2 H_{A_4} = A \sum_{v \in V} \sum_{i = 1}^{N} \left( x_{v,i} - \sum_{u: uv \in E} x_{uv,i} \right)^2 HA5=AvV(j=1Δzv,ju:uvEi=1N/2(xuv,i+xvu,i))2 H_{A_5} = A \sum_{v \in V} \left( \sum_{j=1}^\Delta z_{v,j} - \sum_{u: uv \in E} \sum_{i=1}^{N/2} (x_{uv,i} + x_{vu,i}) \right)^2 HA6=AuvEi=1N/2xuv,i(2xu,i1xv,i) H_{A_6} = A \sum_{uv \in E} \sum_{i=1}^{N/2} x_{uv,i} (2 - x_{u, i-1} - x_{v, i}) HA=HA1+HA2+HA3+HA4+HA5+HA6 H_A = H_{A_1} + H_{A_2} + H_{A_3} + H_{A_4} + H_{A_5} + H_{A_6} HB=BuvEi=1N/2cuvxuv,i H_B = B \sum_{uv \in E} \sum_{i=1}^{N/2} c_{uv} x_{uv,i} H=HA+HB H = H_A + H_B

エネルギー関数の最小化により問題が解けるのは何故か #

上で定義したエネルギー関数はどのような制約を課しているのでしょうか。

まず、 HA1H_{A_1} について考えます。 これは、ただ一つの頂点が存在して深さが0となる、すなわち根であるという制約を課しています。

次に、 HA2H_{A_2} について考えます。 これは各頂点について、ちょうど一つの根からの深さ ii が存在するという制約を課しています。

次に、 HA3H_{A_3} について考えます。 これは各辺について、その辺を木 TT が含む時に、ある一つの根からの深さ ii が存在して、その辺がちょうど深さ ii の位置にあるという制約を課しています。

次に、 HA4H_{A_4} について考えます。 これは各頂点 vv と根からの各深さ ii について、その頂点が木においてその深さを持つ時に、ちょうど一つの vv を含む辺 uvEuv \in E が存在して、その辺が深さ ii を持って、 uu の方が vv よりも根に近いということが成り立つ、という制約を課しています。

次に、 HA5H_{A_5} について考えます。 これは各頂点 vv について、その頂点が木において次数 jj を持つ時に、その頂点を含む辺のうち木に含まれるものがちょうど jj 本であるという制約を課しています。

次に、 HA6H_{A_6} について考えます。 これは各辺 uvuv \in と根からの各深さ ii に対して、辺 uvuv uu の方が vv よりも根に近く深さ ii の状況で木に含まれているときに、 頂点 uu が木に深さ i1i-1 で含まれ、頂点 vv が木に深さ ii で含まれるという制約を課しています。

以上の6つの制約により、得られる木が次数制約を満たす全域木であることが保証されます。

最後に、 HBH_B について考えます。 これは、得られる木の重さを最小にするためのものです。 実際、これは得られる木の重みの総和になっています。 したがって、これを最小化することで重みを最小化していることになります。

もちろん、次数制約を満たす全域木であることは、最小性よりも優先する条件ですから、 A>Bmax(cuv)A > Bmax(c_{uv}) でなくてはいけません。

参考文献 #

A. Lucas, “Ising formulations of many NP problems” (open access)

Calendar 2019-01-07