■連分数展開の問題

 連分数展開が有限で終わることと有理数であることは同値です.そこで,2次方程式の解となる√nの連分数展開を求めると,たとえば

  √2=[1:2,2,2,2,・・・]

  √3=[1:1,2,1,2,1,2,1,2,・・・]

  √7=[2:1,1,1,4,1,1,1,4,・・・]

のように循環型の単純連分数に展開されることが知られています.一般に,2次の無理数(整数係数の2次方程式の解)は周期的な連分数展開をもちます(ラグランジュの定理).

 平方根を無限連分数に表す手順はわかりやすく,たとえば,1<√2<2であるから

  √2=1+(√2−1)

    =1+1/(√2+1)    2<√2+1<3

    =1+1/{2+(√2−1)}

    =1+1/{2+1/(√2+1)}

    =1+1/{2+1/(2+(√2−1)}

    =1+1/{2+1/(2+1/(√2+1)}

    =1+1/{2+1/{2+1/{2+1/{2+・・・

の手順を何度も繰り返すことにより,

  √2=[1:2,2,2,2,・・・]

ができあがります.また,黄金比φ=(1+√5)/2は,

  φ=[1:1,1,1,,1,・・・]

で表されます.黄金比φ=(1+√5)/2が,無限連分数

  φ=[1:1,1,1,,1,・・・]

や無限の入れ子の根号

  φ=√(1+√(1+√(1+√(1+・・・

で3通りにも表されるという事実は魔法のようにさえ思えます.

 ここでは,連分数展開を用いて数の集合を定義してみますが,たとえば,正の実数が無限連分数展開され,そのすべての部分商が1または2であるような実数の集合のハウスドルフ次元は0.531280506・・・であることが計算されています.

 3次以上の方程式の解,たとえば3√2の連分数展開を求めると,

  3√2=[1:3,1,5,1,1,4,1,1,8,1,14,1,10,2,1,4,・・・]

の一般項は求めることができません.この展開に現れる整数に最大値があることも示すこともできないのです.

 有理数は有限連分数,無理数で代数的数の場合は無限循環連分数,超越数は無限非循環連分数になります.たとえば,超越数eの連分数展開は,

  e=[2;1,2,1,1,4,1,1,6,1,1,8,1,1,10,1,1,12,1,1,14,1,1,16,・・・]

と書け,数字の出方が自然数順になっていることがわかります.すなわち,

  e=[2;1,2,1,1,4,1,1,6,1,・・・,1,2n,1,・・・]

 πの連分数展開

  π=[3;7,15,1,292,1,1,1,2,1,3,1,14,2,1,1,2,2,2,2,1,84,2,1,1,15,3,13,1,4,2,6,6,99,1,2,2,6,3,5,1,1,6,・・・]

にはなんの規則性も見あたらないようにみえます.もちろん,一般項は見つかっていません.10進数表現しても

  e=2.718281827459045・・・

 π=3.141592653589793・・・

eには何かパターンがありそうに見えますが,πの数の並び方には何のパターンもありません.しかし,単純連分数(分子がすべて1)に限らなければ,

  π/4=1/{1+1^2/{2+3^2/{2+5^2/{2+7^2/{2+9^2/{2+・・・}

分子には奇数の平方が並んでいるというパターンを見つけることができます.

===================================

 単純循環連分数

  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)

 同様に,2項が循環する連分数は

  L=[a:b,c,b,c,・・・]

  L=ab−bc/2+((bc)^2/4+bc)^(1/2)

 ところで,数を連分数で表示すると数字1が大量に出現することに気づきます.そこで,連分数の部分商の分布について考えてみます.

  [参]Havil著,新妻弘監訳「オイラーの定数ガンマ」共立出版

によると,整数部を除いた[0:a1,a2,a3,・・・,an]がxより小さい小数となる確率は

 P([0:a1,a2,a3,・・・,an]<x)=log2(1+x)+εn

で与えられますが,1928年にクズミンはほとんどすべての連分数に対して,

  εn=O(q^√n)  0<q<1

1929年にレヴィは

  εn=O(q^n)  q=0.7

であることを示しました.どちらも誤差項εnは漸近的に0になることを示しています.

 連分数の部分商の確率密度関数は

  P(an=k)=P(k<εn<k+1)=P(εn<k+1)−P(εn<k)

→log2(1+1/k)−log2(1+1/(k+1))=log2(1+1/k(k+2))

 したがって,十分大きなnに対する部分商の起こる確率Pは

k        1  2  3  4  5  6  7  8 9+

P(an =k)  .41 .17  .09  .06  .04 .03 .02 .02 .16

となることがわかります.

===================================

 次にanの平均値を求めてみます. ヒンチンは一般の連分数

  [a0:a1,a2,a3,・・・,an,・・・]

の大多数についてあてはまる法則を発見しています.

 ヒンチンの定理とは,幾何平均(a1a2・・・an)^1/nの値がn→∞のとき,ある無限乗積から定まる定数

  (a1a2・・・an)^1/n→Π(1+1/k(k+2))^logk/log2=2.685452001・・・

に収束するというものです.κ=2.68545・・・はヒンチンの定数として知られています.

 ただし,分母に明確なパターンのある代数的数やeをはじめとするいくつかの超越数は例外になります.

  (eの場合,(a1a2・・・an)^1/n→0.6259・・・)

 算術平均は発散するのに対し幾何平均は収束するというわけですが,ほとんどすべての連分数の場合,調和平均も収束し,その極限値は

  n/(1/a1+1/a2+・・・+1/an)→1.74540568・・・

===================================