跳到主要内容

東京大学 情報理工学研究科 数理情報学 2018年8月実施 第3問

Author

hari64boli64

Description

Rn\mathbb{R}^nnn 次元ユークリッド空間とする。 x,yRnx, y \in \mathbb{R}^n の内積を x,y\langle x, y \rangle で表し、xx のノルムを x:=x,x\lVert x \rVert := \sqrt{\langle x, x \rangle} とする。 以下の設問に答えよ。

(1) 非空な閉集合 KRnK \subseteq \mathbb{R}^n と点 xRnx \in \mathbb{R}^n に対して、以下を示せ:

  • (1-1) xy=infzKxz\lVert x - y \rVert = \inf_{z \in K} \lVert x - z \rVert を満たす yKy \in K が存在する。
  • (1-2) KK が凸なら、そのような点 yy は一意に定まる。
  • (1-3) KK が凸で、xxKK に含まれないなら、ある cRnc \in \mathbb{R}^ndRd \in \mathbb{R} が存在して、
(c,x)>d,(c,z)d(zK)\begin{aligned} &(c, x) > d, \\ &(c, z) \leq d \quad (z \in K) \end{aligned}

   となる。

(2) Rn\mathbb{R}^n の部分集合 A\mathcal{A} に対して C(A)RnC(\mathcal{A}) \subseteq \mathbb{R}^nA\mathcal{A} の元たちの非負結合全体とする。 すなわち、A\mathcal{A} の元 a1,a2,,am (m1)a_1, a_2, \ldots, a_m \ (m \geq 1) と非負実数 λ1,λ2,,λm\lambda_1, \lambda_2, \ldots, \lambda_m によって、 x=λ1a1+λ2a2++λmamx = \lambda_1 a_1 + \lambda_2 a_2 + \cdots + \lambda_m a_m とかける点 xx からなる集合が C(A)C(\mathcal{A}) である。

  • (2-1) ARn\mathcal{A} \subseteq \mathbb{R}^n に対して、以下が成立することを示せ:
C(A)=BC(B).C(\mathcal{A}) = \bigcup_B C(\mathcal{B}).

   ここで、和は、A\mathcal{A} のすべての一次独立な部分集合 B\mathcal{B} にわたってとる。

  • (2-2) ARn\mathcal{A} \subseteq \mathbb{R}^n が有限集合なら C(A)C(\mathcal{A}) は閉集合となることを示せ。

(3) n×mn \times m 実行列 ARn×mA \in \mathbb{R}^{n \times m}nn 次元ベクトル xRnx \in \mathbb{R}^n に対して、以下の 2 つの性質 (P)(P), (Q)(Q) を考える:

  • (P)(P) Aλ=x,λ0A \lambda = x, \lambda \geq 0 を満たす λRm\lambda \in \mathbb{R}^m が存在する。
  • (Q)(Q) cA0,cx>0c^\top A \leq 0, c^\top x > 0 を満たす cRnc \in \mathbb{R}^n が存在する。

ここで、\top は転置を表し、ベクトル uu に対して記法 u0 (u0)u \geq 0 \ (u \leq 0)uu の各成分が非負(非正)であることを意味する。以下を示せ:

  • (3-1) (P)(P)(Q)(Q) は同時に成立することはない。
  • (3-2) (P)(P)(Q)(Q) のどちらかは成立する。

Kai

(1)

これは分離定理を示す問題である。

(1-1)

Rn\mathbb{R}^n における有界な閉集合 KK はコンパクト集合であること、及び、コンパクト集合上で定義される連続関数 f:KR,f(z)=xzf:K \to \mathbb{R},f(z)=\lVert x-z \rVert は、最大値と最小値を持つことから、xy=infzKf(z)\lVert x-y \rVert=\inf_{z\in K}f(z) を達成する yKy \in K が取れる。

なお、KK は有界とは限らないが、本問に関しては、xx からの距離で適当に区切っても題意に影響がない為、有界な閉集合に限定出来る。

(1-2)

凸性より明らか。省略する。

(1-3)

d=0d=0 として良い。

イメージとしては、xxyy の間の垂直二等分線 (分離超平面) で、xxyy が別々の領域に分かれる為、内積の正負が反転するという感じである。

しかし、厳密にこれを記述することはかなり難しいと思われる。

ここでは省略する。

(「分離超平面定理」を参照)

(2)

(2-1)

C(A)BC(B)C(\mathcal{A}) \supseteq \bigcup_{\mathcal{B}}C(\mathcal{B}) は明らか。

C(A)BC(B)C(\mathcal{A}) \subseteq \bigcup_{\mathcal{B}}C(\mathcal{B}) を示す。

C(A)C(\mathcal{A}) の元 xx を取る。

x=λ1a1++λmamx=\lambda_1 a_1+ \cdots + \lambda_m a_m に関して、am=μ1a1++μm1am1a_m=\mu_1 a_1 + \cdots + \mu_{m-1} a_{m-1} と書けたと仮定する。

全ての 1i<m1 \leq i < m に対して、λi+λmμi0\lambda_i+\lambda_m \mu_i \geq 0 が成立すれば、xC({ai1i<m})x \in C(\{a_i | 1 \leq i < m \}) となる。

そうでない場合、λi+kμi\lambda_i+k \mu_i が最小の kk00 になるようなインデックスを ii に取ると、λm=λm+k\lambda_m=\lambda_m'+k として、

x=λ1a1++λiai++λmam=λ1a1++λiai++(λm+k)am=λ1a1++λiai++λmam+k(μ1a1++μm1am1)=(kμ1+λ1)a1++(kμi+λi)ai++(kμm1+λm1)am1+λmam=(kμ1+λ1)a1++0++(kμm1+λm1)am1+λmamC({aj1jm,ji})\begin{aligned} x & =\lambda_1 a_1 + \cdots + \lambda_i a_i + \cdots + \lambda_m a_m \\ & =\lambda_1 a_1 + \cdots + \lambda_i a_i + \cdots +(\lambda_m'+k) a_m \\ & =\lambda_1 a_1 + \cdots + \lambda_i a_i + \cdots + \lambda_m' a_m + k(\mu_1a_1+\cdots+\mu_{m-1}a_{m-1}) \\ & =(k\mu_1+\lambda_1) a_1 + \cdots + (k\mu_i+\lambda_i)a_i + \cdots + (k \mu_{m-1}+\lambda_{m-1}) a_{m-1}+\lambda_m' a_m \\ & =(k\mu_1+\lambda_1) a_1 + \cdots + 0 + \cdots + (k \mu_{m-1}+\lambda_{m-1}) a_{m-1}+\lambda_m' a_m \\ & \in C(\{a_j | 1 \leq j \leq m, j\neq i \}) \end{aligned}

となり、要素数を減らすことが出来る。

以下、帰納的に一次独立になるまで、この議論を繰り返せばよい。

よって、示された。

(2-2)

C(B)C(\mathcal{B}) が閉集合であることが言えれば、有限個の閉集合の和集合は閉集合であることから、C(A)C(\mathcal{A}) も閉集合であると言える。

C(B)C(\mathcal{B}) が閉集合であることを示す。

これは、有限生成錐の閉性を示せば良い。

これもかなり記述は難しい気がする。

ここでは省略する。

(「有限生成錐が閉集合になることについて」を参照)

(3)

以上を基に、ファルカスの補題を示す。

(3-1)

(P)λRm s.t. Aλ=x,λ0\begin{aligned} (P) & \Rightarrow \exists \lambda \in \mathbb{R}^m \text{ s.t. } A\lambda =x,\lambda \geq 0 \end{aligned}

この時、cTA0cTAλ=cTx0c^TA \leq 0 \Rightarrow c^TA\lambda =c^Tx \leq 0 となり、(P)(Q)(P)\Rightarrow \overline{(Q)} が示された。

(3-2)

AA の各列ベクトルが生成する有限生成錐 C(A={ai}1im)C(\mathcal{A}=\{a_i\}_{1\leq i \leq m}) は非空な凸閉集合である。

(P)\overline{(P)} ならば、xC(A)x \notin C(\mathcal{A}) であり、分離定理より、c,x>0\langle c,x\rangle >0 かつ、c,ai0\langle c,a_i\rangle \leq 0 となる cc が存在する。これは (Q)(Q) に他ならない。

よって、(P)(Q)\overline{(P)}\Rightarrow (Q) が示された。

Knowledge

ファルカスの補題は、弱双対定理からも示せる。

minxcTx s.t. Ax=b,x0\begin{aligned} \min_{x} c^Tx \text{ s.t. } Ax = b, x \geq \boldsymbol{0} \end{aligned}
maxybTy s.t. ATyc\begin{aligned} \max_{y} b^Ty \text{ s.t. } A^Ty \leq c \end{aligned}

記号を入れ替えて、

minλ0Tλ s.t. Aλ=x,λ0minλ0 s.t. Aλ=x,λ0\begin{aligned} & \min_{\lambda} \boldsymbol{0}^T\lambda \text{ s.t. } A\lambda = x, \lambda \geq \boldsymbol{0} \\ \Rightarrow & \min_{\lambda} 0 \text{ s.t. } A\lambda = x, \lambda \geq \boldsymbol{0} \end{aligned}
maxcxTc s.t. ATc0maxccTx s.t. cTA0\begin{aligned} & \max_{c} x^Tc \text{ s.t. } A^Tc \leq \boldsymbol{0} \\ \Rightarrow & \max_{c} c^Tx \text{ s.t. } c^TA \leq \boldsymbol{0} \end{aligned}

となる。

(Q)\overline{(Q)} ならば、弱双対定理より、(P)(P) となる。