1 定義と解の合成
定義 1.1 (ペル方程式).Dを平方数でない正の整数とする。整数解(x,y)を求める方程式
x2−Dy2=1をペル方程式 (Pell equation) という。(x,y)=(±1,0)は常に解であり、これを自明解という。
Dが平方数でD=k2と書くことができる場合には、
x2−Dy2=(x−ky)(x+ky)=1
となります。二つの整数の積が1になるのは両方が1、または両方が−1の場合だけなので、この場合には自明解しかありません。したがって、Dが平方数でないという仮定が必要です。
解の合成を記述するため、α=x+yDに対してαˉ=x−yD、N(α)=ααˉとおきます。このときN(α)=x2−Dy2です。
定理 1.2.Dを平方数でない正整数とする。(x1,y1)と(x2,y2)がともにx2−Dy2=1の整数解ならば、
(x1x2+Dy1y2, x1y2+x2y1)もペル方程式の解である。
証明.α1=x1+y1D、α2=x2+y2Dとおく。積を展開すると
α1α2=(x1+y1D)(x2+y2D)=(x1x2+Dy1y2)+(x1y2+x2y1)D,となる。共役をとると
α1α2=(x1x2+Dy1y2)−(x1y2+x2y1)D=(x1−y1D)(x2−y2D)=αˉ1αˉ2である。したがって
N(α1α2)=α1α2⋅α1α2=α1α2⋅αˉ1αˉ2=(α1αˉ1)(α2αˉ2)=N(α1)N(α2).が成り立つ。仮定よりN(α1)=N(α2)=1であるから、N(α1α2)=1である。これは主張の整数の組がx2−Dy2=1を満たすことを意味する。▨
非自明な正の解を一つ得ると、定理 1.2を繰り返し適用して無限個の解を作ることができます。残る問題は、最初の正の解を連分数から構成することと、その構成ですべての解を得ることができることを示すことです。
2 平方根の周期連分数
定理 2.1.Dを平方数でない正整数とし、a0=⌊D⌋とする。整数の組(Pk,Qk)と整数akを
P0=0,Q0=1,ak=⌊QkD+Pk⌋,Pk+1=akQk−Pk,Qk+1=QkD−Pk+12によって定める。この計算では常に
0≤Pk<D,Qk>0,Qk∣D−Pk2が成り立つ。また、(Pℓ,Qℓ)=(a0,1)となる最小の整数ℓ≥1が存在し、
D=[a0;a1,a2,…,aℓ],aℓ=2a0である。したがって、(Pk,Qk)=(a0,1)を初めて得た時点が一周期の計算の停止条件になる。
証明.P0=0、Q0=1とし、表示された漸化式によって各項を順に定める。Q0はD−P02を割る。Qk∣D−Pk2が成り立つとき、Pk+1≡−Pk(modQk)であるから、Qk∣D−Pk+12である。したがってQk+1は整数である。
αk=(D+Pk)/Qkとおく。ak=⌊αk⌋とPk+1=akQk−Pkを用いて分母を有理化すると
αk−ak1=D−Pk+1Qk=Qk+1D+Pk+1=αk+1となる。k=1ではP1=a0、Q1=D−a02>0であり、
Q1P1−D=−D+a01は−1より大きく0より小さい。またα1=1/(D−a0)>1である。αk>1かつ−1<(Pk−D)/Qk<0と仮定する。αk−akは0より大きく1より小さいので、
0<D−Pk+1<Qkである。また(Pk−D)/Qk−ak<−1であるからD+Pk+1>Qk>0である。よって∣Pk+1∣<Dであり、Qk+1=(D−Pk+12)/Qk>0である。さらに
αk+1=αk−ak1>1,Qk+1Pk+1−D=(Pk−D)/Qk−ak1である。右辺の分母は−1より小さいので
−1<Qk+1Pk+1−D<0となる。αk+1>1と最後の不等式を加え、Qk+1>0を用いるとPk+1>0が得られる。帰納法により、主張した整数性、正値、整除関係がすべてのkで成り立つ。
k≥1では、整数Pkは0<Pk<Dを満たし、正整数QkはD−Pk2の約数である。したがって、現れる組(Pk,Qk)の候補は有限個しかない。さらに−1<(Pk−D)/QkからPk+Qk>Dである。次の状態(P′,Q′)が分かれば、直前の状態の第二成分は
Q=Q′D−(P′)2によって定まる。直前の第一成分PはP≡−P′(modQ)とD−Q<P<Dを満たす。この長さQの区間には、同じ剰余をもつ整数は高々一つしかない。したがって、状態の移り方には逆向きにも分岐がない。有限個の状態を進む列は、最初の状態(P1,Q1)へ戻ってから同じ列を繰り返す。
(P1,Q1)=(a0,D−a02)の直前の状態ではQ=(D−a02)/Q1=1である。またD−1<P<Dを満たす整数はP=a0だけである。よって、ある最小のℓ≥1で(Pℓ,Qℓ)=(a0,1)となり、その次に(P1,Q1)が再び現れる。akは状態から一意に定まるので、a1,…,aℓが繰り返される。最後にaℓ=⌊D+a0⌋=2a0である。▨
定理 2.2.Dを平方数でない正整数とし、定理 2.1の周期長をℓとする。近似分数pk/qk=[a0;a1,…,ak]の分子と分母を§A4.11 命題 3.2の漸化式によって定めると、すべてのk≥0について
pk2−Dqk2=(−1)k+1Qk+1が成り立つ。したがって、pk2−Dqk2=1が初めて成り立つ添字は、ℓが偶数ならk=ℓ−1、ℓが奇数ならk=2ℓ−1である。
証明. 完全商αk+1=(D+Pk+1)/Qk+1に§A4.11 命題 3.2を適用すると
D=qkαk+1+qk−1pkαk+1+pk−1となる。αk+1を代入し、有理数部分とDの係数を比較すると
pk=Pk+1qk+Qk+1qk−1,Dqk=Pk+1pk+Qk+1pk−1を得る。したがって
pk2−Dqk2=pk(pk−Pk+1qk)−Qk+1pk−1qk=Qk+1(pkqk−1−pk−1qk)=(−1)k+1Qk+1.Qk+1=1なら、定理 2.1の証明で得たD−Qk+1<Pk+1<DによりPk+1=a0である。したがってQk+1=1となるのはk+1がℓの正の倍数である場合に限る。(−1)k+1=1も必要なので、条件を満たす最小のk+1は、ℓが偶数ならℓ、ℓが奇数なら2ℓである。
D=2では(P1,Q1)=(1,1)でℓ=1であり、指定された添字k=1ではp1/q1=[1;2]=3/2、32−2⋅22=1となる。D=3では(P1,Q1)=(1,2)、(P2,Q2)=(1,1)でℓ=2であり、指定された添字k=1ではp1/q1=[1;1]=2/1、22−3⋅12=1となる。二つの小さい場合でも添字は一致する。▨
3 基本解と全解の生成
補題 3.1.θを無理数、cを整数、dを正整数とし、c/dは既約であるとする。
θ−dc<2d21ならば、c/dはθの単純連分数の近似分数である。
証明.θ=[a0;a1,a2,…]の近似分数の分子と分母を、§A4.11 命題 3.2の漸化式によってpk,qkと定める。q0=1、q1=a1≥1であり、k≥0ではqk+2≥qk+1+qkである。したがって、qkは広義に増加し、上に有界でない。qj≤dを満たす最大の添字j≥0を取ると、qj≤d<qj+1となる。
Q=qj、R=qj−1、t=[aj+1;aj+2,…]、a=⌊t⌋とおく。§A4.11 命題 3.2により
θ=Qt+Rpjt+pj−1,qj+1=aQ+R,であり、0≤R≤Qである。pjR−pj−1Q=±1なので、ある整数U,Vが存在して
(c,d)=U(pj,Q)+V(pj−1,R)と書くことができる。したがって
θ−dc=d(Qt+R)∣U−Vt∣である。V=0なら、c/dが既約でd>0であることからU=1となり、c/d=pj/qjである。
V=0とする。V<0ならd=UQ+VR≥QからU>0であり、2d∣U−Vt∣≥2Q(1+t)>Qt+Rである。V>0かつU≤0なら2d∣U−Vt∣≥2Qt>Qt+Rである。U,V>0なら、d<aQ+Rから1≤U≤a−1であり、a≥2である。V≥1ではUQ+VR>0、Vt−U≥t−U>0であり、前者はVとともに減少せず、後者は狭義に増加する。したがって、2(UQ+VR)(Vt−U)はVとともに増加し、V=1の場合を調べれば足りる。r=R/Qとおくと、必要な不等式は
2(U+r)(t−U)≥t+rである。左辺と右辺の差はUについて下向きに開く二次式なので、1≤U≤a−1における最小値は端点で得られる。t=a+f、0<f<1と書くと、U=1での差は
t−2+r(2t−3)≥0,U=a−1での差は
a−2+(2a−3)f+r(1+2f)≥0である。よってV=0のすべての場合に
θ−dc≥2d21となり、仮定に反する。したがってV=0であり、c/dは近似分数である。▨
定理 3.2 (ラグランジュの定理).Dを平方数でない正整数とする。ペル方程式は非自明な正整数解をもち、x+yD>1が最小となる正整数解(x,y)が存在する。この解を基本解と呼ぶ。定理 2.2で指定される近似分数から基本解が得られる。基本解をε=x1+y1Dと書けば、すべての正整数解は
εn=xn+ynD(n=1,2,3,…)によって尽くされる。整数解全体は(±1,0)と、n≥1に対する(±xn,±yn)であり、二つの符号は独立に選ぶことができる。
証明.定理 2.2で指定される近似分数をp/qとすれば、p2−Dq2=1であり、p,qは正整数である。したがって非自明な正整数解が存在する。
正整数解(x,y)に対してx2−Dy2=1ならgcd(x,y)=1であり、
0<D−yx=y2(D+x/y)1<2y21である。補題 3.1により、ペル方程式のすべての正整数解x/yはDの近似分数に現れる。定理 2.2では値1が初めて得られる添字が特定されているので、そこで得る解をε=x1+y1Dとすれば、x1+y1D>1は正整数解の中で最小である。実際、t>1に対して(t−t−1)/(2D)は狭義に増加し、解ではy=(t−t−1)/(2D)であるから、より小さい解はより小さい分母をもつ。近似分数の分母は漸化式により添字とともに広義に増加し、指定された近似分数より小さい分母はそれより前の添字にしか現れない。したがって、そのような解は存在しない。
定理 1.2の合成式により、εのすべての正整数乗は正整数係数をもち、N(εn)=1(n≥1)を満たす。α=x+yDを任意の正整数解とする。ε>1であり、εn≥1+n(ε−1)であるから、ある整数r≥0が存在して
εr≤α<εr+1となる。εr=A+BDと書けばε−r=A−BDであり、
β=αε−r=(xA−DyB)+(yA−xB)Dは整数係数をもち、1≤β<εかつββˉ=1を満たす。β=X+YDと書くと
X=2β+β−1≥1,Y=2Dβ−β−1≥0である。β>1なら(X,Y)は基本解より小さい非自明な正整数解となるので、基本解の最小性に反する。よってβ=1であり、α=εrである。α>1なので、r≥1である。
最後に、非自明な整数解(x,y)ではxとyはともに0でない。方程式はxとyのそれぞれの符号を変えても保たれるので、(∣x∣,∣y∣)に正整数解の分類を適用することができる。これに自明解(±1,0)を加えると、主張した整数解全体が得られる。▨
4 具体例
例 4.1.D=2では2=[1;2]=[1;2,2,2,…]である。近似分数を計算すると
[1]=1,[1;2]=1+21=23.3/2に対応する(x,y)=(3,2)について
32−2×22=9−8=1であるから、この組が基本解である。定理 1.2により基本解を合成すると
(3+22)2=9+122+8=17+122であるから、次の解は(x,y)=(17,12)である。元の方程式へ代入すると
172−2×122=289−288=1となる。以後も(3+22)3,(3+22)4,…によって解が得られる。
問題 4.2. 正整数nに対し、Tn=n(n+1)/2を三角数と呼ぶ。Tn=m2を満たす正整数の組(n,m)をすべて求めよ。また、nが小さい順に三つの組を書け。
解答.
n(n+1)/2=m2の両辺を8倍して1を加えると
(2n+1)2−8m2=1となる。したがって、x=2n+1、y=mとおけば、(x,y)はペル方程式x2−8y2=1の正整数解である。逆に、このペル方程式の正整数解ではx2≡1(mod8)なのでxは奇数である。またy≥1からx≥3である。よってn=(x−1)/2、m=yは正整数であり、Tn=m2を満たす。二つの方程式の正整数解は、この対応によって一対一に対応する。
8の連分数を求めると、a0=2であり、
(P1,Q1)=(2,4),a1=1,(P2,Q2)=(2,1),a2=4となる。したがって8=[2;1,4]で、周期長は2である。定理 2.2により、近似分数[2;1]=3/1から解(x1,y1)=(3,1)が得られる。定理 3.2により、この解は基本解であり、すべての正整数解は
xr+yr8=(3+8)r(r=1,2,3,…)で尽くされる。よって、求める組の全体は
(n,m)=(2xr−1,yr)(r=1,2,3,…)である。
最初の三つのべきは
3+8,(3+8)2=17+68,(3+8)3=99+358である。したがって、最初の三つの組は
(n,m)=(1,1), (8,6), (49,35)となる。解の合成式からxr+1=3xr+8yr>xrであるため、nもこの順に大きくなる。▨
5 例題
次のペル方程式 x² − D y²=1(D は平方数でない正整数)について答えよ。基本解とは、自明解 (±1, 0) 以外で x, y がともに最小の正の解のことである。
解法の型解は (x + y√D) のノルム 1 の元。基本解α1=x1 + y1√D のべきα1ⁿ を展開すると、すべての正の解が得られる
例 5.1. ペル方程式 x² − 2 y²=1 の基本解は (x, y)=(3, 2) である。(3 + 2√2)² を展開して、次の解を求めよ。
(3+22)2=x+y2⇒(x,y)=?答(x,y)=(17, 12)
注意 展開すると (3 + 2√2)²=3² + 2×2² + 2×3×2√2=17 + 12√2。ノルムN(α)=x² − 2y² は乗法的なのでN(α²)=N(α)²=1、つまりこれも解。検算: 17² − 2×12²=289 − 288=1。
例 5.2. ペル方程式 x² − 5 y²=1 の基本解(自明解以外で最小の正の解)を求めよ。
x2−5y2=1答(x,y)=(9, 4)
注意 検算: 9² − 5×4²=81 − 80=1。基本解は √5 の連分数展開の近似分数から得られる: √5=[2; 4, 4, 4, …](4 の繰り返し)。周期の長さは 1 で、基本解の大きさは周期の長さと密接に関係している(D=13 のように周期が長いと基本解は一気に大きくなる)。
例 5.3. ペル方程式 x² − 6 y²=1 の基本解は (x, y)=(5, 2) である。(5 + 2√6)² を展開して、次の解を求めよ。
(5+26)2=x+y6⇒(x,y)=?答(x,y)=(49, 20)
注意 展開すると (5 + 2√6)²=5² + 6×2² + 2×5×2√6=49 + 20√6。ノルムN(α)=x² − 6y² は乗法的なのでN(α²)=N(α)²=1、つまりこれも解。検算: 49² − 6×20²=2401 − 2400=1。
例 5.4. ペル方程式 x² − 5 y²=1 の基本解は (x, y)=(9, 4) である。(9 + 4√5)² を展開して、次の解を求めよ。
(9+45)2=x+y5⇒(x,y)=?答(x,y)=(161, 72)
注意 展開すると (9 + 4√5)²=9² + 5×4² + 2×9×4√5=161 + 72√5。ノルムN(α)=x² − 5y² は乗法的なのでN(α²)=N(α)²=1、つまりこれも解。検算: 161² − 5×72²=25921 − 25920=1。
例 5.5. ペル方程式 x² − 2 y²=1 について、基本解 (3, 2) と 2 番目の解 (17, 12) が分かっている。3 番目の解を求めよ。
(3+22)3=(3+22)(17+122)=x+y2⇒(x,y)=?答(x,y)=(99, 70)
注意 解の合成則 (x1, y1)·(x2, y2)=(x1x2 + D y1y2, x1y2 + x2y1) を使う: x=3×17 + 2×2×12=99、y=3×12 + 2×17=70。検算: 99² − 2×70²=9801 − 9800=1。基本解のべきα1ⁿ を取り続ければ解は無限に生成できる。
例 5.6. ペル方程式 x² − 8 y²=1 について、基本解 (3, 1) と 2 番目の解 (17, 6) が分かっている。3 番目の解を求めよ。
(3+8)3=(3+8)(17+68)=x+y8⇒(x,y)=?答(x,y)=(99, 35)
注意 解の合成則 (x1, y1)·(x2, y2)=(x1x2 + D y1y2, x1y2 + x2y1) を使う: x=3×17 + 8×1×6=99、y=3×6 + 1×17=35。検算: 99² − 8×35²=9801 − 9800=1。基本解のべきα1ⁿ を取り続ければ解は無限に生成できる。
例 5.7. ペル方程式 x² − 6 y²=1 について、基本解 (5, 2) と 2 番目の解 (49, 20) が分かっている。3 番目の解を求めよ。
(5+26)3=(5+26)(49+206)=x+y6⇒(x,y)=?答(x,y)=(485, 198)
注意 解の合成則 (x1, y1)·(x2, y2)=(x1x2 + D y1y2, x1y2 + x2y1) を使う: x=5×49 + 6×2×20=485、y=5×20 + 2×49=198。検算: 485² − 6×198²=235225 − 235224=1。基本解のべきα1ⁿ を取り続ければ解は無限に生成できる。
例 5.8. ペル方程式 x² − 5 y²=1 について、基本解 (9, 4) と 2 番目の解 (161, 72) が分かっている。3 番目の解を求めよ。
(9+45)3=(9+45)(161+725)=x+y5⇒(x,y)=?答(x,y)=(2889, 1292)
注意 解の合成則 (x1, y1)·(x2, y2)=(x1x2 + D y1y2, x1y2 + x2y1) を使う: x=9×161 + 5×4×72=2889、y=9×72 + 4×161=1292。検算: 2889² − 5×1292²=8346321 − 8346320=1。基本解のべきα1ⁿ を取り続ければ解は無限に生成できる。
例 5.9. ペル方程式 x² − 2 y²=1 の基本解(自明解以外で最小の正の解)を求めよ。
x2−2y2=1答(x,y)=(3, 2)
注意 検算: 3² − 2×2²=9 − 8=1。基本解は √2 の連分数展開の近似分数から得られる: √2=[1; 2, 2, 2, …](2 の繰り返し)。周期の長さは 1 で、基本解の大きさは周期の長さと密接に関係している(D=13 のように周期が長いと基本解は一気に大きくなる)。
例 5.10. ペル方程式 x² − 13 y²=1 の基本解(自明解以外で最小の正の解)を求めよ。
x2−13y2=1答(x,y)=(649, 180)
注意 検算: 649² − 13×180²=421201 − 421200=1。基本解は √13 の連分数展開の近似分数から得られる: √13=[3; 1, 1, 1, 1, 6, …](1, 1, 1, 1, 6 の繰り返し)。周期の長さは 5 で、基本解の大きさは周期の長さと密接に関係している(D=13 のように周期が長いと基本解は一気に大きくなる)。
6 演習
次のペル方程式 x² − D y²=1(D は平方数でない正整数)について答えよ。基本解とは、自明解 (±1, 0) 以外で x, y がともに最小の正の解のことである。
解法の型解は (x + y√D) のノルム 1 の元。基本解α1=x1 + y1√D のべきα1ⁿ を展開すると、すべての正の解が得られる
問題 6.1. ペル方程式 x² − 5 y²=1 の基本解(自明解以外で最小の正の解)を求めよ。
x2−5y2=1解答.
答(x,y)=(9, 4)
注意 検算: 9² − 5×4²=81 − 80=1。基本解は √5 の連分数展開の近似分数から得られる: √5=[2; 4, 4, 4, …](4 の繰り返し)。周期の長さは 1 で、基本解の大きさは周期の長さと密接に関係している(D=13 のように周期が長いと基本解は一気に大きくなる)。▨
問題 6.2. ペル方程式 x² − 3 y²=1 の基本解は (x, y)=(2, 1) である。(2 + √3)² を展開して、次の解を求めよ。
(2+3)2=x+y3⇒(x,y)=?解答.
答(x,y)=(7, 4)
注意 展開すると (2 + √3)²=2² + 3 + 2×2√3=7 + 4√3。ノルムN(α)=x² − 3y² は乗法的なのでN(α²)=N(α)²=1、つまりこれも解。検算: 7² − 3×4²=49 − 48=1。▨
問題 6.3. ペル方程式 x² − 6 y²=1 の基本解は (x, y)=(5, 2) である。(5 + 2√6)² を展開して、次の解を求めよ。
(5+26)2=x+y6⇒(x,y)=?解答.
答(x,y)=(49, 20)
注意 展開すると (5 + 2√6)²=5² + 6×2² + 2×5×2√6=49 + 20√6。ノルムN(α)=x² − 6y² は乗法的なのでN(α²)=N(α)²=1、つまりこれも解。検算: 49² − 6×20²=2401 − 2400=1。▨
問題 6.4. ペル方程式 x² − 7 y²=1 の基本解は (x, y)=(8, 3) である。(8 + 3√7)² を展開して、次の解を求めよ。
(8+37)2=x+y7⇒(x,y)=?解答.
答(x,y)=(127, 48)
注意 展開すると (8 + 3√7)²=8² + 7×3² + 2×8×3√7=127 + 48√7。ノルムN(α)=x² − 7y² は乗法的なのでN(α²)=N(α)²=1、つまりこれも解。検算: 127² − 7×48²=16129 − 16128=1。▨
問題 6.5. ペル方程式 x² − 6 y²=1 について、基本解 (5, 2) と 2 番目の解 (49, 20) が分かっている。3 番目の解を求めよ。
(5+26)3=(5+26)(49+206)=x+y6⇒(x,y)=?解答.
答(x,y)=(485, 198)
注意 解の合成則 (x1, y1)·(x2, y2)=(x1x2 + D y1y2, x1y2 + x2y1) を使う: x=5×49 + 6×2×20=485、y=5×20 + 2×49=198。検算: 485² − 6×198²=235225 − 235224=1。基本解のべきα1ⁿ を取り続ければ解は無限に生成できる。▨
問題 6.6. ペル方程式 x² − 8 y²=1 の基本解は (x, y)=(3, 1) である。(3 + √8)² を展開して、次の解を求めよ。
(3+8)2=x+y8⇒(x,y)=?解答.
答(x,y)=(17, 6)
注意 展開すると (3 + √8)²=3² + 8 + 2×3√8=17 + 6√8。ノルムN(α)=x² − 8y² は乗法的なのでN(α²)=N(α)²=1、つまりこれも解。検算: 17² − 8×6²=289 − 288=1。▨
問題 6.7. ペル方程式 x² − 6 y²=1 の基本解(自明解以外で最小の正の解)を求めよ。
x2−6y2=1解答.
答(x,y)=(5, 2)
注意 検算: 5² − 6×2²=25 − 24=1。基本解は √6 の連分数展開の近似分数から得られる: √6=[2; 2, 4, 2, 4, …](2, 4 の繰り返し)。周期の長さは 2 で、基本解の大きさは周期の長さと密接に関係している(D=13 のように周期が長いと基本解は一気に大きくなる)。▨
問題 6.8. ペル方程式 x² − 8 y²=1 について、基本解 (3, 1) と 2 番目の解 (17, 6) が分かっている。3 番目の解を求めよ。
(3+8)3=(3+8)(17+68)=x+y8⇒(x,y)=?解答.
答(x,y)=(99, 35)
注意 解の合成則 (x1, y1)·(x2, y2)=(x1x2 + D y1y2, x1y2 + x2y1) を使う: x=3×17 + 8×1×6=99、y=3×6 + 1×17=35。検算: 99² − 8×35²=9801 − 9800=1。基本解のべきα1ⁿ を取り続ければ解は無限に生成できる。▨
問題 6.9. ペル方程式 x² − 2 y²=1 の基本解(自明解以外で最小の正の解)を求めよ。
x2−2y2=1解答.
答(x,y)=(3, 2)
注意 検算: 3² − 2×2²=9 − 8=1。基本解は √2 の連分数展開の近似分数から得られる: √2=[1; 2, 2, 2, …](2 の繰り返し)。周期の長さは 1 で、基本解の大きさは周期の長さと密接に関係している(D=13 のように周期が長いと基本解は一気に大きくなる)。▨
問題 6.10. ペル方程式 x² − 13 y²=1 の基本解(自明解以外で最小の正の解)を求めよ。
x2−13y2=1解答.
答(x,y)=(649, 180)
注意 検算: 649² − 13×180²=421201 − 421200=1。基本解は √13 の連分数展開の近似分数から得られる: √13=[3; 1, 1, 1, 1, 6, …](1, 1, 1, 1, 6 の繰り返し)。周期の長さは 5 で、基本解の大きさは周期の長さと密接に関係している(D=13 のように周期が長いと基本解は一気に大きくなる)。▨
問題 6.11. ペル方程式 x² − 8 y²=1 の基本解(自明解以外で最小の正の解)を求めよ。
x2−8y2=1解答.
答(x,y)=(3, 1)
注意 検算: 3² − 8×1²=9 − 8=1。基本解は √8 の連分数展開の近似分数から得られる: √8=[2; 1, 4, 1, 4, …](1, 4 の繰り返し)。周期の長さは 2 で、基本解の大きさは周期の長さと密接に関係している(D=13 のように周期が長いと基本解は一気に大きくなる)。▨
問題 6.12. ペル方程式 x² − 3 y²=1 の基本解(自明解以外で最小の正の解)を求めよ。
x2−3y2=1解答.
答(x,y)=(2, 1)
注意 検算: 2² − 3×1²=4 − 3=1。基本解は √3 の連分数展開の近似分数から得られる: √3=[1; 1, 2, 1, 2, …](1, 2 の繰り返し)。周期の長さは 2 で、基本解の大きさは周期の長さと密接に関係している(D=13 のように周期が長いと基本解は一気に大きくなる)。▨
問題 6.13. ペル方程式 x² − 7 y²=1 について、基本解 (8, 3) と 2 番目の解 (127, 48) が分かっている。3 番目の解を求めよ。
(8+37)3=(8+37)(127+487)=x+y7⇒(x,y)=?解答.
答(x,y)=(2024, 765)
注意 解の合成則 (x1, y1)·(x2, y2)=(x1x2 + D y1y2, x1y2 + x2y1) を使う: x=8×127 + 7×3×48=2024、y=8×48 + 3×127=765。検算: 2024² − 7×765²=4096576 − 4096575=1。基本解のべきα1ⁿ を取り続ければ解は無限に生成できる。▨
問題 6.14. ペル方程式 x² − 2 y²=1 について、基本解 (3, 2) と 2 番目の解 (17, 12) が分かっている。3 番目の解を求めよ。
(3+22)3=(3+22)(17+122)=x+y2⇒(x,y)=?解答.
答(x,y)=(99, 70)
注意 解の合成則 (x1, y1)·(x2, y2)=(x1x2 + D y1y2, x1y2 + x2y1) を使う: x=3×17 + 2×2×12=99、y=3×12 + 2×17=70。検算: 99² − 2×70²=9801 − 9800=1。基本解のべきα1ⁿ を取り続ければ解は無限に生成できる。▨
問題 6.15. ペル方程式 x² − 5 y²=1 の基本解は (x, y)=(9, 4) である。(9 + 4√5)² を展開して、次の解を求めよ。
(9+45)2=x+y5⇒(x,y)=?解答.
答(x,y)=(161, 72)
注意 展開すると (9 + 4√5)²=9² + 5×4² + 2×9×4√5=161 + 72√5。ノルムN(α)=x² − 5y² は乗法的なのでN(α²)=N(α)²=1、つまりこれも解。検算: 161² − 5×72²=25921 − 25920=1。▨