メルセンヌ数
メルセンヌ数 (メルセンヌすう、英 : Mersenne number )とは、2の冪 よりも 1 小さい自然数 、すなわち 2n − 1 (n は自然数)の形の自然数のことである。これを Mn で表すことが多い。メルセンヌ数を小さい順に列挙すると
1 ,
3 ,
7 ,
15 ,
31 ,
63 ,
127 ,
255 ,
511 ,
1023 ,
2047 , 4095,
8191 , ... (
オンライン整数列大辞典 の数列
A000225 )
となる。メルセンヌ数は2進法 表記で n 桁の 11⋯11 、すなわちレピュニット となる。
「メルセンヌ数」の名は、この形の素数に関する予想を発表したマラン・メルセンヌ に由来する。
基本的な性質
Mn が素数ならば n もまた素数である。このことは、nが合成数のメルセンヌ数が次のように因数分解できることから分かる[ 1] :
2ab − 1 = (2a − 1)(1 + 2a + 22a + ⋯ + 2(b −1)a ).
また、この等式より、m | n のとき M m | M n である。一方、 p が素数でも Mp が素数とは限らない。最小の反例は p = 11 の場合であり、M 11 = 2047 = 23 × 89 が成り立つ。
奇素数 p に対して Mp が素数であるかどうかは、リュカ-レーマー・テストによって判定できる (#素数判定法 節を参照)。
完全数
Mp = 2p − 1 が素数ならば、2p −1 (2p − 1) は完全数 である[ 1] [ 3] 。この定理はすでに紀元前3世紀 頃のユークリッド原論 で証明されていた。したがって、完全数の探索はメルセンヌ素数の探索に終始された。
2p −1 (2p − 1) は明らかに偶数 であるが、偶数の完全数でこの生成式から得られるもの以外はないのか2000年間にわたって未解決であったが、18世紀 にオイラー によりこの形に限ることが証明された[ 3] 。
メルセンヌ数の素因数
p を素数とする。
Mp の素因数は 2p を法として 1 と合同、かつ 8 を法として 1 または −1 と合同である[ 6] 。
p ≡ 3 (mod 4) のとき、Mp が 2p + 1 で割れることと、2p + 1 が素数であることは同値である[ 6] 。
ある計算可能な正定数 c が存在して、Mp の最大素因数を q について、q ≥ cp log p [ 7] 。
メルセンヌ素数
メルセンヌ素数 (メルセンヌそすう、Mersenne prime)とは、素数であるメルセンヌ数のことである。
2024年10月現在発見されているメルセンヌ素数は全部で52個ある。その中で最大のものは2024年 10月 に発見された2136279841 − 1 であり、十進法 で表記したときの桁数は4102万4320桁に及ぶ[ GIMPS 1] 。
これより大きい素数は、2024年10月現在メルセンヌ素数以外でも発見されていない。
メルセンヌ素数の発見の歴史
古代~中世
メルセンヌ素数の探求は紀元前3世紀 ごろに端を発する。古代エジプト の数学者エウクレイデス は『原論 』の中で、「2n − 1 が素数ならば、2n − 1 (2n − 1) は完全数 である」ことを証明した[ 8] 。ここから、メルセンヌ素数の探索は完全数の探索にも繋がることとなる[ 9] [ 注釈 1] 。
小さいメルセンヌ素数がいつから知られているかは定かではないが、少なくとも最初の4つの完全数はゲラサのニコマコス の『算術入門 』ですでに言及されている[ 10] [ 11] 。5番目から7番目の完全数は、13世紀イスラム の数学者イブン・ファッルース (fr:Ibn Fallus ) が論文に記している[ 12] 。ヨーロッパでは、5番目の完全数が1456年と1461年の日付が付された古い写本に記されて[ 13] [ 14] おり、6番目と7番目のメルセンヌ素数および完全数も1603年 にピエトロ・カタルディ (英 : Pietro Cataldi ) によって発見されている[ 12] 。
メルセンヌの予想
メルセンヌの予想の表: p ≦ 263
〇:Mp が素数の場合/×:Mp が合成数の場合水色が正解 、ピンク色が間違い を示す[ 15] 。
p
2
3
5
7
11
13
17
19
Mp
〇
〇
〇
〇
×
〇
〇
〇
p
23
29
31
37
41
43
47
53
Mp
×
×
〇
×
×
×
×
×
p
59
61
67
71
73
79
83
89
Mp
×
〇
×
×
×
×
×
〇
p
97
101
103
107
109
113
127
131
Mp
×
×
×
〇
×
×
〇
×
p
137
139
149
151
157
163
167
173
Mp
×
×
×
×
×
×
×
×
p
179
181
191
193
197
199
211
223
Mp
×
×
×
×
×
×
×
×
p
227
229
233
239
241
251
257
263
Mp
×
×
×
×
×
×
×
×
1644年 、マラン・メルセンヌ は「素数 p で 2p − 1 が素数になるのは、p ≦ 257 では p = 2, 3, 5, 7, 13, 17, 19, 31, 67, 127, 257 の11個の場合だけである」という予想を公表した[ 12] [ 16] 。しかしメルセンヌ自身はその予想を証明することができず、しかもその予想の一部は誤っていた。
成果を見るのはメルセンヌが予想を公表してから128年後、1772年 、オイラー (p = 31 では素数)[ 1] 。その次の成果はさらに104年後、1876年 、リュカ (効率的な素数判定法 リュカ・テスト (英語版 ) を考案、p = 67 では素数でない、p = 127 では素数[ 1] )であった。その後リュカ・テストは改良が加えられ、メルセンヌが予想した範囲にない3個が付け加えられた(p = 61 (1883年 )、p = 89 (1911年 )、p = 107 (1914年 ))。メルセンヌが予想した最後の数 p = 257 について決着がついたのは1922年 のことであり、 p = 257 も合成数だった[ 1] 。
結局メルセンヌの11個の予想のうち2つは外れた。なおかつ、間に予想できなかった3つが含まれていたことを考えれば予想は正しかったとはいえないが、その後の歴史を見ても大きな原動力となり先駆的であったことに敬意を表し、素数であるメルセンヌ数をメルセンヌ素数 という[要出典 ] 。
1903年 10月 、アメリカの数学者フランク・ネルソン・コール は実際の素因数分解を探し求め、ニューヨーク で開かれたアメリカ数学会 の会議で 193707721 × 761838257287 を黒板に計算し、M 67 と一致することを証明した。この間一言もしゃべらず、席に戻った後、少し間を置いて拍手が沸き起こったと伝えられている。
1952年 、ラファエル・M・ロビンソン が SWAC を利用して M 521 から M 2281 まで、5つのメルセンヌ素数を発見[ 1] して以降、発見にはコンピュータ が使用されており、コンピュータの進歩と共に新たなメルセンヌ素数が発見されつつある。
GIMPSによる発見
1996年 、メルセンヌ素数を発見することを目的として作られた分散コンピューティング によるプロジェクト GIMPS が発足し、35番目のメルセンヌ素数 M 1,398,269 (1996年11月13日 、Joel Armengaud[ GIMPS 2] )以来、GIMPSによるメルセンヌ素数の発見が続いている。
2008年 8月23日 、GIMPS は46番目の素数候補が、カリフォルニア大学ロサンゼルス校 の数学部のコンピュータによって発見されたと報じた[要出典 ] 。この素数は電子フロンティア財団 が賞金を懸けた1000万桁以上の最初の素数となるため、GIMPS によって同校数学部に50,000ドル 、慈善事業に25,000ドル、残りを前の6つのメルセンヌ素数の発見者へ分配することになった[ GIMPS 3] 。
2008年9月6日 、GIMPS は45番目の素数候補が、ドイツで発見されたと報じた[要出典 ] 。これは、GIMPS によって発見された中では、発見順序と桁数が逆転した初めてのケースである。
素数判定法
知られている素数の中で最大のものが1876年以降ほぼ一貫してメルセンヌ素数である理由は、この判定法にある[要出典 ] 。
リュカ・テスト
p が (4j + 3) 型の素数のとき、S 0 = 3 , Sn = S n −12 − 2 (n ≥ 1) で {Sn } を定義すると、
M
p
∤
S
k
(
0
≦ ≦ -->
k
≦ ≦ -->
p
− − -->
2
)
{\displaystyle M_{p}\not \mid S_{k}\ (0\leqq k\leqq p-2)}
ならば、Mp は合成数である
M
p
∣ ∣ -->
S
p
− − -->
2
{\displaystyle M_{p}\mid S_{p-2}}
ならば、Mp は素数である
アルゴリズム
アルゴリズムは以下の擬似コードで表される。
入力 : p : (4j + 3) 型の素数であるテスト対象の整数
出力 : PRIME :素数の場合, COMPOSIT :合成数の場合
Lucas_Test (p ):
var s = 3
var MP = (1 << p) − 1
for n in range (2, p):
s = (s2 − 2) % MP
if s == 0 then :
return PRIME
else :
return COMPOSIT
リュカ–レーマー・テスト
p が奇素数のとき、S 0 = 4 , Sn = S n −12 − 2 (n ≥ 1) で{Sn } を定義すると、
M
p
∤
S
k
(
0
≦ ≦ -->
k
≦ ≦ -->
p
− − -->
2
)
{\displaystyle M_{p}\not \mid S_{k}\ (0\leqq k\leqq p-2)}
ならば、Mp は合成数である
M
p
∣ ∣ -->
S
p
− − -->
2
{\displaystyle M_{p}\mid S_{p-2}}
ならば、Mp は素数である
リュカ–レーマー・テストは二進 計算機 用のアルゴリズム に向いており、コンピュータによるメルセンヌ素数の発見には、この判定法が用いられてきた。例えば、2p ≡ 1 (mod Mp ) より、A ·2p + B ≡ A + B (mod Mp ) が成り立つので、Mp で割る割り算の代わりに、二進法で p 桁のシフト演算 と足し算だけで計算できる。
アルゴリズム
アルゴリズムは以下の擬似コードで表される。
入力 : p :奇素数であるテスト対象の整数
出力 : PRIME :素数の場合, COMPOSIT :合成数の場合
Lucas_Lehmer_Test (p ):
var s = 4
var MP = (1 << p) − 1
for n in range (2, p):
s = (s2 − 2) % MP
if s == 0 then :
return PRIME
else :
return COMPOSIT
入力 : p :奇素数であるテスト対象の整数
出力 : PRIME :素数の場合, COMPOSIT :合成数の場合
Lucas_Lehmer_Test_FAST (p ):
var s = 4
var m = 2p − 1
for n in range (2, p):
var s2 = s × s
s = (s2 & m) + (s2 >> p)
if s >= m then
s = s − m
s = s − 2
if s == 0 then
return PRIME
else
return COMPOSIT
メルセンヌ素数の一覧
2024年10月現在知られているメルセンヌ素数は下記の表の52個である。ただし、メルセンヌ素数としての番号が確定しているものは48番目までであり、これより大きいメルセンヌ素数については、間に未発見のメルセンヌ素数がないかどうか検証中である。
メルセンヌ素数は小さい順番から並べると
となる。
#
p
Mp の 桁数
発見日
発見者
1
2
1
紀元前500年?[ 27]
2
3
1
紀元前500年?[ 27]
3
5
2
紀元前275年?[ 27]
4
7
3
紀元前275年?[ 27]
5
13
4
1456年[ 9]
不明[ 9]
6
17
6
1603年[ 12]
ピエトロ・カタルディ (英語版 )
7
19
6
1603年[ 12]
ピエトロ・カタルディ
8
31
10
1772年
レオンハルト・オイラー
9
61
19
1883年
イヴァン・パヴシン (英語版 )
10
89
27
1911年
R・E・パワーズ (英語版 )
11
107
33
1914年
R・E・パワーズ[ 28]
12
127
39
1876年
エドゥアール・リュカ
13
521
157
1952年1月30日
ラファエル・M・ロビンソン (英語版 ) , 使用:SWAC
14
607
183
1952年1月30日
ラファエル・M・ロビンソン
15
1,279
386
1952年6月25日
ラファエル・M・ロビンソン
16
2,203
664
1952年10月7日
ラファエル・M・ロビンソン
17
2,281
687
1952年10月9日
ラファエル・M・ロビンソン
18
3,217
969
1957年9月8日
ハンス・リーゼル , 使用:BESK
19
4,253
1,281
1961年11月3日
アレクサンダー・フルウィッツ , 使用:IBM 7090
20
4,423
1,332
1961年11月3日
アレクサンダー・フルウィッツ
21
9,689
2,917
1963年5月11日
ドナルド・ギリース , 使用:ILLIAC II
22
9,941
2,993
1963年5月16日
ドナルド・ギリース
23
11,213
3,376
1963年6月2日
ドナルド・ギリース
24
19,937
6,002
1971年3月4日
ブライアント・タッカーマン (英語版 ) , 使用:IBM 360 /91
25
21,701
6,533
1978年10月30日
ランドン・カート・ノル (英語版 ) & ローラ・ニッケル , 使用:CDC Cyber 174
26
23,209
6,987
1979年2月9日
ランドン・カート・ノル
27
44,497
13,395
1979年4月8日
ハリー・ネルソン & デイヴィッド・スローウィンスキー (英語版 )
28
86,243
25,962
1982年9月25日
デイヴィッド・スローウィンスキー
29
110,503
33,265
1988年1月28日
ウォルター・コルキット & ルーク・ウェルシュ
30
132,049
39,751
1983年9月19日[ 27]
デイヴィッド・スローウィンスキー
31
216,091
65,050
1985年9月1日[ 27]
デイヴィッド・スローウィンスキー
32
756,839
227,832
1992年2月19日
デイヴィッド・スローウィンスキー & ポール・ゲイジ (英語版 ) 使用:Harwell Lab Cray-2 [ 29]
33
859,433
258,716
1994年1月4日[ 30]
デイヴィッド・スローウィンスキー & ポール・ゲイジ
34
1,257,787
378,632
1996年9月3日
デイヴィッド・スローウィンスキー & ポール・ゲイジ[ 31]
35
1,398,269
420,921
1996年11月13日
GIMPS / Joel Armengaud[ GIMPS 2]
36
2,976,221
895,932
1997年8月24日
GIMPS / Gordon Spence[ GIMPS 4]
37
3,021,377
909,526
1998年1月27日
GIMPS / Roland Clarkson[ GIMPS 5]
38
6,972,593
2,098,960
1999年6月1日
GIMPS / Nayan Hajratwala[ GIMPS 6]
39
13,466,917
4,053,946
2001年11月14日
GIMPS / Michael Cameron[ GIMPS 7]
40
20,996,011
6,320,430
2003年11月17日
GIMPS / Michael Shafer[ GIMPS 8]
41
24,036,583
7,235,733
2004年5月15日
GIMPS / Josh Findley[ GIMPS 9]
42
25,964,951
7,816,230
2005年2月18日
GIMPS / Martin Nowak et al.[ GIMPS 10]
43
30,402,457
9,152,052
2005年12月15日
GIMPS / カーティス・クーパー (英語版 ) , Steven Boone[ GIMPS 11]
44
32,582,657
9,808,358
2006年9月4日
GIMPS / カーティス・クーパー, Steven Boone[ GIMPS 12]
45
37,156,667
11,185,272
2008年9月6日
GIMPS / Hans-Michael Elvenich[ GIMPS 3]
46
42,643,801
12,837,064
2009年4月12日
GIMPS / Odd Magnar Strindmo
47
43,112,609
12,978,189
2008年8月23日
GIMPS / エドソン・スミス[ GIMPS 3]
48
57,885,161
17,425,170
2013年1月25日
GIMPS / カーティス・クーパー[ 32] [ GIMPS 13]
[ * 1]
74,207,281
22,338,618
2016年1月7日
GIMPS / カーティス・クーパー[ GIMPS 14]
[ * 1]
77,232,917
23,249,425
2017年12月26日
GIMPS / Jonathan Pace[ GIMPS 15]
[ * 1]
82,589,933
24,862,048
2018年12月7日
GIMPS / Patrick Laroche[ GIMPS 16]
[ * 1]
136,279,841
41,024,320
2024年10月21日
GIMPS / Luke Durant[ GIMPS 1]
^ a b c d 49番目以降はメルセンヌ素数としての順番が確定していない。
未解決問題
メルセンヌ素数は無数に存在するか?
素数 p に対して Mp が合成数であるとき、これをメルセンヌ合成数 と呼ぶことにして、それは無数に存在するか?
平方因子を持つメルセンヌ数 Mp (p は素数)が存在するか?
n を奇数とするとき、次の3つの条件のうち2つが満たされれば、残りの1つも満足されると予想されており、n < 105 に対してこの予想は正しいと確認されている[ 33] 。
Mn が素数
n = 2k ± 1 または 4k ± 3
(2n + 1)/3 が素数
脚注
注釈
^ 2n − 1 (2n − 1) は偶数であるため、この式は奇数の完全数 について何も言及しない。また、偶数の完全数がこの形に限られることは18世紀にレオンハルト・オイラー が証明するまで未解決であった。
出典
^ a b c d e f 淡中 1982 , pp. 65–67.
^ a b 和田 1981 , pp. 59–61
^ a b 和田 1981 , p. 193
^ Theorem 1, Erdős & Shoray 1976
^ 『原論』 第9巻, 命題36. (pp. 225–226, ユークリッド 1971 )
^ a b c “Mersenne Primes: History, Theorems and Lists ” (英語). PrimePages . 2022年2月21日時点のオリジナル よりアーカイブ。2022年2月22日 閲覧。
^ Voight, John (1998年5月31日). “PERFECT NUMBERS: AN ELEMENTARY INTRODUCTION ” (PDF) (英語). 2021年6月28日時点のオリジナル よりアーカイブ。2022年2月21日 閲覧。
^ Nicomachus of Gerasa (1926). Introduction to Arithmetic . Martin Luther D'Ooge (trans). The Macmillan Company. pp. 207–212. https://archive.org/details/NicomachusIntroToArithmetic
^ a b c d e O'Conner, John; Edmund F. Robertson . “Perfect numbers ” (英語). MacTutor History of Mathematics . 2022年1月11日時点のオリジナル よりアーカイブ。2022年2月22日 閲覧。
^ Dickson, Leonard E. (1919) (英語). History of the Theory of Numbers . 1 . Carnegie Institution of Washington. p. 6. https://archive.org/details/historyoftheoryo01dick/page/n5/mode/2up
^ Anonymous (1456), Calendarium ecclesiasticum - BSB Clm 14908 (Anonymous Manuscript), https://iiif.biblissima.fr/collections/manifest/4dca8a1c8b0be3df1ce6954d8d1ba60590cd04cc 2022年2月21日 閲覧。
^ オンライン整数列大辞典 の数列 A000043
^ Mersenne, Marin (1644) (ラテン語). Cogitata Physico Mathematica . Paris: Antonii Bertier. doi :10.3931/e-rara-11036
^ a b c d e f Landon Curt Noll , Mersenne Prime Digits and Names .
^ The Prime Pages, M107 : Fauquembergue or Powers? .
^ The Prime Pages, The finding of the 32nd Mersenne .
^ Chris Caldwell, The Largest Known Primes .
^ The Prime Pages, A Prime of Record Size! 21257787 − 1 .
^ CASEY JOHNSTON (2013年2月7日). “「これまでで最大の素数」を発見” . WIRED (WIRED.jp). http://wired.jp/2013/02/07/volunteer-discovers-a-new-17-million-digit-prime-number/ 2013年2月10日 閲覧。
^ P. T. Bateman, J. L. Selfridge, S. S. Wagstaff Jr., "The new Mersenne conjecture", American Mathematical Monthly, 96 (1989), 125–128.
GIMPS
^ a b "GIMPS Discovers Largest Known Prime Number: 2136,279,841 -1" (Press release) (英語). GIMPS. 21 October 2024. 2024年10月21日閲覧 。
^ a b "GIMPS Discovers 35th Mersenne Prime, 21,398,269 -1 is now the Largest Known Prime" (Press release) (英語). GIMPS. 23 November 1996.
^ a b c "GIMPS Discovers 45th and 46th Mersenne Primes, 243,112,609 -1 is now the Largest Known Prime" (Press release) (英語). GIMPS. 15 September 2008. 2022年2月19日閲覧 。
^ "GIMPS Discovers 36th Mersenne Prime, 22,976,221 -1 is now the Largest Known Prime" (Press release) (英語). GIMPS. 1 September 1997. 2022年2月19日閲覧 。
^ "GIMPS Discovers 37th Mersenne Prime, 23,021,377 -1 is now the Largest Known Prime" (Press release) (英語). GIMPS. 2 February 1998.
^ "GIMPS Discovers 38th Mersenne Prime 26,972,593 -1 is now the Largest Known Prime. Stakes Claim to $50,000 EFF Award" (Press release) (英語). GIMPS. 30 June 1999. 2022年2月19日閲覧 。
^ "GIMPS Discovers 39th Mersenne Prime, 213,466,917 -1 is now the Largest Known Prime. Researchers Discover Largest Multi-Million-Digit Prime Using Entropia Distributed Computing Grid" (Press release) (英語). GIMPS. 6 December 2001. 2022年2月19日閲覧 。
^ "GIMPS Discovers 40th Mersenne Prime, 220,996,011 -1 is now the Largest Known Prime. Mersenne Project Discovers Largest Known Prime Number on World-Wide Volunteer Computer Grid" (Press release) (英語). GIMPS. 2 December 2003. 2022年2月19日閲覧 。
^ "GIMPS Discovers 41st Mersenne Prime, 224,036,583 -1 is now the Largest Known Prime. Project Leaders Believe $100,000 Award Within Reach" (Press release) (英語). GIMPS. 28 May 2004. 2022年2月19日閲覧 。
^ "GIMPS Discovers 42nd Mersenne Prime, 225,964,951 -1 is now the Largest Known Prime" (Press release) (英語). GIMPS. 27 February 2005. 2022年2月19日閲覧 。
^ "GIMPS Discovers 43rd Mersenne Prime, 230,402,457 -1 is now the Largest Known Prime" (Press release) (英語). GIMPS. 24 December 2005. 2022年2月19日閲覧 。
^ "GIMPS Discovers 44th Mersenne Prime, 232,582,657 -1 is now the Largest Known Prime" (Press release) (英語). GIMPS. 11 September 2006. 2022年2月19日閲覧 。
^ "GIMPS Discovers 48th Mersenne Prime, 257,885,161 -1 is now the Largest Known Prime" (Press release) (英語). GIMPS. 5 February 2013. 2022年2月19日閲覧 。
^ "GIMPS Project Discovers Largest Known Prime Number: 274,207,281 -1" (Press release) (英語). GIMPS. 19 January 2016. 2022年3月30日閲覧 。
^ "GIMPS Project Discovers Largest Known Prime Number: 277,232,917 -1" (Press release) (英語). GIMPS. 3 January 2018. 2022年2月19日閲覧 。
^ "GIMPS Discovers Largest Known Prime Number: 282,589,933 -1" (Press release) (英語). GIMPS. 21 December 2018. 2022年2月19日閲覧 。
参考文献
関連項目
外部リンク
生成式 漸化式 (英語版 ) 各種の性質 基数依存 組
互いに素
双子 (p , p + 2 )
Bi-twin chain (n − 1, n + 1, 2n − 1, 2n + 1, … )
三つ子 (p , p + 2 or p + 4, p + 6 )
四つ子 (p , p + 2, p + 6, p + 8 )
k −Tuple
いとこ (p , p + 4 )
セクシー (p , p + 6 )
陳
ソフィー・ジェルマン (p , 2p + 1 )
カニンガム鎖 (p , 2p ± 1, … )
安全 (p , (p − 1)/2 )
算術数列 (英語版 ) (p + an ; n = 0, 1, … )
平衡 (p − n , p , p + n )
桁数 複素数 合成数 関連する話題 最初の50個
素数の一覧