# 平成9年度 東大院 物理学専攻 修士 数学 解答

> 平成9年度 修士課程 数学（全3問）の解答。2次形式の対角化と回転行列、上半平面の曲線の微分方程式（双曲幾何の測地線）、ピタゴラス数の一般解を扱います。
> https://rikai.mugen-giken.com/exams/utokyo-physics/master-h9-math

:::caution[出典と、この解答の位置づけ]
本記事が扱う問題は、**東京大学大学院理学系研究科物理学専攻の大学院入試問題**からの引用です（引用元: [https://www.phys.s.u-tokyo.ac.jp/](https://www.phys.s.u-tokyo.ac.jp/)）。

- 掲載している**解答・解説は東京大学が公表したものではなく**、夢現技研合同会社が独自に作成したものです。誤りが含まれる可能性があります。
- 本記事に**問題文は転載していません**。問題は上記の専攻公式サイトでご確認ください。
- 本記事の内容について、東京大学および同専攻へのお問い合わせはご遠慮ください。誤りのご指摘は[お問い合わせ](/contact)までお願いします。
:::

平成9年度（試験実施は平成8年8月）の数学は全3問で、3問すべてに解答する形式です。第1問は2次形式の対角化から始めて、直交行列を回転として読む幾何の問題に発展します。第2問は連立2階微分方程式の求積で、実体はポアンカレ上半平面の測地線です。第3問はピタゴラス数の一般解を段階的に証明させる整数問題で、素因数分解の一意性を使う (c) の証明をきちんと書けるかが差になります。なお、この問題冊子には数学と物理学の両方が収録されていますが、この記事では数学の3問のみを扱います。

| 問題 | 分野 | 主題 |
|---|---|---|
| 第1問 | 線形代数 | 2次形式の対角化と回転行列の位数 |
| 第2問 | 微分方程式 | 上半平面の曲線族（半直線と半円） |
| 第3問 | 整数論 | ピタゴラス数の一般解 |

## 第1問 2次形式の対角化と回転行列の位数

### 設定

実変数 $x,y,z$ の2次形式

$$
f(x,y,z)=2x^2+3y^2+4z^2-\sqrt{6}\,xy+\sqrt{6}\,yz
$$

を考えます。$\vec r=(x,y,z)^T$ として $f=\vec r^{\,T}A\vec r$ と書いたときの対称行列 $A$ の固有値・固有ベクトルを求め、固有ベクトルを並べた直交行列 $P$ を座標回転として調べる問題です。

### 設問(a)

対称行列で書くには交差項の係数を半分ずつ振り分けます。

$$
A=\begin{pmatrix}2 & -\dfrac{\sqrt6}{2} & 0\\[4pt] -\dfrac{\sqrt6}{2} & 3 & \dfrac{\sqrt6}{2}\\[4pt] 0 & \dfrac{\sqrt6}{2} & 4\end{pmatrix}
$$

特性方程式 $\det(A-\lambda E)=0$ を第1行で展開すると

$$
\begin{aligned}
\det(A-\lambda E)&=(2-\lambda)\left[(3-\lambda)(4-\lambda)-\frac32\right]-\frac{\sqrt6}{2}\cdot\frac{\sqrt6}{2}(4-\lambda)\\
&=-\lambda^3+9\lambda^2-23\lambda+15=-(\lambda-1)(\lambda-3)(\lambda-5)
\end{aligned}
$$

となるので、固有値は $\lambda=1,3,5$ です。

$\lambda=1$: $(A-E)\vec p=0$ の第1行から $x=\frac{\sqrt6}{2}y$、第3行から $z=-\frac{\sqrt6}{6}y$ となり、$y=\sqrt6$ と取ると $\vec p\propto(3,\sqrt6,-1)^T$。ノルムの2乗は $9+6+1=16$ です。

$\lambda=3$: 第1行から $x=-\frac{\sqrt6}{2}y$、第3行から $z=-\frac{\sqrt6}{2}y$ となり、$\vec p\propto(-\sqrt6,2,-\sqrt6)^T$。ノルムの2乗は $6+4+6=16$ です。

$\lambda=5$: 第1行から $x=-\frac{\sqrt6}{6}y$、第3行から $z=\frac{\sqrt6}{2}y$ となり、$\vec p\propto(-1,\sqrt6,3)^T$。ノルムの2乗は $1+6+9=16$ です。

第2成分が非負になる符号を選んで規格化すると、答えは

$$
\lambda_1=1:\ \vec p_1=\frac14\begin{pmatrix}3\\ \sqrt6\\ -1\end{pmatrix},\qquad
\lambda_2=3:\ \vec p_2=\frac14\begin{pmatrix}-\sqrt6\\ 2\\ -\sqrt6\end{pmatrix},\qquad
\lambda_3=5:\ \vec p_3=\frac14\begin{pmatrix}-1\\ \sqrt6\\ 3\end{pmatrix}
$$

です。検算として $A\vec p_1$ の第2成分は $-\frac{\sqrt6}{2}\cdot 3+3\sqrt6-\frac{\sqrt6}{2}=\sqrt6$ となり、たしかに $\vec p_1$ の第2成分の1倍です（他成分・他固有値も同様に確認できます）。

### 設問(b)

$P=(\vec p_1,\vec p_2,\vec p_3)$ は実対称行列の異なる固有値に属する単位固有ベクトルを並べたものなので、列が正規直交系をなし、$P^TP=E$、すなわち $P$ は直交行列です（実際 $\vec p_1\cdot\vec p_2=\frac{1}{16}(-3\sqrt6+2\sqrt6+\sqrt6)=0$ などで直接確認できます）。$\vec r=P\vec r{\,}'$ を代入すると

$$
f=\vec r^{\,T}A\vec r=\vec r{\,}'^T P^TAP\,\vec r{\,}'=\vec r{\,}'^T\begin{pmatrix}1&0&0\\0&3&0\\0&0&5\end{pmatrix}\vec r{\,}'
$$

となります。$P^TAP$ が固有値を並べた対角行列になるのは $AP=P\,\mathrm{diag}(1,3,5)$ と $P^T=P^{-1}$ によります。答えは

$$
f=x'^2+3y'^2+5z'^2
$$

で、交差項のない標準形に変換されます。また、任意のベクトル $\vec v$ に対して

$$
\lVert P\vec v\rVert^2=(P\vec v)^T(P\vec v)=\vec v^{\,T}P^TP\vec v=\vec v^{\,T}\vec v=\lVert\vec v\rVert^2
$$

なので、$P$ による変換はベクトルの大きさを変えません。これが示すべきことでした。

### 設問(c)

まず $P$ が回転であることを確認します。行列式を計算すると

$$
\det P=\frac{1}{64}\left[3(6+6)+\sqrt6(3\sqrt6+\sqrt6)-(-6+2)\right]=\frac{36+24+4}{64}=1
$$

なので、$P$ は行列式 $+1$ の直交行列、すなわち回転です。回転軸の方向ベクトル $\vec n$ は $P\vec n=\vec n$ を満たすベクトルです。$(P-E)\vec n=0$ を成分で書くと、第1行と第3行はともに $-x-\sqrt6 y-z=0$、第2行は $\sqrt6 x-2y+\sqrt6 z=0$ を与えます。第2行から $y=\frac{\sqrt6}{2}(x+z)$ で、これを第1行に入れると $-(x+z)-3(x+z)=0$、よって $z=-x$、$y=0$ です。答えは

$$
\vec n=\frac{1}{\sqrt2}\begin{pmatrix}1\\0\\-1\end{pmatrix}
$$

（およびその逆向き）です。実際 $P(1,0,-1)^T=\frac14(3+1,\ \sqrt6-\sqrt6,\ -1-3)^T=(1,0,-1)^T$ となり、不変であることが確かめられます。

### 設問(d)

回転角 $\theta$ はトレースから求まります。回転行列では $\operatorname{Tr}P=1+2\cos\theta$ であり、

$$
\operatorname{Tr}P=\frac{3+2+3}{4}=2\quad\Longrightarrow\quad \cos\theta=\frac12,\qquad \theta=\frac{\pi}{3}
$$

です（回転の向きは軸の向きの取り方に依存しますが、以下の議論に影響しません）。$P^n$ は同じ軸まわりの角 $n\theta=n\pi/3$ の回転なので、$P^n=E$ となるのは $n\pi/3$ が $2\pi$ の整数倍のとき、すなわち $n$ が $6$ の倍数のときに限ります。$n=1,\dots,5$ では回転角が $2\pi$ の整数倍にならないため $P^n\neq E$ です。よって最小の自然数は $n=6$ です。

## 第2問 上半平面の曲線族（半直線と半円）

### 設定

上半平面 $-\infty<x<\infty,\ y>0$ 内のなめらかな曲線を実パラメータ $s$ で $(x(s),y(s))$ と表し、連立微分方程式

$$
\begin{cases}
y\,x''-2\,x'y'=0\\
y\,y''+x'^2-y'^2=0
\end{cases}
$$

を満たすとします。$'$ は $s$ 微分です。これはポアンカレ上半平面の測地線の方程式で、解曲線が $x$ 軸に垂直な半直線と $x$ 軸上に中心を持つ半円になることを、変数変換 $X=x'/y$、$Y=y'/y$ を通して導く問題です。

### 設問(a)

$X=x'/y$ を $s$ で微分し、第1式 $x''=2x'y'/y$ を使うと

$$
X'=\frac{x''}{y}-\frac{x'y'}{y^2}=\frac{2x'y'}{y^2}-\frac{x'y'}{y^2}=\frac{x'y'}{y^2}=XY
$$

となります。同様に $Y=y'/y$ を微分し、第2式 $y''=(y'^2-x'^2)/y$ を使うと

$$
Y'=\frac{y''}{y}-\frac{y'^2}{y^2}=\frac{y'^2-x'^2}{y^2}-\frac{y'^2}{y^2}=-\frac{x'^2}{y^2}=-X^2
$$

です。答えは

$$
X'=XY,\qquad Y'=-X^2
$$

です。この2式から

$$
\frac{d}{ds}\left(X^2+Y^2\right)=2XX'+2YY'=2X^2Y-2X^2Y=0
$$

となるので、$X^2+Y^2$ は $s$ によらない定数 $C$ に等しく、$X^2+Y^2=C$ が成り立ちます。これが示すべきことでした。

### 設問(b)

$C=1$ とします。このとき $x'^2+y'^2=y^2$ であり、$s$ は計量 $ds^2=(dx^2+dy^2)/y^2$ に関する弧長パラメータです。

$X=0$ の場合。$X'=XY$ から $X$ は恒等的に $0$ で（ある $s$ で $X=0$ なら初期値 $0$ の線形方程式の解として常に $0$）、$x'=0$、すなわち $x=x_0$（定数）です。$Y^2=1$ から $y'/y=\pm1$ で、$y=y_0e^{\pm s}$（$y_0>0$）となります。曲線は

$$
x=x_0,\qquad y>0
$$

の、$x$ 軸に垂直な半直線の族です。$y=y_0e^{\pm s}$ は $y\,y''-y'^2=y^2-y^2=0$ と第1式（両辺 $0$）を満たすことが直接確認できます。

$X\neq0$ の場合。$Y'=-X^2=Y^2-1$ を解きます。$\lvert Y\rvert<1$（$X\neq0$ より）なので変数分離して

$$
\int\frac{dY}{Y^2-1}=s+\text{const.}\quad\Longrightarrow\quad Y=-\tanh(s-s_0)
$$

を得ます（$s_0$ は積分定数）。すると $X^2=1-Y^2=\operatorname{sech}^2(s-s_0)$ より $X=\pm\operatorname{sech}(s-s_0)$ です（符号は $s$ によらず一定。$X$ は連続で $0$ にならないためです）。$u=s-s_0$ と書きます。$y'/y=Y=-\tanh u$ を積分して

$$
y=\frac{R}{\cosh u}\qquad(R>0)
$$

となり、$x'=Xy=\pm R\operatorname{sech}^2u$ を積分して

$$
x=x_0\pm R\tanh u
$$

を得ます。この2式から $u$ を消去すると

$$
(x-x_0)^2+y^2=R^2\tanh^2u+R^2\operatorname{sech}^2u=R^2
$$

です。よって曲線は、中心 $(x_0,0)$ を $x$ 軸上に持つ半径 $R$ の半円（上半平面にある部分）の族です。中心の位置 $x_0$ と半径 $R>0$ は任意です。逆にこの解が元の方程式を満たすことは、$x'=R\operatorname{sech}^2u$、$y'=-R\operatorname{sech}u\tanh u$、$x''=-2R\operatorname{sech}^2u\tanh u$、$y''=R\operatorname{sech}u(\tanh^2u-\operatorname{sech}^2u)$ を代入すると、第1式は $-2R^2\operatorname{sech}^3u\tanh u+2R^2\operatorname{sech}^3u\tanh u=0$、第2式は $R^2\operatorname{sech}^2u(\tanh^2u-\operatorname{sech}^2u)+R^2\operatorname{sech}^4u-R^2\operatorname{sech}^2u\tanh^2u=0$ となって確かめられます。

まとめると、答えは次のとおりです。$X=0$ のとき、曲線群は $x$ 軸に垂直な半直線 $x=x_0$（$y>0$）。$X\neq0$ のとき、曲線群は $x$ 軸上に中心を持つ半円 $(x-x_0)^2+y^2=R^2$（$y>0$）。これらはポアンカレ上半平面模型の測地線にほかなりません。

## 第3問 ピタゴラス数の一般解

### 設定

$x^2+y^2=z^2$ を満たす正の整数 $x,y,z$ で、どの2つも互いに素（最大公約数が 1）であるものをすべて求める方法を、小問を積み重ねて構成します。以下、$\gcd$ で最大公約数を表します。

### 設問(a)

$x$ と $y$ の一方が偶数、他方が奇数であることを示します。

まず、$x,y$ がともに偶数だとすると $\gcd(x,y)\geq2$ となり、互いに素という仮定に反します。

次に、$x,y$ がともに奇数だとします。奇数の2乗は $(2k+1)^2=4k(k+1)+1\equiv1\pmod 4$ なので、$x^2+y^2\equiv2\pmod 4$ です。一方 $z^2$ は、$z$ が偶数なら $z^2\equiv0\pmod4$、奇数なら $z^2\equiv1\pmod4$ であり、いずれにせよ $2\pmod4$ にはなりません。よって $x^2+y^2=z^2$ と矛盾します。

したがって残る可能性は、一方が偶数でもう一方が奇数の場合だけです。これが示すべきことでした。なおこのとき $z^2=x^2+y^2$ は奇数なので $z$ も奇数です。

### 設問(b)

偶数の方を $x$ とすると、(a) の帰結から $y$ と $z$ はともに奇数です。よって $z+y$ と $z-y$ はともに偶数で、

$$
A=\frac{z+y}{2},\qquad B=\frac{z-y}{2}
$$

はどちらも整数です。また $z>y>0$（$z^2=x^2+y^2>y^2$ と正値性から）なので $A>B>0$ であり、$A,B$ は正の整数です。

互いに素であることを示します。$d$ を $A$ と $B$ の公約数とすると、$d$ は $A+B=z$ と $A-B=y$ の両方を割ります。仮定より $\gcd(y,z)=1$ なので $d=1$、すなわち $\gcd(A,B)=1$ です。これが示すべきことでした。

あとで使うので付記すると、$A$ と $B$ の積は

$$
AB=\frac{z^2-y^2}{4}=\frac{x^2}{4}=\left(\frac{x}{2}\right)^2
$$

となり、$x$ が偶数なので $x/2$ は正の整数です。つまり互いに素な $A,B$ の積が平方数になっています。

### 設問(c)

互いに素な正の整数 $A,B$ が $AB=C^2$（$C$ は正の整数）を満たすとき、$A=\alpha^2$、$B=\beta^2$ となる整数 $\alpha,\beta$ が存在することを、素因数分解の一意性を用いて証明します。

$p$ を任意の素数とし、整数 $m$ の素因数分解における $p$ の指数を $\operatorname{ord}_p m$ と書きます。素因数分解の一意性から、指数は乗法的で

$$
\operatorname{ord}_p(AB)=\operatorname{ord}_p A+\operatorname{ord}_p B=\operatorname{ord}_p(C^2)=2\operatorname{ord}_p C
$$

が成り立ちます。ここで $\gcd(A,B)=1$ なので、$p$ が $A$ と $B$ の両方を割ることはなく、$\operatorname{ord}_p A$ と $\operatorname{ord}_p B$ の少なくとも一方は $0$ です。したがって $\operatorname{ord}_p A\neq0$ ならば $\operatorname{ord}_p B=0$ であり、上式から $\operatorname{ord}_p A=2\operatorname{ord}_p C$ は偶数です。すなわち $A$ の素因数分解に現れるすべての素数の指数は偶数であり、

$$
\alpha=\prod_p p^{\operatorname{ord}_p A/2}
$$

とおけば $A=\alpha^2$ です。まったく同様に $B$ のすべての素因数の指数も偶数で、$B=\beta^2$ となる整数 $\beta$ が存在します。これが示すべきことでした。

### 設問(d)

(b) の $A,B$ は互いに素な正の整数で $AB=(x/2)^2$ を満たすので、(c) より $A=\alpha^2$、$B=\beta^2$ となる正の整数 $\alpha,\beta$ が取れます。$A>B>0$ から $\alpha>\beta>0$ です。このとき

$$
z=A+B=\alpha^2+\beta^2,\qquad y=A-B=\alpha^2-\beta^2,\qquad \left(\frac{x}{2}\right)^2=AB=\alpha^2\beta^2
$$

なので、答えは

$$
x=2\alpha\beta,\qquad y=\alpha^2-\beta^2,\qquad z=\alpha^2+\beta^2\qquad(\alpha>\beta>0)
$$

です。実際 $x^2+y^2=4\alpha^2\beta^2+(\alpha^2-\beta^2)^2=(\alpha^2+\beta^2)^2=z^2$ が恒等式として成り立ちます。なお $x,y,z$ がどの2つも互いに素になるためには、$\gcd(\alpha,\beta)=1$ かつ $\alpha,\beta$ の偶奇が異なることが必要十分です。$\gcd(\alpha,\beta)=1$ は $\gcd(A,B)=1$ から従い、$\alpha,\beta$ がともに奇数だと $y,z$ がともに偶数になって互いに素でなくなるため、偶奇は異なります。

### 設問(e)

$\alpha\leq6$ の範囲で、(d) の式に $\alpha>\beta>0$、$\gcd(\alpha,\beta)=1$、偶奇が異なる、という条件を満たす組 $(\alpha,\beta)$ をすべて代入します。除外されるのは、偶奇が一致する $(3,1)$、$(5,1)$、$(5,3)$（ともに奇数。たとえば $(3,1)$ は $(x,y,z)=(6,8,10)$ となり互いに素でない）と、公約数を持つ $(4,2)$、$(6,2)$、$(6,3)$、$(6,4)$ です。残る組と得られる $(x,y,z)$ は次のとおりです。

| $(\alpha,\beta)$ | $x=2\alpha\beta$ | $y=\alpha^2-\beta^2$ | $z=\alpha^2+\beta^2$ |
|---|---|---|---|
| $(2,1)$ | $4$ | $3$ | $5$ |
| $(3,2)$ | $12$ | $5$ | $13$ |
| $(4,1)$ | $8$ | $15$ | $17$ |
| $(4,3)$ | $24$ | $7$ | $25$ |
| $(5,2)$ | $20$ | $21$ | $29$ |
| $(5,4)$ | $40$ | $9$ | $41$ |
| $(6,1)$ | $12$ | $35$ | $37$ |
| $(6,5)$ | $60$ | $11$ | $61$ |

いずれも $x^2+y^2=z^2$ を満たし（たとえば $60^2+11^2=3600+121=3721=61^2$）、どの2つも互いに素です。答えはこの8組です。

出典: 東京大学大学院理学系研究科 物理学専攻 平成9年度 修士課程 入学試験問題 数学。問題文は要約して引用しています。