跳到主要内容

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

Author​

hari64boli64, 祭音Myyura

Description​

nn 次実正方行列 A=(aij)∈Rn×nA = (a_{ij}) \in \mathbb{R}^{n \times n}, B=(bij)∈Rn×nB = (b_{ij}) \in \mathbb{R}^{n \times n} は

aij≥0(i,j=1,2,…,n),∑i=1naij=1(j=1,2,…,n),B=αA+1−αn11⊤a_{ij} \geq 0 \quad (i, j = 1, 2, \ldots, n), \quad \sum_{i=1}^{n} a_{ij} = 1 \quad (j = 1, 2, \ldots, n), \quad B = \alpha A + \frac{1 - \alpha}{n} \boldsymbol{1} \boldsymbol{1}^\top

を満たすものとする。 ただし、⊤\top は転置を表し、α\alpha は 0<α<10 < \alpha < 1 を満たす実数、1\boldsymbol{1} はすべての成分が 11 であるベクトル 1=(1,…,1)⊤∈Rn\boldsymbol{1} = (1, \ldots, 1)^\top \in \mathbb{R}^n とする。 また、v=(v1,…,vn)⊤∈Rnv = (v_1, \ldots, v_n)^\top \in \mathbb{R}^n について ∥v∥1:=∑i=1n∣vi∣\lVert v \rVert_1 := \sum_{i=1}^{n} |v_i| とする。以下の設問に答えよ。

(1) 行列 AA の固有値の絶対値で最大となるものを求めよ。

(2) 各成分が非負のベクトル x∈Rnx \in \mathbb{R}^n で、

Bx=x,1⊤x=1Bx = x, \quad \boldsymbol{1}^\top x = 1

となるものが存在することを示せ。なお、次の事実を用いてよい:

「Rn\mathbb{R}^n の非空なコンパクト凸集合からそれ自身への連続写像には不動点が存在する。」

(3) 1⊤q=0\boldsymbol{1}^\top q = 0 を満たすベクトル q=(q1,…,qn)⊤∈Rnq = (q_1, \ldots, q_n)^\top \in \mathbb{R}^n に対して、

∣∑j=1nbijqj∣≤∑j=1nbij∣qj∣−1−αn∥q∥1(i=1,2,…,n)\left| \sum_{j=1}^{n} b_{ij} q_j \right| \leq \sum_{j=1}^{n} b_{ij} |q_j| - \frac{1 - \alpha}{n} \|q\|_1 \qquad (i = 1, 2, \ldots, n)

となることを示せ。

(4) (2) の条件を満たすベクトル x∈Rnx \in \mathbb{R}^n は、正の整数 NN に対して、

∥BN1n−x∥1≤αN∥1n−x∥1\left \lVert B^N \frac{\boldsymbol{1}}{n} - x \right \rVert_1 \leq \alpha^N \left \lVert \frac{\boldsymbol{1}}{n} - x \right \rVert_1

を満たすことを示せ。

题目描述​

设实矩阵 A=(aij),B=(bij)∈Rn×nA=(a_{ij}),B=(b_{ij})\in\mathbb R^{n\times n} 满足

aij≥0,∑i=1naij=1(j=1,…,n),a_{ij}\ge0,\qquad \sum_{i=1}^na_{ij}=1\quad(j=1,\ldots,n),
B=αA+1−αn11⊤,0<α<1,B=\alpha A+\frac{1-\alpha}{n}\mathbf1\mathbf1^\top, \qquad0<\alpha<1,

其中 1=(1,…,1)⊤\mathbf1=(1,\ldots,1)^\top。向量的 1-范数定义为

∥v∥1=∑i=1n∣vi∣.\|v\|_1=\sum_{i=1}^n|v_i|.
  1. 求 AA 的特征值模的最大值。

  2. 证明存在分量非负的概率向量 xx 满足

    Bx=x,1⊤x=1.Bx=x,\qquad\mathbf1^\top x=1.

    可使用“非空紧凸集到自身的连续映射存在不动点”这一事实。

  3. 对任意满足 1⊤q=0\mathbf1^\top q=0 的向量 q=(q1,…,qn)⊤q=(q_1,\ldots,q_n)^\top,证明对每个 ii,

    ∣∑j=1nbijqj∣≤∑j=1nbij∣qj∣−1−αn∥q∥1.\left|\sum_{j=1}^nb_{ij}q_j\right| \le\sum_{j=1}^nb_{ij}|q_j| -\frac{1-\alpha}{n}\|q\|_1.
  4. 对第 2 问的任一不动点 xx 和任意正整数 NN,证明

    ∥BN1n−x∥1≤αN∥1n−x∥1.\left\|B^N\frac{\mathbf1}{n}-x\right\|_1 \le\alpha^N\left\|\frac{\mathbf1}{n}-x\right\|_1.

Kai​

(1)​

一般に、転置行列の固有値は、固有方程式が同じになることから、元の行列と同じ固有値を取る。

ATA^T は、1\boldsymbol{1} を固有ベクトルとして、11 を固有値に持つ。

(AT1)i=∑jnAijT1j=∑jnaji=1⇒AT1=11\begin{aligned} & (A^T\boldsymbol{1})_i=\sum_{j}^n A^T_{ij}\boldsymbol{1}_j=\sum_{j}^na_{ji}=1 \\ \Rightarrow & A^T\boldsymbol{1}=1\boldsymbol{1} \end{aligned}

そして、11 より絶対値が大きな固有値を持つことはない。

固有ベクトル x\boldsymbol{x} の、絶対値に関する最大値を xix_i で取るとすると、

(ATx)i=λxi⇔∑jnajixj=λxi\begin{aligned} & (A^T\boldsymbol{x})_i =\lambda x_i \\ \Leftrightarrow & \sum_{j}^n a_{ji}x_j =\lambda x_i \\ \end{aligned}

一方、

∣∑j=1najixj∣≤∑j=1naji∣xj∣≤∣xi∣∑j=1naji=∣xi∣<∣λ∣ ∣xi∣=∣λxi∣\begin{aligned} \left|\sum_{j=1}^n a_{ji}x_j\right| \leq \sum_{j=1}^n a_{ji}|x_j| \leq |x_i|\sum_{j=1}^n a_{ji} =|x_i|<|\lambda|\,|x_i|=|\lambda x_i| \end{aligned}

つまり、∣λ∣>1|\lambda|>1 なら、この不等式に反し矛盾。

よって、AA の固有値の絶対値で最大となるものは 11 である。

(2)​

関数 f:S→Sf:S \to S を f(x)=Bxf(\boldsymbol{x})=B\boldsymbol{x} で定義する。

ただし、S={x∈R≥0n∣1Tx=1}S=\{\boldsymbol{x}\in \mathbb{R}_{\geq 0}^n \mid \boldsymbol{1}^T\boldsymbol{x}=1 \} である。

まず Bx∈SB\boldsymbol{x}\in S を示す。

Bx=αAx+1−αn11Tx\begin{aligned} B\boldsymbol{x}=\alpha A\boldsymbol{x}+\frac{1-\alpha}{n}\boldsymbol{1}\boldsymbol{1}^T\boldsymbol{x} \end{aligned}

という表式と、aij,αa_{ij},\alpha の値域より、Bx∈R≥0nB\boldsymbol{x}\in \mathbb{R}_{\geq 0}^n は明らか。

1TBx=1\boldsymbol{1}^TB\boldsymbol{x}=1 を示す。

1TBx=1T(αAx+1−αn11Tx)=α(1TA)x+1−αn1T11Tx=α1Tx+(1−α)1Tx=α+(1−α)=1\begin{aligned} \boldsymbol{1}^TB\boldsymbol{x} & =\boldsymbol{1}^T\left(\alpha A\boldsymbol{x}+\frac{1-\alpha}{n}\boldsymbol{1}\boldsymbol{1}^T\boldsymbol{x}\right) \\ & =\alpha (\boldsymbol{1}^TA)\boldsymbol{x}+\frac{1-\alpha}{n}\boldsymbol{1}^T\boldsymbol{1}\boldsymbol{1}^T\boldsymbol{x} \\ & =\alpha \boldsymbol{1}^T\boldsymbol{x}+(1-\alpha)\boldsymbol{1}^T\boldsymbol{x} \\ & =\alpha+(1-\alpha) \\ & =1 \end{aligned}

よって、f(S)⊂Sf(S)\subset S である。

線形写像 ff は連続である。SS は非空な凸集合であり、また有界閉集合なのでコンパクトである。

以上より、ヒントで与えられているブラウワーの不動点定理より、

x∈Rn,Bx=x,1Tx=1\begin{aligned} \boldsymbol{x}\in \mathbb{R}^n, B\boldsymbol{x}=\boldsymbol{x}, \boldsymbol{1}^T\boldsymbol{x}=1 \end{aligned}

の条件を満たすような x\boldsymbol{x} が存在する。

(3)​

∣∑j=1nbijqj∣=∣(Bq)i∣=∣((αA+1−αn11T)q)i∣=∣(αAq)i∣(∵1Tq=0)=∣∑j=1nαaijqj∣≤∑j=1n(αaij)∣qj∣=∑j=1n(B−1−αn11T)ij∣qj∣=∑j=1n(bij−1−αn)∣qj∣=∑j=1nbij∣qj∣−1−αn∥q∥1\begin{aligned} & \left|\sum_{j=1}^nb_{ij}q_j\right| \\ = & |(B\boldsymbol{q})_i| \\ = & \left|\left((\alpha A+\frac{1-\alpha}{n}\boldsymbol{1}\boldsymbol{1}^T)\boldsymbol{q}\right)_i\right| \\ = & \left|(\alpha A\boldsymbol{q})_i\right| \quad (\because \boldsymbol{1}^T\boldsymbol{q}=0) \\ = & \left|\sum_{j=1}^n\alpha a_{ij}q_j\right| \\ \leq & \sum_{j=1}^n(\alpha a_{ij})|q_j| \\ = & \sum_{j=1}^n \left(B-\frac{1-\alpha}{n}\boldsymbol{1}\boldsymbol{1}^T\right)_{ij}|q_j| \\ = & \sum_{j=1}^n \left(b_{ij}-\frac{1-\alpha}{n}\right)|q_j| \\ = & \sum_{j=1}^n b_{ij}|q_j|-\frac{1-\alpha}{n}\|\boldsymbol{q}\|_1 \\ \end{aligned}

(4)​

まず、q=1n−x\boldsymbol{q}=\frac{\boldsymbol{1}}{n}-\boldsymbol{x} と定義される q\boldsymbol{q} を用いると、題意は、

∥BN1n−x∥1≤αN∥1n−x∥1⇔∥BN1n−BNx∥1≤αN∥q∥1⇔∥BNq∥1≤αN∥q∥1\begin{aligned} & \|B^N\frac{\boldsymbol{1}}{n}-\boldsymbol{x}\|_1 \leq \alpha^N \|\frac{\boldsymbol{1}}{n}-\boldsymbol{x}\|_1 \\ \Leftrightarrow & \|B^N\frac{\boldsymbol{1}}{n}-B^N\boldsymbol{x}\|_1 \leq \alpha^N \|\boldsymbol{q}\|_1 \\ \Leftrightarrow & \|B^N\boldsymbol{q}\|_1 \leq \alpha^N \|\boldsymbol{q}\|_1 \\ \end{aligned}

になる。

特に、1Tq=1T1n−1Tx=nn−1=0\boldsymbol{1}^T\boldsymbol{q}=\boldsymbol{1}^T\frac{\boldsymbol{1}}{n}-\boldsymbol{1}^T\boldsymbol{x}=\frac{n}{n}-1=0 であるが、この性質は、以下に示すように、BB を乗じても変わらない。

1T(Bq)=1T(αA+1−αn11T)q=α1TAq+1−αn1T11Tq=α1Tq+(1−α)1Tq=1Tq=0(∵1Tq=0)\begin{aligned} \boldsymbol{1}^T(B\boldsymbol{q}) & =\boldsymbol{1}^T \left(\alpha A+\frac{1-\alpha}{n}\boldsymbol{1}\boldsymbol{1}^T\right)\boldsymbol{q} \\ & =\alpha \boldsymbol{1}^TA\boldsymbol{q}+\frac{1-\alpha}{n}\boldsymbol{1}^T\boldsymbol{1}\boldsymbol{1}^T\boldsymbol{q} \\ & =\alpha \boldsymbol{1}^T\boldsymbol{q}+(1-\alpha)\boldsymbol{1}^T\boldsymbol{q} \\ & =\boldsymbol{1}^T\boldsymbol{q} \\ & =0 \quad (\because \boldsymbol{1}^T\boldsymbol{q}=0) \end{aligned}

よって、1Tq=0\boldsymbol{1}^T\boldsymbol{q}=0 を満たす q\boldsymbol{q} に関して、∥Bq∥1≤α∥q∥1\|B\boldsymbol{q}\|_1\leq \alpha\|\boldsymbol{q}\|_1 であることを示せばよい。

これは、(3) で示したことから、

∥Bq∥1≤α∑i=1n∑j=1naij∣qj∣=α∑j=1n(∑i=1naij)∣qj∣=α∥q∥1.\begin{aligned} \|B\boldsymbol{q}\|_1 &\leq \alpha\sum_{i=1}^n\sum_{j=1}^n a_{ij}|q_j| \\ &=\alpha\sum_{j=1}^n\left(\sum_{i=1}^na_{ij}\right)|q_j| \\ &=\alpha\|\boldsymbol{q}\|_1. \end{aligned}

となり、直ちに従う。

Knowledge​

コンパクト集合​

Rn\mathbb{R}^n において、有界閉集合がコンパクト集合であって、閉集合がコンパクト集合である訳では無いことに注意。

例として、R∖(−1,1)\mathbb{R} \setminus (-1,1) は閉集合であるが、当然コンパクト集合ではない。

ブラウワーの不動点定理​

(2) のヒントは、ブラウワーの不動点定理と呼ばれている。