MathLabs

算術と数論

ペル方程式

整数解が連分数と結びつく方程式 x2−dy2=1x^2-dy^2=1。

直観直感:双曲線上で d\sqrt{d} を近似する整数点

x2−2y2=1x^2-2y^2=1 を満たす整数 x,yx,y を探す。小さい値を試すと (x,y)=(3,2)(x,y)=(3,2) が成り立つ、なぜなら 9−8=19-8=1 だからである。整理すると (xy)2=2+1y2\left(\frac{x}{y}\right)^2=2+\frac{1}{y^2} となり、x/y=3/2=1.5x/y=3/2=1.5 はすでに 2≈1.41421\sqrt{2}\approx 1.41421 への驚くほど良い近似である。x2−dy2=1x^2 - d y^2 = 1 の形の方程式——ペル方程式と呼ばれる——は連分数と最良有理近似の理論と切り離せないことが分かる。

連分数の近似分数のパスグラフで、ペル方程式の基本解として一つの頂点が強調されている。
d\sqrt{d} の連分数の近似分数をパスグラフとして描く:各頂点は近似分数 pk/qkp_k/q_k であり、強調された頂点が基本解である。

発展定義、基本解、構造

定義: ペル方程式と基本解

完全平方数でない正の整数 dd に対し、x2−dy2=1x^2 - d y^2 = 1 をペル方程式と呼ぶ。その正の整数解 (x,y)(x,y) のうち xx が最小のもの(同値に yy が最小のもの)を基本解 (x1,y1)(x_1,y_1) と呼ぶ。

x2−dy2=1,d∈Z>0 not a perfect squarex^2 - d y^2 = 1,\qquad d\in\mathbb{Z}_{>0}\ \text{not a perfect square}

左辺を (x−yd)(x+yd)=1(x-y\sqrt d)(x+y\sqrt d)=1 と因数分解すると、ペル方程式は環 Z[d]\mathbb{Z}[\sqrt d] に関する主張になる:解はノルム 11 の元に対応し、そのような二つの元を掛けると別のものが得られる。特に (x1,y1)(x_1,y_1) が分かれば、それ以降のすべての解はべき乗を取ることで得られる:xn+ynd=(x1+y1d)nx_n+y_n\sqrt{d} = (x_1+y_1\sqrt{d})^n。

(x1+y1d)(x1−y1d)=1 ⟹ xn+ynd=(x1+y1d)n(x_1+y_1\sqrt d)(x_1-y_1\sqrt d)=1 \ \Longrightarrow\ x_n+y_n\sqrt d = (x_1+y_1\sqrt d)^n
小さい dd に対する基本解
ddd\sqrt d の連分数基本解 (x1,y1)(x_1,y_1)
22[1;2‾][1;\overline{2}](3,2)(3,2)
33[1;1,2‾][1;\overline{1,2}](2,1)(2,1)
55[2;4‾][2;\overline{4}](9,4)(9,4)
77[2;1,1,1,4‾][2;\overline{1,1,1,4}](8,3)(8,3)

発展定理と証明

完全平方数でないすべての正の整数 dd に対し、x2−dy2=1x^2 - d y^2 = 1 は正の整数 x,yx,y による解を持つ。

なぜ正しいのか?

双曲線 x2−dy2=1x^2-dy^2=1 が実数点を持つことは確かだとしても、それが格子点を通ることは全く自明ではない——この定理は、有理近似に関する巧妙な鳩の巣論法により、平方数でないすべての dd に対して常にそうなることを保証する。

証明

ディリクレの近似定理により、任意の整数 Q>0Q>0 に対し 1≤q≤Q1\le q\le Q かつ ∣d−pq∣<1qQ\left|\sqrt d - \frac{p}{q}\right|<\frac{1}{qQ} を満たす整数 p,qp,q が存在する。Q→∞Q\to\infty とすると ∣d−pq∣<1q2\left|\sqrt d-\frac{p}{q}\right|<\frac{1}{q^2} を満たす組 (p,q)(p,q) が無限に得られる。

そのような各組に対し ∣p−qd∣<1q|p-q\sqrt d|<\frac{1}{q} なので ∣p2−dq2∣=∣p−qd∣⋅∣p+qd∣<1q(2qd+1q)<2d+1|p^2-dq^2|=|p-q\sqrt d|\cdot|p+q\sqrt d|<\frac{1}{q}\left(2q\sqrt d+\frac{1}{q}\right)<2\sqrt d+1 となる。よって p2−dq2p^2-dq^2 は区間 (−2d−1, 2d+1)(-2\sqrt d-1,\,2\sqrt d+1) 内の有限個の整数値しか取らないが、組 (p,q)(p,q) は無限に存在する。

鳩の巣原理により、その有限範囲内のある固定された非零整数 kk が、無限個の組 (p,q)(p,q) に対し p2−dq2=kp^2-dq^2=k を満たす。これら無限個の組の中で、再び鳩の巣原理により、無限個が同じ剰余類 p≡p0, q≡q0(mod∣k∣)p\equiv p_0,\ q\equiv q_0\pmod{|k|} を持つ。

そのような相異なる二組 (p1,q1)≠(p2,q2)(p_1,q_1)\neq(p_2,q_2) で p12−dq12=p22−dq22=kp_1^2-dq_1^2=p_2^2-dq_2^2=k かつ ∣k∣|k| を法として剰余が一致するものを取る。x+yd=(p1+q1d)(p2−q2d)kx+y\sqrt d = \frac{(p_1+q_1\sqrt d)(p_2-q_2\sqrt d)}{k} とおく;展開すると x=p1p2−dq1q2kx=\frac{p_1p_2-dq_1q_2}{k} と y=p1q2−p2q1ky=\frac{p_1q_2-p_2q_1}{k} は、∣k∣|k| を法とする剰余が一致することにより正真正銘の整数であることが分かり、ノルム N(a+bd)=a2−db2N(a+b\sqrt d)=a^2-db^2 の乗法性より x2−dy2=k⋅kk2=1x^2-dy^2=\frac{k\cdot k}{k^2}=1 が得られる。(p1,q1)≠(p2,q2)(p_1,q_1)\neq(p_2,q_2) だが極限で同じ比を与えるので y≠0y\neq 0 が確認でき、(x,y)(x,y) を (∣x∣,∣y∣)(|x|,|y|) に置き換えても(平方しか現れないので依然解である)正の整数解が得られる。

(x1,y1)(x_1,y_1) が x2−dy2=1x^2 - d y^2 = 1 の基本解であるとき、すべての正整数解 (x,y)(x,y) はある n≥1n\ge 1 に対する (xn,yn)(x_n,y_n) に等しい。ここで xn+ynd=(x1+y1d)nx_n+y_n\sqrt{d} = (x_1+y_1\sqrt{d})^n である。

なぜ正しいのか?

これは、無限に多い解が謎めいた散在集合ではなく、最小の解を繰り返し「掛け合わせる」ことで生成される完全に予測可能な等比数列的な列であることを述べている——無限の探索を一つの数を見つけることに帰着させる。

証明

まず xn+ynd=(x1+y1d)nx_n+y_n\sqrt d=(x_1+y_1\sqrt d)^n で定義される (xn,yn)(x_n,y_n) がすべての nn で確かに解であることを確認する:共役を取ると xn−ynd=(x1−y1d)nx_n-y_n\sqrt d=(x_1-y_1\sqrt d)^n なので xn2−dyn2=(xn+ynd)(xn−ynd)=[(x1+y1d)(x1−y1d)]n=(x12−dy12)n=1n=1x_n^2-dy_n^2=(x_n+y_n\sqrt d)(x_n-y_n\sqrt d)=\left[(x_1+y_1\sqrt d)(x_1-y_1\sqrt d)\right]^n=(x_1^2-dy_1^2)^n=1^n=1 となる。

次に、この形でない任意の正整数解 (x,y)(x,y) があると仮定する;n→∞n\to\infty で xn→∞x_n\to\infty なので、xn+ynd≤x+yd<xn+1+yn+1d=(xn+ynd)(x1+y1d)x_n+y_n\sqrt d \le x+y\sqrt d < x_{n+1}+y_{n+1}\sqrt d = (x_n+y_n\sqrt d)(x_1+y_1\sqrt d) を満たす一意な nn が存在する。

両辺を (xn+ynd)(x_n+y_n\sqrt d) で割る、すなわちその逆元 (xn−ynd)(x_n-y_n\sqrt d)(xn2−dyn2=1x_n^2-dy_n^2=1 より有効)を掛ける:x′+y′d=(x+yd)(xn−ynd)x'+y'\sqrt d = (x+y\sqrt d)(x_n-y_n\sqrt d) とおく。すると 1≤x′+y′d<x1+y1d1\le x'+y'\sqrt d < x_1+y_1\sqrt d であり、x′2−dy′2=(x2−dy2)(xn2−dyn2)=1⋅1=1x'^2-dy'^2=(x^2-dy^2)(x_n^2-dy_n^2)=1\cdot 1=1 となるので (x′,y′)(x',y') もペル方程式の解である。

x′+y′d≥1x'+y'\sqrt d\ge 1 と x′2−dy′2=1x'^2-dy'^2=1 を用いた簡単な計算により x′≥1x'\ge 1、y′≥0y'\ge 0 が分かる(x′+y′d≥1x'+y'\sqrt d\ge 1 だが y′<0y'<0 の解は x′>x1x'>x_1 を強制し、x′+y′d<x1+y1dx'+y'\sqrt d<x_1+y_1\sqrt d とノルムの式に矛盾する)。もし y′>0y'>0 なら、(x′,y′)(x',y') は x′+y′d<x1+y1dx'+y'\sqrt d<x_1+y_1\sqrt d を満たす正の解となり、基本解の最小性に矛盾する。よって y′=0y'=0、したがって x′=1x'=1 となり、x+yd=xn+yndx+y\sqrt d=x_n+y_n\sqrt d、すなわち結局 (x,y)=(xn,yn)(x,y)=(x_n,y_n) となる——これは仮定に矛盾する。

したがってすべての正の解はちょうどある (xn,yn)(x_n,y_n) であり、基本解が解全体の集合を生成することが証明された。

発展実世界での応用と具体例

ペル方程式は単なるパズルではない:機械式歯車設計で使われる最良有理近似を支配し、その数論は古典的な整数分解アルゴリズムの基礎になっており、大きな dd に対する基本解の計算は、まさにHallgrenの量子アルゴリズムが多項式時間で解く計算問題である——効率的な古典的方法は知られていない。

例: d=2d=2 の基本解を求める

2\sqrt{2} の連分数を用いて x2−2y2=1x^2-2y^2=1 の基本解を求め、次の解を生成せよ。

解答

2\sqrt2 の連分数は [1;2‾]=1+12+12+⋯[1;\overline{2}]=1+\cfrac{1}{2+\cfrac{1}{2+\cdots}} であり、近似分数は 1, 32, 75, 1712, 4129,…1,\ \frac{3}{2},\ \frac{7}{5},\ \frac{17}{12},\ \frac{41}{29},\dots である。

近似分数 32\frac{3}{2} を確認する:32−2⋅22=9−8=13^2-2\cdot 2^2=9-8=1。これが成り立つ最小の近似分数なので、基本解は (x1,y1)=(3,2)(x_1,y_1)=(3,2) である。

漸化式 xn+ynd=(x1+y1d)nx_n+y_n\sqrt{d} = (x_1+y_1\sqrt{d})^n を n=2n=2 で使う:x2+y22=(3+22)2=9+122+8=17+122x_2+y_2\sqrt2=(3+2\sqrt2)^2=9+12\sqrt2+8=17+12\sqrt2、よって (x2,y2)=(17,12)(x_2,y_2)=(17,12) である。

検算:172−2⋅122=289−288=117^2-2\cdot 12^2=289-288=1、これで次の解が確認でき、上で見つけた近似分数 1712\frac{17}{12} と一致する——理論の予測どおりである。

例: 工学的近似:2\sqrt{2} に近い歯車比の設計

ある機構は、小さな整数の歯数を用いて歯数比が 2\sqrt2 にできるだけ近い、噛み合う二つの歯車を必要とし、多くの回転の後でも誤差が小さいままであるようにしたい。ペル方程式の近似分数を使って歯数を選び、誤差を評価せよ。

解答

x2−2y2=1x^2-2y^2=1 の基本解 (3,2)(3,2) から、比 3/2=1.53/2=1.5 は ∣3/2−2∣≈0.0858|3/2-\sqrt2|\approx 0.0858 を与える——使えるが、精密機械には粗い。

3+223+2\sqrt2 を二乗して得られる次の近似分数 (17,12)(17,12) を使うと:歯数 1717 と 1212 の歯車対は比 17/12≈1.4166717/12\approx 1.41667 を与え、∣17/12−2∣≈0.00245|17/12-\sqrt2|\approx 0.00245 となり、歯数がわずかに増えるだけで 3/23/2 設計よりおよそ 3535 倍正確になる。

連分数の任意の近似分数 pk/qkp_k/q_k に対する一般的な評価は ∣pkqk−2∣<1qk2\left|\frac{p_k}{q_k}-\sqrt2\right|<\frac{1}{q_k^2} であり、技術者が次のペル解 (x3,y3)(x_3,y_3)((3+22)3=99+702(3+2\sqrt2)^3=99+70\sqrt2、すなわち 99/7099/70 から得られる)に移ると歯数は 7070 に増えるが、誤差は 1/702≈0.00021/70^2\approx 0.0002 未満に縮小する。

これはペル方程式に明示される工学上の基本的なトレードオフを示している:連続する各解 (xn,yn)(x_n,y_n) は、製造コストの増加(歯数の増加、すなわちより大きな yny_n)と引き換えに、定量化可能で急速に縮小する近似誤差を得られ、設計者は予算と精度要件に合うこの曲線上の点を正確に選ぶことができる。

x2−2y2=1x^2-2y^2=1 の基本解は何か?

(x1,y1)(x_1,y_1) が基本解のとき、解 (x2,y2)(x_2,y_2) はどのように得られるか?

x2−3y2=−1x^2-3y^2=-1 に整数解はあるか?

Hallgrenのアルゴリズムはどのような計算を用いてペル方程式を多項式時間で解くか?

参考文献

  1. Sean Hallgren (2007). Polynomial-time quantum algorithms for Pell's equation and the principal ideal problem · DOI:10.1145/1206035.1206039
  2. Hendrik W. Lenstra Jr. (2002). Solving the Pell Equation