✌🏽 😺 👩🏿‍🎤 メルセンヌ素数とルーク・レマー検定 👩🏿‍🤝‍👩🏼 👩🏿‍⚕️ 🚶🏿

ジョン・マギーの投稿「メルセンヌ・プライムスとルーカス-レーマー・テスト」の翻訳。

記事のコードはここからダウンロードできます。

出版の翻訳と準備の手助けをしてくれたPolina Sologubに深く感謝します。

内容

- はじめに。

- オイラーおよびメルセンヌの乗数定理

- ルークとレーマー

- から

M_{13}

${M_ {13}}$ 前に

M_{20}

${M_ {20}}$

- 完璧な数字

-21、22、23番目のメルセンヌ数

-24、25、26番のメルセンヌ番号。

- 27日と28日のメルセンヌ

- 第29回メルセンヌ

- 30日と31日のメルセンヌ

- メルセンヌ番号の優れたインターネット検索

- メルセンヌ数の因数分解

はじめに

メルセンヌ素数は次の形式の素数です

M_{p} = 2^{p} - 1

${M_p} = {2 ^ p}-1$ （次数pの値も単純でなければなりません）。これらの素数は、17世紀前半にこの形式の素数のリストを作成したフランスの数学者および宗教学者メルセンヌに由来しています。それらの最初の4つは長い間知られています。

M_{2} = 3

${M_2} = 3$ 、

M_{3} = 7

${M_3} = 7$ 、

M_{5} = 31

${M_5} = 31$ そして

M_{7} = $ 12

${M_7} = $ 12$ 。

メルセンヌは

2^{p} - 1

${2 ^ p}-1$ 素数に対しては単純になります

p l e q s l a n t 257

$p \ leqslant 257$ 多くに属する

p i n l e f t t e x t 2 t e x t 、 3 t e x t 、 5 t e x t 、 7 t e x t 、 13 \テ キ ス ト 、 17 \テ キ ス ト 、 19 \テ キ ス ト 、 31 \テ キ ス ト 、 67 \テ キ ス ト 、 127 \テ キ ス ト 、 257 \右

$p \ in \ left \ {{{\ text {2}} {\ text {、3}} {\ text {、5}} {\ text {、7}} {\ text {、13}} {\テキスト{、17}} {\テキスト{、19}} {\テキスト{、31}} {\テキスト{、67}} {\テキスト{、127}} {\テキスト{、257}}} \右\ }$ 。彼がすべてについて正しかったかどうか、 Wolfram言語 -PrimeQ関数を使用して確認できます。これは、数字を合成することを証明するために特定の要因を検索する必要がない、簡単にするために数字をテストする最新の方法を使用します。

彼の声明は

M_{67}

${M_ {67}}$ -素数、単なるタイプミス、そして実際には彼は

M_{61}

${M_ {61}}$ 。核分裂テストは利用可能な数少ないツールの1つであるため、メルセンヌの生涯の間に単純性のチェックが困難であったことは容易に理解できます。たとえば、

M_{257}

${M_ {257}}$ 最小のファクターは15桁の数値であるため、最新の因数分解法を使用してもそれほど簡単に見つけることはできません。 FactorInteger関数は、分解法を大きな整数に適用できる最も高度な方法を使用します。

オイラーおよびメルセンヌの乗数定理

シンプルテストの分野における最初の重要な成果は、偉大な数学者レナード・オイラーに属します。

M_{31}

${M_ {31}}$ 簡単です。彼は、単純な除数

M_{31}

${M_ {31}}$ 1または62（mod 248）に等しくなければなりません。

このような比較的短いリストは、オイラーの時点でさえ、（手動で）試行分割を使用して確認できます。彼はメルセンヌの定理の適用に属します。

q

$q$ 仕切りです

M_{p}

${M_p}$ それから

q e q u i v 1 v e e - 1 \左 （ b m o d 8 \右 ）

$q \ equiv 1 \ vee-1 \左（{\ bmod 8} \右）$ 、

q e q u i v 1 \左 （ b m o d p \右 ）

$q \ equiv 1 \左（{\ bmod p} \右）$ そして

q 2 k p + 1 $ 相

$q \ 2kp + 1 $相$ 正の整数の場合

k

$k$ 。これらの事実により、可能な除数の数が大幅に制限されます。

M_{p}

${M_p}$ 。以下の関数を使用して、この定理を使用して可能な除数のリストを提供する方法を示します

M_{p}

${M_p}$ より小さい

s q r t 2^{p} - 1

$\ sqrt {{2 ^ p}-1}$ 。

これらの機能を使用して、除数をすばやく見つけます。

2^{41} - 1

${2 ^ {41}}-1$ 。ご注意ください

q

$q$ 仕切りです

2^{p} - 1

${2 ^ {p}}-1$ 場合にのみ

2^{p} e q u i v 1 \左 （ b m o d q \右 ）

${2 ^ p} \ equiv 1 \左（{\ bmod q} \右）$ 。これにより、効果的な累乗法を提供するPowerMod関数を使用できるようになります。

以下は、161649文字のメルセンヌ数です。

ルークとレーマー

次の重要なステップは、 Eduard Lukeによるこの形式の数字の単純さをチェックする方法の発見でした。彼は1876年に彼の方法を使用して、

M_{127}

${M_ {127}}$ （「プレコンピューター」時代の最大のメルセンヌ数）は単純です。 20世紀の初め、二項算術と代数の基礎が広く知られるようになったとき、デレクヘンリーレマーはルークの方法を完成させました。結果として得られたLuc-Lemer数の単純性テストは、この形式の数が素数であるかどうかの効果的なチェックを提供しました。チェックは、モジュロの比較を使用して実行されました。

つまり、

k

$k$ によって表される数と同一

p

$p$ 下位ビットに加えて、残りのビットで表される数値。この比率は、次の限り再帰的に適用できます。

k < 2^{p} - 1

$k <{2 ^ p}-1$ 。

次の例を考えてみましょう。ここでそれを示します

のために

k = t e x t 1234567891

$k = {\ text {1234567891}}$ 。ご注意ください

（23下位ビット）および

-残りのビットが最下位にシフトされることを意味します。

以下の関数は、この計算方法を定義します

k b m o d \左 （ 2^{p} - 1 \右 ）

$k \ bmod \左（{{2 ^ p}-1} \右）$ ビット演算のみを使用します（除算なし）。ご注意ください

2^{n} - 1

${2 ^ n}-1$ バイナリ形式を持っています

111 . . . 111_{2}

${111 ... 111_2}$ 、単一の0はありません。したがって、

p

$p$ 数値の下位ビット

k

$k$ 。

次の関数は、Luc-Lemer単純性テスト（LLT）を実装します。シーケンスを定義する

s_{0} = 4

${s_0} = 4$ 、

s_{i} = s_{i - 1}^{2} - 2; i > 0

${s_i} = s_ {i-1} ^ 2-2; i> 0$ 。それから

M_{p} = 2^{p} - 1

${M_p} = {2 ^ p}-1$ 場合にのみ

s_{p - 1} e q u i v 0 \左 （ b m o d M_{p} \右 ）

${s_ {p-1}} \ equiv 0 \左（{\ bmod {M_p}} \右）$ 。

経験上、これらの関数にかかる時間のほとんどは整数演算に費やされることが示されています。

かどうかを確認するには

2^{p} - 1

${2 ^ p}-1$ シンプルな場合、最初に小さなシンプルな仕切りでシンプルさをテストし、他のシンプルさチェックを実行するのが最善です。まずcheckMPDivisorsでエンコードされたメルセンヌの素数の除数定理を使用し、次にPrimeQ関数を使用します。これが機能しない場合は、Luc-Lemerテストを適用します。

ここでは、 PrimeQ関数の拡張バージョンを紹介します。これは、多数のフォームに対してLuc-Lemerテストを適用します

2^{p} - 1

${2 ^ p}-1$ 。

から ${M_ {13}}$ 前に ${M_ {20}}$

Luc-Lemerテストを使用してコンピューターで見つかった最初のメルセンヌ素数は

M_{521}

${M_ {521}}$ 、ランプコンピューターSWAC（Standards Western Automatic Computer）に基づいて1952年1月30日にRafael Robinsonによって発見されました。以下に、それぞれ37ビットの256ワードを含むこのコンピューターのメモリブロックを示します。

20番目のメルセンヌ素数は、IBM 7090での50分間のLuc-Lemerテストの結果として、1961年11月にAlexander Hurwitzによって発見されました。

この種の作業に適したWolfram言語の機能の1つは、高速整数演算です。メルセンヌ素数の検索は、コンピューター化された検索の時代のd明期に大きな課題になりました。研究者は、高速フーリエ変換法を採用して、2つの大きな整数を乗算するタスクを、変換された数字の単純な要素単位の積に変換しました。整数の素早い乗算は、Luc-Lemerテストの2乗ステップに必要です。 Wolfram言語は数十億文字の正確な整数で動作するように最適化された最新のアルゴリズムを使用しています。たとえば、最後のものが

M_{4423}

${M_ {4423}}$ 、実際には素数のメルセンヌ数であり、そのすべての数を示しています。

完璧な数字

メルセンヌ素数と完全数の間には興味深い関係があります。完全数とは、そのすべての約数の合計に等しい数です（数自体を除く）。ユークリッドは仮定し、オイラーはすべての完全な数が次の形式であることを証明しました

P = 2^{p - 1} \左 （ 2^{p} - 1 \右 ） = 2^{p - 1} M_{p}

$P = {2 ^ {p-1}} \左（{{2 ^ p}-1} \右）= {2 ^ {p-1}} {M_p}$ 。 Wolfram言語のPerfectNumberQ関数は、数値が完全かどうかをチェックします。このプロパティのデモを行います

M_{31}

${M_ {31}}$ 。

21、22、および23メルセンヌ

21番、22番、および23番のメルセンヌ数の再発見に進みます（さらに、次の形式で示します）。

\＃ 21

$\＃21$ ）：

\＃ 21 = M_{9689}

$\＃21 = {M_ {9689}}$ 、

\＃ 22 = M_{9941}

$\＃22 = {M_ {9941}}$ 、

\＃ 23 = M_{11213}

$\＃23 = {M_ {11213}}$ 。それらはすべて、ドナルドギリスによって発見されました。ドナルドギリスは、1963年春にイリアックIIでLLTを開始しました（こちらを参照）。フォームのすべての番号を確認するには

2^{p} - 1

${2 ^ p}-1$ 間の素数の場合

7927 l e q s l a n t p l e q s l a n t 17389

$7927 \ leqslant p \ leqslant 17389$ 約6分かかりました。

24、25、26番目のメルセンヌ数

次に、検索を展開して見つけます

\＃ 24 = M_{19937}

$\＃24 = {M_ {19937}}$ 、

\＃ 25 = M_{21701}

$\＃25 = {M_ {21701}}$ 、

\＃ 26 = M_{23209}

$\＃26 = {M_ {23209}}$ 。これらの最後のものは1979年2月にランドン・カート・ノールとローラ・ニッケルによって発見されました。彼らは

M_{21001}

${M_ {21001}}$ 前に

M_{24499}

${M_ {24499}}$ CDC Cyber 174スーパーコンピューター上で（これについてはこちらをご覧ください）。計算が長くなっているため、並列処理を使用し始めています。テストは独立しているため、 ParallelMapを使用して作業を高速化できます。範囲チェック

17393 l e q s l a n t p l e q s l a n t $ 27, 44

$17393 \ leqslant p \ leqslant $ 27,44$ 約3分半かかります（4つのコアが使用されます）。

特化されたLuc-Lemerテストは、より一般的なPrimeQ関数（指定されたメルセンヌ数）よりもはるかに高速であることに注意してください。

27日と28日、メルセンヌ

次に、範囲を確認しました

27457 l e q s l a n t p l e q s l a n t 48611

$27457 \ leqslant p \ leqslant 48611$ 番号を見つける

\＃ 27 = M_{44497}

$\＃27 = {M_ {44497}}$ 。それは1979年4月にハリー・ネルソンと彼のチームによって発見されました（彼らはCray-1スーパーコンピューターを使用しました）。検索は15分で完了しました。

次のメルセンヌ数、

\＃ 28 = M_{86243}

$\＃28 = {M_ {86243}}$ 。 1982年9月にDavid Slowinskiによって-これもCray-1で開かれました。このスーパーコンピューターの重量は約5トンで、約115キロワットの電力を消費し、そのコンピューティングパフォーマンスは160 メガフロップに達しました。 100万個の64ビットメモリワード（8メガバイト）が付属し、現在の価格では約1,600,000ドルかかりました。以下は、その冷却システムの詳細です。比較のために： Raspberry Piの重さは数オンスで、4ワットで動作し、約410メガフロップスを提供し、1GBのRAMを搭載しています。これはすべて40ドルです。 Mathematicaにもすぐに付属しています。

数

\＃ 28 = M_{86243}

$\＃28 = {M_ {86243}}$ 25962桁が含まれます。この値は1時間14分で（Raspberry Piではなくラップトップで）検出され、範囲でテストされました

48619 l e q s l a n t p l e q s l a n t 87533

$48619 \ leqslant p \ leqslant 87533$ 。

第29回メルセンヌ

現在、より深刻なコンピューター時間が必要なので、実行ごとにタイムスタンプも生成します。次に範囲を確認します

87557 l e q s l a n t p l e q s l a n t 110597

$87557 \ leqslant p \ leqslant 110597$ 。 1時間44分後に、コンピューターは以下を示しました。

\＃ 29 = M_{110503}

$\＃29 = {M_ {110503}}$ 。この番号は、1988年1月29日にWalker ColquittとLuke Welsh（NEC DX-2スーパーコンピューター。記事はここにあります）によって最初に発見されました。

30日と31日のメルセンヌ

次の2つのメルセンヌ素数：

M_{132049}

${M_ {132049}}$ そして

M_{216091}

${M_ {216091}}$ 、-は実際に29日以前に発見されました（28日を発見した同じチームによって）。彼らはCray X-MPを使用して、1983年9月に30番目のメルセンヌ数を、1985年9月に31番目を見つけました。確認する

\＃ 30

$\＃30$ 範囲検索機能を使用する

110603 l e q s l a n t p l e q s l a n t 139901

$110603 \ leqslant p \ leqslant 139901$ 。それぞれを確認するには

M_{p}

${M_p}$ この範囲では、4時間8分かかりました。

グレートインターネットメルセンヌ数検索

34番目のメルセンヌ素数の発見により、

M_{1} 257787

${M_1257787}$ -1996年9月、メルセンヌ素数を検索するためのスーパーコンピューターの時代は終わりました。次の15個は、グレートインターネットメルセンヌシンプルナンバーサーチ（ GIMPS ）のボランティアによって発見されました。これは、パーソナルコンピュータで（バックグラウンドプロセスとして）Luc-Lemerテストのバリアントを実行します。この大規模な分散コンピューティングプロジェクトは、現在、1秒あたり約300テラフロップスに相当するパフォーマンスレベルを達成しており、アイドル時間は130万台を超えています。

Luc-Lemerテストを使用して、34番目のメルセンヌ数を確認しましょう。この段階で、パソコンの機能の限界に達しました。この範囲で数千のメルセンヌ数をテストするには、何日もかかります。 1つの大きな素数をテストするために必要な数十億の計算の中で1つの算術エラーが誤った結論につながり、損失を意味するため、Luc-Lemerテストはコンピューターハードウェアおよびソフトウェアの信頼性のストレステストとしてよく使用されることに注意してくださいメルセンヌ数または間違った答え。それぞれを確認したという事実

M_{p}

${M_p}$ 2から139901までの素数の場合、MathematicaシステムとWolfram言語での多数の算術演算と2項演算の信頼性を説得力をもって証明します。

メルセンヌ数の因数分解

すでに見てきたように、フォームの数の拡大で考えられる要因の数

2^{p} - 1

${2 ^ p}-1$ メルセンヌ乗数の定理によれば、それは有界です。これにより、この形式の大きな除数の効率的なコンピューター検索が可能になりました。この記事の執筆時点では、すべてのメルセンヌ数（

2^{1201} - 1

${2 ^ {1201}}-1$ ）完全に因数分解されました（図はここにあります）。 GIMPSプロジェクトは、単純な数字の検証の副産物として、多くのメルセンヌ数の大きな乗数の発見にもつながりました。この問題を解決するための最新のアプローチである新しい記事（および17の新しい因子分解）は、ここにあります。最大の因数分解された数は