跳到主要内容

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

Author

hari64boli64

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}

aij0(i,j=1,2,,n),i=1naij=1(j=1,2,,n),B=αA+1αn11a_{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 は転置を表し、α\alpha0<α<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 について v1:=i=1nvi\lVert v \rVert_1 := \sum_{i=1}^{n} |v_i| とする。以下の設問に答えよ。

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

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

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

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

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

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

j=1nbijqjj=1nbijqj1αnq1(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) の条件を満たすベクトル xRnx \in \mathbb{R}^n は、正の整数 NN に対して、

BN1nx1αN1nx1\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

を満たすことを示せ。

Kai

(1)

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

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

(AT1)i=jnAijT1j=jnaji=1AT1=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=λxijnajixj=λxi\begin{aligned} & (A^T\boldsymbol{x})_i =\lambda x_i \\ \Leftrightarrow & \sum_{j}^n a_{ji}x_j =\lambda x_i \\ \end{aligned}

一方、

jnajixjjnajixjjnajixixijnaji=xi<λxi=λxi\begin{aligned} \left|\sum_{j}^n a_{ji}x_j\right| \leq \sum_{j}^n |a_{ji}\|x_j| \leq \sum_{j}^n |a_{ji}\|x_i| \leq |x_i|\sum_{j}^n a_{ji} = |x_i| < |\lambda\|x_i|=|\lambda x_i| \end{aligned}

つまり、λ\lambda11 より大きな値を取ると、この不等式に反し矛盾。

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

(2)

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

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

ここで、S=TS=T を証明する。

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 の値域より、TR0nT \subset \mathbb{R}_{\geq 0}^n は明らか。

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

1TBx=1T(αAx+1αn11Tx)=α(1TA)x+1αn1T11Tx=α1Tx+1αn1Tx=α+1αn=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}+\frac{1-\alpha}{n}\boldsymbol{1}^T\boldsymbol{x} \\ & =\alpha+\frac{1-\alpha}{n} \\ & =1 \end{aligned}

よって、T=ST=S である。

SSRn\mathbb{R}^n の非空なコンパクト集合であることは、SSRn\mathbb{R}^n の有界閉集合であることから、これはコンパクト。

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

xRn,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αaijqjj=1n(αaij)qj=j=1n(B1αn11T)ijqj=j=1n(bij1αn)qj=j=1nbijqj1αnq1\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=1nx\boldsymbol{q}=\frac{\boldsymbol{1}}{n}-\boldsymbol{x} と定義される q\boldsymbol{q} を用いると、題意は、

BN1nx1αN1nx1BN1nBNx1αNq1BNq1αNq1\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=1T1n1Tx=nn1=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αn1Tq=(α+1αn)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}+\frac{1-\alpha}{n}\boldsymbol{1}^T\boldsymbol{q} \\ & =\left(\alpha+\frac{1-\alpha}{n}\right)\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} に関して、Bq1αq1\|B\boldsymbol{q}\|_1\leq \alpha\|\boldsymbol{q}\|_1 であることを示せばよい。

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

(Bq)ij=1n(αaij)qjαj=1nqj=αq1\begin{aligned} |(B\boldsymbol{q})_i| & \leq \sum_{j=1}^n(\alpha a_{ij})|q_j| \\ & \leq \alpha \sum_{j=1}^n |q_j| \\ & =\alpha \|\boldsymbol{q}\|_1 \end{aligned}

となり、直ちに従う。

Knowledge

コンパクト集合

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

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

ブラウワーの不動点定理

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