跳到主要内容

早稲田大学 創造理工学研究科 経営システム工学専攻 2017年7月実施 オペレーションズリサーチ 問題7

Author

祭音Myyura

Description

次の線形計画問題 (P)(P) について答えよ。

minimizez=3x12x2subject tox1+x2+x3=9,3x1+x2+x4=18,x1+x5=7,x1,x2,x3,x4,x50.\begin{array}{ll} \text{minimize}&z=-3x_1-2x_2\\ \text{subject to} &x_1+x_2+x_3=9,\\ &3x_1+x_2+x_4=18,\\ &x_1+x_5=7,\\ &x_1,x_2,x_3,x_4,x_5\geq0. \end{array}
  1. 単体法の計算過程と pivot の列・行を示して解け。
  2. 双対問題を示せ。
  3. 主問題と双対問題の最適目的関数値が一致することを示せ。

Kai

[小問 1]

w=z=3x1+2x2w=-z=3x_1+2x_2 の最大化問題とみなす。x3,x4,x5x_3,x_4,x_5 はスラック変数であり、表の最下段は cjzjc_j-z_j である。

初期表

基底cBc_B右辺x1x_1x2x_2x3x_3x4x_4x5x_5
x3x_30911100
x4x_401831010
x5x_50710001
cjzjc_j-z_j32000

x1x_1 列を選び、比 9/1,18/3,7/19/1,18/3,7/1 の最小値から x4x_4 行を pivot 行とする。

第1 pivot 後

基底cBc_B右辺x1x_1x2x_2x3x_3x4x_4x5x_5
x3x_30302/32/311/3-1/30
x1x_13611/31/301/31/30
x5x_50101/3-1/301/3-1/31
cjzjc_j-z_j010-10

次に x2x_2 列を選ぶ。正の係数に対する比は 3/(2/3)=9/23/(2/3)=9/26/(1/3)=186/(1/3)=18 なので、x3x_3 行を pivot 行とする。

第2 pivot 後

基底cBc_B右辺x1x_1x2x_2x3x_3x4x_4x5x_5
x2x_229/29/2013/23/21/2-1/20
x1x_139/29/2101/2-1/21/21/20
x5x_505/25/2001/21/21/2-1/21
cjzjc_j-z_j003/2-3/21/2-1/20

すべての cjzj0c_j-z_j\leq0 なので最適である。したがって

(x1,x2,x3,x4,x5)=(92,92,0,0,52)\boxed{ (x_1,x_2,x_3,x_4,x_5) =\left(\frac92,\frac92,0,0,\frac52\right) }

であり、

w=452,z=452.\boxed{w^*=\frac{45}{2}},\qquad \boxed{z^*=-\frac{45}{2}}.

[小問 2]

等式制約に対する双対変数を y1,y2,y3y_1,y_2,y_3 とすると、本来これらは自由変数である。各スラック変数 x3,x4,x50x_3,x_4,x_5\geq0 に対応する双対制約から、さらに y1,y2,y30y_1,y_2,y_3\leq0 が課される。双対問題は

maximize9y1+18y2+7y3subject toy1+3y2+y33,y1+y22,y1,y2,y30.\begin{array}{ll} \text{maximize} &9y_1+18y_2+7y_3\\ \text{subject to} &y_1+3y_2+y_3\leq-3,\\ &y_1+y_2\leq-2,\\ &y_1,y_2,y_3\leq0. \end{array}

同値に ui=yi0u_i=-y_i\geq0 とおけば、ww の最大化問題に対する双対

minimize9u1+18u2+7u3subject tou1+3u2+u33,u1+u22,u1,u2,u30\begin{array}{ll} \text{minimize} &9u_1+18u_2+7u_3\\ \text{subject to} &u_1+3u_2+u_3\geq3,\\ &u_1+u_2\geq2,\\ &u_1,u_2,u_3\geq0 \end{array}

を得る。

[小問 3]

双対の実行可能解

(u1,u2,u3)=(32,12,0)(u_1,u_2,u_3)=\left(\frac32,\frac12,0\right)

では、2本の制約がともに等号で成立し、目的値は

932+1812+70=452=w.9\cdot\frac32+18\cdot\frac12+7\cdot0 =\frac{45}{2}=w^*.

元の双対変数では y=(3/2,1/2,0)y=(-3/2,-1/2,0) であり、

9y1+18y2+7y3=452=z.9y_1+18y_2+7y_3=-\frac{45}{2}=z^*.

主・双対の実行可能解の目的値が一致したため、弱双対性から両者は最適である。