跳到主要内容

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

Author

祭音Myyura

Description

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

minimizez=3x12x24x3subject tox1+x2+2x34,2x1+2x35,2x1+x2+3x37,x1,x2,x30.\begin{array}{ll} \text{minimize}&z=-3x_1-2x_2-4x_3\\ \text{subject to} &x_1+x_2+2x_3\leq4,\\ &2x_1+2x_3\leq5,\\ &2x_1+x_2+3x_3\leq7,\\ &x_1,x_2,x_3\geq0. \end{array}
  1. 単体法の計算過程と pivot の列・行を示して解け。
  2. 双対問題を示せ。
  3. 主問題と双対問題の最適目的関数値が一致することを示せ。

Kai

[小問 1]

w=z=3x1+2x2+4x3w=-z=3x_1+2x_2+4x_3 の最大化問題へ変換し、スラック変数 s1,s2,s3s_1,s_2,s_3 を加える。表の最下段は cjzjc_j-z_j である。

初期表

基底cBc_B右辺x1x_1x2x_2x3x_3s1s_1s2s_2s3s_3
s1s_104112100
s2s_205202010
s3s_307213001
cjzjc_j-z_j324000

x3x_3 列を選び、比 4/2,5/2,7/34/2,5/2,7/3 の最小値から s1s_1 行を pivot 行とする。

第1 pivot 後

基底cBc_B右辺x1x_1x2x_2x3x_3s1s_1s2s_2s3s_3
x3x_3421/21/21/21/211/21/200
s2s_2011-10-110
s3s_3011/21/21/2-1/203/2-3/201
cjzjc_j-z_j100-200

x1x_1 列を選び、比 2/(1/2),1/1,1/(1/2)2/(1/2),1/1,1/(1/2) の最小値から s2s_2 行を pivot 行とする。

第2 pivot 後

基底cBc_B右辺x1x_1x2x_2x3x_3s1s_1s2s_2s3s_3
x3x_343/23/201111/2-1/20
x1x_1311-10-110
s3s_301/21/2000-11/2-1/21
cjzjc_j-z_j010-1-10

x2x_2 列を選ぶ。正の係数を持つのは x3x_3 行だけなので、そこを pivot 行とする。

第3 pivot 後

基底cBc_B右辺x1x_1x2x_2x3x_3s1s_1s2s_2s3s_3
x2x_223/23/201111/2-1/20
x1x_135/25/210101/21/20
s3s_301/21/2000-11/2-1/21
cjzjc_j-z_j00-1-21/2-1/20

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

(x1,x2,x3)=(52,32,0),w=212,z=212.\boxed{ (x_1,x_2,x_3)=\left(\frac52,\frac32,0\right),\qquad w^*=\frac{21}{2},\qquad z^*=-\frac{21}{2} }.

[小問 2]

元の最小化問題 (P)(P) の双対変数を y1,y2,y30y_1,y_2,y_3\leq0 とすると、双対問題 (D)(D)

maximize4y1+5y2+7y3subject toy1+2y2+2y33,y1+y32,2y1+2y2+3y34,y1,y2,y30.\begin{array}{ll} \text{maximize} &4y_1+5y_2+7y_3\\ \text{subject to} &y_1+2y_2+2y_3\leq-3,\\ &y_1+y_3\leq-2,\\ &2y_1+2y_2+3y_3\leq-4,\\ &y_1,y_2,y_3\leq0. \end{array}

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

minimize4u1+5u2+7u3subject tou1+2u2+2u33,u1+u32,2u1+2u2+3u34,u1,u2,u30\begin{array}{ll} \text{minimize} &4u_1+5u_2+7u_3\\ \text{subject to} &u_1+2u_2+2u_3\geq3,\\ &u_1+u_3\geq2,\\ &2u_1+2u_2+3u_3\geq4,\\ &u_1,u_2,u_3\geq0 \end{array}

となる。

[小問 3]

双対の実行可能解

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

を取ると、左辺は

ATu=(3,2,5)T(3,2,4)TA^{\mathsf T}u=(3,2,5)^{\mathsf T}\geq(3,2,4)^{\mathsf T}

であり、目的値は

42+512+70=212=w.4\cdot2+5\cdot\frac12+7\cdot0 =\frac{21}{2}=w^*.

すなわち元の双対変数では y=(2,1/2,0)y=(-2,-1/2,0)

4y1+5y2+7y3=212=z.4y_1+5y_2+7y_3=-\frac{21}{2}=z^*.

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