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.
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.
n n n 次二重確率行列A A A と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 )} と置く。∣ N ( S ) ∣ − ∣ S ∣ = ∑ i ∉ I ∑ j ∈ J a i j \lvert N(S)\rvert-\lvert S\rvert=\sum_{i\notin I}\sum_{j\in J}a_{ij} ∣ N ( S )∣ − ∣ S ∣ = ∑ i ∈ / I ∑ j ∈ J a ij を示せ。さらに、空でなくX X X と異なるS ⊆ X S\subseteq X S ⊆ X であって∣ N ( S ) ∣ = ∣ S ∣ \lvert N(S)\rvert=\lvert S\rvert ∣ N ( S )∣ = ∣ S ∣ を満たすものが存在することと、ある整数1 ≤ k ≤ n − 1 1\le k\le n-1 1 ≤ k ≤ n − 1 と置換α , β \alpha,\beta α , β について、b i j = a α ( i ) β ( j ) b_{ij}=a_{\alpha(i)\beta(j)} b ij = a α ( i ) β ( j ) で定まる行列B = ( b i j ) B=(b_{ij}) B = ( b ij ) が「i ≤ k < j i\le k<j i ≤ k < j またはj ≤ k < i j\le k<i j ≤ k < i ならばb i j = 0 b_{ij}=0 b ij = 0 」を満たすことが同値であることを示せ。
n n n 次二重確率行列A A A について、定理 3.3 の表示を、置換の個数m m m がp ( A ) − n + 1 p(A)-n+1 p ( A ) − n + 1 以下であるようにとることができることを示せ。これを用いて、すべてのn n n 次二重確率行列はn 2 − n + 1 n^{2}-n+1 n 2 − n + 1 個以下の置換行列の凸結合として表されることを示せ。
p ( A ) > n p(A)>n p ( A ) > n であるn n n 次二重確率行列A A A と、補題 2.3 による置換σ \sigma σ をとり、0 < λ ′ < min i a i σ ( i ) 0<\lambda'<\min_ia_{i\sigma(i)} 0 < λ ′ < min i a iσ ( i ) を満たす実数λ ′ \lambda' λ ′ に対してB ′ = ( A − λ ′ P σ ) / ( 1 − λ ′ ) B'=(A-\lambda'P_\sigma)/(1-\lambda') B ′ = ( A − λ ′ P σ ) / ( 1 − λ ′ ) と置く。B ′ B' B ′ が二重確率行列でありp ( B ′ ) = p ( A ) p(B')=p(A) p ( B ′ ) = p ( A ) を満たすことを示せ。さらに、n = 2 n=2 n = 2 の二重確率行列A 0 A_0 A 0 から出発し、各回でこの形の置き換えA k ↦ A k + 1 A_k\mapsto A_{k+1} A k ↦ A k + 1 を行って得られる列A 0 , A 1 , A 2 , … A_0,A_1,A_2,\dots A 0 , A 1 , A 2 , … であって、どのA k A_k A k も置換行列でないものを作れ。
n ≥ 2 n\ge2 n ≥ 2 とする。成分が非負で行和がすべて1 1 1 であるn n n 次正方行列であって、p ( A ) = n p(A)=n p ( A ) = n を満たすが置換行列でないものの例を挙げよ。さらに、n n n 次二重確率行列A A A が置換行列でないならばp ( A ) ≥ n + 2 p(A)\ge n+2 p ( A ) ≥ n + 2 であることを示し、p ( A ) = n + 2 p(A)=n+2 p ( A ) = n + 2 を満たすn n 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 が置換にわたる最大値以下であることを、同じ線形性の議論で示せ。
解答. (1) を示す。補題 2.2 の証明のとおり、i ∈ I i\in I i ∈ I かつj ∉ J j\notin J j ∈ / J ならばa i j = 0 a_{ij}=0 a ij = 0 であり、行和がすべて1 1 1 であるから∣ I ∣ = ∑ i ∈ I ∑ j ∈ J a i j \lvert I\rvert=\sum_{i\in I}\sum_{j\in J}a_{ij} ∣ I ∣ = ∑ i ∈ I ∑ j ∈ J a ij である。列和がすべて1 1 1 であるから
∣ J ∣ = ∑ j ∈ J ∑ i = 1 n a i j = ∑ i ∈ I ∑ j ∈ J a i j + ∑ i ∉ I ∑ j ∈ J a i j \lvert J\rvert=\sum_{j\in J}\sum_{i=1}^{n}a_{ij}=\sum_{i\in I}\sum_{j\in J}a_{ij}+\sum_{i\notin I}\sum_{j\in J}a_{ij} ∣ J ∣ = j ∈ J ∑ i = 1 ∑ n a ij = i ∈ I ∑ j ∈ J ∑ a ij + i ∈ / I ∑ j ∈ J ∑ a ij である。∣ S ∣ = ∣ I ∣ \lvert S\rvert=\lvert I\rvert ∣ S ∣ = ∣ I ∣ 、∣ N ( S ) ∣ = ∣ J ∣ \lvert N(S)\rvert=\lvert J\rvert ∣ N ( S )∣ = ∣ J ∣ であるから、二つの等式の差をとって∣ N ( S ) ∣ − ∣ S ∣ = ∑ i ∉ I ∑ j ∈ J a i j \lvert N(S)\rvert-\lvert S\rvert=\sum_{i\notin I}\sum_{j\in J}a_{ij} ∣ N ( S )∣ − ∣ S ∣ = ∑ i ∈ / I ∑ j ∈ J a ij を得る。
空でなくX X X と異なるS S S が∣ N ( S ) ∣ = ∣ S ∣ \lvert N(S)\rvert=\lvert S\rvert ∣ N ( S )∣ = ∣ S ∣ を満たすとし、k = ∣ S ∣ k=\lvert S\rvert k = ∣ S ∣ と置くと1 ≤ k ≤ n − 1 1\le k\le n-1 1 ≤ k ≤ n − 1 である。前段の等式の右辺は非負の数の和であり0 0 0 に等しいので、i ∉ I i\notin I i ∈ / I かつj ∈ J j\in J j ∈ J ならばa i j = 0 a_{ij}=0 a ij = 0 である。またi ∈ I i\in I i ∈ I かつj ∉ J j\notin J j ∈ / J ならばa i j = 0 a_{ij}=0 a ij = 0 である。∣ I ∣ = ∣ J ∣ = k \lvert I\rvert=\lvert J\rvert=k ∣ I ∣ = ∣ J ∣ = k であるから、α ( { 1 , … , k } ) = I \alpha(\{1,\dots,k\})=I α ({ 1 , … , k }) = I を満たす置換α \alpha α とβ ( { 1 , … , k } ) = J \beta(\{1,\dots,k\})=J β ({ 1 , … , k }) = J を満たす置換β \beta β が存在する。このときα ( { k + 1 , … , n } ) \alpha(\{k+1,\dots,n\}) α ({ k + 1 , … , n }) はI I I の補集合であり、β ( { k + 1 , … , n } ) \beta(\{k+1,\dots,n\}) β ({ k + 1 , … , n }) はJ J J の補集合である。i ≤ k < j i\le k<j i ≤ k < j ならばα ( i ) ∈ I \alpha(i)\in I α ( i ) ∈ I 、β ( j ) ∉ J \beta(j)\notin J β ( j ) ∈ / J であるからb i j = a α ( i ) β ( j ) = 0 b_{ij}=a_{\alpha(i)\beta(j)}=0 b ij = a α ( i ) β ( j ) = 0 である。j ≤ k < i j\le k<i j ≤ k < i ならばα ( i ) ∉ I \alpha(i)\notin I α ( i ) ∈ / I 、β ( j ) ∈ J \beta(j)\in J β ( j ) ∈ J であるからb i j = 0 b_{ij}=0 b ij = 0 である。
逆に、1 ≤ k ≤ n − 1 1\le k\le n-1 1 ≤ k ≤ n − 1 と置換α , β \alpha,\beta α , β についてB B B が条件を満たすとする。I = α ( { 1 , … , k } ) I=\alpha(\{1,\dots,k\}) I = α ({ 1 , … , k }) 、S = { x i : i ∈ I } S=\{x_i:\ i\in I\} S = { x i : i ∈ I } と置くと、∣ S ∣ = k \lvert S\rvert=k ∣ S ∣ = k であるからS S S は空でなくX X X と異なる。i ∈ I i\in I i ∈ I とj j j について、i = α ( i ′ ) i=\alpha(i') i = α ( i ′ ) (i ′ ≤ k i'\le k i ′ ≤ k )と書きj ′ = β − 1 ( j ) j'=\beta^{-1}(j) j ′ = β − 1 ( j ) と置くとa i j = b i ′ j ′ a_{ij}=b_{i'j'} a ij = b i ′ j ′ であり、j ′ > k j'>k j ′ > k ならばこれは0 0 0 である。よってy j ∈ N ( S ) y_j\in N(S) y j ∈ N ( S ) ならばj ∈ β ( { 1 , … , k } ) j\in\beta(\{1,\dots,k\}) j ∈ β ({ 1 , … , k }) であり、∣ N ( S ) ∣ ≤ k = ∣ S ∣ \lvert N(S)\rvert\le k=\lvert S\rvert ∣ N ( S )∣ ≤ k = ∣ S ∣ である。補題 2.2 により∣ N ( S ) ∣ ≥ ∣ S ∣ \lvert N(S)\rvert\ge\lvert S\rvert ∣ N ( S )∣ ≥ ∣ S ∣ であるから∣ N ( S ) ∣ = ∣ S ∣ \lvert N(S)\rvert=\lvert S\rvert ∣ N ( S )∣ = ∣ S ∣ である。
(2) を、p ( A ) p(A) p ( A ) についての累積帰納法で示す。p ( A ) = n p(A)=n p ( A ) = n のとき、命題 3.2 によりA A A はある置換τ \tau τ の置換行列P τ P_\tau P τ であり、m = 1 = p ( A ) − n + 1 m=1=p(A)-n+1 m = 1 = p ( A ) − n + 1 の表示A = 1 ⋅ P τ A=1\cdot P_\tau A = 1 ⋅ P τ がある。p ( A ) > n p(A)>n p ( A ) > n とし、正の成分の個数がp ( A ) p(A) p ( A ) より少ないすべての二重確率行列B B B について、置換の個数がp ( B ) − n + 1 p(B)-n+1 p ( B ) − n + 1 以下の表示が存在すると仮定する。定理 3.3 の証明で示したとおり、補題 2.3 による置換σ \sigma σ とλ = min i a i σ ( i ) \lambda=\min_ia_{i\sigma(i)} λ = min i a iσ ( i ) について0 < λ < 1 0<\lambda<1 0 < λ < 1 であり、B = ( A − λ P σ ) / ( 1 − λ ) B=(A-\lambda P_\sigma)/(1-\lambda) B = ( A − λ P σ ) / ( 1 − λ ) は二重確率行列であってp ( B ) ≤ p ( A ) − 1 p(B)\le p(A)-1 p ( B ) ≤ p ( A ) − 1 を満たす。帰納法の仮定により、置換τ 1 , … , τ l \tau_1,\dots,\tau_l τ 1 , … , τ l (l ≤ p ( B ) − n + 1 ≤ p ( A ) − n l\le p(B)-n+1\le p(A)-n l ≤ p ( B ) − n + 1 ≤ p ( A ) − n )と、総和が1 1 1 である正の実数μ 1 , … , μ l \mu_1,\dots,\mu_l μ 1 , … , μ l が存在してB = ∑ k μ k P τ k B=\sum_k\mu_kP_{\tau_k} B = ∑ k μ k P τ k である。同じ証明のとおり
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 は係数が正で総和が1 1 1 の表示であり、置換の個数はl + 1 ≤ p ( A ) − n + 1 l+1\le p(A)-n+1 l + 1 ≤ p ( A ) − n + 1 である。よって主張が成り立つ。
成分の位置はn 2 n^{2} n 2 個であるからp ( A ) ≤ n 2 p(A)\le n^{2} p ( A ) ≤ n 2 であり、すべてのn n n 次二重確率行列はn 2 − n + 1 n^{2}-n+1 n 2 − n + 1 個以下の置換行列P σ k P_{\sigma_k} P σ k と正の係数λ k \lambda_k λ k (∑ k λ k = 1 \sum_k\lambda_k=1 ∑ k λ k = 1 )によってA = ∑ k λ k P σ k A=\sum_k\lambda_kP_{\sigma_k} A = ∑ k λ k P σ k と表される。たとえば例 5.1 の行列ではp ( A ) = 6 p(A)=6 p ( A ) = 6 、n = 3 n=3 n = 3 であり、上界は4 4 4 である。実際に得られた表示の置換の個数は2 2 2 であり、上界の範囲に収まる。
(3) を示す。λ = min i a i σ ( i ) \lambda=\min_ia_{i\sigma(i)} λ = min i a iσ ( i ) と置く。定理 3.3 の証明により、p ( A ) > n p(A)>n p ( A ) > n のときλ < 1 \lambda<1 λ < 1 であるから0 < λ ′ < λ < 1 0<\lambda'<\lambda<1 0 < λ ′ < λ < 1 であり、1 − λ ′ > 0 1-\lambda'>0 1 − λ ′ > 0 である。b i j ′ = ( a i j − λ ′ ( P σ ) i j ) / ( 1 − λ ′ ) b'_{ij}=(a_{ij}-\lambda'(P_\sigma)_{ij})/(1-\lambda') b ij ′ = ( a ij − λ ′ ( P σ ) ij ) / ( 1 − λ ′ ) と置く。j = σ ( i ) j=\sigma(i) j = σ ( i ) のときa i σ ( i ) − λ ′ ≥ λ − λ ′ > 0 a_{i\sigma(i)}-\lambda'\ge\lambda-\lambda'>0 a iσ ( i ) − λ ′ ≥ λ − λ ′ > 0 であるからb i σ ( i ) ′ > 0 b'_{i\sigma(i)}>0 b iσ ( i ) ′ > 0 であり、a i σ ( i ) > 0 a_{i\sigma(i)}>0 a iσ ( i ) > 0 でもある。j ≠ σ ( i ) j\ne\sigma(i) j = σ ( i ) のときb i j ′ = a i j / ( 1 − λ ′ ) b'_{ij}=a_{ij}/(1-\lambda') b ij ′ = a ij / ( 1 − λ ′ ) であるから、b i j ′ ≥ 0 b'_{ij}\ge0 b ij ′ ≥ 0 であり、b i j ′ > 0 b'_{ij}>0 b ij ′ > 0 とa i j > 0 a_{ij}>0 a ij > 0 は同値である。よってB ′ B' B ′ の成分は非負であり、B ′ B' B ′ の正の成分の位置はA A A の正の成分の位置に一致してp ( B ′ ) = p ( A ) p(B')=p(A) p ( B ′ ) = p ( A ) である。命題 1.3 によりP σ P_\sigma P σ の行和と列和は1 1 1 であるから、第i i i 行の和は( 1 − λ ′ ) / ( 1 − λ ′ ) = 1 (1-\lambda')/(1-\lambda')=1 ( 1 − λ ′ ) / ( 1 − λ ′ ) = 1 であり、列和も同様に1 1 1 である。よってB ′ B' B ′ は二重確率行列である。
列を作る。0 < a < 1 0<a<1 0 < a < 1 に対して
M ( a ) = ( a 1 − a 1 − a a ) M(a)=\begin{pmatrix}a&1-a\\ 1-a&a\end{pmatrix} M ( a ) = ( a 1 − a 1 − a a ) と置くと、M ( a ) M(a) M ( a ) は二重確率行列でありp ( M ( a ) ) = 4 > 2 p(M(a))=4>2 p ( M ( a )) = 4 > 2 であるから、命題 3.2 により置換行列ではない。a 0 = 1 / 2 a_0=1/2 a 0 = 1/2 、A 0 = M ( a 0 ) A_0=M(a_0) A 0 = M ( a 0 ) とする。A k = M ( a k ) A_k=M(a_k) A k = M ( a k ) (0 < a k < 1 0<a_k<1 0 < a k < 1 )に対し、σ \sigma σ を恒等置換にとるとa 11 = a 22 = a k > 0 a_{11}=a_{22}=a_k>0 a 11 = a 22 = a k > 0 であるからσ \sigma σ は補題 2.3 の条件を満たし、min i a i σ ( i ) = a k \min_ia_{i\sigma(i)}=a_k min i a iσ ( i ) = a k である。λ ′ = a k / 2 \lambda'=a_k/2 λ ′ = a k /2 と置くと0 < λ ′ < a k 0<\lambda'<a_k 0 < λ ′ < a k であり、置き換えの結果の対角成分は( a k − a k / 2 ) / ( 1 − a k / 2 ) = a k / ( 2 − a k ) (a_k-a_k/2)/(1-a_k/2)=a_k/(2-a_k) ( a k − a k /2 ) / ( 1 − a k /2 ) = a k / ( 2 − a k ) 、非対角成分は( 1 − a k ) / ( 1 − a k / 2 ) = 1 − a k / ( 2 − a k ) (1-a_k)/(1-a_k/2)=1-a_k/(2-a_k) ( 1 − a k ) / ( 1 − a k /2 ) = 1 − a k / ( 2 − a k ) である。よってa k + 1 = a k / ( 2 − a k ) a_{k+1}=a_k/(2-a_k) a k + 1 = a k / ( 2 − a k ) と置くとA k + 1 = M ( a k + 1 ) A_{k+1}=M(a_{k+1}) A k + 1 = M ( a k + 1 ) である。a k = 1 / ( 2 k + 1 ) a_k=1/(2^{k}+1) a k = 1/ ( 2 k + 1 ) ならば
a k + 1 = 1 / ( 2 k + 1 ) 2 − 1 / ( 2 k + 1 ) = 1 2 ( 2 k + 1 ) − 1 = 1 2 k + 1 + 1 a_{k+1}=\frac{1/(2^{k}+1)}{2-1/(2^{k}+1)}=\frac{1}{2(2^{k}+1)-1}=\frac{1}{2^{k+1}+1} a k + 1 = 2 − 1/ ( 2 k + 1 ) 1/ ( 2 k + 1 ) = 2 ( 2 k + 1 ) − 1 1 = 2 k + 1 + 1 1 であるから、帰納法によりすべてのk k k についてa k = 1 / ( 2 k + 1 ) a_k=1/(2^{k}+1) a k = 1/ ( 2 k + 1 ) であり、0 < a k < 1 0<a_k<1 0 < a k < 1 が保たれる。したがって列A 0 , A 1 , … A_0,A_1,\dots A 0 , A 1 , … はすべての段で定まり、どのA k A_k A k も置換行列ではない。
(4) を示す。第1 1 1 列の成分がすべて1 1 1 で、他の成分がすべて0 0 0 であるn n n 次正方行列をA A A とする。成分は非負であり、各行の和は1 1 1 であり、p ( A ) = n p(A)=n p ( A ) = n である。n ≥ 2 n\ge2 n ≥ 2 であるから第2 2 2 列の和は0 0 0 であり、命題 1.3 により置換行列の列和はすべて1 1 1 であるから、A A A は置換行列ではない。
次に、A A A を置換行列でないn n n 次二重確率行列とする。命題 3.2 によりp ( A ) ≥ n p(A)\ge n p ( A ) ≥ n であり、p ( A ) = n p(A)=n p ( A ) = n ならばA A A は置換行列であるからp ( A ) ≥ n + 1 p(A)\ge n+1 p ( A ) ≥ n + 1 である。p ( A ) = n + 1 p(A)=n+1 p ( A ) = n + 1 と仮定する。各行は少なくとも一つの正の成分をもつので、ある一つの行i 0 i_0 i 0 はちょうど二つの正の成分をもち、他の各行i ≠ i 0 i\ne i_0 i = i 0 はちょうど一つの正の成分a i τ ( i ) a_{i\tau(i)} a i τ ( i ) をもつ。i ≠ i 0 i\ne i_0 i = i 0 のとき、行和が1 1 1 であるからa i τ ( i ) = 1 a_{i\tau(i)}=1 a i τ ( i ) = 1 である。相異なるi , i ′ ≠ i 0 i,i'\ne i_0 i , i ′ = i 0 についてτ ( i ) = τ ( i ′ ) \tau(i)=\tau(i') τ ( i ) = τ ( i ′ ) とすると第τ ( i ) \tau(i) τ ( i ) 列の和が2 2 2 以上になるので、τ ( i ) ≠ τ ( i ′ ) \tau(i)\ne\tau(i') τ ( i ) = τ ( i ′ ) である。よって集合K = { τ ( i ) : i ≠ i 0 } K=\{\tau(i):\ i\ne i_0\} K = { τ ( i ) : i = i 0 } はn − 1 n-1 n − 1 個の列からなる。j = τ ( i ) ∈ K j=\tau(i)\in K j = τ ( i ) ∈ K について第j j j 列の和は1 = a i j 1=a_{ij} 1 = a ij であり成分は非負であるから、a i 0 j = 0 a_{i_0j}=0 a i 0 j = 0 である。したがって行i 0 i_0 i 0 の正の成分はK K K に属さない列にあり、そのような列はn − ( n − 1 ) = 1 n-(n-1)=1 n − ( n − 1 ) = 1 個しかない。これは行i 0 i_0 i 0 が相異なる二つの列に正の成分をもつことに反する。よってp ( A ) ≥ n + 2 p(A)\ge n+2 p ( A ) ≥ n + 2 である。
等号の例として、a 11 = a 12 = a 21 = a 22 = 1 / 2 a_{11}=a_{12}=a_{21}=a_{22}=1/2 a 11 = a 12 = a 21 = a 22 = 1/2 、3 ≤ i ≤ n 3\le i\le n 3 ≤ i ≤ n についてa i i = 1 a_{ii}=1 a ii = 1 とし、他の成分を0 0 0 とするn n n 次正方行列A A A をとる。成分は非負であり、第1 1 1 、2 2 2 行と第1 1 1 、2 2 2 列の和は1 / 2 + 1 / 2 = 1 1/2+1/2=1 1/2 + 1/2 = 1 、第i i i 行と第i i i 列(i ≥ 3 i\ge3 i ≥ 3 )の和は1 1 1 であるから、A A A は二重確率行列である。正の成分は左上の四つとa i i a_{ii} a ii (3 ≤ i ≤ n 3\le i\le n 3 ≤ i ≤ n )であるからp ( A ) = 4 + ( n − 2 ) = n + 2 p(A)=4+(n-2)=n+2 p ( A ) = 4 + ( n − 2 ) = n + 2 である。
(5) を示す。C = ( c i j ) C=(c_{ij}) C = ( c ij ) を成分が実数であるn n n 次正方行列とし、
γ ′ = max σ ∑ i = 1 n c i σ ( i ) \gamma'=\max_{\sigma}\ \sum_{i=1}^{n}c_{i\sigma(i)} γ ′ = σ max i = 1 ∑ n c iσ ( i ) と置く。ここで最大は置換σ \sigma σ の全体をわたる。置換の個数はn ! n! n ! 個で有限であるからγ ′ \gamma' γ ′ は定まり、γ ′ \gamma' γ ′ を与える置換σ ♯ \sigma^{\sharp} σ ♯ が存在する。定義 1.2 によりP σ P_\sigma P σ の第i i i 行で1 1 1 に等しい成分は第σ ( i ) \sigma(i) σ ( i ) 列の成分だけであるから、任意の置換σ \sigma σ について
∑ 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 ) が成り立つ。任意の二重確率行列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 ) を得る。γ ′ \gamma' γ ′ の定義により各k k k について∑ i c i σ k ( i ) ≤ γ ′ \sum_ic_{i\sigma_k(i)}\le\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}\ \le\ \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 = γ ′ が成り立つ。さらに命題 1.3 によりP σ ♯ P_{\sigma^{\sharp}} P σ ♯ は二重確率行列であり、上の等式により∑ i , j c i j ( P σ ♯ ) i j = γ ′ \sum_{i,j}c_{ij}(P_{\sigma^{\sharp}})_{ij}=\gamma' ∑ i , j c ij ( P σ ♯ ) ij = γ ′ である。よって二重確率行列の全体における∑ i , j c i j a i j \sum_{i,j}c_{ij}a_{ij} ∑ i , j c ij a ij の最大値は存在してγ ′ \gamma' γ ′ に等しく、置換行列で達成される。最小化の場合の証明と異なるのは、各項の評価∑ i c i σ k ( i ) ≤ γ ′ \sum_ic_{i\sigma_k(i)}\le\gamma' ∑ i c i σ k ( i ) ≤ γ ′ の不等号の向きだけである。正の係数λ k \lambda_k λ k を掛けて足し合わせる段は、どちらの向きの不等号も保つ。▨