跳到主要内容

東京工業大学 情報理工学院 数理・計算科学系 2022年8月実施 运筹学

Author

思齐塾, 祭音Myyura

Description

問 4

nn を 2 以上の整数とし, R++\mathbb{R}_{++} を正の実数の集合とする. nn 次元実ベクトル c,aR++n\boldsymbol{c}, \boldsymbol{a} \in \mathbb{R}_{++}^n と正の実数 bR++b \in \mathbb{R}_{++} に対して,線形計画問題 (P) を以下で定める:

(P)maximizecTxsubject toaTxb,0xi1(i=1,,n).\begin{aligned} \text{(P)} \quad \text{maximize} \quad & \boldsymbol{c}^T \boldsymbol{x} \\ \text{subject to} \quad & \boldsymbol{a}^T \boldsymbol{x} \le b, \\ & 0 \le x_i \le 1 \quad (i = 1, \dots, n). \end{aligned}

ここで, yT\boldsymbol{y}^T はベクトル y\boldsymbol{y} の転置を表す.ベクトル c\boldsymbol{c}a\boldsymbol{a} の要素は

c1a1c2a2cnan,aib(i=1,,n)\begin{gathered} \frac{c_1}{a_1} \ge \frac{c_2}{a_2} \ge \dots \ge \frac{c_n}{a_n}, \\ a_i \le b \quad (i = 1, \dots, n) \end{gathered}

を満たすとし, i=1kaib<i=1k+1ai\sum_{i=1}^k a_i \le b < \sum_{i=1}^{k+1} a_i となる添字 k{1,,n1}k \in \{1, \dots, n-1\} が存在するとする.このとき,以下の問いに答えよ.

(1) 線形計画問題 (P) の双対問題 (D) を書け.

(2) 線形計画問題 (P) と (D) それぞれの最適解を求め,これらが最適解であることを双対定理に基づいて確認せよ.

(3) 線形計画問題 (P) の最適値を α\alpha とし,整数計画問題

maximizecTxsubject toaTxb,x{0,1}n\begin{aligned} \text{maximize} \quad & \boldsymbol{c}^T \boldsymbol{x} \\ \text{subject to} \quad & \boldsymbol{a}^T \boldsymbol{x} \le b, \\ & \boldsymbol{x} \in \{0, 1\}^n \end{aligned}

の最適値を αˉ\bar{\alpha} とする.このとき αˉ12α\bar{\alpha} \ge \frac{1}{2} \alpha を示せ.

题目描述

n2n\geq2R++\mathbb R_{++} 表示正实数集合。给定

c,aR++n,bR++,\boldsymbol c,\boldsymbol a\in\mathbb R_{++}^n, \qquad b\in\mathbb R_{++},

考虑线性规划

(P)最大化cx,约束条件axb,0xi1(i=1,,n).\begin{aligned} \text{(P)}\qquad \text{最大化}\quad&\boldsymbol c^\top\boldsymbol x,\\ \text{约束条件}\quad&\boldsymbol a^\top\boldsymbol x\leq b,\\ &0\leq x_i\leq1\quad(i=1,\ldots,n). \end{aligned}

其中上标 \top 表示转置。假设

c1a1c2a2cnan,aib(i=1,,n),\frac{c_1}{a_1}\geq\frac{c_2}{a_2}\geq\cdots\geq\frac{c_n}{a_n}, \qquad a_i\leq b\quad(i=1,\ldots,n),

并且存在 k{1,,n1}k\in\{1,\ldots,n-1\} 满足

i=1kaib<i=1k+1ai.\sum_{i=1}^k a_i\leq b <\sum_{i=1}^{k+1}a_i.
  1. 写出 (P) 的对偶问题 (D)。
  2. 分别求 (P) 与 (D) 的最优解,并通过对偶定理验证它们确实最优。
  3. 记 (P) 的最优值为 α\alpha,记相应的 0011 整数规划
最大化cx,约束条件axb,x{0,1}n\begin{aligned} \text{最大化}\quad&\boldsymbol c^\top\boldsymbol x,\\ \text{约束条件}\quad&\boldsymbol a^\top\boldsymbol x\leq b,\\ &\boldsymbol x\in\{0,1\}^n \end{aligned}

的最优值为 αˉ\bar\alpha。证明

αˉ12α.\bar\alpha\geq\frac12\alpha.

Kai

解答

Sk=i=1kaiS_k=\sum_{i=1}^k a_ir=bSkr=b-S_k とおく。仮定から 0r<ak+10\le r<a_{k+1} である。

(1) 双対問題

容量制約に双対変数 y0y\ge0 、上限制約 xi1x_i\le1zi0z_i\ge0 を対応させると

(D)minimizeby+i=1nzisubject toaiy+zici(i=1,,n),y0,zi0(i=1,,n)\begin{aligned} \text{(D)}\quad \text{minimize}\quad &by+\sum_{i=1}^n z_i\\ \text{subject to}\quad &a_i y+z_i\ge c_i\quad(i=1,\ldots,n),\\ &y\ge0,\quad z_i\ge0\quad(i=1,\ldots,n) \end{aligned}

となる。

(2) 最適解

(P) に対して

xi={1,1ik,r/ak+1,i=k+1,0,k+2inx_i^*=\begin{cases} 1,&1\le i\le k,\\ r/a_{k+1},&i=k+1,\\ 0,&k+2\le i\le n \end{cases}

とする。 0r/ak+1<10\le r/a_{k+1}<1 かつ aTx=Sk+r=b\boldsymbol a^T\boldsymbol x^*=S_k+r=b なので実行可能である。

(D) に対して

y=ck+1ak+1,zi={ciaiy,1ik,0,k+1iny^*=\frac{c_{k+1}}{a_{k+1}},\qquad z_i^*=\begin{cases}c_i-a_i y^*,&1\le i\le k,\\0,&k+1\le i\le n\end{cases}

とする。 ci/aic_i/a_i の降順性から、 iki\le k では zi0z_i^*\ge0ik+1i\ge k+1 では aiycia_i y^*\ge c_i であり、双対実行可能である。

両目的値は

cTx=i=1kci+rak+1ck+1\boldsymbol c^T\boldsymbol x^*=\sum_{i=1}^k c_i+\frac{r}{a_{k+1}}c_{k+1}

および

by+izi=by+i=1k(ciaiy)=i=1kci+rck+1ak+1by^*+\sum_i z_i^*=by^*+\sum_{i=1}^k(c_i-a_i y^*) =\sum_{i=1}^k c_i+r\frac{c_{k+1}}{a_{k+1}}

で一致する。弱双対性により両者は最適解であり、

α=i=1kci+bSkak+1ck+1.\boxed{\alpha=\sum_{i=1}^k c_i+\frac{b-S_k}{a_{k+1}}c_{k+1}}.

(3)

最初の kk 個だけを選ぶ整数解は実行可能で価値 Ck=i=1kciC_k=\sum_{i=1}^k c_i をもつ。また ak+1ba_{k+1}\le b より、第 k+1k+1 項だけを選ぶ整数解も実行可能で価値 ck+1c_{k+1} をもつ。よって

αˉmax{Ck,ck+1}.\bar\alpha\ge\max\{C_k,c_{k+1}\}.

一方、 0r/ak+1<10\le r/a_{k+1}<1 なので

α=Ck+rak+1ck+1Ck+ck+12max{Ck,ck+1}2αˉ.\alpha=C_k+\frac{r}{a_{k+1}}c_{k+1} \le C_k+c_{k+1} \le2\max\{C_k,c_{k+1}\} \le2\bar\alpha.

従って

αˉ12α.\boxed{\bar\alpha\ge\frac12\alpha}.