院試hub

東京大学 院試 過去問 解答例

東大 数理科学研究科 数理科学専攻 専門科目A 2014年度 院試 解答例・解説

東京大学 数理科学研究科 数理科学専攻 専門科目A 2014年度の院試 過去問について、設問ごとの解法方針・部分点の置き所を解説。全7問収録の解答・解説PDFと併用できます。問題本文は含みません。

最終更新:

1 — 線形写像と正定値性

差を取って階段化する

この行列は同じ値がブロック状に並ぶので,掃き出しを機械的に行うよりも, 行同士の差を見る方が速い。第1行と第2行,第3行と第2行を比べるだけで x2+x3+x4=0x_2+x_3+x_4=0x3+x4=0x_3+x_4=0 が出る。

正定値性の見せ方

単に固有値を計算する必要はない。直交補空間を (x,y,z,z)(x,y,z,z) とパラメータ表示し,平方完成で (x+y+2z)2+(a1)(y+2z)2+4(ba)z2 (x+y+2z)^2+(a-1)(y+2z)^2+4(b-a)z^2 と書けば,仮定 1<a<b1<a<b がどこで効くかが明確になる。

2014年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1) 問題の行列を AA とおく。v=(x1,x2,x3,x4)Tv=(x_1,x_2,x_3,x_4)^T に対して Av=0Av=0 を考える。第1行と第2行の差から (1a)(x2+x3+x4)=0 (1-a)(x_2+x_3+x_4)=0 であるから,x2+x3+x4=0x_2+x_3+x_4=0 である。これを第1行に戻すと x1=0x_1=0 である。さらに第3行と第2行の差から (ba)(x3+x4)=0 (b-a)(x_3+x_4)=0 となり,1<a<b1<a<b より x3+x4=0x_3+x_4=0 である。したがって x2=0x_2=0 も従う。

よって kerf={(0,0,t,t)TtR}=R(0,0,1,1)T. \ker f=\{(0,0,t,-t)^T\mid t\in\mathbb R\} =\mathbb R(0,0,1,-1)^T . この直交補空間は V={(x,y,z,z)Tx,y,zR} V=\{(x,y,z,z)^T\mid x,y,z\in\mathbb R\} である。

(2) いま v=(x,y,z,z)TVv=(x,y,z,z)^T\in V とおくと,直接計算により vTAv=x2+2xy+4xz+ay2+4ayz+4bz2=(x+y+2z)2+(a1)(y+2z)2+4(ba)z2. \begin{aligned} v^TAv &=x^2+2xy+4xz+ay^2+4ayz+4bz^2\\ &=(x+y+2z)^2+(a-1)(y+2z)^2+4(b-a)z^2 . \end{aligned} a1>0, ba>0a-1>0,\ b-a>0 であるから右辺は非負である。さらに右辺が 00 ならば z=0,y+2z=0,x+y+2z=0 z=0,\qquad y+2z=0,\qquad x+y+2z=0 となり,x=y=z=0x=y=z=0 である。したがって VV 上の二次形式は正定値である。

最終答

(1)kerf=R(0,0,1,1)T,V={(x,y,z,z)Tx,y,zR}. \text{(1)}\quad \ker f=\mathbb R(0,0,1,-1)^T,\qquad V=\{(x,y,z,z)^T\mid x,y,z\in\mathbb R\}. また (2)vTAv=(x+y+2z)2+(a1)(y+2z)2+4(ba)z2 \text{(2)}\quad v^TAv=(x+y+2z)^2+(a-1)(y+2z)^2+4(b-a)z^2 となるので,VV 上で正定値である。

2 — 乗法的Cauchy方程式

割って加法形にする

式の右辺に xf(y)x f(y)yf(x)y f(x) が出るときは,f(x)/xf(x)/x を作るのが定石である。 これにより,正の実数の乗法群から実数の加法群への連続準同型になる。

連続性の役割

有理数乗の公式だけなら代数的に出せるが,微分可能性には連続性が本質的である。 連続な加法的関数が一次関数であることを使うと,解の形が cxlogxcx\log x に一意に決まる。

2014年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1) x>0x>0 に対して g(x)=f(x)x g(x)=\frac{f(x)}{x} とおく。与えられた関係式を xyxy で割ると g(xy)=g(x)+g(y)(x,y>0) g(xy)=g(x)+g(y)\qquad (x,y>0) を得る。ff は連続であるから gg も連続である。

αQ\alpha\in\mathbb Q とする。まず正の整数 mm について g(xm)=mg(x) g(x^m)=m g(x) が帰納法で従う。また x=(x1/m)mx=(x^{1/m})^m より g(x1/m)=1mg(x) g(x^{1/m})=\frac1m g(x) である。従って α=m/n\alpha=m/n と書けば g(xα)=αg(x) g(x^\alpha)=\alpha g(x) である。これを f(x)=xg(x)f(x)=xg(x) に戻すと f(xα)=xαg(xα)=αxαg(x)=αxα1f(x) f(x^\alpha)=x^\alpha g(x^\alpha) =\alpha x^\alpha g(x) =\alpha x^{\alpha-1}f(x) となる。

(2) 次に微分可能性を示す。x=etx=e^t と置き, h(t)=g(et) h(t)=g(e^t) と定めると,hh は連続で h(s+t)=h(s)+h(t) h(s+t)=h(s)+h(t) を満たす。連続な加法的関数は一次関数であるから,ある定数 cc が存在して h(t)=ct h(t)=ct である。よって g(x)=clogx,f(x)=cxlogx. g(x)=c\log x,\qquad f(x)=cx\log x . 右辺は (0,)(0,\infty) 上で微分可能である。

最終答

(1)f(xα)=αxα1f(x)(αQ, x>0) \text{(1)}\quad \boxed{f(x^\alpha)=\alpha x^{\alpha-1}f(x)\qquad(\alpha\in\mathbb Q,\ x>0)} (2)ある定数 cR によりf(x)=cxlogx(x>0) \text{(2)}\quad \text{ある定数 }c\in\mathbb R\text{ により}\quad \boxed{ f(x)=cx\log x\qquad (x>0) } である。特に ff は微分可能である。

3 — 距離関数の微分

下限を取る順序

距離関数の Lipschitz 性は,三角不等式を書いたあとに AA 上の下限を取れば終わる。 閉集合であることはこの部分では不要で,最近点の存在を使う後半で効いている。

勾配の大きさ

1-Lipschitz 性から勾配の大きさは高々 11 である。一方,最近点へ向かう直線上では 距離が速度 11 で減るため,ある方向微分の絶対値が 11 になる。両側から挟むのが要点である。

2014年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1) 任意の x,yRnx,y\in\mathbb R^naAa\in A に対して三角不等式より xaxy+ya \|x-a\|\le \|x-y\|+\|y-a\| である。aAa\in A について下限を取ると d(x)xy+d(y) d(x)\le \|x-y\|+d(y) となる。x,yx,y を入れ替えれば d(y)xy+d(x) d(y)\le \|x-y\|+d(x) も成り立つ。したがって d(x)d(y)xy |d(x)-d(y)|\le \|x-y\| であり,dd は Lipschitz 連続である。

(2) 次に x0Ax_0\notin Add が全微分可能であるとする。r=d(x0)>0r=d(x_0)>0 とおく。 AA は閉集合なので,x0x_0 からの距離を実現する点 a0Aa_0\in A が存在する。 u=a0x0a0x0 u=\frac{a_0-x_0}{\|a_0-x_0\|} とおくと,0tr0\le t\le r について d(x0+tu)=rt d(x_0+tu)=r-t である。実際,右辺以下であることは a0Aa_0\in A から従い,右辺以上であることは 三角不等式から従う。

したがって方向微分は Dud(x0)=limt+0d(x0+tu)d(x0)t=1. D_ud(x_0)=\lim_{t\to+0}\frac{d(x_0+tu)-d(x_0)}{t}=-1 . 一方,全微分可能性より Dud(x0)=d(x0)u D_ud(x_0)=\nabla d(x_0)\cdot u である。よって Cauchy--Schwarz の不等式から 1=d(x0)ud(x0) 1=|\nabla d(x_0)\cdot u|\le \|\nabla d(x_0)\| を得る。すでに dd は 1-Lipschitz なので,任意の単位ベクトル ww に対して d(x0)w1 |\nabla d(x_0)\cdot w|\le 1 であり,d(x0)1\|\nabla d(x_0)\|\le 1 である。従って d(x0)=1. \|\nabla d(x_0)\|=1 .

最終答

(1)d(x)d(y)xy \text{(1)}\quad |d(x)-d(y)|\le \|x-y\| である。さらに x0Ax_0\notin Add が全微分可能なら (2)d(x0)=1 \text{(2)}\quad \|\nabla d(x_0)\|=1 である。

4 — Poisson核型積分

符号付きの ana^n

aa は負でもよい。展開は a<1|a|<1 で絶対収束するため,係数はそのまま ana^n であり,nn の偶奇による符号も自然に含まれる。

a>1|a|>1 の処理

a>1|a|>1 の場合は逆数 r=1/ar=1/a に直す。ここで ana^{-n} が現れるため, 単に an|a|^{-n} と書くと負の aa の場合に誤る。

2014年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

まず a<1|a|<1 とする。このとき 112acosθ+a2=1(1aeiθ)(1aeiθ) \frac{1}{1-2a\cos\theta+a^2} =\frac{1}{(1-ae^{i\theta})(1-ae^{-i\theta})} であり,絶対収束する幾何級数を掛け合わせて 112acosθ+a2=11a2(1+2k=1akcoskθ) \frac{1}{1-2a\cos\theta+a^2} =\frac{1}{1-a^2}\left(1+2\sum_{k=1}^{\infty}a^k\cos k\theta\right) を得る。三角関数の直交性より,正の整数 nn について In=11a22an02πcos2(nθ)dθ=2πan1a2. I_n=\frac{1}{1-a^2}\cdot 2a^n\int_0^{2\pi}\cos^2(n\theta)\,d\theta =\frac{2\pi a^n}{1-a^2}.

次に a>1|a|>1 とする。r=1/ar=1/a とおけば r<1|r|<1 で, 12acosθ+a2=a2(12rcosθ+r2) 1-2a\cos\theta+a^2 =a^2(1-2r\cos\theta+r^2) である。したがって先ほどの結果を rr に適用して In=1a22πrn1r2=2πana21. I_n =\frac{1}{a^2}\frac{2\pi r^n}{1-r^2} =\frac{2\pi a^{-n}}{a^2-1}.

最終答

正の整数 nn に対して In={2πan1a2,a<1,2πana21,a>1. \boxed{ I_n= \begin{cases} \dfrac{2\pi a^n}{1-a^2}, & |a|<1,\\[1.1em] \dfrac{2\pi a^{-n}}{a^2-1}, & |a|>1. \end{cases}}

5 — 回転面の面積条件

全区間で面積が一致する意味

S(p,q)=2πT(p,q)S(p,q)=2\pi T(p,q) が任意の区間で成り立つので,積分値だけでなく被積分関数が一致する。 ここを明記しないと,微分方程式への移行が不十分になる。

単調性の使い方

ff が単調増加であるため,平方根を取ったあと f0f'\ge 0 の枝を選ぶ。 この符号選択がないと,同じ式から減少解も形式的には出てしまう。

2014年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

曲線を zz 軸のまわりに回転したとき,区間 [p,q][p,q] に対応する回転面の面積は S(p,q)=2πpqf(x)1+f(x)2dx S(p,q)=2\pi\int_p^q f(x)\sqrt{1+f'(x)^2}\,dx である。一方,対応する平面領域の面積は T(p,q)=pq(f(x)+1)dx T(p,q)=\int_p^q (f(x)+1)\,dx である。条件はすべての 0p<q10\le p<q\le 1 で成り立つので,被積分関数が等しく f(x)1+f(x)2=f(x)+1 f(x)\sqrt{1+f'(x)^2}=f(x)+1 でなければならない。

ff は正値かつ単調増加であるから f(x)=2f(x)+1f(x) f'(x)=\frac{\sqrt{2f(x)+1}}{f(x)} である。y=f(x)y=f(x) と見て逆関数を用いると dxdy=y2y+1. \frac{dx}{dy}=\frac{y}{\sqrt{2y+1}} . 初期条件 f(0)=1f(0)=1 から y=1y=1 のとき x=0x=0 である。したがって x=1yu2u+1du. x=\int_1^y \frac{u}{\sqrt{2u+1}}\,du . s=2u+1s=\sqrt{2u+1} と置くと u2u+1du=(2u+1)3/262u+12 \int \frac{u}{\sqrt{2u+1}}\,du =\frac{(2u+1)^{3/2}}{6}-\frac{\sqrt{2u+1}}{2} であり,下端 u=1u=1 ではこの値は 00 である。よって x=(2y+1)3/232y+16=(y1)2y+13. x=\frac{(2y+1)^{3/2}-3\sqrt{2y+1}}{6} =\frac{(y-1)\sqrt{2y+1}}{3}.

最終答

曲線は xxyy の関数として x=(y1)2y+13 \boxed{x=\frac{(y-1)\sqrt{2y+1}}{3}} と表される。

6 — 可換部分空間の次元

条件 (b) の読み替え

条件 (a) は VV 内の可換性,条件 (b) はそれ以上大きな可換部分空間に拡大できないことを 述べている。つまり VV は極大可換部分空間である。

33 次行列での上限

33 次行列では,可換部分空間の最大次元は 33 である。典型例は span{I,A,A2} \operatorname{span}\{I,A,A^2\} や対角行列全体である。試験ではこの上限と,極大性によってその上限まで拡大されることを 分けて書くと論理が明確になる。

2014年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

まず,VV は包含に関して極大な可換部分ベクトル空間である。実際,条件 (b) は VV を真に含む任意の部分空間 WW が可換でないことを意味している。 また,単位行列 IIVV に含まれていなければ V+CIV+\mathbb C I はなお可換で VV を真に含むので,極大性に反する。従って IV I\in V である。

次の事実を用いる。複素 33 次正方行列全体の中の可換な線形部分空間の次元は高々 33 である。確認しておく。可換な行列族は同時三角化できる。もしその中に最小多項式の次数が 33 の行列 AA があれば,その中心化代数は C[A]=span{I,A,A2} \mathbb C[A]=\operatorname{span}\{I,A,A^2\} であり,その可換部分空間は高々 33 次元である。そうでない場合は,各非スカラー元は 一次元部分空間と二次元部分空間からなる共通の分解または旗を保存する。二次元ブロック内の 可換部分空間は高々 22 次元であり,残りの一次元ブロックの自由度を足しても高々 33 次元である。従って dimV3. \dim V\le 3 .

一方,極大性から dimV\dim V33 未満ではあり得ない。実際,非スカラー元 AVA\in V を取る。最小多項式の次数が 33 なら,VC[A]V\subset \mathbb C[A] であり, C[A]\mathbb C[A]33 次元可換部分空間であるから,極大性より V=C[A]V=\mathbb C[A] である。 最小多項式の次数が 22 以下の場合も,Jordan標準形で見ると AA と可換な行列の中に I,AI,A と独立な行列を選び,33 次元の可換部分空間まで拡大できる。したがって dimV<3\dim V<3 なら VV は真に大きい可換部分空間に含まれ,条件 (b) に反する。

以上より dimV=3 \dim V=3 である。

最終答

dimCV=3 \boxed{\dim_{\mathbb C}V=3}

7 — Lorenz型ベクトル場と体積収縮

吸収領域の確認

境界で Φ\Phi の軌道微分が負なら,流れは外へ出ない。領域内部で微分が正になる点が あっても問題はなく,境界での向きだけを確認すればよい。

体積は発散で決まる

流れ写像のJacobianは exp(01divV(X(t))dt)=e5 \exp\left(\int_0^1 \operatorname{div}V(X(t))\,dt\right)=e^{-5} である。非線形項があっても,発散が定数であるため体積収縮率は一発で決まる。

2014年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1) ベクトル場を V=(V1,V2,V3)V=(V_1,V_2,V_3) と書くと, divV=V1x+V2y+V3z=311=5. \operatorname{div}V =\frac{\partial V_1}{\partial x} +\frac{\partial V_2}{\partial y} +\frac{\partial V_3}{\partial z} =-3-1-1=-5 .

(2) 次に Φ(x,y,z)=2x2+3y2+3z224z2 \Phi(x,y,z)=2x^2+3y^2+3z^2-24z-2 とおく。領域は Φ0\Phi\le 0 で表される。軌道に沿った微分は ddtΦ(X(t))=Φ(X(t))V(X(t))=12x26y26z2+24z=6{2x2+y2+(z2)24}. \begin{aligned} \frac{d}{dt}\Phi(X(t)) &=\nabla\Phi(X(t))\cdot V(X(t))\\ &=-12x^2-6y^2-6z^2+24z\\ &=-6\{2x^2+y^2+(z-2)^2-4\}. \end{aligned} Ω\partial\Omega 上では 2x2+3y2+3(z4)2=50 2x^2+3y^2+3(z-4)^2=50 である。この制約のもとで 2x2+y2+(z2)2 2x^2+y^2+(z-2)^2 の最小値は,x=0x=0 として yzyz 平面内の円周 y2+(z4)2=503 y^2+(z-4)^2=\frac{50}{3} 上で点 (0,2)(0,2) への距離を最小にすると得られ, (5032)2>4 \left(\sqrt{\frac{50}{3}}-2\right)^2>4 である。ここで xx 成分を増やすと制約の残りが減る一方で目的関数に 正の 2x22x^2 が加わるため,最小は x=0x=0 で実現される。従って Ω\partial\Omega 上で ddtΦ(X(t))<0 \frac{d}{dt}\Phi(X(t))<0 であり,ベクトル場は境界で内向きである。よって X(1,Ω)Ω X(1,\Omega)\subset \Omega が成り立つ。

(3) 最後に体積を求める。発散が定数 5-5 なので,Liouville の公式より時間 11 の流れは 体積を e5e^{-5} 倍する。したがって Vol(X(1,Ω))=e5Vol(Ω). \operatorname{Vol}(X(1,\Omega))=e^{-5}\operatorname{Vol}(\Omega). Ω\Omega2x2+3y2+3(z4)250 2x^2+3y^2+3(z-4)^2\le 50 であるから,半径は 5,503,503 5,\qquad \sqrt{\frac{50}{3}},\qquad \sqrt{\frac{50}{3}} である。従って Vol(Ω)=4π35503=1000π9. \operatorname{Vol}(\Omega) =\frac{4\pi}{3}\cdot 5\cdot\frac{50}{3} =\frac{1000\pi}{9}. 求める体積は Vol(ΩX(1,Ω))=(1e5)1000π9. \operatorname{Vol}(\Omega\setminus X(1,\Omega)) =(1-e^{-5})\frac{1000\pi}{9}.

最終答

(1)divV=5, \text{(1)}\quad \operatorname{div}V=-5, (2)X(1,Ω)Ω, \text{(2)}\quad X(1,\Omega)\subset\Omega, かつ (3)Vol(ΩX(1,Ω))=1000π9(1e5). \text{(3)}\quad \boxed{ \operatorname{Vol}\bigl(\Omega\setminus X(1,\Omega)\bigr) =\frac{1000\pi}{9}(1-e^{-5}) } .

東京大学 専門科目A — 他の年度