■いわゆるペル方程式(その20)
b・243-a・256=1なるディオファントス解を求めたい。
(256,243)=1となることを。ユークリッドの互除法で示したい。
256=1・243+13
243=18・13+9
13=1・9+4
9=2・4+1
4=4・1
===================================
互除法に現れる計算を、行列を使って書き直すと
[256]=[1,1][243]
[243] {1,0][ 13]
[256]=[1,1][18,1][13]
[243] {1,0][ 1,0][ 9]
[256]=[1,1][18,1][1,1][9]
[243] {1,0][ 1,0][1,0][4]
[256]=[1,1][18,1][1,1][2,1][4]
[243] {1,0][ 1,0][1,0][1,0][1]
[256]=[1,1][18,1][1,1][2,1][4,1][1]
[243] {1,0][ 1,0][1,0][1,0][1,0][0]
===================================
有限連分数77/23=[3:2,1,7]の近似分数の計算を、行列を使って書き直すと
[p1,p0]=[3,1][2,1]=[7,3]
[q1,q0] [1,0][1,0] [2,1], 3/1→7/2
[p2,p1]=[3,1][2,1][1,1]=[7,3][1,1]=[10,7]
[q2,q1] [1,0][1,0][1,0] [2,1][1,0] [ 3,2], →10/3
[p3,p2]=[3,1][2,1][1,1][7,1]=[10,7][7,1]=[77,10]
[q3,q2] [1,0][1,0][1,0][1,0] [ 3,2][1,0] [23, 3] →77/23
===================================
無理数αと連分数の関係は
{a0:a1,a2,・・・,ak,α]=[a0,1][a1,1]・・・[ak,1]α
[ 1,0][ 1,0] [ 1,0]
で与えられる。たとえば、
√3=[1:1,2,1,2,・・・]は,α=[1,2]とおくと
√3=[1,1]α
[1,0]
===================================
α=[1,3,5]は簡約2次無理数である。
α=[1,1][3,1][5,1]α
[1,0][1,0][1,0]
α=[21,4]α=(21α+4)/(16α+3)
[16,3]
16α^2-18α-4=0
α>1より、α=(9+√145)/16
===================================