改訂単体法のソースを表示

'''改訂単体法'''（かいていたんたいほう、'''改訂シンプレックス法'''、{{lang-en-short|Revised simplex method}}）とは、数理最適化において、[[ジョージ・ダンツィーグ]]によって考案された[[線形計画問題]]に対する[[単体法]]に工夫を施した解法である。

改訂単体法は（シンプレックス表を用いた）単体法と同様の手順で実行可能基底解と呼ばれる線形計画問題の最適解となり得る解を求めていくが、実行可能基底解の生成方法が異なっている。基底変数の組み合わせを表現する辞書（基底形式表現）を表したシンプレックス表を用いる代わりに、[[基底 (線型代数学)|基底]]変数に対応する制約式の係数を要素とする[[行列]]を生成していく。線形計画問題を行列形式として表現することで行列の疎な構造を利用して効率良く解くことができる{{sfn|Morgan|1997|loc=§2}}{{Sfn|今野浩|1987|pp=60-63}}。

== 定式化 ==
これ以降では、線形計画問題が下記のような標準形として与えられているものとして説明を行う。
:<math>
\begin{array}{rl}
\text{Minimize} & \boldsymbol{c}^{\mathrm{T}} \boldsymbol{x} \\
\text{s.t.} & \boldsymbol{Ax} = \boldsymbol{b}, \boldsymbol{x} \ge \boldsymbol{0}
\end{array}
</math>
ただし {{math|'''''A''''' ∈ ℝ<sup>''m''×''n''</sup>}} は制約行列であり、{{math|'''''x''''' ≥ '''0'''}} の下で {{math|'''''Ax''''' {{=}} '''''b'''''}} を満たすような解が少なくとも一つ存在し、{{math|'''''A'''''}} が行フルランクであると仮定する。もし、{{math|'''''A'''''}} がフルランクでない場合は、制約条件に冗長な制約式が存在するか、実行不可能な制約が含まれている。{{math|'''''A'''''}} がこのような場合であったとしても、前処理の段階で対処することができる。

== アルゴリズム ==
=== 最適性の条件 ===
線形計画問題では、[[KKT条件]]を満たす解が最適解となるための[[必要十分条件]]となる。線形計画問題に対するKKT条件は以下の通りである: 

:<math>
\begin{align}
\boldsymbol{Ax} & = \boldsymbol{b}, \\
\boldsymbol{A}^{\mathrm{T}} \boldsymbol{\lambda} + \boldsymbol{s} & = \boldsymbol{c}, \\
\boldsymbol{x} & \ge \boldsymbol{0}, \\
\boldsymbol{s} & \ge \boldsymbol{0}, \\
\boldsymbol{s}^{\mathrm{T}} \boldsymbol{x} & = 0
\end{align}
</math>

ただし {{math|'''''λ'''''}} および {{math|'''''s'''''}} は、それぞれ {{math|'''''Ax''''' {{=}} '''''b'''''}}、{{math|'''''x''''' ≥ '''0'''}} に対応する[[ラグランジュの未定乗数法|ラグランジュ乗数]]（KKT乗数）である{{sfn|Nocedal|Wright|2006|p=358|loc=Eq.&nbsp;13.4}}。また最下段の条件 {{math|''s''<sup>T</sup> ''x'' {{=}} ''0''}} は、{{math|''s<sub>i</sub>x<sub>i</sub>'' {{=}} 0}}, {{math|''i'' {{=}} 1, ..., ''n''}} と等価な条件である。この条件は一般的に相補性条件（相補スラック条件）と呼ばれる。

{{仮リンク|線形計画の基本定理|en|Fundamental theorem of linear programming}}から、実行可能多面体の頂点 {{math|'''''x'''''}} は、行列 {{math|'''''A'''''}} の列から選ばれた基底 {{math|'''''B'''''}} によって導かれる{{efn|またこの定理は、実行可能な多面体が構成されるときは、少なくとも一つの頂点が生成され、最適解となりうる頂点も存在することを主張している{{sfn|Nocedal|Wright|2006|p=363|loc=Theorem&nbsp;13.2}}。}}。{{math|'''''A'''''}}が行フルランクであるとき、{{math|'''''B'''''}} は正則である。そして、{{math|'''''A'''''}} は一般性を失うことなく {{math|'''''A''''' {{=}} &#91;'''''B'''''&ensp;'''''N'''''&#93;}} と分割できたと仮定すると、{{math|'''''x'''''}} は

:<math>
\boldsymbol{x} =
\begin{bmatrix}
\boldsymbol{x_B} \\
\boldsymbol{x_N}
\end{bmatrix} =
\begin{bmatrix}
\boldsymbol{B}^{-1} \boldsymbol{b} \\
\boldsymbol{0}
\end{bmatrix}
</math>

と表される。ただし、{{math|'''''x<sub>B</sub>''''' ≥ '''0'''}} である。また {{math|'''''c'''''}} および {{math|'''''s'''''}} はそれぞれ

:<math>
\begin{align}
\boldsymbol{c} & =
\begin{bmatrix}
\boldsymbol{c_B} \\
\boldsymbol{c_N}
\end{bmatrix}, \\
\boldsymbol{s} & =
\begin{bmatrix}
\boldsymbol{s_B} \\
\boldsymbol{s_N}
\end{bmatrix}.
\end{align}
</math>

のように分割される。{{math|'''''x<sub>B</sub>''''' ≥ '''0'''}} であることから相補性条件より {{math|'''''s<sub>B</sub>''''' {{=}} '''0'''}} を満たさなければならない。このことから、最適性の条件は

:<math>
\begin{align}
\boldsymbol{B}^{\mathrm{T}} \boldsymbol{\lambda} & = \boldsymbol{c_B}, \\
\boldsymbol{N}^{\mathrm{T}} \boldsymbol{\lambda} + \boldsymbol{s_N} & = \boldsymbol{c_N},
\end{align}
</math>

と書き直される。これらの条件の下で

:<math>
\begin{align}
\boldsymbol{\lambda} & = (\boldsymbol{B}^{\mathrm{T}})^{-1} \boldsymbol{c_B}, \\
\boldsymbol{s_N} & = \boldsymbol{c_N} - \boldsymbol{N}^{\mathrm{T}} \boldsymbol{\lambda}.
\end{align}
</math>

ここで {{math|'''''s<sub>N</sub>''''' ≥ '''0'''}} を満たしていれば、{{math|'''''x'''''}} は最適性の条件であるKKT条件を満たし、最適解であるといえる。

=== ピボット操作 ===
もし現在の解がKKT条件を満たしていない場合には、非基底変数の中で {{math|s<sub>N</sub> < 0}} に対応する {{math|x <sub>N</sub>}} を制約条件（主実行可能性）が満たされるように基底変数と入れ替えるピボット操作を行う。[[退化 (数学)|退化]]していない場合は、ピボット操作によって目的関数値 {{math|'''''c'''''<sup>T</sup>'''''x'''''}} は厳密に減少する。したがって、線形計画問題が非退化仮定などの仮定の下では、制約条件が表す多面体の頂点が有限個であることから、改訂単体法はピボット操作の反復により最適解の頂点に到達する{{sfn|Nocedal|Wright|2006|p=370|loc=Theorem&nbsp;13.4}}。

{{math|''m'' < ''q'' ≤ ''n''}} を満たし、{{math|''s<sub>q</sub>'' < 0}} となる添字 q を選び、これを基底に入る変数の添字と呼ぶ。対応する行列 {{math|'''''A'''''}} の列 {{math|'''''A'''<sub>q</sub>''}} が基底に移動され、{{math|''x<sub>q</sub>''}} は 0 から増大される。このときの目的関数値の単位増加量は以下の通りである。

:<math>\frac{\partial (\boldsymbol{c}^{\mathrm{T}} \boldsymbol{x})}{\partial x_q} = s_q,</math>

すなわち、{{math|''x<sub>q</sub>''}} を 1 増加させると、目的関数値 {{math|'''''c'''''<sup>T</sup>'''''x'''''}} は {{math|−''s<sub>q</sub>''}} 増加する{{sfn|Nocedal|Wright|2006|p=369|loc=Eq.&nbsp;13.24}}。また、

:<math>\boldsymbol{B x_B} + \boldsymbol{A}_q x_q = \boldsymbol{b},</math>

であることから、{{math|'''''x<sub>B</sub>'''''}} は、{{math|Δ'''''x<sub>B</sub>''''' {{=}} '''''B'''''<sup>−1</sup>'''''A'''<sub>q</sub>x<sub>q</sub>''}} によって対応するように減少させる必要があり、これは変数の非負条件から {{math|'''''x<sub>B</sub>''''' − Δ'''''x<sub>B</sub>''''' ≥ '''0'''}} を満たさなければならない。ここで、{{math|'''''d''''' {{=}} '''''B'''''<sup>−1</sup>'''''A'''<sub>q</sub>''}} とおく。

もし {{math|'''''d''''' ≤ '''0'''}} ならば {{math|''x<sub>q</sub>''}} をどれだけ増加させても {{math|'''''x<sub>B</sub>''''' − Δ'''''x<sub>B</sub>''''' ≥ '''0'''}} は非負のままである。この場合、{{math|'''''c'''''<sup>T</sup>'''''x'''''}} は任意に減少可能であることから、与えられた線形計画問題の[[実行可能領域]]は非有界{{Efn|{{Lang-en-short|unbounded}}}}である（最適解をもたない）。

そうでない場合は、基底から出る変数の添字として
{{math|''p'' {{=}} argmin<sub>1≤''i''≤''m''</sub> {{(}}''x<sub>i</sub>''/''d<sub>i</sub>'' {{!}} ''d<sub>i</sub>'' > 0{{)}}}}

を選択する。この選択により、{{math|''x<sub>q</sub>''}} を 0 から増加させながら、制約条件（主実行可能性）を満たしつつ {{math|'''''A'''<sub>p</sub>''}} を 0 まで減少させることができる。結果としてピボット操作後の基底は {{math|'''''A'''<sub>q</sub>''}} から {{math|'''''A'''<sub>p</sub>''}} に置き換わる。

== 具体例 ==
<!-- {{seealso|Simplex method#Example}} -->
下記の線形計画問題

:<math>
\begin{align}
\boldsymbol{c} & =
\begin{bmatrix}
-2 & -3 & -4 & 0 & 0
\end{bmatrix}^{\mathrm{T}}, \\
\boldsymbol{A} & =
\begin{bmatrix}
3 & 2 & 1 & 1 & 0 \\
2 & 5 & 3 & 0 & 1
\end{bmatrix}, \\
\boldsymbol{b} & =
\begin{bmatrix}
10 \\
15
\end{bmatrix}.
\end{align}
</math>

が与えられたときに改訂単体法により最適解を求める。基底行列・非基底行列を

:<math>
\begin{align}
\boldsymbol{B} & =
\begin{bmatrix}
\boldsymbol{A}_4 & \boldsymbol{A}_5
\end{bmatrix}, \\
\boldsymbol{N} & =
\begin{bmatrix}
\boldsymbol{A}_1 & \boldsymbol{A}_2 & \boldsymbol{A}_3
\end{bmatrix}
\end{align}
</math>

としたとき、初期実行可能基底解は {{math|'''''x''''' {{=}} &#91;0&ensp;0&ensp;0&ensp;10&ensp;15&#93;<sup>T</sup>}} であり、{{math|'''''λ'''''}}、{{math|'''''S<sub>N</sub>'''''}} は

:<math>
\begin{align}
\boldsymbol{\lambda} & =
\begin{bmatrix}
0 & 0
\end{bmatrix}^{\mathrm{T}}, \\
\boldsymbol{s_N} & =
\begin{bmatrix}
-2 & -3 & -4
\end{bmatrix}^{\mathrm{T}}.
\end{align}
</math>

と求まる。

{{math|''q'' {{=}} 3}} を基底に入る変数の添字として選択する。すると、{{math|'''''d''''' {{=}} &#91;1&ensp;3&#93;<sup>T</sup>}} となり、これは {{math|''x''<sub>3</sub>}} を 1 単位増加させると {{math|''x''<sub>4</sub>}} と {{math|''x''<sub>5</sub>}} がそれぞれ 1 と 3 減少することを意味する。したがって、{{math|''x''<sub>3</sub>}} を 5 まで増加させると、その時点で {{math|''x''<sub>5</sub>}} が 0 まで減少し、{{math|''p'' {{=}} 5}} が基底から出る変数の添字となる。

ピボット操作後では

:<math>
\begin{align}
\boldsymbol{B} & =
\begin{bmatrix}
\boldsymbol{A}_3 & \boldsymbol{A}_4
\end{bmatrix}, \\
\boldsymbol{N} & =
\begin{bmatrix}
\boldsymbol{A}_1 & \boldsymbol{A}_2 & \boldsymbol{A}_5
\end{bmatrix}.
\end{align}
</math>

となり、

:<math>
\begin{align}
\boldsymbol{x} & =
\begin{bmatrix}
0 & 0 & 5 & 5 & 0
\end{bmatrix}^{\mathrm{T}}, \\
\boldsymbol{\lambda} & =
\begin{bmatrix}
0 & -4/3
\end{bmatrix}^{\mathrm{T}}, \\
\boldsymbol{s_N} & =
\begin{bmatrix}
2/3 & 11/3 & 4/3
\end{bmatrix}^{\mathrm{T}}.
\end{align}
</math>

が求まる。このとき {{math|'''''s<sub>N</sub>'''''}} が非負であることから、{{math|'''''x'''''}} は最適性の条件（KKT条件）を満たし、最適解 {{math|'''''x'''''}} が求まった。

== 実用上で起こり得る問題点 ==
=== 退化 ===
{{seealso|:en:Simplex method#Degeneracy: stalling and cycling|label1=単体法#退化:停滞と巡回}}
改訂単体法は単体法と同様に最適解を求めるため、入力の {{math|''''' A ''''',''''' b ''''', ''''' c '''''}} によってはピボット操作を行ったときに目的関数値 {{math|'''''c'''''<sup>T</sup>'''''x'''''}} を改善することができない退化やピボット操作に求まる辞書が退化により何度も同一のものが表れる巡回現象が生じ得る。巡回が発生する場合では、{{math|''''' b '''''}} の各成分に互いに異なる微小な値 {{math|'''''ε'''''}} を加える辞書式摂動法やBlandの最小添字規則に従ったピボット操作を行うことで巡回を防止し、有限回のピボット操作で最適解を求めることができる{{sfn|Nocedal|Wright|2006|p=381|loc=§13.5}}。

=== 基底表現 ===
改訂単体法において、基底{{math|'''''B'''''}}が2種類の[[線型方程式系|方程式系]]を満たすものとする。

:<math>
\begin{align}
\boldsymbol{B z} & = \boldsymbol{y}, \\
\boldsymbol{B}^{\mathrm{T}} \boldsymbol{z} & = \boldsymbol{y}.
\end{align}
</math>

通常は {{math|'''''B'''''}} を反復ごとに因子分解するのではなく、[[LU分解]]された行列が直接更新される。その戦略として、Forrest−Tomlin 法や Bartels−Golub 法などが挙げられる。しかし、行列の更新にかかるデータの量や数値誤差がピボット操作を繰り返すごとに蓄積されるため、定期的な因子分解が必要となる{{sfn|Morgan|1997|loc=§2}}{{sfn|Nocedal|Wright|2006|p=372|loc=§13.4}}。

==脚注==
{{脚注ヘルプ}}
=== 注釈 ===
{{Notelist}}

=== 出典 ===
{{Reflist|30em}}

=== 参考文献 ===
{{refbegin}}
* {{Cite book |和書 |author=今野浩 |author-link=今野浩 |year=1987 |title=線形計画法 |publisher=日科技連 |isbn=4-8171-5014-9 |ref=harv }}
* {{cite thesis |last=Morgan |first=S. S. |degree=MSc |title=A Comparison of Simplex Method Algorithms |publisher=[[フロリダ大学]] |year=1997 |url=http://www.cise.ufl.edu/research/sparse/Morgan/index.htm |archive-url=https://web.archive.org/web/20110807134509/http://www.cise.ufl.edu/research/sparse/Morgan/index.htm |archive-date=7 August 2011 }}
* {{cite book |last1=Nocedal |first1=J.|last2=Wright |first2=S. J. |editor1-last=Mikosch |editor1-first=T. V. |editor2-last=Resnick |editor2-first=S. I. |editor3-last=Robinson |editor3-first=S. M. |title=Numerical Optimization |edition=2nd |year=2006 |series=Springer Series in Operations Research and Financial Engineering |publisher=[[Springer Science+Business Media|Springer]] |location=New York, NY, USA |isbn=978-0-387-30303-1 |url=https://www.springer.com/mathematics/book/978-0-387-30303-1 |ref=harv }}
{{refend}}

{{アルゴリズム}}
{{最適化アルゴリズム}}
{{Normdaten}}
{{Comp-stub}}
{{Math-stub}}

{{DEFAULTSORT:かいていたんたいほう}}
[[Category:線型計画法]]
[[Category:数学に関する記事]]