連分数とは,
a1/(b1+a2/(b2+a3/(b3+a4/(b4+a5/b5+・・・)
のような分数を続けた式で,実用上は最初にa0+をつけた形が使われます.整数論で使われる連分数は普通,ak=1,bkが正の整数である標準連分数です.
(その1)では超黄金比を連分数を使って説明しましたが,今回のコラムではそれを補足してみたいと思います.
===================================
[Q]大長方形から正方形をa個切り取った後に残る中長方形から,正方形をb個切り取った後に残る小長方形が中長方形と相似になるのは?
[A]単純循環連分数
L=[a:b,b,b,b,・・・]
で表される数Lを求めてみることにしましょう.
L−a=R=[0:b,b,b,b,・・・]=1/(b+R)
R^2+bR−1=0 → R=(−b+(b^2+4)^(1/2))/2
L=a+R=a−b/2+(b^2/4+1)^(1/2)
a=1,b=1 → L=(1+√5)/2=[1;1,1,1,1,1,・・・]
a=1,b=2 → L=√2=[1;2,2,2,2,2,・・・]
a=2,b=4 → L=√5=[2;4,4,4,・・・]
[Q]大長方形から正方形をa個切り取った後に残る中長方形から,正方形をb個切り取り,さらに正方形をc個切り取った後に残る小長方形が中長方形と相似になるのは?
[A]同様に,2項が循環する連分数は
L=[a:b,c,b,c,・・・]
L−a=R=[0:b,c,b,c,・・・]=1/(b+1/(c+R))
bR^2+bcR−c=0 → R=(−c+(c^2+4c/b)^(1/2))/2
L=a−c/2+(c^2/4+c/b)^(1/2)
a=1,b=1,c=2 → L=√3=[1;1,2,1,2,1,2,・・・]
a=2,b=2,c=4 → L=√6=[2;2,4,2,4,2,・・・]
−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
連分数展開によって
(1+√5)/2=[1;1,1,1,1,1,・・・]
√2=[1;2,2,2,2,2,・・・]
のように,1や2が無限に繰り返されるという規則性を見ることができます.
√3=[1;1,2,1,2,1,2,・・・]
では交互に1,2が現れる循環連分数となります.以下,
√5=[2;4,4,4,・・・]
√6=[2;2,4,2,4,2,・・・]
√7=[2;1,1,1,4,1,1,1,4,・・・]
一般に,√mの連分数展開は循環連分数となり周期性が証明されます.これは既約分数の小数展開が循環小数になることと対比するとおもしろい事実です.
その際,
√m=[q0;q1,q2,・・・,qn-1,2q0,・・・]
という周期nの連分数展開が得られます.
√2=[1;2,・・・]
√3=[1;1,2,・・・]
√5=[2;4,・・・]
√6=[2;2,4,・・・]
√7=[2;1,1,1,4,・・・]
すなわち,どの循環節もqn=2q0=[2√m]で終わっています.
たとえば,√199の展開
√199=[14;9,2,1,2,2,5,4,1,1,13,1,1,4,5,2,2,1,2,9,28,・・・]
をみると,14で始まり28で終わるというのもこの理由によります.
このように,標準無限連分数のうち,部分分母列のあるところから先が巡回的になる循環連分数は2次の無理数(整数係数の2次方程式の解として表される数)に収束します.この性質により,整数項の標準連分数はいわゆるペル方程式:x^2−my^2=d(多くは±1,±4)の解法など整数論の分野で活用されます.
また,√199の循環節の最後の28を除くと13を中心として対称になっていることにも気付かされます.
√43=[6;1,1,3,1,5,1,3,1,1,12,・・・]
√54=[7;2,1,6,1,2,14,・・・]
√76=[8;1,2,1,1,5,4,5,1,1,2,1,16,・・・]
√94=[9;1,2,3,1,1,5,1,8,1,5,1,1,3,2,1,18,・・・]
√1000=[31;1,1,1,1,1,6,2,2,15,2,2,6,1,1,1,1,1,62,・・・]
循環部の最後の項を除いた部分は回文(前から読んでも後から読んでも同じ)になっているという事実も,199のみならず,2次の無理数√mに共通していえる性質です.
√m=[q0;q1,q2,・・,q2,q1,2q0,・・・]
なお,2次の無理数には循環連分数が対応しますが,連分数による実数の最良近似は解を下方と上方から近似していく方法であって,ユークリッドの互除法に直結しています.
−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
以上のことから,最も素朴な循環連分数は
√m=[q0;2q0,2q0,2q0,・・・]
で表されるものと考えられます.
このとき,
P=2q0^2+1,Q=2q0
より,mは
(2q0^2+1)^2−m・4q0^2=±1
を満たす整数となるのですが,結局,このようなmは
m=q0^2+1=2,5,10,・・・
となることが導き出されます.
√2=[1;2,2,2,・・・]
√5=[2;4,4,4,・・・]
===================================