リンド数学パピルス
エジプト式分数 (エジプトしきぶんすう、単にエジプト分数 とも、英 : Egyptian fraction )とは、いくつかの異なる単位分数 (分子が 1 の分数)の和、あるいは分数 をそのように表す方式を意味する。例えば、通常 5 / 6 で表す分数を 1 / 2 + 1 / 3 などと表す。任意の正 の有理数 はこの形式で表すことができるが、表し方は一意ではない。この形式で分数を扱う方法は、古くは古代エジプト のリンド・パピルス に見られ、ヨーロッパ では中世 まで広く用いられた。現代でも数論 の分野において、エジプト式分数に端を発する数学上の未解決問題 が多く残されている。
単位分数展開
以下、特に断らない限り、単に「分数」といった場合、正の真分数、すなわち 0 より大きく 1 より小さな分数のみを考えているものとする。
例えば 2 / 5 は単位分数の和として 1 / 5 + 1 / 5 と表せるが、エジプト式分数では同じ単位分数を繰り返し用いることはせず、2 / 5 = 1 / 3 + 1 / 15 のように表す。いかなる分数に対してもこのような単位分数展開が必ず存在することは自明ではないが、後述するように今日ではあらゆる分数が無数に多くの単位分数展開を持つことが証明されている(#強欲算法 の節参照)。さらに例を挙げると、3 / 7 = 1 / 4 + 1 / 7 + 1 / 28 = 1 / 6 + 1 / 7 + 1 / 14 + 1 / 21 であって、前者の展開は項数が最小であり、後者の展開は最大分母の値が最小である[ 1] 。このように、どのような単位分数展開が最も「単純」であるか、は明らかではない。
古代エジプト
ホルスの目
エジプト中王国 では、ホルスの目 を用いたそれ以前の不完全な分数体系(1 / 2k (k = 1, 2, …, 6) の和で表す)に替わって、エジプト式分数による方法が発達した。エジプト式分数が見られる古い文献としては、エジプト数学革巻き 、モスクワ・パピルス 、レイズナー・パピルス (英語版 ) 、カフン・パピルス (英語版 ) 、アクミム木刻版 (英語版 ) がある。特に有名なリンド・パピルスは、紀元前1650年頃に書かれたものであり、5 以上 101 以下の奇数 n に対して 2 / n を単位分数の和で表している(#リンド・パピルスの展開一覧 の節参照)。
古代エジプト人が、いちいちこのように単位分数の和で表した理由については、よく分かっていない。ただ、リンド・パピルスにはパンを分け合う問題がいくつもあって、実際にパンを分け合うにはエジプト式の表示が理に適っている場合がある。例えば、リンド・パピルスの問題3は、6斤のパンを10人で分け合うとき、1人分は 1 / 2 + 1 / 10 であることを答とする。6斤のパンをそれぞれ5等分するよりも、5斤を1斤づつ2等分して1片ずつ取り、残りの1斤を10等分する方が簡単である[ 2] 。一方では、合理的とは思えない表示を選ぶ場合もある。リンド・パピルスの問題4は、7斤のパンを10人で分け合う問題であるが、1 / 2 + 1 / 5 ではなく、2 / 3 + 1 / 30 を答としている[ 3] 。2 / 3 は単位分数ではないから、この表示は狭い意味でエジプト式ではないが、古代エジプト人にとって 2 / 3 は特別な数であったらしい。2 / 3 = 1 / 2 + 1 / 6 であることを知っていたにもかかわらず、好んでこの数を用いている。
リンド・パピルスにおける 2 / n の表を参照すれば、分母が 100 以下の奇数である多くの分数が、機械的に単位分数の和で表せる。例えば、表より 2 / 21 = 1 / 14 + 1 / 42 であるから、
5 / 21 = 1 / 21 + (1 / 14 + 1 / 42 ) + (1 / 14 + 1 / 42 ) = 1 / 21 + 1 / 7 + 1 / 21 = 1 / 7 + 1 / 14 + 1 / 42
と計算できる[ 4] 。リンド・パピルスにおいて、2 / n に特に注意が払われているのは、古代エジプトの乗法 (英語版 ) アルゴリズムが2倍を基礎においているためであろう、とも考えられている[ 5] 。
表記
古代エジプト人たちは、2 / 3 を唯一の例外として、単位分数のみを表記した。単位分数 1 / n を表すために、神官文字 では点を、神聖文字 では
を n を表す記号の上に置いた。例えば
=
1
3
{\displaystyle ={\frac {1}{3}}}
=
1
10
{\displaystyle ={\frac {1}{10}}}
といった具合である。1 / 2 と 2 / 3 のみ、特別なグリフ
=
1
2
{\displaystyle ={\frac {1}{2}}}
=
2
3
{\displaystyle ={\frac {2}{3}}}
を持つ。2 / 3 のグリフは、正確には右の縦線が若干長い。長い方が 1 を、短い方が 1 / 2 を表し、全体としてはその和 3 / 2 の逆数を意味している[ 2] 。
計算方法
現代の数学史家は、リンド・パピルスなどの古文書を調べ、古代エジプト人のエジプト式分数による計算方法がどのようなものであったかを研究した。特に、リンド・パピルスに書かれた 2 / n の表現がどのように得られたのかに注目し、様々な説を立てている。古代エジプト人が、分数を単位分数の和に表す系統的な方法を知っていたかどうかは不明であるが、少なくとも単一の方法のみを用いたのではなさそうである。恒等式 2 / 2m + 1 = 1 / m + 1 + 1 / (m + 1)(2m + 1) を用いれば、単一の方法で2つの単位分数の和に表せるにもかかわらず、分母が大きくなるのを嫌ってか、リンド・パピルスでは3項あるいは4項の和に表しているものもある。数学史家たちの分析によれば、分母が素数 の場合と合成数 の場合で、リンド・パピルスの著者は異なる方法を用いており、それぞれの場合においても複数の方法を用いている。
分母が奇素数の場合(1)
小さな奇素数 p = 2m + 1 (3, 5, 7, 11, 23) に対しては、恒等式 2 / 2m + 1 = 1 / m + 1 + 1 / (m + 1)(2m + 1) が用いられている。この方法は奇素数に限らず、任意の奇数に対して使用できる。
分母が奇素数の場合(2)
大きめの奇素数 p (13, 17, 19, 29, …) に対しては、恒等式 2 / p = 1 / A + 2A − p / Ap が用いられている。ここで、A は p / 2 < A < p を満たし、約数 を多く持つ数が選ばれる。2A − p / Ap について、分子が A のいくつかの約数の和に表すことができれば、約分して単位分数の和を得る。例えば、p = 37 に対して A = 24 とすると、2A − p = 11 = 3 + 8 で 3 と 8 は 24 の約数であるから、リンド・パピルスの展開 2 / 37 = 1 / 24 + 1 / 111 + 1 / 296 を得る。A を取り替えたり、約数の和に分解する方法を変えたりすると別の展開を得る。
分母が半素数の場合(1)
分母が2つの奇素数の積として pq であるとき、a = p + 1/ 2 として恒等式 2 / pq = 1 / aq + 1 / apq を用いることができる。例えば、p = 3, q = 7 のとき、a = 2 より 2 / 21 = 1 / 14 + 1 / 42 を得る。この方法で、リンド・パピルスの単位分数展開のうち、分母が半素数 であるものの多くは説明が付く。
分母が半素数の場合(2)
分母が半素数の場合、r = p + q / 2 として恒等式 2 / pq = 1 / pr + 1 / qr を用いることもできる[ 6] 。例えば、p = 5, q = 7 とすると、リンド・パピルスの表示 2 / 35 = 1 / 30 + 1 / 42 を得る。2 / 91 についても同様である。
分母がその他の合成数の場合
その他の合成数 n については、n の約数 m に対する 1 / m の単位分数展開から得られる。例えば、2 / 19 = 1 / 12 + 1 / 76 + 1 / 114 を 5 で割ることにより、2 / 95 = 1 / 60 + 1 / 380 + 1 / 570 を得る。実際は、1 / 380 + 1 / 570 = 1 / 228 であるから、より簡単な展開を得るが、リンド・パピルスでは簡約化されていないものが記されている。3つ以上の素数の積、27, 45, 63, 75, 81, 99 に対してもこの方法で説明が付く。
分母が101の場合
リンド・パピルスの最後の単位分数展開 2 / 101 = 1 / 101 + 1 / 202 + 1 / 303 + 1 / 606 は、以上のどれにも当てはまらないが、恒等式 2 / p = 1 / p + 1 / 2p + 1 / 3p + 1 / 6p に p = 101 を代入して得られる。これと同等の等式は『エジプト数学羊皮紙巻子本』でも用いられている。
リンド・パピルスの展開一覧
リンド・パピルスの最初に記された単位分数展開の一覧[ 7] を下記の表に記す。2 / 3 は別格として特別の注意が払われている。単位分数展開は一意ではないが、リンド・パピルスでは、1つの分数に対して1つの展開だけが記されており、それは必ずしも最も単純な展開ではない。例えば、2 / 13 = 1 / 7 + 1 / 91 であるが、なぜかこれよりも項数が多く、分母も大きなものが記されている。
背景が水色のセルはリンド・パピルスに記されている展開方法を示す。
リンド・パピルスに記された 2 / n の単位分数展開一覧
分母
種類
奇数
奇素数
半素数
半素数
合成数
3
素数
(2 / 3 = 1 / 2 + 1 / 6 )
5
素数
2 / 5 = 1 / 3 + 1 / 15
7
素数
2 / 7 = 1 / 4 + 1 / 28
9
半素数
2 / 9 = 1 / 5 + 1 / 45
2 / 9 = 1 / 6 + 1 / 18 (a =2, p =3, q =3)
2 / 9 = 1 / 6 + 1 / 18 (2 / 3 = 1 / 2 + 1 / 6 を 3 で割る。)
11
素数
2 / 11 = 1 / 6 + 1 / 66
13
素数
2 / 13 = 1 / 7 + 1 / 91
2 / 13 = 1 / 8 + 1 / 52 + 1 / 104 (A =8, p =13, 2A -p =3=2+1)
15
半素数
2 / 15 = 1 / 8 + 1 / 120
2 / 15 = 1 / 10 + 1 / 30 (a =2, p =3, q =5)
2 / 15 = 1 / 12 + 1 / 20 (r =4, p =3, q =5)
2 / 15 = 1 / 10 + 1 / 30 (2 / 3 = 1 / 2 + 1 / 6 を 5 で割る。)
17
素数
2 / 17 = 1 / 9 + 1 / 153
2 / 17 = 1 / 12 + 1 / 51 + 1 / 68 (A =12, p =17, 2A -p =7=4+3)
19
素数
2 / 19 = 1 / 10 + 1 / 190
2 / 19 = 1 / 12 + 1 / 76 + 1 / 114 (A =12, p =19, 2A -p =5=3+2)
21
半素数
2 / 21 = 1 / 11 + 1 / 231
2 / 21 = 1 / 14 + 1 / 42 (a =2, p =3, q =7)
2 / 21 = 1 / 15 + 1 / 35 (r =5, p =3, q =7)
2 / 21 = 1 / 14 + 1 / 42 (2 / 3 = 1 / 2 + 1 / 6 を 7 で割る。)
23
素数
2 / 23 = 1 / 12 + 1 / 276
25
半素数
2 / 25 = 1 / 13 + 1 / 325
2 / 25 = 1 / 15 + 1 / 75 (a =3, p =5, q =5)
2 / 25 = 1 / 15 + 1 / 75 (2 / 5 = 1 / 3 + 1 / 15 を 5 で割る。)
27
合成数
2 / 27 = 1 / 14 + 1 / 378
2 / 27 = 1 / 18 + 1 / 54 (2 / 3 = 1 / 2 + 1 / 6 を 9 で割る。)
29
素数
2 / 29 = 1 / 15 + 1 / 435
2 / 29 = 1 / 24 + 1 / 58 + 1 / 174 + 1 / 232 (A =24, p =29, 2A -p =19=12+4+3)
31
素数
2 / 31 = 1 / 16 +1 / 496
2 / 31 = 1 / 20 +1 / 124 + 1 / 155 (A =20, p =31, 2A -p =9=5+4)
33
半素数
2 / 33 = 1 / 17 + 1 / 561
2 / 33 = 1 / 22 + 1 / 66 (a =2, p =3, q =11)
2 / 33 = 1 / 21 + 1 / 77 (r =7, p =3, q =11)
2 / 33 = 1 / 22 + 1 / 66 (2 / 3 = 1 / 2 + 1 / 6 を 11 で割る。)
35
半素数
2 / 35 = 1 / 18 + 1 / 630
2 / 35 = 1 / 21 + 1 / 105 (a =3, p =5, q =7)
2 / 35 = 1 / 30 + 1 / 42 (r =6, p =5, q =7)
2 / 35 = 1 / 21 + 1 / 105 (2 / 5 = 1 / 3 + 1 / 15 を 7 で割る。)
37
素数
2 / 37 = 1 / 19 + 1 / 703
2 / 37 = 1 / 24 + 1 / 111 + 1 / 296 (A =24, p =37, 2A -p =11=8+3)
39
半素数
2 / 39 = 1 / 20 + 1 / 780
2 / 39 = 1 / 26 + 1 / 78 (a =2, p =3, q =13)
2 / 39 = 1 / 24 + 1 / 104 (r =8, p =3, q =13)
2 / 39 = 1 / 26 + 1 / 78 (2 / 3 = 1 / 2 + 1 / 6 を 13 で割る。)
41
素数
2 / 41 = 1 / 21 + 1 / 861
2 / 41 = 1 / 24 + 1 / 246 + 1 / 328 (A =24, p =41, 2A -p =7=4+3)
43
素数
2 / 43 = 1 / 22 + 1 / 946
2 / 43 = 1 / 42 + 1 / 86 + 1 / 129 + 1 / 301 (A =42, p =43, 2A -p =41=21+14+6)
45
合成数
2 / 45 = 1 / 23 + 1 / 1035
2 / 45 = 1 / 30 + 1 / 90 (2 / 3 = 1 / 2 + 1 / 6 を 15 で割る。)
47
素数
2 / 47 = 1 / 24 + 1 / 1128
2 / 47 = 1 / 30 + 1 / 141 + 1 / 470 (A =30, p =47, 2A -p =13=10+3)
49
半素数
2 / 49 = 1 / 25 + 1 / 1225
2 / 49 = 1 / 28 + 1 / 196 (a =4, p =7, q =7)
2 / 49 = 1 / 28 + 1 / 196 (2 / 7 = 1 / 4 + 1 / 28 を 7 で割る。)
51
半素数
2 / 51 = 1 / 26 + 1 / 1326
2 / 51 = 1 / 34 + 1 / 102 (a =2, p =3, q =17)
2 / 51 = 1 / 30 + 1 / 170 (r =10, p =3, q =17)
2 / 51 = 1 / 34 + 1 / 102 (2 / 3 = 1 / 2 + 1 / 6 を 17 で割る。)
53
素数
2 / 53 = 1 / 27 + 1 / 1431
2 / 53 = 1 / 30 + 1 / 318 + 1 / 795 (A =30, p =53, 2A -p =7=5+2)
55
半素数
2 / 55 = 1 / 28 + 1 / 1540
2 / 55 = 1 / 33 + 1 / 165 (a =3, p =5, q =11)
2 / 55 = 1 / 40 + 1 / 88 (r =8, p =5, q =11)
2 / 55 = 1 / 30 + 1 / 330 (2 / 11 = 1 / 6 + 1 / 66 を 5 で割る。)
57
半素数
2 / 57 = 1 / 29 + 1 / 1653
2 / 57 = 1 / 38 + 1 / 114 (a =2, p =3, q =19)
2 / 57 = 1 / 33 + 1 / 209 (r =11, p =3, q =19)
2 / 57 = 1 / 38 + 1 / 114 (2 / 3 = 1 / 2 + 1 / 6 を 19 で割る。)
59
素数
2 / 59 = 1 / 30 + 1 / 1770
2 / 59 = 1 / 36 + 1 / 236 + 1 / 531 (A =36, p =59, 2A -p =13=9+4)
61
素数
2 / 61 = 1 / 31 + 1 / 1891
2 / 61 = 1 / 40 + 1 / 244 + 1 / 488 + 1 / 610 (A =40, p =61, 2A -p =19=10+5+4)
63
合成数
2 / 63 = 1 / 32 + 1 / 2016
2 / 63 = 1 / 42 + 1 / 126 (2 / 3 = 1 / 2 + 1 / 6 を 21 で割る。)
65
半素数
2 / 65 = 1 / 33 + 1 / 2145
2 / 65 = 1 / 39 + 1 / 195 (a =3, p =5, q =13)
2 / 65 = 1 / 45 + 1 / 117 (r =9, p =5, q =13)
2 / 65 = 1 / 39 + 1 / 195 (2 / 5 = 1 / 3 + 1 / 15 を 13 で割る。)
67
素数
2 / 67 = 1 / 34 + 1 / 2278
2 / 67 = 1 / 40 + 1 / 335 + 1 / 536 (A =40, p =67, 2A -p =13=8+5)
69
半素数
2 / 69 = 1 / 35 + 1 / 2415
2 / 69 = 1 / 46 + 1 / 138 (a =2, p =3, q =23)
2 / 69 = 1 / 39 + 1 / 299 (r =13, p =3, q =23)
2 / 69 = 1 / 46 + 1 / 138 (2 / 3 = 1 / 2 + 1 / 6 を 23 で割る。)
71
素数
2 / 71 = 1 / 36 + 1 / 2556
2 / 71 = 1 / 40 + 1 / 568 + 1 / 710 (A =40, p =71, 2A -p =9=5+4)
73
素数
2 / 73 = 1 / 37 + 1 / 2701
2 / 73 = 1 / 60 + 1 / 219 + 1 / 292 + 1 / 365 (A =60, p =73, 2A -p =47=20+15+12)
75
合成数
2 / 75 = 1 / 38 + 1 / 2850
2 / 75 = 1 / 50 + 1 / 150 (2 / 3 = 1 / 2 + 1 / 6 を 25 で割る。)
77
半素数
2 / 77 = 1 / 39 + 1 / 3003
2 / 77 = 1 / 44 + 1 / 308 (a =4, p =7, q =11)
2 / 77 = 1 / 63 + 1 / 99 (r =9, p =7, q =11)
2 / 77 = 1 / 44 + 1 / 308 (2 / 7 = 1 / 4 + 1 / 28 を 11 で割る。)
79
素数
2 / 79 = 1 / 40 + 1 / 3160
2 / 79 = 1 / 60 + 1 / 237 + 1 / 316 + 1 / 790 (A =60, p =79, 2A -p =41=20+15+6)
81
合成数
2 / 81 = 1 / 41 + 1 / 3321
2 / 81 = 1 / 54 + 1 / 162 (2 / 3 = 1 / 2 + 1 / 6 を 27 で割る。)
83
素数
2 / 83 = 1 / 42 + 1 / 3486
2 / 83 = 1 / 60 + 1 / 332 + 1 / 415 + 1 / 498 (A =60, p =83, 2A -p =37=15+12+10)
85
半素数
2 / 85 = 1 / 43 + 1 / 3655
2 / 85 = 1 / 51 + 1 / 255 (a =3, p =5, q =17)
2 / 85 = 1 / 55 + 1 / 187 (r =11, p =5, q =17)
2 / 85 = 1 / 51 + 1 / 255 (2 / 5 = 1 / 3 + 1 / 15 を 17 で割る。)
87
半素数
2 / 87 = 1 / 44 + 1 / 3828
2 / 87 = 1 / 58 + 1 / 174 (a =2, p =3, q =29)
2 / 87 = 1 / 78 + 1 / 964 (r =16, p =3, q =29)
2 / 87 = 1 / 58 + 1 / 174 (2 / 3 = 1 / 2 + 1 / 6 を 29 で割る。)
89
素数
2 / 89 = 1 / 45 + 1 / 4005
2 / 89 = 1 / 60 + 1 / 356 + 1 / 534 + 1 / 890 (A =60, p =89, 2A -p =31=15+10+6)
91
半素数
2 / 91 = 1 / 46 + 1 / 4416
2 / 91 = 1 / 52 + 1 / 364 (a =4, p =7, q =13)
2 / 91 = 1 / 70 + 1 / 130 (r =10, p =7, q =13)
2 / 91 = 1 / 52 + 1 / 364 (2 / 7 = 1 / 4 + 1 / 28 を 13 で割る。)
93
半素数
2 / 93 = 1 / 47 + 1 / 4371
2 / 93 = 1 / 62 + 1 / 186 (a =2, p =3, q =31)
2 / 93 = 1 / 51 + 1 / 527 (r =17, p =3, q =31)
2 / 93 = 1 / 62 + 1 / 186 (2 / 3 = 1 / 2 + 1 / 6 を 31 で割る。)
95
半素数
2 / 95 = 1 / 48 + 1 / 4560
2 / 95 = 1 / 57 + 1 / 285 (a =3, p =5, q =19)
2 / 95 = 1 / 60 + 1 / 228 (r =12, p =5, q =19)
2 / 95 = 1 / 60 + 1 / 380 + 1 / 570 (2 / 19 = 1 / 12 + 1 / 76 + 1 / 114 を 5 で割る。)
97
素数
2 / 97 = 1 / 49 + 1 / 4753
2 / 97 = 1 / 56 + 1 / 679 + 1 / 776 (A =56, p =97, 2A -p =15=8+7)
99
合成数
2 / 99 = 1 / 50 + 1 / 4950
2 / 99 = 1 / 66 + 1 / 198 (2 / 3 = 1 / 2 + 1 / 6 を 33 で割る。)
101の場合
101
素数
2 / 101 = 1 / 52 + 1 / 5252
2 / 101 = 1 / 101 + 1 / 202 + 1 / 303 + 1 / 606
中世
中国やインドでは、古くから分子に任意の自然数 を許す今日の分数表現を用いたが、ヨーロッパでは17世紀頃までエジプト式が用いられていた。エジプト式分数は実用的な計算には向いておらず、これに固執したことが数学の発展を遅らせたと主張する歴史家もいる[ 1] 。一方で六十進法 で数字を表記したバビロニア では、早い段階から1未満の数を表すのに小数 を導入していた。古代ローマ の天文学者 プトレマイオス は、著書『アルマゲスト 』において「複雑な計算にはエジプト式分数ではなく六十進法を用いる」という趣旨の言を残している[ 8] 。これはプトレマイオスに限った話ではなく、多くの学者が天文計算に六十進法を用いており、角度を度数法で表す際の1度未満の度数単位や、1時間未満の時間の単位が六十進法であるのは、これに由来する。ロナルド・グラハム によると、20世紀を代表する数学者の一人アンドレ・ヴェイユ は、古代エジプト人がエジプト式分数を用いたことについて「間違った方向へ進んだのだ」と語った[ 9] 。
1202年、フィボナッチ は『算盤の書 』において、任意の分数を単位分数の和に表すアルゴリズム をいくつか発表した。まず、分母 n が性質「n 未満の任意の自然数は、いくつかの n の約数の和で表せる」を持つとき、分子を n の約数の和で表して約分することにより、単位分数の和に表せる。そのような性質を持つ n は
1, 2, 4, 6, 8, 12, 16, 18, 20, 24, 28, 30, 32, 36, 40, 42, 48, …(オンライン整数列大辞典 の数列 A5153 )
と続くが、『算盤の書』では例として分母が 6, 8, 12, 20, 24, 60, 100 であるものについて、単位分数展開のリストを与えている。例えば、5 / 12 の分子 5 は、分母 12 の約数の和として 4 + 1 と表せるので、5 / 12 = 1 / 3 + 1 / 12 である。
分母がそのような性質を持たない場合について、フィボナッチは次のような方法を提示している。a / b に対して、b / 2 < c < b を満たし、多くの約数を持つ c を取る。ac / bc の分子を bc の約数の和に表すことができれば、約分して単位分数の和となる。この方法は、リンド・パピルスの表に対して現代数学史家が推測した方法の一つと似ている。
その他の方法のいくつかは
a / ab − 1 = 1 / b + 1 / b (ab − 1)
のような恒等式を用いる。フィボナッチは、例として 8 / 11 を挙げた。まず、分母に 1 を加えた 12 を分子が割るように、2 / 11 + 6 / 11 と分解し、それから恒等式を適用して
8 / 11 = 1 / 2 + 1 / 22 + 1 / 6 + 1 / 66
を導いている。
強欲算法
以上のいずれの方法の通用しない場合に対して、フィボナッチは強欲算法 [ 注釈 1] (greedy algorithm) と呼ばれる方法を提案した。単位分数の和に展開しようとする分数に対して、それ以下の最大の単位分数を取る。それを引いた残りに対しても、繰り返し最大の単位分数を取る。式で書けば、分数 x / y を
x
y
=
1
⌈
y
/
x
⌉
+
x
−
(
y
mod
x
)
y
⌈
y
/
x
⌉
{\displaystyle {\frac {x}{y}}={\frac {1}{\lceil y/x\rceil }}+{\frac {x-(y\,{\bmod {\,}}x)}{y\lceil y/x\rceil }}}
と、次々に置き換える方法である。ここに、括弧は天井関数 である。例えば、4 / 13 に強欲算法を適用すると、
4 / 13 = 1 / 4 + 1 / 18 + 1 / 468
となる。
フィボナッチは、強欲算法の手続きが有限回で終了することの証明を与えてはいない。後にシルベスター はこの方法を再発見し、有限回で終了することの証明も与え、1880年に発表した[ 1] 。実際、一度の手続きで分子は少なくとも 1 小さくなるので、a / b は多くとも a 個の単位分数の和で表せる。2人の名を取って、強欲算法は「フィボナッチ=シルベスターのアルゴリズム」とも呼ばれる。
フィボナッチ自身も注意したように、強欲算法はときに複雑な単位分数展開を与える。例えば、5 / 121 に強欲算法を適用すると、
5 / 121 = 1 / 25 + 1 / 757 + 1 / 763309 + 1 / 873960180913 + 1 / 1527612795642093418846225
となるが、ずっと簡潔な単位分数展開
5 / 121 = 1 / 33 + 1 / 121 + 1 / 363
を持つ。強欲算法は単純で分かりやすく、任意の分数が異なる単位分数の和で表せることの易しい証明も与えるが、このように複雑な展開になる場合もあるため、フィボナッチ自身は、最初の分解の後は他の方法を適用することを勧めている。
単位分数 1 / n に強欲算法を適用すると、恒等式
1 / n = 1 / n + 1 + 1 / n (n + 1)
を得る。これより、単位分数は2つの単位分数の和に表せるので、任意の分数は無数に多くの単位分数展開を持つ。
現代
現代の数論 の研究者は、エジプト式分数に関する多くの問題について研究している。例えば、単位分数展開の項数や分母の大きさを評価すること、ある性質を持った単位分数に限った展開を与えること、単位分数展開のアルゴリズムを与えること、などである。リチャード・ガイ (英語版 ) の本『数論未解決問題の事典』第3版 D11 に、これまでの研究成果と未解決問題の概略がある。
研究成果
任意に与えられた分数 x / y は、単位分数展開として、その最大分母が高々
O
(
y
log
2
y
log
log
y
)
{\displaystyle O\left({\frac {y\log ^{2}y}{\log \log y}}\right)}
であるものを持つ[ 10] 。また、項数が高々
O
(
log
y
)
{\displaystyle O\left({\sqrt {\log y}}\right)}
であるものを持つ[ 11] 。ここに、O はランダウの記号 である。
グラハムは、2以上の任意の自然数 n に対し、分母を n 乗数に限った場合にエジプト式分数として表せるような有理数を特徴付けた[ 12] 。例えば、有理数 q がいくつかの平方数 (1 を含める)の逆数の和として表せるための必要十分条件は、q が2つの半開区間 の和集合
[
0
,
π
2
6
−
1
)
∪
[
1
,
π
2
6
)
{\displaystyle \left[0,{\frac {\pi ^{2}}{6}}-1\right)\cup \left[1,{\frac {\pi ^{2}}{6}}\right)}
に含まれることである。ここに現れる
π
2
6
{\displaystyle {\frac {\pi ^{2}}{6}}}
は、リーマンゼータ関数 ζ の特殊値 ζ (2) である。
エルデシュ とグラハムは、2 以上の整数の集合を有限個の集合に分割 した場合、それがどのような分割であっても、そのうちの一つの集合の有限部分集合 S を取って
∑
n
∈
S
1
n
=
1
{\displaystyle \sum _{n\in S}{\frac {1}{n}}=1}
とできると予想 した。予想の内容は、よく次のように言い換えられる。「単位分数を有限個の色でどのように色分けしても、そのうちの単色のみを用いて 1 の単位分数展開が得られる。」この予想は2003年に証明された[ 13] 。
未解決問題
エジプト式分数に関するオープンプロブレム を以下に挙げる。
複雑性クラス
任意の分数に対し、項数や最大分母が最小の単位分数展開を総当たり法で見つけることはできるが、この問題が計算複雑性理論 においてどの複雑性クラス に属するのか、例えば多項式時間 で見つけられるかは知られていない。
エルデシュ=シュトラウス予想
エルデシュ=シュトラウス予想 は、すべての整数 n ≥ 2 に対し
4
n
=
1
x
+
1
y
+
1
z
{\displaystyle {\frac {4}{n}}={\frac {1}{x}}+{\frac {1}{y}}+{\frac {1}{z}}}
は正の整数解を持つ、という予想である[ 14] 。エルデシュ=シュトラウス予想が成り立つことは n < 1014 まで確かめられている[ 15] 。また、n が 840 を法として 12 , 112 , 132 , 172 , 192 , 232 に合同 な場合を除き、予想が成り立つことが示されている[ 16] 。
シェルピンスキー予想
ヴァツワフ・シェルピニスキ は、2以上の任意の整数 n に対し、
5 / n = 1 / x + 1 / y + 1 / z
は正の整数解を持つと予想した[ 14] 。さらに「任意の整数 k に対し、n が十分大きければ、真分数 k / n は3つ以下の単位分数の和で表せる」とも予想した[ 1] 。
分母が奇数である単位分数に限れば、その和も分母が奇数になる。逆に、分母が奇数である分数は、分母が奇数である単位分数展開が可能であることが知られている[ 1] 。しかし、分母を奇数に限った場合に、強欲算法が有限回で終了するかどうかは知られていない。
数学パズル
単位分数展開は、しばしば数学パズル の題材にもなる。特に人気があるのは、1 の単位分数展開である。1 に対して強欲算法を用いると、有名な表示
1 = 1 / 2 + 1 / 3 + 1 / 6
を得る。
用いることができる単位分数に制限をかけると、より深みのある類似の問題が多く考えられる。例えば、分母を奇数に限ると、強欲算法によって 1 は13個の単位分数の和に表され、最後の項の分母は
209525411280522638000804396401925664136495425904830384693383280180439963265695525939102230139815
となる[ 17] 。項数が最小なのは9項のものであり、そのようなものは全部で5通りあることが知られている[ 1] [ 18] 。その中でも、最大分母が最小であるものは
1 = 1 / 3 + 1 / 5 + 1 / 7 + 1 / 9 + 1 / 11 + 1 / 15 + 1 / 35 + 1 / 45 + 1 / 231
である。項数を度外視した場合、最大分母が最小であるものは、11項の和
1 = 1 / 3 + 1 / 5 + 1 / 7 + 1 / 9 + 1 / 11 + 1 / 33 + 1 / 35 + 1 / 45 + 1 / 55 + 1 / 77 + 1 / 105
である[ 1] 。
1の分割表
分母の和 n
分割数 p (n )
1の分割
0
0
―
1
0
―
2
0
―
3
1
―
4
1
―
5
2
―
6
3
―
7
4
―
8
5
―
9
7
―
10
9
―
11
11
1 = 1 / 2 + 1 / 3 + 1 / 6 .
12
14
―
13
17
―
14
21
―
15
26
―
16
31
―
17
37
―
18
45
―
19
53
―
20
63
―
21
75
―
22
88
―
23
103
―
24
121
1 = 1 / 2 + 1 / 4 + 1 / 6 + 1 / 12 .
25
141
―
26
164
―
27
191
―
28
221
―
29
255
―
30
295
1 = 1 / 2 + 1 / 3 + 1 / 10 + 1 / 15 .
31
339
1 = 1 / 2 + 1 / 4 + 1 / 5 + 1 / 20 .
32
389
1 = 1 / 2 + 1 / 3 + 1 / 9 + 1 / 18 .
33
447
―
34
511
―
35
584
―
36
667
―
37
759
1 = 1 / 2 + 1 / 3 + 1 / 8 + 1 / 24 .
38
863
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 6 + 1 / 20 .
39
981
―
40
1112
―
41
1259
―
42
1425
―
43
1609
1 = 1 / 2 + 1 / 4 + 1 / 10 + 1 / 12 + 1 / 15 .
44
1815
―
45
2047
1 = 1 / 2 + 1 / 4 + 1 / 9 + 1 / 12 + 1 / 18 ,
1 = 1 / 2 + 1 / 5 + 1 / 6 + 1 / 12 + 1 / 20 .
46
2303
―
47
2589
―
48
2909
―
49
3263
―
50
3657
1 = 1 / 2 + 1 / 4 + 1 / 8 + 1 / 12 + 1 / 24 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 10 + 1 / 12 + 1 / 15 .
51
4096
―
52
4581
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 9 + 1 / 12 + 1 / 18 .
53
5119
1 = 1 / 2 + 1 / 5 + 1 / 6 + 1 / 10 + 1 / 30 .
54
5717
1 = 1 / 2 + 1 / 3 + 1 / 7 + 1 / 42 .
55
6377
1 = 1 / 2 + 1 / 4 + 1 / 7 + 1 / 14 + 1 / 28 .
56
7107
―
57
7916
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 10 + 1 / 15 + 1 / 20 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 8 + 1 / 12 + 1 / 24 .
58
8807
―
59
9791
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 9 + 1 / 18 + 1 / 20 .
60
10879
1 = 1 / 2 + 1 / 6 + 1 / 9 + 1 / 10 + 1 / 15 + 1 / 18 .
61
12075
1 = 1 / 2 + 1 / 4 + 1 / 6 + 1 / 21 + 1 / 28 .
62
13393
1 = 1 / 2 + 1 / 4 + 1 / 6 + 1 / 20 + 1 / 30 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 7 + 1 / 14 + 1 / 28 .
63
14847
―
64
16443
1 = 1 / 2 + 1 / 4 + 1 / 8 + 1 / 10 + 1 / 40 ,
1 = 1 / 2 + 1 / 5 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 20 ,
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 8 + 1 / 20 + 1 / 24 ,
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 10 + 1 / 12 + 1 / 30 .
65
18199
1 = 1 / 2 + 1 / 6 + 1 / 8 + 1 / 10 + 1 / 15 + 1 / 24 .
66
20131
1 = 1 / 2 + 1 / 3 + 1 / 12 + 1 / 21 + 1 / 28 ,
1 = 1 / 2 + 1 / 4 + 1 / 6 + 1 / 18 + 1 / 36 ,
1 = 1 / 2 + 1 / 5 + 1 / 9 + 1 / 12 + 1 / 18 + 1 / 20 .
67
22249
1 = 1 / 2 + 1 / 3 + 1 / 12 + 1 / 20 + 1 / 30 ,
1 = 1 / 2 + 1 / 4 + 1 / 7 + 1 / 12 + 1 / 42 ,
1 = 1 / 2 + 1 / 5 + 1 / 6 + 1 / 9 + 1 / 45 ,
1 = 1 / 2 + 1 / 6 + 1 / 8 + 1 / 9 + 1 / 18 + 1 / 24 .
68
24575
―
69
27129
1 = 1 / 2 + 1 / 3 + 1 / 14 + 1 / 15 + 1 / 35 ,
1 = 1 / 2 + 1 / 6 + 1 / 7 + 1 / 12 + 1 / 14 + 1 / 28 .
70
29926
―
71
32991
1 = 1 / 2 + 1 / 3 + 1 / 11 + 1 / 22 + 1 / 33 ,
1 = 1 / 2 + 1 / 3 + 1 / 12 + 1 / 18 + 1 / 36 ,
1 = 1 / 2 + 1 / 5 + 1 / 8 + 1 / 12 + 1 / 20 + 1 / 24 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 8 + 1 / 10 + 1 / 40 ,
1 = 1 / 3 + 1 / 4 + 1 / 9 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 18 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 20 .
72
36351
―
73
40025
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 9 + 1 / 12 + 1 / 18 + 1 / 20 .
74
44045
1 = 1 / 2 + 1 / 5 + 1 / 9 + 1 / 10 + 1 / 18 + 1 / 30 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 7 + 1 / 12 + 1 / 42 .
75
48445
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 8 + 1 / 15 + 1 / 40 .
76
53249
1 = 1 / 2 + 1 / 4 + 1 / 6 + 1 / 16 + 1 / 48 ,
1 = 1 / 2 + 1 / 5 + 1 / 7 + 1 / 14 + 1 / 20 + 1 / 28 ,
1 = 1 / 3 + 1 / 4 + 1 / 8 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 24 .
77
58498
―
78
64233
1 = 1 / 2 + 1 / 6 + 1 / 8 + 1 / 10 + 1 / 12 + 1 / 40 ,
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 9 + 1 / 12 + 1 / 45 ,
1 = 1 / 3 + 1 / 4 + 1 / 8 + 1 / 9 + 1 / 12 + 1 / 18 + 1 / 24 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 8 + 1 / 12 + 1 / 20 + 1 / 24 .
79
70487
1 = 1 / 2 + 1 / 3 + 1 / 10 + 1 / 24 + 1 / 40 ,
1 = 1 / 2 + 1 / 5 + 1 / 8 + 1 / 10 + 1 / 24 + 1 / 30 .
80
77311
1 = 1 / 2 + 1 / 4 + 1 / 10 + 1 / 15 + 1 / 21 + 1 / 28 .
81
84755
1 = 1 / 2 + 1 / 3 + 1 / 12 + 1 / 16 + 1 / 48 ,
1 = 1 / 2 + 1 / 4 + 1 / 10 + 1 / 15 + 1 / 20 + 1 / 30 ,
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 7 + 1 / 20 + 1 / 42 ,
1 = 1 / 3 + 1 / 4 + 1 / 7 + 1 / 10 + 1 / 14 + 1 / 15 + 1 / 28 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 9 + 1 / 10 + 1 / 18 + 1 / 30 .
82
92863
1 = 1 / 2 + 1 / 4 + 1 / 9 + 1 / 18 + 1 / 21 + 1 / 28 ,
1 = 1 / 2 + 1 / 4 + 1 / 12 + 1 / 14 + 1 / 15 + 1 / 35 ,
1 = 1 / 2 + 1 / 5 + 1 / 6 + 1 / 20 + 1 / 21 + 1 / 28 ,
1 = 1 / 2 + 1 / 5 + 1 / 8 + 1 / 12 + 1 / 15 + 1 / 40 ,
1 = 1 / 2 + 1 / 6 + 1 / 7 + 1 / 10 + 1 / 15 + 1 / 42 .
83
101697
1 = 1 / 2 + 1 / 4 + 1 / 9 + 1 / 18 + 1 / 20 + 1 / 30 ,
1 = 1 / 3 + 1 / 4 + 1 / 7 + 1 / 9 + 1 / 14 + 1 / 18 + 1 / 28 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 7 + 1 / 14 + 1 / 20 + 1 / 28 .
84
111321
1 = 1 / 2 + 1 / 4 + 1 / 11 + 1 / 12 + 1 / 22 + 1 / 33 ,
1 = 1 / 2 + 1 / 6 + 1 / 7 + 1 / 9 + 1 / 18 + 1 / 42 .
85
121791
1 = 1 / 2 + 1 / 4 + 1 / 10 + 1 / 14 + 1 / 20 + 1 / 35 ,
1 = 1 / 2 + 1 / 4 + 1 / 10 + 1 / 15 + 1 / 18 + 1 / 36 ,
1 = 1 / 2 + 1 / 5 + 1 / 8 + 1 / 10 + 1 / 20 + 1 / 40 .
86
133183
1 = 1 / 2 + 1 / 5 + 1 / 9 + 1 / 10 + 1 / 15 + 1 / 45 ,
1 = 1 / 2 + 1 / 8 + 1 / 9 + 1 / 10 + 1 / 15 + 1 / 18 + 1 / 24 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 8 + 1 / 10 + 1 / 24 + 1 / 30 .
87
145577
1 = 1 / 2 + 1 / 4 + 1 / 6 + 1 / 15 + 1 / 60 ,
1 = 1 / 2 + 1 / 4 + 1 / 8 + 1 / 21 + 1 / 24 + 1 / 28 ,
1 = 1 / 2 + 1 / 5 + 1 / 6 + 1 / 18 + 1 / 20 + 1 / 36 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 10 + 1 / 15 + 1 / 21 + 1 / 28 ,
1 = 1 / 4 + 1 / 5 + 1 / 6 + 1 / 9 + 1 / 10 + 1 / 15 + 1 / 18 + 1 / 20 .
88
159045
1 = 1 / 2 + 1 / 4 + 1 / 8 + 1 / 20 + 1 / 24 + 1 / 30 ,
1 = 1 / 2 + 1 / 5 + 1 / 7 + 1 / 12 + 1 / 20 + 1 / 42 ,
1 = 1 / 2 + 1 / 7 + 1 / 10 + 1 / 12 + 1 / 14 + 1 / 15 + 1 / 28 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 10 + 1 / 15 + 1 / 20 + 1 / 30 ,
1 = 1 / 3 + 1 / 4 + 1 / 7 + 1 / 8 + 1 / 14 + 1 / 24 + 1 / 28 .
89
173681
1 = 1 / 2 + 1 / 3 + 1 / 9 + 1 / 30 + 1 / 45 ,
1 = 1 / 2 + 1 / 6 + 1 / 7 + 1 / 8 + 1 / 24 + 1 / 42 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 9 + 1 / 18 + 1 / 21 + 1 / 28 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 12 + 1 / 14 + 1 / 15 + 1 / 35 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 8 + 1 / 12 + 1 / 15 + 1 / 40 .
90
189585
1 = 1 / 2 + 1 / 7 + 1 / 9 + 1 / 12 + 1 / 14 + 1 / 18 + 1 / 28 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 9 + 1 / 18 + 1 / 20 + 1 / 30 .
91
206847
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 11 + 1 / 12 + 1 / 22 + 1 / 33 .
92
225584
1 = 1 / 2 + 1 / 3 + 1 / 12 + 1 / 15 + 1 / 60 ,
1 = 1 / 2 + 1 / 4 + 1 / 5 + 1 / 36 + 1 / 45 ,
1 = 1 / 2 + 1 / 4 + 1 / 8 + 1 / 18 + 1 / 24 + 1 / 36 ,
1 = 1 / 2 + 1 / 4 + 1 / 10 + 1 / 12 + 1 / 24 + 1 / 40 ,
1 = 1 / 2 + 1 / 5 + 1 / 6 + 1 / 14 + 1 / 30 + 1 / 35 ,
1 = 1 / 2 + 1 / 5 + 1 / 6 + 1 / 15 + 1 / 24 + 1 / 40 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 10 + 1 / 14 + 1 / 20 + 1 / 35 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 10 + 1 / 15 + 1 / 18 + 1 / 36 ,
1 = 1 / 3 + 1 / 4 + 1 / 8 + 1 / 9 + 1 / 10 + 1 / 18 + 1 / 40 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 8 + 1 / 10 + 1 / 20 + 1 / 40 ,
1 = 1 / 3 + 1 / 5 + 1 / 9 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 18 + 1 / 20 ,
1 = 1 / 4 + 1 / 5 + 1 / 6 + 1 / 8 + 1 / 10 + 1 / 15 + 1 / 20 + 1 / 24 .
93
245919
1 = 1 / 2 + 1 / 5 + 1 / 8 + 1 / 9 + 1 / 24 + 1 / 45 ,
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 12 + 1 / 20 + 1 / 21 + 1 / 28 ,
1 = 1 / 3 + 1 / 4 + 1 / 7 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 42 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 9 + 1 / 10 + 1 / 15 + 1 / 45 ,
1 = 1 / 3 + 1 / 6 + 1 / 8 + 1 / 9 + 1 / 10 + 1 / 15 + 1 / 18 + 1 / 24 .
94
267967
1 = 1 / 2 + 1 / 6 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 21 + 1 / 28 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 8 + 1 / 21 + 1 / 24 + 1 / 28 ,
1 = 1 / 4 + 1 / 5 + 1 / 6 + 1 / 8 + 1 / 9 + 1 / 18 + 1 / 20 + 1 / 24 ,
1 = 1 / 4 + 1 / 5 + 1 / 6 + 1 / 9 + 1 / 10 + 1 / 12 + 1 / 18 + 1 / 30 .
95
291873
1 = 1 / 2 + 1 / 3 + 1 / 9 + 1 / 27 + 1 / 54 ,
1 = 1 / 2 + 1 / 3 + 1 / 10 + 1 / 20 + 1 / 60 ,
1 = 1 / 2 + 1 / 4 + 1 / 8 + 1 / 9 + 1 / 72 ,
1 = 1 / 2 + 1 / 4 + 1 / 9 + 1 / 15 + 1 / 20 + 1 / 45 ,
1 = 1 / 2 + 1 / 4 + 1 / 10 + 1 / 15 + 1 / 16 + 1 / 48 ,
1 = 1 / 2 + 1 / 6 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 20 + 1 / 30 ,
1 = 1 / 2 + 1 / 7 + 1 / 8 + 1 / 12 + 1 / 14 + 1 / 24 + 1 / 28 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 8 + 1 / 20 + 1 / 24 + 1 / 30 ,
1 = 1 / 3 + 1 / 4 + 1 / 7 + 1 / 9 + 1 / 12 + 1 / 18 + 1 / 42 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 7 + 1 / 12 + 1 / 20 + 1 / 42 ,
1 = 1 / 3 + 1 / 6 + 1 / 7 + 1 / 10 + 1 / 12 + 1 / 14 + 1 / 15 + 1 / 28 .
96
317787
1 = 1 / 2 + 1 / 5 + 1 / 7 + 1 / 10 + 1 / 30 + 1 / 42 ,
1 = 1 / 2 + 1 / 6 + 1 / 9 + 1 / 12 + 1 / 18 + 1 / 21 + 1 / 28 ,
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 14 + 1 / 15 + 1 / 20 + 1 / 35 ,
1 = 1 / 4 + 1 / 5 + 1 / 6 + 1 / 7 + 1 / 12 + 1 / 14 + 1 / 20 + 1 / 28 .
97
345855
1 = 1 / 2 + 1 / 4 + 1 / 9 + 1 / 16 + 1 / 18 + 1 / 48 ,
1 = 1 / 2 + 1 / 5 + 1 / 6 + 1 / 16 + 1 / 20 + 1 / 48 ,
1 = 1 / 2 + 1 / 6 + 1 / 9 + 1 / 12 + 1 / 18 + 1 / 20 + 1 / 30 ,
1 = 1 / 3 + 1 / 5 + 1 / 8 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 20 + 1 / 24 ,
1 = 1 / 3 + 1 / 6 + 1 / 7 + 1 / 9 + 1 / 12 + 1 / 14 + 1 / 18 + 1 / 28 .
98
376255
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 11 + 1 / 20 + 1 / 22 + 1 / 33 ,
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 12 + 1 / 18 + 1 / 20 + 1 / 36 .
99
409173
1 = 1 / 2 + 1 / 4 + 1 / 8 + 1 / 15 + 1 / 30 + 1 / 40 ,
1 = 1 / 2 + 1 / 6 + 1 / 9 + 1 / 14 + 1 / 15 + 1 / 18 + 1 / 35 ,
1 = 1 / 2 + 1 / 6 + 1 / 10 + 1 / 11 + 1 / 15 + 1 / 22 + 1 / 33 ,
1 = 1 / 2 + 1 / 6 + 1 / 10 + 1 / 12 + 1 / 14 + 1 / 20 + 1 / 35 ,
1 = 1 / 2 + 1 / 6 + 1 / 10 + 1 / 12 + 1 / 15 + 1 / 18 + 1 / 36 ,
1 = 1 / 2 + 1 / 8 + 1 / 9 + 1 / 10 + 1 / 12 + 1 / 18 + 1 / 40 ,
1 = 1 / 3 + 1 / 4 + 1 / 5 + 1 / 6 + 1 / 36 + 1 / 45 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 8 + 1 / 18 + 1 / 24 + 1 / 36 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 10 + 1 / 12 + 1 / 24 + 1 / 40 ,
1 = 1 / 3 + 1 / 5 + 1 / 8 + 1 / 9 + 1 / 12 + 1 / 18 + 1 / 20 + 1 / 24 ,
1 = 1 / 4 + 1 / 5 + 1 / 6 + 1 / 8 + 1 / 10 + 1 / 12 + 1 / 24 + 1 / 30 .
100
444792
1 = 1 / 2 + 1 / 6 + 1 / 7 + 1 / 8 + 1 / 21 + 1 / 56 ,
1 = 1 / 3 + 1 / 4 + 1 / 7 + 1 / 8 + 1 / 12 + 1 / 24 + 1 / 42 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 8 + 1 / 9 + 1 / 24 + 1 / 45 .
1 の単位分数展開で分母の和が 50 および 100 となるものは、それぞれ
1 = 1 / 2 + 1 / 4 + 1 / 8 + 1 / 12 + 1 / 24 ,
1 = 1 / 3 + 1 / 4 + 1 / 6 + 1 / 10 + 1 / 12 + 1 / 15 .
および
1 = 1 / 2 + 1 / 6 + 1 / 7 + 1 / 8 + 1 / 21 + 1 / 56 ,
1 = 1 / 3 + 1 / 4 + 1 / 7 + 1 / 8 + 1 / 12 + 1 / 24 + 1 / 42 ,
1 = 1 / 3 + 1 / 5 + 1 / 6 + 1 / 8 + 1 / 9 + 1 / 24 + 1 / 45 .
である[ 注釈 2] 。
なお、正整数の逆数の和である調和級数 は無限大に発散するから、(真分数に限らず)任意の正の有理数は単位分数の和で表すことができる。ただし、調和級数の発散は非常にゆっくりであるため、比較的小さな有理数であっても多くの単位分数が必要になる。例えば 10 を単位分数の和で表すには、2万個を超える項が必要である[ 1] 。
計算法の例
a
n
=
1
x
+
1
y
+
1
z
{\displaystyle {\frac {a}{\;n\;}}={\frac {1}{\;x\;}}+{\frac {1}{\;y\;}}+{\frac {1}{\;z\;}}}
に対して,次の計算法が報告されている[ 19] 。
x
=
⌊
n
a
⌋
+
b
,
{\displaystyle x={\bigg \lfloor }{\frac {\;n\;}{a}}{\bigg \rfloor }+b,}
y
=
2
c
n
x
d
,
{\displaystyle y={\frac {\;2\,c\,n\,x\;}{d}},}
z
=
2
c
n
x
2
c
(
a
x
−
n
)
−
d
,
{\displaystyle z={\frac {\;2\,c\,n\,x\;}{\;2\,c\,(a\,x-n)-d\;\;}},}
b
=
1
,
2
,
3
,
⋯
.
c
=
1
,
2
,
3
,
⋯
.
d
=
1
,
2
,
3
,
⋯
,
c
(
a
x
−
n
)
.
{\displaystyle b=1,2,3,\cdots .\;\;\;\;c=1,2,3,\cdots .\;\;\;\;d=1,2,3,\cdots ,c\,(a\,x-n).}
上記の記号で,
⌊
⋆
⌋
{\displaystyle \lfloor \star \rfloor }
は,床関数 である。
y の右辺 2cnx / d および,
z の右辺 2cnx / 2c (ax − n ) − d が整除されれば解を得る。
脚注
注釈
^ ガードナーの記事の一松による訳[ 1] 。他に、ホフマンの本の平石による訳では「欲張り展開法」[ 9] 。計算機科学ではgreedy algorithmは貪欲算法 と訳される。
^ 例えば分母の和が 100 の解を得るには、100 を相異なる自然数の和に分割したリストを生成して、リストの各自然数の逆数の和が 1 となるものを選べばよい。
出典
^ a b c d e f g h i ガードナー 2010 , pp. 54–58
^ a b カッツ 2005 , p. 13
^ ボイヤー 2009 , p. 19
^ カジョリ 1997 , p. 31
^ 上垣 2006 , pp. 21–23
^ ボイヤー 2009 , p. 17
^ 上垣 2006 , p. 24
^ 一松 2007 , p. 33
^ a b ホフマン 2000 , pp. 169–170
^ Tenenbaum, G.; Yokota, H. (1990), "Length and denominators of Egyptian fractions", Journal of Number Theory 35: 150--156, doi :10.1016/0022-314X(90)90109-5 , MR 1057319 .
^ Vose, M. (1985), "Egyptian fractions", Bulletin of the London Mathematical Society 17: 21, doi :10.1112/blms/17.1.21 , MR 0766441 .
^ Graham, R. L. (1964), "On finite sums of reciprocals of distinct nth powers" , Pacific Journal of Mathematics 14 (1): 85--92, MR 0159788 .
^ Croot, Ernest S., III (2003). "On a coloring conjecture about unit fractions". Annals of Mathematics 157 (2): 545--556. doi :10.4007/annals.2003.157.545 . arXiv :math.NT/0311421 .
^ a b “2011 日本数学コンクールのまとめ ” (PDF). 日本数学コンクール委員会. p. 9. 2018年9月8日 閲覧。
^ Allan Swett, The Erdos-Straus Conjecture
^ ガイ 2010 , p. 244
^ オンライン整数列大辞典 の数列 A130738
^ 『数学セミナー 』1971年12月号
^ 長島 隆廣 『数学セミナー』第22巻,第11号,通巻264号,1983年11月発行,日本評論社,pp.92-93。
参考文献
関連文献
Anshel, Michael M.; Goldfeld, Dorian (1991), “Partitions, Egyptian fractions, and free products of finite abelian groups”, Proceedings of the American Mathematical Society 111 (4): 889--899, doi :10.1090/S0002-9939-1991-1065083-1 , MR 1065083 .
Beeckmans, L. (1993), “The splitting algorithm for Egyptian fractions”, Journal of Number Theory 43 : 173--185, doi :10.1006/jnth.1993.1015 , MR 1207497 .
Botts, Truman (1967), “A chain reaction process in number theory” , Mathematics Magazine 40 (2): 55--65, doi :10.2307/2688508 , MR 0209217 , http://www.jstor.org/stable/2688508 .
Breusch, R. (1954), “A special case of Egyptian fractions, solution to advanced problem 4512”, American Mathematical Monthly 61 : 200--201 .
Bruins, Evert M. (1957), “Platon et la table 2/n e'gyptiennes [Plato and the Egyptian table 2/n]” (French), Janus 46 : 253--263 .
Eves, Howard (1953), An Introduction to the History of Mathematics , Holt, Reinhard, and Winston, ISBN 0030295580 .
Gardner, Milo (2002), “The Egyptian Mathematical Leather Roll, attested short term and long term”, in Gratton-Guiness, Ivor, History of the Mathematical Sciences , Hindustan Book Co, pp. 119--134, ISBN 8185931453
Gillings, Richard J. (1982), Mathematics in the Time of the Pharaohs , Dover, ISBN 048624315X .
Hultsch, Friedrich (1895), “Die Elemente der a"gyptischen Theilungsrechnung: Erste Anhandlung” (German), Abhandlungen der philologisch-historischen Classe der Ko"niglich-Sa"chsischen Gesellschaft der Wissenschaften, Sa"chsische Akademie der Wissenschaften zu Leipzig Philologisch-Historische Klasse (Leipzig: S. Hirzel) 17 (1) .
Knorr, Wilbur R. (1982), “Techniques of fractions in ancient Egypt and Greece”, Historia Mathematica 9 : 133--171, doi :10.1016/0315-0860(82)90001-5 , MR 0662138 .
Martin, G. (1999), “Dense Egyptian fractions”, Transactions of the American Mathematical Society 351 : 3641--3657, doi :10.1090/S0002-9947-99-02327-2 , MR 1608486 .
Sigler, Laurence E. (trans.) (2002), Fibonacci's Liber Abaci , Springer-Verlag, ISBN 0387954198
Stewart, B. M. (1954), “Sums of distinct divisors” , American Journal of Mathematics 76 (4): 779--785, doi :10.2307/2372651 , MR 0064800 , http://jstor.org/stable/2372651 .
Stewart, I. (1992), “The riddle of the vanishing camel”, Scientific American (June): 122--124 .
Struik, Dirk J. (1967), A Concise History of Mathematics , Dover, pp. 20--25, ISBN 0486602559 .
Takenouchi, T. (1921), “On an indeterminate equation”, Proceedings of the Physico-Mathematical Society of Japan [Nippon Suugaku-Buturigakkwai Kizi], 3rd ser. 3 : 78--92 .
Wagon, Stan (1999), Mathematica in Action , Springer, pp. 321--329, ISBN 0387986847 .
関連項目
外部リンク