高校数学解説

Pell方程式(D=2023)を解く

302

今年の一問

$x^{2} - 2023 y^{2} = 1$ を満たす最小の自然数 $(x, y)$ の組を求めよ

はじめに

まず $x^{2} - D y^{2} = 1$
このような形であらわされる不定方程式をPell方程式と呼びます
今回は連分数展開を使った解法を紹介します

今回は先人たちの偉大な研究結果を用います
これらの事実はこの問題を解くにあたって強力な飛び道具となります
厳密な証明は大学数学の「代数的整数論・二次体の整数論」で示されます

連分数展開

まずは問題を解くののカギとなる連分数について軽く紹介します

連分数展開

ある実数 $a$ が
$a = a_{0} + \frac{1}{a_{1} + \frac{1}{a_{2} + \frac{1}{a_{3} + \frac{1}{⋱}}}}$
表されるとき $a = [a_{0}; a_{1}, a_{2}, a_{3}, \dots]$ と表現する
これを $a$ の連分数展開という

実際にこれを計算するにはその数の整数部分を $a_{0}$ として小数部分の逆数をとって、
新しく分母に来た数に同じ操作を…と繰り返していくと求められます
もちろん有理数に対してこの操作を行うと有限回で終了してしまい、面白くありません
しかし、無理数に対して行うと、どうやら面白いことがわかります

$\sqrt{n}$ の連分数展開は $a_{1}$ 以降で循環する
( $n$ は平方数でないとき)

実際に例を見てみましょう
$\sqrt{2} = [1; 2, 2, 2, 2, 2, \dots]$
$\sqrt{3} = [1; 1, 2, 1, 2, 1, \dots]$
周期は $n$ によってバラバラで、 $n = 61$ のときは周期が $11$ と非常に長くなっています

無理数の連分数近似を途中で打ち切って、普通の分数に戻すと
打ち切りが遅かった分だけ精度のいい分数の近似が得られる

どういうことでしょう
無理数の典型例 $π$ で試してみましょう
$π = [3; 7, 15, 1, \dots]$
ですので $a_{3}$ で打ち切って計算すると
$π \approx 3 + \frac{1}{7 + \frac{1}{15 + \frac{1}{1}}} = \frac{355}{113} = 3.1415929 . . .$
たった4個目までの計算で小数第6位まで一致しています!
これだけでも面白いんですが、肝心のPell方程式との関係はどうなっているんでしょうか

Pell方程式の解き方

まずは $\sqrt{2023}$ を連分数展開します
$\sqrt{2023} = [44; 1, 43, 1, 88, \dots]$ です( $1, 3, 1, 88$ の周期4で繰り返します)
$a_{3}$ (周期の最後だけ含めない)で区切って近似値を求めると $\frac{2024}{45}$
そしてなんと $(2024, 45)$ はこの方程式の最小自然数解(基本解)になっています
なのでこれが答えです