1 二重確率行列と置換行列
定義 1.1. n n n 次実正方行列A = ( a i j ) A=(a_{ij}) A = ( a ij ) が二重確率行列 (doubly stochastic matrix ) であるとは、次の三条件を満たすことをいう。
すべてのi , j i,j i , j についてa i j ≥ 0 a_{ij}\ge0 a ij ≥ 0 である。
すべてのi i i について∑ j = 1 n a i j = 1 \sum_{j=1}^{n}a_{ij}=1 ∑ j = 1 n a ij = 1 である。
すべてのj j j について∑ i = 1 n a i j = 1 \sum_{i=1}^{n}a_{ij}=1 ∑ i = 1 n a ij = 1 である。
定義 1.2. 置換σ \sigma σ に対し、n n n 次実正方行列P σ = ( p i j ) P_\sigma=(p_{ij}) P σ = ( p ij ) を
p i j = { 1 ( j = σ ( i ) ) , 0 ( j ≠ σ ( i ) ) p_{ij}=\begin{cases}1 & (j=\sigma(i)),\\ 0 & (j\ne\sigma(i))\end{cases} p ij = { 1 0 ( j = σ ( i )) , ( j = σ ( i )) で定め、置換行列 (permutation matrix ) という。
命題 1.3. 任意の置換σ \sigma σ に対しP σ P_\sigma P σ は二重確率行列である。
証明. 成分は0 0 0 または1 1 1 であるから非負である。第i i i 行ではj = σ ( i ) j=\sigma(i) j = σ ( i ) のときだけ成分が1 1 1 であり、他は0 0 0 であるから行の和は1 1 1 である。第j j j 列では、p i j = 1 p_{ij}=1 p ij = 1 となるi i i はσ ( i ) = j \sigma(i)=j σ ( i ) = j を満たすもの、すなわちi = σ − 1 ( j ) i=\sigma^{-1}(j) i = σ − 1 ( j ) ただ一つであるから列の和も1 1 1 である。▨
2 正成分の定める二部グラフと Hall の条件
定義 2.1. n n n 次二重確率行列A = ( a i j ) A=(a_{ij}) A = ( a ij ) に対し、二部グラフG A = ( X ⊔ Y , E A ) G_A=(X\sqcup Y,E_A) G A = ( X ⊔ Y , E A ) を次で定める。X = { x 1 , … , x n } X=\{x_1,\dots,x_n\} X = { x 1 , … , x n } を行に対応する頂点集合、Y = { y 1 , … , y n } Y=\{y_1,\dots,y_n\} Y = { y 1 , … , y n } を列に対応する頂点集合とし、
E A = { x i y j : a i j > 0 } E_A=\{\,x_iy_j:\ a_{ij}>0\,\} E A = { x i y j : a ij > 0 } とする。G A G_A G A をA A A の正成分の二部グラフ (positive-entry bipartite graph ) という。
補題 2.2. n n n 次二重確率行列A A A に対し、任意のS ⊆ X S\subseteq X S ⊆ X について∣ N ( S ) ∣ ≥ ∣ S ∣ \lvert N(S)\rvert\ge\lvert S\rvert ∣ N ( S )∣ ≥ ∣ S ∣ が成り立つ。ここでN ( S ) N(S) N ( S ) はS S S の頂点に隣接するY Y Y の頂点の全体である。
証明. S ⊆ X S\subseteq X S ⊆ X をとり、I = { i : x i ∈ S } I=\{i:\ x_i\in S\} I = { i : x i ∈ S } 、J = { j : y j ∈ N ( S ) } J=\{j:\ y_j\in N(S)\} J = { j : y j ∈ N ( S )} と置く。∣ I ∣ = ∣ S ∣ \lvert I\rvert=\lvert S\rvert ∣ I ∣ = ∣ S ∣ かつ∣ J ∣ = ∣ N ( S ) ∣ \lvert J\rvert=\lvert N(S)\rvert ∣ J ∣ = ∣ N ( S )∣ である。
i ∈ I i\in I i ∈ I かつj ∉ J j\notin J j ∈ / J ならばa i j = 0 a_{ij}=0 a ij = 0 である。実際、a i j > 0 a_{ij}>0 a ij > 0 とするとx i y j ∈ E A x_iy_j\in E_A x i y j ∈ E A でありy j ∈ N ( S ) y_j\in N(S) y j ∈ N ( S ) となってj ∈ J j\in J j ∈ J に反する。したがって、I I I に属する行についての行和を足し合わせると
∣ I ∣ = ∑ i ∈ I ∑ j = 1 n a i j = ∑ i ∈ I ∑ j ∈ J a i j \lvert I\rvert=\sum_{i\in I}\sum_{j=1}^{n}a_{ij}=\sum_{i\in I}\sum_{j\in J}a_{ij} ∣ I ∣ = i ∈ I ∑ j = 1 ∑ n a ij = i ∈ I ∑ j ∈ J ∑ a ij となる。ここで最初の等号は行和がすべて1 1 1 であることによる。次に、成分がすべて非負であることから、J J J に属する列についてI I I に属する行だけの和を全体の列和で上から評価することができ、
∑ i ∈ I ∑ j ∈ J a i j = ∑ j ∈ J ∑ i ∈ I a i j ≤ ∑ j ∈ J ∑ i = 1 n a i j = ∣ J ∣ \sum_{i\in I}\sum_{j\in J}a_{ij}=\sum_{j\in J}\sum_{i\in I}a_{ij}\le\sum_{j\in J}\sum_{i=1}^{n}a_{ij}=\lvert J\rvert i ∈ I ∑ j ∈ J ∑ a ij = j ∈ J ∑ i ∈ I ∑ a ij ≤ j ∈ J ∑ i = 1 ∑ n a ij = ∣ J ∣ となる。最後の等号は列和がすべて1 1 1 であることによる。二つを合わせて∣ S ∣ = ∣ I ∣ ≤ ∣ J ∣ = ∣ N ( S ) ∣ \lvert S\rvert=\lvert I\rvert\le\lvert J\rvert=\lvert N(S)\rvert ∣ S ∣ = ∣ I ∣ ≤ ∣ J ∣ = ∣ N ( S )∣ を得る。▨
証明. 補題 2.2 によりG A G_A G A は Hall の条件を満たす。したがって§D2.12 定理 1.3 により、X X X のすべての頂点を被覆するマッチングM M M が存在する。M M M の各辺はちょうど一つのX X X の頂点と一つのY Y Y の頂点を端点にもち、M M M はマッチングであるから、x i x_i x i を被覆するM M M の辺はただ一本である。その辺をx i y σ ( i ) x_iy_{\sigma(i)} x i y σ ( i ) と書いて、写像σ \sigma σ を定める。
σ \sigma σ が単射であることを示す。σ ( i ) = σ ( i ′ ) = j \sigma(i)=\sigma(i')=j σ ( i ) = σ ( i ′ ) = j かつi ≠ i ′ i\ne i' i = i ′ とすると、M M M の二辺x i y j x_iy_j x i y j とx i ′ y j x_{i'}y_j x i ′ y j がともにy j y_j y j を端点にもち、M M M がマッチングであることに反する。{ 1 , … , n } \{1,\dots,n\} { 1 , … , n } から自分自身への単射は全単射であるから、σ \sigma σ は置換である。x i y σ ( i ) ∈ E A x_iy_{\sigma(i)}\in E_A x i y σ ( i ) ∈ E A であるからa i σ ( i ) > 0 a_{i\sigma(i)}>0 a iσ ( i ) > 0 が成り立つ。▨
3 置換行列の凸結合表示
定義 3.1. n n n 次実正方行列A = ( a i j ) A=(a_{ij}) A = ( a ij ) に対し、a i j > 0 a_{ij}>0 a ij > 0 を満たす添字の対( i , j ) (i,j) ( i , j ) の個数をp ( A ) p(A) p ( A ) と書き、A A A の正の成分の個数 (number of positive entries ) という。
命題 3.2. n n n 次二重確率行列A A A に対しp ( A ) ≥ n p(A)\ge n p ( A ) ≥ n が成り立つ。さらにp ( A ) = n p(A)=n p ( A ) = n であることと、A A A が置換行列であることは同値である。
証明. 各i i i について、第i i i 行の成分はすべて非負でありその和は1 1 1 であるから、a i j > 0 a_{ij}>0 a ij > 0 を満たすj j j が少なくとも一つ存在する。行ごとにそのような対を一つずつ選ぶと相異なるn n n 個の対が得られるのでp ( A ) ≥ n p(A)\ge n p ( A ) ≥ n である。
p ( A ) = n p(A)=n p ( A ) = n とする。上の評価で等号が成り立つのは、各行に正の成分がちょうど一つある場合に限る。第i i i 行の正の成分をa i τ ( i ) a_{i\tau(i)} a i τ ( i ) とすると、同じ行の他の成分は0 0 0 であり行和が1 1 1 であるからa i τ ( i ) = 1 a_{i\tau(i)}=1 a i τ ( i ) = 1 である。第j j j 列の和はτ ( i ) = j \tau(i)=j τ ( i ) = j を満たすi i i の個数に等しく、これが1 1 1 であるから、各j j j に対してτ ( i ) = j \tau(i)=j τ ( i ) = j を満たすi i i はちょうど一つである。よってτ \tau τ は置換であり、定義 1.2 によりA = P τ A=P_\tau A = P τ である。
逆にA = P τ A=P_\tau A = P τ とすると、正の成分は各行の( i , τ ( i ) ) (i,\tau(i)) ( i , τ ( i )) ただ一つずつであるからp ( A ) = n p(A)=n p ( A ) = n である。▨
3.1 証明方針
p ( A ) p(A) p ( A ) についての累積帰納法で示す。
p ( A ) p(A) p ( A ) の最小値は命題 3.2 によりn n n であり、そのときA A A 自身が置換行列である。これが基底の場合である。
p ( A ) > n p(A)>n p ( A ) > n のときは、補題 2.3 が与える置換σ \sigma σ をとり、係数を
λ = min 1 ≤ i ≤ n a i σ ( i ) \lambda=\min_{1\le i\le n}a_{i\sigma(i)} λ = 1 ≤ i ≤ n min a iσ ( i )
と定める。λ > 0 \lambda>0 λ > 0 である。まずλ < 1 \lambda<1 λ < 1 を示す。λ = 1 \lambda=1 λ = 1 とすると、各行でa i σ ( i ) = 1 a_{i\sigma(i)}=1 a iσ ( i ) = 1 となり、行和が1 1 1 で成分が非負であることから同じ行の他の成分はすべて0 0 0 となってp ( A ) = n p(A)=n p ( A ) = n になり、仮定に反するからである。
次にB = ( A − λ P σ ) / ( 1 − λ ) B=(A-\lambda P_\sigma)/(1-\lambda) B = ( A − λ P σ ) / ( 1 − λ ) が再び二重確率行列であることを、成分の非負性、行和、列和の順に確かめる。λ \lambda λ を最小値にとったことが非負性を保証する。さらに、最小値を与える位置の成分がB B B では0 0 0 になるのでp ( B ) ≤ p ( A ) − 1 p(B)\le p(A)-1 p ( B ) ≤ p ( A ) − 1 であり、帰納法の仮定をB B B へ適用することができる。最後に、B B B の表示に( 1 − λ ) (1-\lambda) ( 1 − λ ) を掛けてλ P σ \lambda P_\sigma λ P σ を加えるとA A A の表示が得られ、係数の総和がλ + ( 1 − λ ) = 1 \lambda+(1-\lambda)=1 λ + ( 1 − λ ) = 1 となることを確かめる。
定理 3.3 (Birkhoff–von Neumann の定理). A A A をn n n 次二重確率行列とする。このとき、有限個の置換σ 1 , … , σ m \sigma_1,\dots,\sigma_m σ 1 , … , σ m と正の実数λ 1 , … , λ m \lambda_1,\dots,\lambda_m λ 1 , … , λ m が存在して
A = ∑ k = 1 m λ k P σ k , ∑ k = 1 m λ k = 1 A=\sum_{k=1}^{m}\lambda_kP_{\sigma_k},
\qquad
\sum_{k=1}^{m}\lambda_k=1 A = k = 1 ∑ m λ k P σ k , k = 1 ∑ m λ k = 1 が成り立つ。
証明. 正の成分の個数p ( A ) p(A) p ( A ) についての累積帰納法で示す。命題 3.2 によりp ( A ) ≥ n p(A)\ge n p ( A ) ≥ n である。
p ( A ) = n p(A)=n p ( A ) = n の場合 。命題 3.2 によりA A A は置換行列P τ P_\tau P τ である。m = 1 m=1 m = 1 、λ 1 = 1 \lambda_1=1 λ 1 = 1 ととれば主張が成り立つ。
p ( A ) > n p(A)>n p ( A ) > n の場合 。正の成分の個数がp ( A ) p(A) p ( A ) より少ないすべての二重確率行列について主張が成り立つと仮定する。
補題 2.3 により、すべてのi i i についてa i σ ( i ) > 0 a_{i\sigma(i)}>0 a iσ ( i ) > 0 を満たす置換σ \sigma σ が存在する。λ = min i a i σ ( i ) \lambda=\min_i a_{i\sigma(i)} λ = min i a iσ ( i ) と置くとλ > 0 \lambda>0 λ > 0 である。また、行和が1 1 1 で成分が非負であるから各成分は1 1 1 以下でありλ ≤ 1 \lambda\le1 λ ≤ 1 である。λ = 1 \lambda=1 λ = 1 と仮定すると、各i i i についてa i σ ( i ) = 1 a_{i\sigma(i)}=1 a iσ ( i ) = 1 となり、第i i i 行の他の成分は非負で和が0 0 0 であるからすべて0 0 0 となる。すると各行の正の成分がちょうど一つとなりp ( A ) = n p(A)=n p ( A ) = n が成り立つので、p ( A ) > n p(A)>n p ( A ) > n という仮定に反する。よって0 < λ < 1 0<\lambda<1 0 < λ < 1 である。
B = ( b i j ) B=(b_{ij}) B = ( b ij ) を
b i j = a i j − λ ( P σ ) i j 1 − λ b_{ij}=\frac{a_{ij}-\lambda\,(P_\sigma)_{ij}}{1-\lambda} b ij = 1 − λ a ij − λ ( P σ ) ij で定める。B B B が二重確率行列であることを示す。j = σ ( i ) j=\sigma(i) j = σ ( i ) のときb i j = ( a i σ ( i ) − λ ) / ( 1 − λ ) b_{ij}=(a_{i\sigma(i)}-\lambda)/(1-\lambda) b ij = ( a iσ ( i ) − λ ) / ( 1 − λ ) であり、λ \lambda λ が最小値であるからa i σ ( i ) − λ ≥ 0 a_{i\sigma(i)}-\lambda\ge0 a iσ ( i ) − λ ≥ 0 となりb i j ≥ 0 b_{ij}\ge0 b ij ≥ 0 である。j ≠ σ ( i ) j\ne\sigma(i) j = σ ( i ) のときb i j = a i j / ( 1 − λ ) ≥ 0 b_{ij}=a_{ij}/(1-\lambda)\ge0 b ij = a ij / ( 1 − λ ) ≥ 0 である。行和については、命題 1.3 によりP σ P_\sigma P σ の第i i i 行の和が1 1 1 であるから
∑ j = 1 n b i j = 1 1 − λ ( ∑ j = 1 n a i j − λ ∑ j = 1 n ( P σ ) i j ) = 1 − λ 1 − λ = 1 \sum_{j=1}^{n}b_{ij}=\frac{1}{1-\lambda}\Bigl(\sum_{j=1}^{n}a_{ij}-\lambda\sum_{j=1}^{n}(P_\sigma)_{ij}\Bigr)=\frac{1-\lambda}{1-\lambda}=1 j = 1 ∑ n b ij = 1 − λ 1 ( j = 1 ∑ n a ij − λ j = 1 ∑ n ( P σ ) ij ) = 1 − λ 1 − λ = 1 となる。列和についても、P σ P_\sigma P σ の各列の和が1 1 1 であるから同じ計算により∑ i b i j = 1 \sum_{i}b_{ij}=1 ∑ i b ij = 1 となる。よってB B B は二重確率行列である。
次にp ( B ) < p ( A ) p(B)<p(A) p ( B ) < p ( A ) を示す。b i j > 0 b_{ij}>0 b ij > 0 ならばa i j − λ ( P σ ) i j > 0 a_{ij}-\lambda(P_\sigma)_{ij}>0 a ij − λ ( P σ ) ij > 0 であり、( P σ ) i j ≥ 0 (P_\sigma)_{ij}\ge0 ( P σ ) ij ≥ 0 かつλ > 0 \lambda>0 λ > 0 であるからa i j > 0 a_{ij}>0 a ij > 0 である。よってB B B の正の成分の位置はA A A の正の成分の位置に含まれる。さらに、最小値を与える添字i 0 i_0 i 0 、すなわちa i 0 σ ( i 0 ) = λ a_{i_0\sigma(i_0)}=\lambda a i 0 σ ( i 0 ) = λ を満たすi 0 i_0 i 0 をとると、b i 0 σ ( i 0 ) = 0 b_{i_0\sigma(i_0)}=0 b i 0 σ ( i 0 ) = 0 である一方a i 0 σ ( i 0 ) = λ > 0 a_{i_0\sigma(i_0)}=\lambda>0 a i 0 σ ( i 0 ) = λ > 0 である。よってp ( B ) ≤ p ( A ) − 1 p(B)\le p(A)-1 p ( B ) ≤ p ( A ) − 1 である。
帰納法の仮定をB B B へ適用すると、置換τ 1 , … , τ l \tau_1,\dots,\tau_l τ 1 , … , τ l と正の実数μ 1 , … , μ l \mu_1,\dots,\mu_l μ 1 , … , μ l が存在して
B = ∑ k = 1 l μ k P τ k , ∑ k = 1 l μ k = 1 B=\sum_{k=1}^{l}\mu_kP_{\tau_k},
\qquad
\sum_{k=1}^{l}\mu_k=1 B = k = 1 ∑ l μ k P τ k , k = 1 ∑ l μ k = 1 が成り立つ。b i j b_{ij} b ij の定め方から、すべてのi i i とj j j についてa i j = λ ( P σ ) i j + ( 1 − λ ) b i j a_{ij}=\lambda\,(P_\sigma)_{ij}+(1-\lambda)b_{ij} a ij = λ ( P σ ) ij + ( 1 − λ ) b ij が成り立つので、§D3.1 定義 2.1 と§D3.1 定義 1.1 の行列の相等によりA = λ P σ + ( 1 − λ ) B A=\lambda P_\sigma+(1-\lambda)B A = λ P σ + ( 1 − λ ) B である。この右辺へB B B の表示を代入し、§D3.1 命題 2.3 の実数倍についての法則によって係数をまとめると
A = λ P σ + ∑ k = 1 l ( 1 − λ ) μ k P τ k A=\lambda P_\sigma+\sum_{k=1}^{l}(1-\lambda)\mu_kP_{\tau_k} A = λ P σ + k = 1 ∑ l ( 1 − λ ) μ k P τ k となる。λ > 0 \lambda>0 λ > 0 かつ( 1 − λ ) μ k > 0 (1-\lambda)\mu_k>0 ( 1 − λ ) μ k > 0 であり、係数の総和は
λ + ∑ k = 1 l ( 1 − λ ) μ k = λ + ( 1 − λ ) ∑ k = 1 l μ k = λ + ( 1 − λ ) = 1 \lambda+\sum_{k=1}^{l}(1-\lambda)\mu_k=\lambda+(1-\lambda)\sum_{k=1}^{l}\mu_k=\lambda+(1-\lambda)=1 λ + k = 1 ∑ l ( 1 − λ ) μ k = λ + ( 1 − λ ) k = 1 ∑ l μ k = λ + ( 1 − λ ) = 1 である。よってm = l + 1 m=l+1 m = l + 1 として主張の形の表示が得られる。▨
4 線形割当問題
系 4.1. C = ( c i j ) C=(c_{ij}) C = ( c ij ) を成分が実数であるn n n 次正方行列とし、
γ = min σ ∑ i = 1 n c i σ ( i ) \gamma=\min_{\sigma}\ \sum_{i=1}^{n}c_{i\sigma(i)} γ = σ min i = 1 ∑ n c iσ ( i ) と置く。ここで最小は置換σ \sigma σ の全体をわたる。このとき、任意のn n n 次二重確率行列A = ( a i j ) A=(a_{ij}) A = ( a ij ) に対し
∑ i = 1 n ∑ j = 1 n c i j a i j ≥ γ \sum_{i=1}^{n}\sum_{j=1}^{n}c_{ij}a_{ij}\ \ge\ \gamma i = 1 ∑ n j = 1 ∑ n c ij a ij ≥ γ が成り立ち、γ \gamma γ を与える置換σ ∗ \sigma^{*} σ ∗ の置換行列A = P σ ∗ A=P_{\sigma^{*}} A = P σ ∗ で等号が成り立つ。したがって、二重確率行列の全体における∑ i , j c i j a i j \sum_{i,j}c_{ij}a_{ij} ∑ i , j c ij a ij の最小値は存在してγ \gamma γ に等しく、置換行列で達成される。
証明. 置換の個数はn ! n! n ! 個で有限であるからγ \gamma γ は定まり、γ \gamma γ を与える置換σ ∗ \sigma^{*} σ ∗ が存在する。任意の置換σ \sigma σ について、定義 1.2 により
∑ i = 1 n ∑ j = 1 n c i j ( P σ ) i j = ∑ i = 1 n c i σ ( i ) \sum_{i=1}^{n}\sum_{j=1}^{n}c_{ij}(P_\sigma)_{ij}=\sum_{i=1}^{n}c_{i\sigma(i)} i = 1 ∑ n j = 1 ∑ n c ij ( P σ ) ij = i = 1 ∑ n c iσ ( i ) が成り立つ。σ = σ ∗ \sigma=\sigma^{*} σ = σ ∗ ととると、この式の左辺はγ \gamma γ に等しい。命題 1.3 によりP σ ∗ P_{\sigma^{*}} P σ ∗ は二重確率行列であるから、等号を与える二重確率行列が存在する。
任意の二重確率行列A A A をとる。定理 3.3 により、置換σ 1 , … , σ m \sigma_1,\dots,\sigma_m σ 1 , … , σ m と正の実数λ 1 , … , λ m \lambda_1,\dots,\lambda_m λ 1 , … , λ m が存在してA = ∑ k λ k P σ k A=\sum_k\lambda_kP_{\sigma_k} A = ∑ k λ k P σ k かつ∑ k λ k = 1 \sum_k\lambda_k=1 ∑ k λ k = 1 である。各成分についてa i j = ∑ k λ k ( P σ k ) i j a_{ij}=\sum_k\lambda_k(P_{\sigma_k})_{ij} a ij = ∑ k λ k ( P σ k ) ij であるから、和の順序を入れ替えて
∑ i , j c i j a i j = ∑ k = 1 m λ k ∑ i , j c i j ( P σ k ) i j = ∑ k = 1 m λ k ∑ i = 1 n c i σ k ( i ) \sum_{i,j}c_{ij}a_{ij}
=\sum_{k=1}^{m}\lambda_k\sum_{i,j}c_{ij}(P_{\sigma_k})_{ij}
=\sum_{k=1}^{m}\lambda_k\sum_{i=1}^{n}c_{i\sigma_k(i)} i , j ∑ c ij a ij = k = 1 ∑ m λ k i , j ∑ c ij ( P σ k ) ij = k = 1 ∑ m λ k i = 1 ∑ n c i σ k ( i ) となる。ここで和はすべて有限和であるから、順序の入れ替えは加法の結合法則と交換法則による。各k k k について∑ i c i σ k ( i ) ≥ γ \sum_i c_{i\sigma_k(i)}\ge\gamma ∑ i c i σ k ( i ) ≥ γ であり、λ k > 0 \lambda_k>0 λ k > 0 であるから
∑ i , j c i j a i j ≥ ∑ k = 1 m λ k γ = γ ∑ k = 1 m λ k = γ \sum_{i,j}c_{ij}a_{ij}\ \ge\ \sum_{k=1}^{m}\lambda_k\gamma=\gamma\sum_{k=1}^{m}\lambda_k=\gamma i , j ∑ c ij a ij ≥ k = 1 ∑ m λ k γ = γ k = 1 ∑ m λ k = γ を得る。▨
5 検算例
例 5.1 (3 次の二重確率行列の分解と割当の最適値).
A = ( 1 / 2 1 / 2 0 0 1 / 2 1 / 2 1 / 2 0 1 / 2 ) A=\begin{pmatrix}
1/2 & 1/2 & 0\\
0 & 1/2 & 1/2\\
1/2 & 0 & 1/2
\end{pmatrix} A = 1/2 0 1/2 1/2 1/2 0 0 1/2 1/2 と置く。行の和は順に1 / 2 + 1 / 2 + 0 = 1 1/2+1/2+0=1 1/2 + 1/2 + 0 = 1 、0 + 1 / 2 + 1 / 2 = 1 0+1/2+1/2=1 0 + 1/2 + 1/2 = 1 、1 / 2 + 0 + 1 / 2 = 1 1/2+0+1/2=1 1/2 + 0 + 1/2 = 1 である。列の和は順に1 / 2 + 0 + 1 / 2 = 1 1/2+0+1/2=1 1/2 + 0 + 1/2 = 1 、1 / 2 + 1 / 2 + 0 = 1 1/2+1/2+0=1 1/2 + 1/2 + 0 = 1 、0 + 1 / 2 + 1 / 2 = 1 0+1/2+1/2=1 0 + 1/2 + 1/2 = 1 である。成分はすべて非負であるからA A A は二重確率行列であり、正の成分の個数はp ( A ) = 6 p(A)=6 p ( A ) = 6 である。
Hall の条件の確認 。正成分の二部グラフG A G_A G A の辺はx 1 y 1 x_1y_1 x 1 y 1 、x 1 y 2 x_1y_2 x 1 y 2 、x 2 y 2 x_2y_2 x 2 y 2 、x 2 y 3 x_2y_3 x 2 y 3 、x 3 y 1 x_3y_1 x 3 y 1 、x 3 y 3 x_3y_3 x 3 y 3 である。一元集合の近傍はいずれも二元であり、∣ N ( S ) ∣ ≥ ∣ S ∣ \lvert N(S)\rvert\ge\lvert S\rvert ∣ N ( S )∣ ≥ ∣ S ∣ が成り立つ。二元集合については、N ( { x 1 , x 2 } ) = { y 1 , y 2 , y 3 } N(\{x_1,x_2\})=\{y_1,y_2,y_3\} N ({ x 1 , x 2 }) = { y 1 , y 2 , y 3 } 、N ( { x 1 , x 3 } ) = { y 1 , y 2 , y 3 } N(\{x_1,x_3\})=\{y_1,y_2,y_3\} N ({ x 1 , x 3 }) = { y 1 , y 2 , y 3 } 、N ( { x 2 , x 3 } ) = { y 1 , y 2 , y 3 } N(\{x_2,x_3\})=\{y_1,y_2,y_3\} N ({ x 2 , x 3 }) = { y 1 , y 2 , y 3 } であり、いずれも三元である。三元集合についてはN ( X ) = { y 1 , y 2 , y 3 } N(X)=\{y_1,y_2,y_3\} N ( X ) = { y 1 , y 2 , y 3 } である。よって Hall の条件が成り立つ。
置換の取り出しと分解 。恒等置換σ = i d \sigma=\mathrm{id} σ = id についてa 11 = a 22 = a 33 = 1 / 2 > 0 a_{11}=a_{22}=a_{33}=1/2>0 a 11 = a 22 = a 33 = 1/2 > 0 であるから、σ \sigma σ は補題 2.3 の条件を満たす。λ = min { 1 / 2 , 1 / 2 , 1 / 2 } = 1 / 2 \lambda=\min\{1/2,1/2,1/2\}=1/2 λ = min { 1/2 , 1/2 , 1/2 } = 1/2 である。恒等置換の置換行列は、§D3.1 定義 3.6 の単位行列にほかならない。3 3 3 次の単位行列をI I I と書くとP i d = I P_{\mathrm{id}}=I P id = I であり、
B = A − 1 2 I 1 − 1 2 = 2 ( 0 1 / 2 0 0 0 1 / 2 1 / 2 0 0 ) = ( 0 1 0 0 0 1 1 0 0 ) B=\frac{A-\tfrac12 I}{1-\tfrac12}
=2\begin{pmatrix}
0 & 1/2 & 0\\
0 & 0 & 1/2\\
1/2 & 0 & 0
\end{pmatrix}
=\begin{pmatrix}
0 & 1 & 0\\
0 & 0 & 1\\
1 & 0 & 0
\end{pmatrix} B = 1 − 2 1 A − 2 1 I = 2 0 0 1/2 1/2 0 0 0 1/2 0 = 0 0 1 1 0 0 0 1 0 となる。B B B の正の成分の個数は3 3 3 であり、p ( A ) = 6 p(A)=6 p ( A ) = 6 より真に小さい。B B B はτ ( 1 ) = 2 \tau(1)=2 τ ( 1 ) = 2 、τ ( 2 ) = 3 \tau(2)=3 τ ( 2 ) = 3 、τ ( 3 ) = 1 \tau(3)=1 τ ( 3 ) = 1 で定まる置換τ \tau τ の置換行列P τ P_\tau P τ である。よって
A = 1 2 P i d + 1 2 P τ A=\frac12 P_{\mathrm{id}}+\frac12 P_\tau A = 2 1 P id + 2 1 P τ であり、係数はいずれも正でその和は1 1 1 である。右辺を成分ごとに計算すると、第1 1 1 行は( 1 / 2 , 1 / 2 , 0 ) (1/2,\,1/2,\,0) ( 1/2 , 1/2 , 0 ) 、第2 2 2 行は( 0 , 1 / 2 , 1 / 2 ) (0,\,1/2,\,1/2) ( 0 , 1/2 , 1/2 ) 、第3 3 3 行は( 1 / 2 , 0 , 1 / 2 ) (1/2,\,0,\,1/2) ( 1/2 , 0 , 1/2 ) となりA A A に一致する。
割当の最適値 。費用行列を
C = ( 1 2 3 2 3 1 3 1 2 ) C=\begin{pmatrix}1&2&3\\ 2&3&1\\ 3&1&2\end{pmatrix} C = 1 2 3 2 3 1 3 1 2 とする。六つの置換について総費用を数え上げる。恒等置換ではc 11 + c 22 + c 33 = 1 + 3 + 2 = 6 c_{11}+c_{22}+c_{33}=1+3+2=6 c 11 + c 22 + c 33 = 1 + 3 + 2 = 6 、τ \tau τ ではc 12 + c 23 + c 31 = 2 + 1 + 3 = 6 c_{12}+c_{23}+c_{31}=2+1+3=6 c 12 + c 23 + c 31 = 2 + 1 + 3 = 6 、1 1 1 と2 2 2 を入れ替える置換ではc 12 + c 21 + c 33 = 2 + 2 + 2 = 6 c_{12}+c_{21}+c_{33}=2+2+2=6 c 12 + c 21 + c 33 = 2 + 2 + 2 = 6 、1 1 1 と3 3 3 を入れ替える置換ではc 13 + c 22 + c 31 = 3 + 3 + 3 = 9 c_{13}+c_{22}+c_{31}=3+3+3=9 c 13 + c 22 + c 31 = 3 + 3 + 3 = 9 、2 2 2 と3 3 3 を入れ替える置換ではc 11 + c 23 + c 32 = 1 + 1 + 1 = 3 c_{11}+c_{23}+c_{32}=1+1+1=3 c 11 + c 23 + c 32 = 1 + 1 + 1 = 3 、残る置換1 ↦ 3 , 2 ↦ 1 , 3 ↦ 2 1\mapsto3,\ 2\mapsto1,\ 3\mapsto2 1 ↦ 3 , 2 ↦ 1 , 3 ↦ 2 ではc 13 + c 21 + c 32 = 3 + 2 + 1 = 6 c_{13}+c_{21}+c_{32}=3+2+1=6 c 13 + c 21 + c 32 = 3 + 2 + 1 = 6 である。よってγ = 3 \gamma=3 γ = 3 である。
上のA A A については
∑ i , j c i j a i j = 1 2 ( c 11 + c 12 + c 22 + c 23 + c 31 + c 33 ) = 1 2 ( 1 + 2 + 3 + 1 + 3 + 2 ) = 1 2 ⋅ 12 = 6 \sum_{i,j}c_{ij}a_{ij}=\tfrac12(c_{11}+c_{12}+c_{22}+c_{23}+c_{31}+c_{33})=\tfrac12(1+2+3+1+3+2)=\tfrac12\cdot12=6 i , j ∑ c ij a ij = 2 1 ( c 11 + c 12 + c 22 + c 23 + c 31 + c 33 ) = 2 1 ( 1 + 2 + 3 + 1 + 3 + 2 ) = 2 1 ⋅ 12 = 6 であり、分解A = 1 2 P i d + 1 2 P τ A=\tfrac12P_{\mathrm{id}}+\tfrac12P_\tau A = 2 1 P id + 2 1 P τ から得られる1 2 ⋅ 6 + 1 2 ⋅ 6 = 6 \tfrac12\cdot6+\tfrac12\cdot6=6 2 1 ⋅ 6 + 2 1 ⋅ 6 = 6 と一致する。6 ≥ 3 = γ 6\ge3=\gamma 6 ≥ 3 = γ であり、系 4.1 の不等式と整合する。
6 演習
問題 6.1.
補題 2.2 の証明で用いた二つの等号、すなわち行和による∣ I ∣ \lvert I\rvert ∣ I ∣ の表示と列和による∣ J ∣ \lvert J\rvert ∣ J ∣ の表示を、n = 3 n=3 n = 3 の具体的な二重確率行列と∣ S ∣ = 2 \lvert S\rvert=2 ∣ S ∣ = 2 の場合について書き下し、不等号がどこで生じるかを確かめよ。
定理 3.3 の証明でλ < 1 \lambda<1 λ < 1 を示す段を省くと、その後の議論のどこが成り立たなくなるかを述べよ。
定理 3.3 の証明で、λ \lambda λ をmin i a i σ ( i ) \min_i a_{i\sigma(i)} min i a iσ ( i ) ではなく、より小さい正の数にとった場合を考える。B B B が二重確率行列であることは保たれるが、証明のどの主張が成り立たなくなるかを指摘し、手続きが有限回で終わらない可能性を説明せよ。
命題 3.2 の証明において、行和がすべて1 1 1 であることと列和がすべて1 1 1 であることを、それぞれどこで用いたかを指摘せよ。さらに、行和の条件だけを課した非負行列で、p ( A ) = n p(A)=n p ( A ) = n を満たすが置換行列ではない例を挙げよ。
系 4.1 の証明を、最小化ではなく最大化について書き直せ。すなわち、任意の二重確率行列A A A について∑ i , j c i j a i j \sum_{i,j}c_{ij}a_{ij} ∑ i , j c ij a ij が置換にわたる最大値以下であることを、同じ線形性の議論で示せ。
7 つまずいたら
二重確率行列の条件は行和と列和の両方についてである(定義 1.1 )。行和だけが1 1 1 である非負行列は、正成分の二部グラフが Hall の条件を満たすとは限らず、補題 2.3 は成り立たない。
証明に用いるのは Hall の定理であって、最大マッチングと最小頂点被覆を結ぶ定理ではない。必要なのは行の全体を被覆するマッチングの存在だけであり、それは§D2.12 定理 1.3 が与える。
係数λ \lambda λ を対応する成分の最小値にとることが、二つの役割を果たす。第一に、引いた後の成分が非負に保たれる。第二に、最小値を与える位置の成分が0 0 0 になり、正成分の個数が真に減る。どちらか一方だけでは帰納法が回らない。
定理 3.3 の表示は一意ではない(注意 3.4 )。ある表示に現れない置換が別の表示に現れることがある。
8 扱った範囲と次の記事
本記事は、有限次の二重確率行列の正成分が定める二部グラフが Hall の条件を満たすことを行和と列和の数え上げによって証明し、置換行列を逐次取り出す手続きが有限回で終わることを正成分の個数についての帰納法で示して、置換行列の凸結合表示を証明した。さらに、線形割当問題の最適値が置換行列で達成されることを、凸結合表示と目的関数の線形性から導いた。表示の一意性、二重確率行列全体が定める多面体の面の構造、および無限次への拡張は扱っていない。次の記事では、有限マトロイドの基底の補集合によって双対マトロイドを定義し、双対階数公式と回路および余回路の対応を証明する。