跳到主要内容

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

Author

hari64boli64

Description

頂点集合 V={v1,v2,,vn}V = \{v_1, v_2, \ldots, v_n\} と枝集合 EE からなる連結無向グラフ G=(V,E)G = (V, E) を考える。 頂点 viv_i に接続する枝の本数を did_i と書く。 ただし、GG は自己閉路や多重枝は無いものとする。n×nn \times n 行列 A=(aij),L=(lij)A = (a_{ij}), L = (l_{ij})

aij={1({vi,vj}E のとき),0(それ以外のとき),lij={1({vi,vj}E のとき),di(i=j のとき),0(それ以外のとき)a_{ij} = \begin{cases} 1 & \text{($ \{v_i, v_j\} \in E$ のとき)}, \\ 0 & \text{(それ以外のとき)} \end{cases}, \quad l_{ij} = \begin{cases} -1 & \text{($ \{v_i, v_j\} \in E$ のとき)}, \\ d_i & (i = j\text{ のとき}), \\ 0 & \text{(それ以外のとき)} \end{cases}

と定義する。以下の設問に答えよ。

(1) 行列 AA のべき乗 AkA^k(i,j)(i, j) 成分は何を表わすか。

(2) 2 頂点間の距離をその 2 頂点を結ぶ経路の最小枝数で定める。任意の 2 頂点間の距離は、AA の相異なる固有値の個数より小さいことを示せ。

(3) 行列 LL の非零固有値に対応する固有ベクトル u=(ui)u = (u_i) について、i=1nui=0\sum_{i=1}^{n} u_i = 0 となることを示せ。

(4) 行列 LL の固有値はすべて非負実数であることを示せ。

(5) 関数 V:RnRV : \mathbb{R}^n \rightarrow \mathbb{R}

V(x)=121i<jnaij(xixj)2(xRn)V(x) = \frac{1}{2} \sum_{1 \leq i < j \leq n} a_{ij}(x_i - x_j)^2 \quad (x \in \mathbb{R}^n)

と定義し、x(t)=(x1(t),x2(t),,xn(t))x(t) = (x_1(t), x_2(t), \ldots, x_n(t)) に関する微分方程式系

dxi(t)dt=V(x)xix=x(t)(i=1,2,,n)\frac{\text{d}x_i(t)}{\text{d}t} = - \frac{\partial V(x)}{\partial x_i} \bigg|_{x = x(t)} \quad (i = 1, 2, \ldots, n)

を考える。初期値 x(0)=(c1,c2,,cn)x(0) = (c_1, c_2, \ldots, c_n) に対する解 x(t)x(t) の極限 x=limtx(t)\overline{x} = \lim_{t \to \infty} x(t) を求め、収束の速さについて論じよ。

Kai

(1)

頂点 ii から頂点 jj へ長さが kk の経路の数

(2)

ケーリーハミルトンの定理より、背理法

(3)

j=1n0=0j=1n0uj=0j=1n(i=1nLij)uj=0i=1nj=1nLijuj=0i=1n(Lu)i=0i=1nλui=0i=1nui=0\begin{aligned} & \sum_{j=1}^{n}{0}=0 \\ \Leftrightarrow & \sum_{j=1}^{n}{0u_j}=0 \\ \Leftrightarrow & \sum_{j=1}^{n}\left(\sum_{i=1}^{n}{L_{ij}} \right)u_j=0 \\ \Leftrightarrow & \sum_{i=1}^{n}\sum_{j=1}^{n}{L_{ij}u_j}=0 \\ \Leftrightarrow & \sum_{i=1}^{n}{(Lu)_i}=0 \\ \Leftrightarrow & \sum_{i=1}^{n}{\lambda u_i}=0 \\ \Leftrightarrow & \sum_{i=1}^{n}{u_i}=0 \\ \end{aligned}

(4)

xTLx=i,jxiLijxj=i,jxi(DijAij)xj=ixi2Diii<jxixjAiji>jxixjAij=i(xi2jAij)i<jxixjAiji<jxjxiAji=i<j(xi22xixj+xj2)Aij=i<jaij(xixj)20\begin{aligned} \boldsymbol{x}^T L \boldsymbol{x} & =\sum_{i,j}x_i L_{ij} x_j \\ & =\sum_{i,j}x_i (D_{ij}-A_{ij}) x_j \\ & =\sum_{i}x_i^2 D_{ii}-\sum_{i<j}x_i x_j A_{ij}- \sum_{i>j}x_i x_j A_{ij} \\ & =\sum_{i}\left(x_i^2 \sum_{j}{A_{ij}}\right)-\sum_{i<j}x_i x_j A_{ij}- \sum_{i<j}x_j x_i A_{ji} \\ & =\sum_{i<j}(x_i^2 -2x_i x_j +x_j^2) A_{ij} \\ & =\sum_{i<j}a_{ij}(x_i -x_j)^2 \\ & \geq 0 \end{aligned}

(5)

dxdt=Lx\frac{\text{d}\boldsymbol{x}}{\text{d}t}=-L\boldsymbol{x} より、x(t)=eLtx(0)\boldsymbol{x}(t)=e^{-Lt}\boldsymbol{x}(0) となる。

LL の固有値が全て非負実数の為、x=limtx(t)=0\boldsymbol{\overline{x}}=\lim_{t \to \infty}\boldsymbol{x}(t) = \boldsymbol{0} となる。

また、収束の速さは LL の固有値に依存する。