Prime Curios!

「キュリオ」（珍品）って何？

2019年7月27日

「Prime Curios!」のホームページには、最初にこう書かれています。

「Prime Curios!」は、素数に関連する珍しいこと・不思議なこと・トリビアを集めた楽しいコレクションです。立ち止まって、野の花の匂いをかぐこと。珍しいコインを集めること。春の雷雨のとき巻き雲を眺めること…。私はこれまで、そういったことの価値が分からない大勢の人たちに出会いました。昔のことわざで「美は、それを見る人の目の中にある」といいます。私たちの目はどうでしょうか。ぜひ私たちのキュリオを試食してみてください。

コンテンツ・エディター: G. L. Honaker, Jr. / テクニカル・エディター: Chris Caldwell

私たちの目標は、面白い性質や形を持つ一つ一つの素数のコレクション（いわば辞書）を作ることです。あなたは何か、面白い数を知っていますか。他の人が「もっと話を聞きたい」と思うような調子で、あなたのキュリオを一般の人に分かるように説明できますか。もしそうなら、知らせてください。

https://primes.utm.edu/curios/

提出したキュリオたち

2019年7月25日

忘れないうちにメモしておきたい。

最初に提出したキュリオは「フェルマーのクリスマス定理で遊ばせて！」（2018年12月）で書いた「79の不思議」。何となく送ってみた。送信のときメアドを記入するので、採用されたらメールで通知でもくるのかな、と思ってた。実際には何の連絡もないし、ユーザー登録もパスワード認証も何もない。今どきのんきな世界だが、まあ他の素数愛好家になりすます詐欺師もいない（なりすますメリットもない）のだろう。知らないうちに（たぶん提出後、翌日くらいに）採用・掲載されていたらしい。記憶が曖昧だが、2019年1月に送って、2月になってそれに気付いたような。送信した方も、のんきだった。

2番目のキュリオは、4187 についての既存のキュリオと関係がある。後から分かった客観的事実としては「間違いに気付いたから連絡した」という形だが、その時点では「自分の側が何か勘違いしてる？」と疑問に思ってた。「4187は逆数の循環節の長さが素数になるような最小の強擬素数（自明な底を除く）」というネタを送り、コメントとして疑問について尋ねてみたが、反応がなかった。

2019-06-21　擬素数の逆数の循環節の長さ　n > 1 を 10 と互いに素な整数とする。1/n を 10 進小数で表し、その循環節の長さを ℓ とする。　(I) n が素数であるか、または底 10 のフェルマー擬素数であるときに限って、n ≡ 1 (mod ℓ) が成り立つ。　(II) 底 10 のフェルマー擬素数は、10 以上なら、全てこの合同式を満たす。

【例】 n = 91 とすると、1/n の循環節の長さ ℓ = 6 は n ≡ 1 (mod ℓ) を満たす。ところが n = 91 = 7 × 13 は素数でないから、フェルマー擬素数になるはず。実際、10^91 − 1 ≡ 1 (mod 91)。

【証】 10 進小数の場合、n を法として ℓ は 10 の位数。n が素数なら、小定理により ℓ は n − 1 の約数だから n − 1 ≡ 0 (mod ℓ) となり (I) の合同式が成立。よく知られているように、逆は不成立。すなわち、フェルマー擬素数も同様の理由から、同じ合同式を満たす。

これだけなら大したことじゃないが、これは話の出発点…

2019-06-24　最小の強擬素数　（正の）合成数 m が与えられたとき、m − 1 = 2^sd を満たす整数 s ≥ 0 と奇数 d ≥ 1 が一意的に定まる。次の条件のいずれかを満たす自然数 b を知りたい。

条件 (I)　　b^d ≡ 1 (mod m)

条件 (II)　　e = 2^rd が b^e ≡ −1 (mod m) を満たすような、整数 r ∈ [0, s−1] が存在

もし m が偶数なら、s = 0, d = m − 1 となり、(I) は普通の擬素数（フェルマー擬素数）と全く同じ条件、(II) は意味を持たない。従って m を奇数の合成数に限ることにする。このとき、(I) または (II) を満たす数（強擬素数と呼ばれる）は必ず擬素数だが、擬素数は必ずしも強擬素数ではない。普通の擬素数の場合と同様の理由から、b ≡ ±1 (mod m) は「自明な底」。底の検索を [2, m−2] の範囲に限定することができる。

そうすると、b ≥ 2 を固定して「b を底とする最小の強擬素数」を考えることができるのはもちろんだが、底を限定せずに「ある性質を持つ最小の強擬素数」を考えることもできる。例えば、100 以下の自然数のうち、非自明な底に対して強擬素数になり得るものは 25, 49, 65, 85, 91 の5個しかない（容易に全数検索可能）。「底の検索範囲を広げたら、350桁の巨大な（非自明な）底に対して 15 も強擬素数になることが発見された！」といったことは、絶対起きない。従って、例えば「49 は、逆数（の10進展開）が循環小数になる最小の強擬素数である」という命題は（数学的な価値はさておき）論理的に正しい。

「強」ではない普通の擬素数についても同様のことが言えるが、3 のべきを除く全ての奇数の合成数は（非自明な底の）擬素数なので、あまり面白みがない。偶数の擬素数（比較的珍しい）についてなら、何か面白いことが言えるかもしれない。

2019-06-23　どの底の擬素数？　（正の）合成数 m が与えられたとき、b^m−1 ≡ 1 (mod m) を満たす自然数 b を知りたい。一般性を失うことなく 1 ≤ b ≤ m と仮定できる。実際、任意の整数 x について、整数 k が存在して、上記範囲の b を使って x = km + b と書ける。このとき x ≡ b (mod m) であり、二項定理から x^m−1 = (km + b)^m−1 ≡ b^m−1 (mod m)。

従って「与えられた m が擬素数になるような底」をリストアップしたいとき、m より大きい底を考える必要はない。とりあえず m 回の試行で全数検索できる。さらに、合同式 b^m−1 ≡ 1 (mod m) は b = m ≡ 0 なら常に不成立、b = 1 なら常に成立。これらは「自明な底」である。

b = m − 1 つまり b ≡ −1 のとき、上記の合同式は m が奇数なら成立、偶数なら不成立（偶数 m = 2 に対しては成立するが、m は合成数なので 4 以上）。これも「自明な底」であり、結局、実質的に問題になるのは 2 ≤ b ≤ m − 2 の範囲の b に限られる。

【例】 合成数 m = 9 は、どのような底に対して擬素数か。これは「8 乗して 9 で割ったとき 1 余るような整数をリストアップせよ」という問題。b = 9 ≡ 0 が条件を満たさないこと、b = 1, 8 ≡ ±1 が条件を満たすことは明白。それ以外の底について、素朴に計算すると:

U = 2^8 = 256 ≡ 4
V = 3^8 = 6561 ≡ 0
W = 4^8 = 65536 ≡ 7
w = 5^8 = 390625 ≡ 7
v = 6^8 = 1679616 ≡ 0
u = 7^8 = 5764801 ≡ 4

実際には、この場合、底が偶数乗されるのだから、b > (m − 1)/2 については、m − b ≡ −b に対する結果を再利用できる（u = U, v = V, w = W）。結論として、9 は自明な底 1, 8（およびそれに 9 の倍数を足した自然数）について擬素数だが、9 が擬素数になるような非自明な底は存在しない。m が大きい場合、m − 1 乗をばか正直に計算すると計算途中の値が大きくなり過ぎるが、その点については抜け道がある。

2019-06-27　ラビン・ミラー擬素数の検索　底 2 ≤ b ≤ m − 2 の強擬素数である奇数 m を「ラビン・ミラー擬素数」と呼ぶことにする。ラビン・ミラー擬素数であるかないかだけを知りたいなら、2 ≤ b ≤ (m − 3)/2 の範囲の底を考えれば十分。なぜなら底 b のラビン・ミラー擬素数は底 m − b ≡ −b (mod m) のラビン・ミラー擬素数でもある。

実際、b^d ≡ ±1 (mod m) は b を −b に置き換えても成立。これは条件 (I) の全部と条件 (II) の一部であり、≡ −1 は、条件 (II) における r = 0 の場合に当たる。条件 (II) の残りの部分は、r ≥ 1 に対し b^e ≡ −1 が成立するというものだが、その場合、b を −b に置き換えても左辺の値は変わらない（e は偶数なので）。

【例1】　n = 121 なら d = 15。b = 3 に対して b^d ≡ 1 (mod n) であり、b を −b に置き換えれば、d は奇数なので右辺は ≡ −1。すなわち、底 b ≡ ±3 の両方に対して 121 はラビン・ミラー擬素数。

【例2】　n = 125 なら d = 31。b = 57, r = 1, e = 62 に対して b^e ≡ −1 (mod n) であり、b を −b に置き換えても、e は偶数なので結果は変わらない。

2019-06-29　4187という数　The Prime Pages の「素数の珍品」（Prime Curios!）コーナーには、ちょっと前まで次の記述があった。4187 is the smallest Rabin-Miller pseudoprime with an odd reciprocal period. [Post]。

「奇数の逆数周期を持つラビン・ミラー擬素数のうち、最小のものは4187」という曖昧な主張だが、「逆数周期」＝「逆数を10進展開したときの循環節の長さ」、「ラビン・ミラー擬素数」＝「奇数の合成数 q のうち、少なくとも一つの底 b ∈ [2, q − 2] について強擬素数となるもの」と解釈できる。

あるきっかけで、6月21日に上記の記述を目にした。「これは絶対間違ってる」という気がした。実際、全数検索によれば、上記命題には51個も反例がある！　上記の性質を持つ本当の最小のラビン・ミラー擬素数は185（底43など）、その逆数周期は奇数3。仮に逆数が純循環小数になるという制限を付けるとしても、最小は961（底229など）、その逆数周期は奇数465。

主張を正しいものにする一つの方法は「逆数周期が奇数」を「純循環かつ逆数周期が素数」に変えること。別の改訂案は「ラビン・ミラー擬素数」を「底10のラビン・ミラー擬素数」に変えること。後者は「任意の底→底10限定」なので命題が弱くなるが、この弱い形にはある種の理論的な興味がある。というのは、ある合同式を考えると「10進数としての逆数周期」は確率的な素数性テストに利用でき、そこにおいて「底10のフェルマー擬素数」は、やはり擬素数となる（2019-06-21　擬素数の逆数の循環節の長さ）。先週末、強弱両方の主張を新しいキュリオとして提案したが、何の反応もなかった。今週末、コメント送信フォームから上記の指摘を送ってみたら、今度は1時間ほどで返事があり、結局、このキュリオは「弱いバージョン」に改訂された。

現時点の本音としては、この弱いバージョンは、つまらない。これがキュリオ（珍品）なら、任意の底 n について、同様のもの（逆数を n 進展開した周期が奇数になる最小のラビン・ミラー擬素数）を考えることができ、それらの一つ一つがキュリオということになってしまう。「どんな底を選んでも、4187より小さいラビン・ミラー擬素数の逆数を素数周期の純循環小数にできない」という強い主張の方が、少なくとも見掛け上は面白い…。

誤りに気付いただけなのに、データベースのシステム上、筆者がこのキュリオの「新オーナー」になってしまった。これも何かの縁、こうなったら「4187のこの性質は、本当に面白いんだよ！」と言えるように、研究してみたい。上記のように、興味深い問題と結び付いていることは確かなので…。

去年まで Prime Curios! に参加すること自体、考えてもなかったけど、今年1月、思い立って「79の不思議」について初投稿（昨年末に「フェルマーのクリスマス定理で遊ばせて！」で書いたもの）。この79の性質は、本当に面白いし、数学的にも奥が深い。一意分解整域になる世界・ならない世界の分かれ目は何か…。答えがあることも、答えを与えてくれる理論の名前も分かっているが、短絡的に答えを見てしまうのも味気ない。もうちょっと山麓をさまよってみたい。

2019-07-01　汝のあるべき姿に戻れ、4187（良い花）！　更新された「素数の珍品たち」。偽から真に変わったのは良いとして、やっぱりしっくりこない…。

あるべき姿（弱いバージョン）　「4187は、逆数を小数で書いたときに循環周期が奇数になるような、10より大きい最小の擬素数（10を底とする）である」 ← この形なら、意味を説明しやすい（下記参照）。

あるべき姿（強いバージョン）　「4187は、逆数の周期が素数になるような、10と互いに素な、最小のラビン・ミラー擬素数である」 ← 今回は立ち入らない。

現在の姿（中途半端）　「4187は、底10において、逆数の周期が奇数になるような、最小のラビン・ミラー擬素数である」 ← 弱いバージョンと同値だが、無駄に難解。「ラビン・ミラー」という条件は冗長。

説明。「10 を n−1 乗して n で割ったら 1 余る」という性質を持つ自然数 n を考える。式で書くと:

10ⁿ⁻¹ ≡ 1 (mod n)

n が「底の10と互いに素な素数」なら（要するに2, 5以外の素数なら）、この関係は必ず成り立つ。例えば7は素数なので、10の6乗を7で割ると1余る。13も素数なので、10の12乗を13で割ると1余る。n が素数でないときは、一般には、この関係は成り立たない。例えば21は素数でない。そして、10の20乗を21で割った余りは1にならない。

ってことは、この関係を満たすか・満たさないかで、素数か素数でないか、簡易的なテストができる。2と5を例外として、この関係を満たさない数は、素数ではない。それはいいんだけど、時々いるんだよね…合成数のくせに、テストを通ってしまう「ニセ素数」が。例えば、9は合成数なのに、10の8乗を9で割ると1余る。ニセ素数。正式には、底10の擬素数と呼ばれる。もっとも 10 ≡ 1 (mod 9) だから、10の何乗を9で割っても1余るのは自明だが…。自明な例も含めて、底10の擬素数を小さい順に列挙すると:

9, 33, 91, 99, 259, 451, 481, …

これらの数（最初の9を除く）の逆数を小数で書くと、1/33 = 0.030303… は循環の周期が2、1/91 = 0.010989 010989… は周期が6、以下順に周期が 2, 6, 10, 6 となる。「あるべき姿（弱）」の意味は「4187も同じ底10の擬素数だが、これは逆数の周期が奇数。ちょっぴり珍しいよ」。丁寧に説明すれば、誰でも理解できそう。

そんなシンプルな話なのに、なぜわざわざ「ラビン・ミラー」という難しいことを言って話を複雑化するのだろう。そもそも「ラビン・ミラー擬素数」とは何なのか。少し考えると「12以上の奇数のうち、底10の強擬素数」という意味になることが分かる。けど、この条件のうち本当に必要なのは「12以上の擬素数」だけ。「奇数」も「強」も余計。「12以上」が必要というのも、その心は、自明な擬素数9を排除したいだけ。9を排除しないと、その逆数 1/9 = 0.111… の循環周期が1（奇数）なので「奇数の周期は珍しい」という話の前提がいきなり崩壊するから…。だったら素直に、「9を除く」とか「10より大きい」って言えばいいじゃん！

「底10の擬素数」というシンプルな条件なら、120万以下の擬素数は300個あって、そのうち「逆数の周期が奇数」は15個だけ（9を除外すれば14個）。まあまあレア。現在の「ラビン・ミラー擬素数」という重い条件を付けると、同じ範囲にラビン・ミラー擬素数は68個しかなく、そのうち14個が「奇数周期」。同じ14個だが、68個中の14個というのは別に珍奇ではない。

条件が必要十分でなく過剰。論理的に間違ってるわけではないが、しっくりこないし、レア感も低下してる。

2019-07-01　追記　↓エディターの先生が「こんな感じならいい？」と再編集してくれました♪　多分プロの数学者から見ると、論理的に同値＝表面的な言い方の問題にすぎないのでしょうが…。もともとの命題（Post、2009年、偽）　4187 is the smallest Rabin-Miller pseudoprime with an odd reciprocal period. → 1回目の更新（2019年6月。真になった）　4187 is the smallest Rabin-Miller pseudoprime in base 10 with an odd reciprocal period. （正確には覚えてないがこんな感じだった） → 今回の更新（2019年7月）　4187 is the smallest pseudoprime > 10 in base-10 with an odd reciprocal period. （この方がシンプルでいいかと）

2019-07-04　強擬素数よりもっと強い（だいたいこんな感じ？）

q を底 b = 10 の擬素数とする。q − 1 = 2^sd と書く。ここで s ≥ 0 は整数、d は奇数。

q − 1 を 2⁰, 2¹, 2², … で割った商 2^sd, 2^s−1d, 2^s−2d, … のうち、b の肩に乗せると ≡ 1 (mod q) となる最小の自然数を e = 2^td, t ∈ [0, s] とする。q が擬素数という仮定から、必ずこのような e が存在（強擬素数でなければ t = s）。

mod q において b^e ≡ 1 なので、b の位数は e の約数。

t = 0 なら e は奇数なので、b の位数は奇数。この条件を満たす q を第1種擬素数と呼ぶことにしよう（強擬素数のうちの、特別な場合に当たる）。一方、t > 0 なら e は偶数であり、しかも b^e/2 ≢ 1 なので、b の位数は、その因子として少なくとも1個の 2 を含み、奇数ではない。

1/q を b 進小数で書いたとき、1/q の循環節の長さ（言い換えれば b mod q の位数）が奇数ということは、q が第1種擬素数ということと同値。だから 4187 の性質について「強擬素数」と断るのは冗長で、ただの「擬素数」の範囲で考えても同じこと。

自明な擬素数 9 も同じ性質を持っていて、それを除外しないと、この文脈では 4187 は輝いてくれない。擬素数の底が b = 10 なのも、地球人がたまたま 10 進法を使っているというだけで、数学的必然性に欠ける。それらの点を別にすれば、これは面白いパズルかもしれない。

3番目のキュリオは、15319333 に関するもの。33進数では「prime/2」がこの素数になる、というネタだが、これはボツになった。同時に、4187 について、問い合わせフォームから質問の形でコメントを送信したところ（自分の中の疑問を解決したかった）、こちらには反応があり、結局、思った通り、掲載されていた内容に不備があったと判明。普通、数学の命題は真偽がハッキリするものだが、既存のキュリオには定義が曖昧な用語が含まれていて、そのせいでクリアカットに判断できなかった。意外な展開として、もともとあった 4187 のキュリオが自分のキュリオという扱いになった。

4番目のキュリオは、65281 に関するもの。4187 の「16進数バージョン」のようなもので「逆数の循環周期が偶数になるような最小の底16強擬素数。120万以下に他にこのような数はない」。これは採用された。

2019-07-07　キュリオ3個目（1個採用・1個ボツ？）

(1) 底16の強擬素数の逆数を16進小数で書くと、ほとんどの場合、周期（循環節の長さ）が奇数になる。65281 = 0xff01 は最小の反例を与える。このような数は、120万以下の範囲には他にない。（今週末に提出、2日くらいで採用された。）

1/ff01 = 0.000100fffeff 000100fffeff... ← 4種の数字だけで書ける

(2) 33進数では (prime)₃₃/2 は素数。(prime)₃₃ − (curio)₃₃ も素数。こちらは2週間くらい前に提出したが未編集のまま。ボツになったらしい。

単純に「自分が面白い」と思うものが見つかったら、またネタを送ってみたい。たとえボツになっても、考えを深めるきっかけになるし、それ自体が楽しいので…。

「The Prime Pages」は素数関係では国際的に有名なサイト（テネシー大学）。既知の巨大素数のランキング、素数に関する広範な用語解説を含む。「Prime Curios!」は、一般愛好家向け（マニアックな素数ハンターに限らない）コーナー。面白いトリビアを発見したら、誰でも投稿可（匿名でもOK、登録不要） → https://primes.utm.edu/curios/index.php　ちなみに 4187 について間違った予想を投稿した Jonathan Vos Post はソフトウェア関係の大物で、採用数3位（500個以上）のキュリオ貢献者らしい。

5番目のキュリオは、656601 に関するもの。三つ子素数の内側にある最小のカーマイケル数。既に10年くらい前にリベラが「三つ子素数の内側にある底2擬素数」と記述しているのであまり新規性がないのだが、カーマイケルであることを含めて、独立に見つけた。

2019-07-11　トリプル・サンドイッチ・カーマイケル

連続する4個の奇数が、全て素数またはカーマイケル数になるケースがある。最初の例は p = 656597 から始まる。ウェブ上を検索した限りでは、このことはまだ明示的には言及されていないようだが、次のことを組み合わせると、同じ結論に達する。

(1) p = 641 から始まる三つ子素数は (p, p+2, p+6) 型で、かつ p+4 が底2の擬素数。これは、このパターンの最初の例であり、2番目・3番目の例も容易に見つかる（p = 656597, 6212357）。Rivera は2002年に既に同じ発見をしていた（Puzzle 170. Pseudoprimemania）。2番目の4つ組に含まれているのは、ただの擬素数ではなく実はカーマイケル。

底2の擬素数を包み込む三つ子素数の最初の数をリベラ素数と名付けよう。

(2) 2017年、Are there further "sandwich"-Carmichael-numbers? において Peter は、カーマイケル数 n = 656601 について n−2, n+2 が素数であることに注目した（素数にサンドイッチされたカーマイケル数）。実はこれ、ただのサンドイッチではなく n−4 も素数（下側のパンが2枚重ねのトリプル・サンド）。それが2番目のリベラ素数というわけ。

2～4番目のカーマイケル・トリプル・サンドは恐らくいずれも17桁で「下側が2重」タイプ。ただし「小さい順に何番目か」という部分は、Pinch のテーブル（Pseudoprimes and Carmichael numbers）に誤りがないことに依存する。

リベラ自身が（2002年のマシンでは）見つけられなかった4番目のリベラ素数は、p = 18958567877。これについては、19*10^9 以下の三つ子素数を全数検索して確認した。

6番目のキュリオは、113 に関するもの。1/2ポンド、1/3ポンド、1/4ポンドがグラムでいうと（端数は四捨五入）どれも素数になるという観察に基づく軽いネタ。これは前から知っていた純粋な発見で、ずっと「そのうち Curios! に送ろうかな」と思ってたもの。

2019-07-13　クォーターパウンダーというサンドイッチ（ハンバーガー）

パズル（一般向け）　100 より大きい3個の素数 q, r, s を見つけてください。ただし 4q, 3r, 2s は連続する3個の整数（例えば 205, 206, 207）です。

どうでもいいような問題だが、面白いのは最小の解 q, r, s がそれぞれ、1/4ポンド、1/3ポンド、1/2ポンドをグラムで表した値（端数は四捨五入）になること。

このパズルを思い付いたきっかけは「クォーターパウンダー」というサンドイッチ（バーガー）。正直、自分でも他人が面白いと感じるか確信がなかったけど、試しに Curios! に送ってみたら、これが通った。

「クォーター・パウンド」「ワンサード・パウンド」「ハーフ・パウンド」のバーガーは、全部、実在する商品。それらのパティの重さは、グラムでいうと全部素数…というわけ。1/4ポンド、1/2ポンドをパッとグラムに換算できると、実用的にも意外と役立つかも？

7番目のキュリオは、13741。底2擬素数のギャップが6の場所を見つけたので「双子擬素数」に倣って「セクシー擬素数」と名付けた。

2019-07-13　セクシー擬素数

「これはやり過ぎかな」とも思ったが「セクシー擬素数」を Curios! に送ったら、あっさり通った。

「双子素数」「双子擬素数」を拡張して、「セクシー素数」があるんだから「セクシー擬素数」を考えてもいいじゃん、という素朴な発想。(13741, 13747) が、底2における最小のセクシー擬素数。

実は2番目の「セクシー擬素数」ペアは「69」から始まる数。そのことについて一言ジョークを入れても良かったのだが、それはさすがに下品なので思いとどまった。個人のブログなら無問題でも、学校のサイトなので…。「セクシー擬素数」が通ることが分かっていれば、ポーカーフェイスで「2番目の例は、69から始まる7桁の数である」とコメントしたんだけどね！

8番目のキュリオは、4369。「16進レプユニットになる最小の底2擬素数」という観察（0x1111）、その形の2番目は 0x1111111 であるという指摘、それらを結合すると3番目（16進数で1を11個並べた数）になる、という観察。このきれいなパターンには、ちゃんと理由があるのだろうけど、いかにもキュリオっぽい。

9番目のネタは、561825265 に関するもの。「561 も 825265 も底2擬素数、それをつなげた 561825265 も底2擬素数」。これも珍妙だが、こちらはボツになってしまった。掲載基準もけっこう気まぐれ？　桁数が多いと審査が厳しいのかもしれない。調子に乗って「561 も 825265 も20進数では回文数」などと書いたのも、くど過ぎたかな？

10番目のキュリオ。第1・第2千年紀の合計日数は 730487日で素数…というもの。これは通った。純粋数学の世界では、暦・歴史が絡むネタは新鮮だったのかな？

11番目のキュリオは 595039。この素数は16進数 0x595039 = 5853241 と解釈しても素数。それをさらに 0x5853241 = 92615233 と変換してもまだ素数。同様に繰り返して、長さ7の連鎖を作れる最小の数。面白いかどうかは微妙。投稿時、長さ8の連鎖も複数得ていたが、あまり大きい数だと無味乾燥なので、長さ7の方を提出。

12番目のキュリオは 24046。{24043, 24046, 24049} の3個は、どれも底3概素数（本物のセクシー素数のペアの真ん中に、偶数の底3擬素数がある）。このパターンの底3擬素数は 10^10 未満に3個しかなく、しかも2個目は1個目（24046）の倍数。「素数ページで偶数のネタは難しいかな」と思いつつ、面白いからいいや、と提出。

13番目のキュリオは q = 9890881。q と 2q+1 がどちらも底2擬素数になること、そのような数は珍しいこと、4進数のレプユニットはしばしばこのパターンになることを書いた。どれも既によく知られている事柄だったが、Feitsma’s table を使って独立に再発見。「ソフィー・ジェルマン擬素数」と名付け、擬素数の長さ3の Cunningham chain は可能だろうか、と尋ねた。自然な問いなので、これもきっと誰かが既に考えてるんだろうけど…。送信のとき base-4 repunit がNGワード扱いされ当惑した（4th base が俗語で「セックス」だからと思われ。repunit in radix 4 と言い換えて投稿）。

14番目のキュリオは素数 50497 に関するトリビアで、再び暦ネタ。毎月の「2日」「3日」「5日」のような「素数の日付」をカウントすると、400年ごとにこの素数になる。

最初のころは、純粋に「たまたま見つけた面白いネタ」「数学的に自然な問題」を散発的に送っていたが、だんだん「ネタのためのネタ」的要素が混じるようになってしまった。「面白いことがあったら Curios! に提出する」から「Curios! に提出するためにネタを探す」にシフト。ゲームのような面白さもあり、つい夢中になって、2週間くらいの間「週7個まで」の制限ぎりぎり近く投稿してしまった。でも、表面的な現象の不思議さ・物珍しさを次々と追い求め、その背後にある数学を考えないのは浅はかで本末転倒。新しいことを考え学ぶきっかけにはなるのだが、「きっかけ」だけで終わらないように、これからは（既に送った上記のキュリオたちも含めて）現象の意味を掘り下げて考えるようにしたい。

Curios! では採用数ランキングが表示されるので「よし、世界ランク100位に入るまで送ろう」などと思ってしまうが、それは意味のない虚栄心。「次の目標は50位だ・20位だ」などと考えてしまうと、きりがない。巨大素数探しの競争に参加する気になれないのも、同じ理由。

「競争に勝つため」みたいな、目的のために努力するのではなく、それ自体が楽しいから夢中になる、ってのが一番いい。外部的な目的・動機付けなんて要らない。教科書を読んできちんと勉強すれば速いことだって、そうしたければ、自分のやり方で試行錯誤を楽しみたい。非効率だっていい。それ自体が目的だから。結果ではなくそれ自体が幸せだから。

世の中の進歩・学問の発展だって、大きなブレークスルーは「命令されてやる人・義務だから働く人・決められたことを忠実にやる犬」からではなく「それ自体を楽しむ人・気まぐれに行動する猫」から生まれるんじゃないだろうか？

資料を調べたければ（教科書を読みたければ）、そのときはそうすればいい。「今はこうしよう」って感じたら、何でもそうすればいい。「こんなことしてる場合じゃない」って分かっていながら、逃避的にゲームに熱中する…みたいなのとは違うよ。後ろめたい逃避的な遊びでは、純粋に楽しめないからね！

2019-07-27　15番目のキュリオとして、62481801147341（テトラコンタンの異性体の数）を出そうとしたら、2018年9月に既に出されていた。数論の人は純粋志向で、有機化学なんて興味ないだろう、というイメージがあるが、そうでもないのね…。せっかくなのでついでに書いておくと、5より大きい素数になるのは、この C40 に続いて C112 と C149。素数愛好家は、この3個目に反応すると思われ。M149 は有名なメルセンヌ素数なので…。これらは純粋な理論上の数で、物理的な意味で、本当にこの数の異性体が存在し得るわけではない（Isomer Count）。

2021-01-27 （参考）上記の件について、日本語のウェブページを発見しました。
https://web.archive.org/web/20180307122946/http://www.org-chem.org/yuuki/mow/0505/crowd.html 「不可能炭化水素」

2019-07-30　堀井憲一郎という人は、読書記録を付けていたがそれをやめて、その理由についてこう言っているそうだ（*）。「本を読みたいんじゃなくて『読んだ本の記録数を伸ばしたくて読む』という気持ちが出てきたから」「数字を数えだすと、数字のほうが大事になって、本体はどうでもよくなってしまう」　上記の現象も同様で、知らないうちに「掲載数を増やしたい気持ち」が出てきて、「純粋に数学を考える気持ち」に干渉するようになった。

ともかく最初は2019年1月ごろ、気まぐれに1回投稿しただけなのに、7月上旬に編集者とメールでやりとりしたのがきっかけとなって、7月には記録上・合計13個も送ってしまった…。計3個くらいボツになったが、11個掲載。10個掲載でランク90位台になって「幼稚な目標」も達成されたのに、そこでやめず、さらにネタを考えてしまう。これはちょっとやばいし（addictive）、あるべき状態ではないと感じた。

本当にそのこと自体を楽しんでいるか？　それとも、それをやったことを人に見てもらいたい・評価してもらいたいだけなのか？　昔、1週間かけてヘキサデカン（C16）の異性体10359個を手動で数えたことがあるが、誰かに自慢するためにやったわけではないし、実際それをやったこと自体、誰にも言わなかった。多分そういうスタンスが、自分には一番合ってる…。「面白いことが見つかったら共有するけど、それは副産物」みたいな。

その後のキュリオたち

2019年10月31日

(#15) 2019年8月。553 は、固定された 1 < b < q−1 に対して q, 2q−1, 4q−3 がいずれも b-擬素数になるような最小。#13 で自問した擬素数チェーンの一種（弱い例だが、もっと強い例を見つけるのは困難）。

(#16) 2019年9月。p = 53176534057553 は sqrt(k)*p, k=1..15 の整数部分が素数または半素数になるような唯一の素数 < 10^11 である。

(#17) 2019年10月。アルファベット順でフランス語最後の素数は
20*10^120 + 20*10^66 + 23*10^18 + 23*10^6 + 23*10^3 + 349 =
20000000000000000000000000000000000000000000000000000020000000000000000000000000000000000000000000000023000000000023023349
である。もともとリベラが考えたネタだが、リベラの結果を改善し、採用された。整数をあえて分数で表現すれば、アルファベット順でさらに後ろにできるかもしれないが、それは数学というより「言葉遊び」だろう。

2019-10-22　「最後の素数」（その1）

素数に「本当の最後」はないこと（いくらでも大きな素数があること）は、よく知られている。けれど例えば「無量大数」までの数詞を普通に使うという条件で、五十音順で日本語最後の素数を考えることはできる。Rivera は数日前、フランス語について、同様のことを投稿した。いわく:

「S = 23000000000023023661 は、アルファベット順において、フランス語の最後の素数である」

上記の主張は、いろいろな意味で正しくない。ABC順でこれより後ろの素数はたくさんある！　数学パズルや素数が好きな方は、挑戦してみてほしい。

〔フランス語の数詞の参考資料〕1000まで1個ずつ書いてあるページ。数値を単語に変換してくれるサイト。

★考えても分からない方は → ヒント

注: 現在（1961年以降）のフランスでは、公式にはロングスケールが使われる（billion=1兆, trillion＝100京）。実際の用例上は、必ずしもそういう言い方にはならないのだが、このパズルでは公式ルールを使う。英語風のショートスケール（billion=10億, trillion＝1兆）ではない。英語と違い、フランス語の million, billion などは、2倍以上だと末尾に -s が付く。cent, mille は英語の hundred, thousand と同じで、2倍以上でも -s が付かない（※）。

（※）例外として、下3桁が端数のない 200, 300, …, 900 のときは cents が使われる（200=deux cents, 201=deux cent et un）。80=quatre-vingts, 81=quatre-vingt-un における vingt と同様。ここでは素数を考えているので、例外ケース（100の倍数）は無関係。

★「最後の素数」（その1）のヒント: 下記のように、リベラの数 S では ????? の部分に six がある。辞書順で six より後ろの単語が ????? に来るようにできれば、アルファベット順で S より後ろになる例が一つ見つかったことになる。

S: Vingt-trois trillions vingt-trois millions vingt-trois mille six cent soixante et un
T: Vingt-trois trillions vingt-trois millions vingt-trois mille ?????

2019-11-01　「最後の素数」（その2）　アン・ドゥー・トロワ、おフランスざーます！

その1の続き。フランス語なんて縁がない人でも「アン・ドゥー・トロワ」くらいは聞いたことあるでしょう。ABC順で trois の方が six より後なので、????? の部分が trois やそれに似た言葉になれば、Sより後ろの素数の一例となる。

S: 23000000000023023661
Vingt-trois trillions vingt-trois millions vingt-trois mille six cent soixante et un
T: 23000000000023023???
 Vingt-trois trillions vingt-trois millions vingt-trois mille ?????

従って、上記の形の数のうち、下3桁が「3」か「30」台か「300」台になる素数が一つでもあれば十分だが、最初の2種類の可能性は、すぐ否定されてしまう。

23000000000023023003 = 3で割り切れる
23000000000023023031 = 409で割り切れる
23000000000023023033 = 3で割り切れる
23000000000023023035 = 5で割り切れる
23000000000023023037 = 19で割り切れる
23000000000023023039 = 3で割り切れる

そこで下3桁が「300」台を考えると、1個だけ「343」が素数となる。

S: 23000000000023023661
Vingt-trois trillions vingt-trois millions vingt-trois mille SIX cent soixante et un
T: 23000000000023023343
 Vingt-trois trillions vingt-trois millions vingt-trois mille TROIS cent quarante-trois

これで「Sはアルファベット順で、フランス語最後の素数」という主張を崩すことができた。この他、もし仮に下3桁を「001」「013」または「020」台または「600」台のある種の数にできれば S より後ろになるが、そのような素数は「647」（Six cent quarante-sept）だけであり、これは「661」（Six cent soixante et un）よりABC順で前なので、この桁数では T が唯一の反例。

注: フランス語の正確なスペルを覚えている必要はない。その1でリンクした参考サイトからコピペすれば十分。「フランス語」という表面に惑わされず「アルゴリズムの問題」と割り切ろう。

しかしまだ問題が解決したわけではない。実は T も「本当の最後」ではなく、アルファベット順でもっと後ろの素数を作る方法がある。その方法は「トリック」ではなく、あくまで「普通の数詞」の範囲だが、「京・垓・𥝱」の類いの非日常的表現（大きな数の名前）が絡んでくる。「ばびっと数え歌　でかい数編」がヒントになるかもしれない。

最終的には、あるタイプのばかでかい数の下3桁を調整して、なるべくABC順で後ろの素数を見つけることに帰着される。一般の整数について、素数性の判定は簡単ではない。上記の S や T にしても「たったの trillion のオーダー」とはいえ、ロングスケールでは2の64乗を超える。素朴な試行除算だけでは厳しい。

ここでは取りあえず PARI の isprime で素数性を証明した。24桁以内の数については「最初の13個の素数を底に SPRP なら、本物の素数」ということが証明されているので、それを自力でやってもいいのだが、まあ PARI を使うのが手っ取り早いでしょう。

2019-11-02　「最後の素数」（その3）　巨獣召喚

「メキニ・メキニ・ヌダラダラ～　いでよっ、アンデシリョン！」

じゃじゃーん　Undécillion　「ぎゃおーん！」

1000000000000000000000000000000000000000000000000000000000000000000

テロップ「巨大怪獣アンデシリョン　10の66乗（100分の1無量大数）」

S や T の頭にある trillions を undécillions に置き換えれば、それだけで「もっと後ろ」の素数が非常に多数考えられる（アルファベット順で u は t より後ろなので）。S を見たとき最初に考えた抜け道は、これだった。だが安心するのはまだ早い。

ラスボス「ワーハッハッ、66乗ごときで巨大とは笑止千万・無知蒙昧。唸れ！　衝撃の！　ヴィジャンティリョン！」

ぞぞっ・ぞぞっ・ぞぞぞっ　Vigintillion　「うをおーん！」

1000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000

テロップ「檄！帝国・最終巨獣ヴィジャンティリョン　10の120乗」

ラスボス「見ろ！　無量大数がゴミのようだ」

ABC順で v は u よりさらに後ろ。ABC順でこれより後ろの単位はない。これら2種類の単位を使って「フランス語最後の素数」が構成される。V1 より V2 の方が後ろであることに注意すれば、この先はほぼ一本道だろう。

V1: Vingt-trois vigintillions 23*10^120
V2: Vingt vigintillions       20*10^120

2019-11-06　「最後の素数」（その4）　10の120乗まで、キッチリ回せ！！

ここまでの話から、
p = 20*10^120 + 20*10^66 + 23*10^18 + 23*10^6 + 23*10^3 + x
0 ≤ x ≤ 1000
の形の素数が存在すれば、その中にフランス語最後の素数（ABC順）があることが分かる。だが100桁を超える巨大な数で「下3桁しかいじれない」というのは、相対的な大きさで言えば「針の穴のような細い隙間で素数を探す」ようなもの。そんな都合のいい素数がそもそもあるのか。あるとしても、どうやって見つければいいのか？

実は PARI の nextprime を使えば、候補はすぐに見つかる。該当するのは x = 349 または 657 または 727 の三つ。いずれの x に対しても p が本物の素数になることは、別途証明される。349 は trois から始まる言葉だが、残り二つは s から始まる言葉。だから、x = 349 のケースが「ABC順で最後」。

million, billion, trillion, ... と上がっていっても、ABC順で trillion より先には行けないように思える。常識の範囲ではそうだが、諦めずに11番目まで上がると undécillion が現れて「Tの束縛」が破れ、さらに20番目まで上がると vigintillion が現れて「Uの束縛」も破れる。もっと上の数詞は考えられるが、W より後ろの文字から始まる数詞はない。

相対的には「針の穴」でも、調べてみると、案外、素数が詰まっていた。100以下の素数は25個（4%）、1000以下の素数は168個（2%未満）…という減り方（どんどん密度が薄くなる）から素朴に考えると、122桁の数に1000分の3（0.3%）も素数が含まれているのは、ちょっと意外な気もする。エラトステネスのふるいによれば、数が大きくなればなるほど、何重にも何重にも「ふるいの弾幕」が張られ「一度も弾に当たらない（何の倍数にもなっていない）数」はまれになるはず…。もっとも素数が薄くなると「新しい弾幕の砲台」も減る。2の倍数・3の倍数・5の倍数・7の倍数…といった「古い弾幕」はどこまでも有効だが、「新しい弾幕」がだんだん薄まる結果、意外と「生き残る素数」も出てくる。

具体的に、
pi(x) = x/log(x)
pi(10^121) / 10^121 ≈ 0.003589
なので、計算上、122桁の数では整数1000個につき素数が3～4個ある。

素数マニアはよくメルセンヌ素数の類いを持ち出して「とんでもない桁数の素数が見つかりました。世界記録更新！」と大騒ぎする。2019年現在、既知の最大素数は約2500万桁、2番目は2300万桁。一見「こんなに大きな素数は非常に珍しい」かのような印象を与える。ところが上記と同様に考えると、同じ桁数の素数はもっともっといっぱいある（大ざっぱに、整数1億個につき1個以上）。巨大素数の新記録はソフト・ハード両面の進歩、工夫と情熱の産物であり、人類の進歩としては素晴らしいことなのだが、人類を外から見ると「2000万桁程度の素数は何億個も何兆個も何京個もあるのに、そのうち1～2個しか見つけられない」というのは自慢にならない。「おれの素数の方がでかいぞ。1000万桁を突破したぞ！」と虚勢を張っても、実際にはたった1000桁の因数分解ができない…。

(#18) 2019年10月。13323377777 は、英語でスペルアウトしたとき100文字を超える最小の素数である。ただし、全て (resp. 最後) の hundred の後ろに and を挿入するなら、1117373377 (resp. 3323373377) がそれに当たる。「数字をスペルアウトする」系のネタで #17 からつながる。採用された。

2019-10-23　10億102万4867

1001024867 は、普通に英語でスペルアウトしたとき、アルファベットの20種類の文字を必要とする最小の素数。「これは面白いものが見つかった！」と思いつつ、検索してみると…。またしてもリベラに、5年くらい先を越されていた（Puzzle 744）！　さすがはカーロス・リベラ。さすがは無冠の帝王（笑）。リベラは、この手のネタにめちゃくちゃ強いだけでなく、有名なパズルサイトを運営していて、世界中からフィードバックを得ている。ちなみに、半年くらい前「あしたのジョーの登場人物と同姓同名ですよ」とメールに書いてみたが、その点については何の反応もなかった。数学遊び一筋で、漫画には興味がないのかもしれない。

リンク先のページでは、アルファベット22種（アクセント記号の有無は区別しない）を使うスペイン語の素数として、10^120 + 10^42 + 10^36 + 10^12 + Cn の形を扱っている。読者が C6 = 5410587 を報告、リベラはさらに小さな C2 = 3710587 を得ているが、これは本当の最小 C1 ではない。C1 の値は昨日リベラ自身が投稿し、筆者も独立に同じ値を得た。小さい順で C2 は2位、C6 は6位になるようだ。パズル好きの方は、C1 を探してみよう！

これらは、あくまで数学パズル。語学はあまり関係ない。スペイン語なんて「ウノ・ドス」くらいしか知らなかったが、数字をスペルアウトするアルゴリズムは資料を見て適当に考えた。「最後の素数」も同様で、フランス語はあんまり関係ない（たとえフランス語の達人でも、それだけでは解きようがない）。問題の核心は「100桁くらいの数の素数性をどう判定するか」。この部分は難しいが（もちろん試行除算だけでは無理）、やればできる範囲だろう。

（あくまで遊び。何かの役に立つわけではない。入試に役立つみたいな世俗的メリットは、何もない。）

リベラの下記のポカは、本人の計算ミスではなく、読者からのフィードバックをノーチェックでそのまま転載したのが原因だろう。少しでも疑う気持ちがあれば、簡単に反例が見つかる。でも肯定的に考えると、リベラがこれを投稿したから、話が広がった！

(#19) 2019年10月。197 について。オイラーの多項式を別にすると |8*k^2 + 8*k − 197| は、この形の最良の素数生成多項式である。k=0..30 に対して、別々の31種類の素数を生成する。絶対面白いはずだが通らなかった。既知と判断されたのだと思われる。このネタも、実はリベラの 2426256797 についてのネタ（数学的に正しくない）について考える過程で発見。

2019-10-25　連続31素数を生成する2次式

8k² + 8k − 197 の絶対値を考えると、それは k = 0, 1, …, 30 に対して相異なる素数。連続40素数を生成するオイラーの式 k² + k + 41 にはかなわないが、それに近いものである。

オイラーの式を別にすると、この型（1次と2次の係数が同じ2次式）の明示的な世界記録は 6k² + 6k + 31 の連続29素数（A060844）だが、こっちは連続31素数。ある意味、世界記録を塗り替えた！

素数生成多項式は検索が盛んであり、8k² + 8k − 197 のような単純なものは、当然とっくに知られているはず。同じ形の式はなぜか見つからないが、本質的に同じ 8k² − 488k + 7243 が2005年に報告され、A272160 にもなっている。冒頭の形の方がシンプルで良いのに、なぜそんな書き方をするのか。想像だが「k = 0, 1, …, 60 の連続61入力に対して素数になる」と言いたいために、複雑な形にしているのでは…。実際に生成されるのは31種であり、同じ素数を2回ずつ生成して見掛けの数を水増しするのは、美しくない。「1次と2次の係数が同じ（小さい自然数）で、定数項が小さい素数」と決めれば（※）、全数検索は容易であり、従ってランキングを確定できる。その結果、連続29素数を生成する別パターンの式も見つかった。

「素数生成は研究され尽くして、遊びで新規参入する余地などない」という気がするが、アイデア次第で意外と隙間もあるのかもしれない。例えば、連続○素数という数は小さくても、それら○個の素数の全て（あるいはほとんど）が何か特別な性質を持つ…というような珍しいパターンを探すと面白いかも。

（※）この前提自体、必ずしも自然ではないが、定数項が正の場合、そうすれば簡単な漸化式でも書ける。オイラーの k² + k + 41 でいえば、生成される 41, 43, 47, 53, 61, … は定数項を初項として、次々に 2, 4, 6, 8, … を足しているとも解釈可能。

(#20) 2019年10月。f(0) = 193 から始めて f(k) = f(k-1) + 90k for k=1..9 は10個の素数を生成し、それらのどの素数も、桁の和が等しい。上記と同系統の素数生成2次式だが「桁の和一定」という条件を付けたもの。このような最小の例。この時点で既に14要素くらいの同様の例を得ていたが、試しに小さいものを送ってみた。9の倍数を足して桁の和が変わりにくいのは当たり前という印象もあるだろうが、実際には桁の和が固定されるわけではなく、しかも素数という条件と重なるので、要素数が多い場合、シンプルな例はかなり少ない。

(#21) 2019年10月。46662004177 = 0xADD45ADD1 は、次の素数を生成する方法を16進で「自己言及」する最小の素数。14進なら ADD9ADD3 があり。ポインターのアドレスのようで面白いと思った。

2019-10-27　誰も考えことのないタイプの素数

p = 46662004177 は、一見、何の変哲もないランダムな素数に見える（反復桁が多いという特徴はあるが）。

16進数に直すと ADD45ADD1 であり、これは「45を足せ。1を足せ」と読める。実際に p+45+1 を計算すると何が起きるか。結果は q = 46662004223 であり、q は p の次の素数。すなわち概念上、素数 p は f(x) = x+45+1 という操作をコーディングした「ゲーデル数（のようなもの）」と解釈でき、その「ゲーデル数」が表す操作をその「ゲーデル数」自身に適用することで次の素数が生成される…という、ある種の自己言及的性質を持っている。ADD は自然言語であり、本物のゲーデル数ではないが、ちょっとワクワクする！

14進数の ADD9ADD3 すなわち10進数の 1159386637 は、さらに小さな同様の例。

36進数なら、アルファベットの全部の文字が使える。「完全自己言及素数」（自分自身を生成する操作を表す文字列。その文字列の指示に従うと、その数自身になる）を作れるかもしれない。

(#22) 2019年10月。129023 は、2進数で1が15個ある最小の素数。下4桁を並び替えた 2039 は、2進数で1が10個ある最小の素数。「2進数で～」の部分は既知だが、並び替えの部分が面白いと思った。#20 の「桁の和」つながりの小ネタ。

(#23) 2019年10月。285281 は、1より大きい連続奇数から作ることができる最小の4つ子素数の先頭。すなわち (285281, 285283, 285287, 285289) は 281, 283, 285, 287, 289 を文字列として連結することで作ることができる。＃19の出発点となった 2426256797 は10要素のベクトルの先頭要素だが、24, 25, 26 を含んでいる（そのような最小例ではない）。もっと小さい規模の組み合わせで同じようなことを考え、4要素のベクトルに行き着いた。「1より大きい連続奇数」の「1より大きい」を外すと、一桁の奇数を並べた形のトリビアルな例が見つかる。

2019-10-26　「対リベラ戦」第3ラウンド（笑）

数日前、Rivera は P(0) = 2426256797 が次の性質を持つことを指摘した。

性質「k=1..9 に対して P(k) = P(k-1) + 2k も素数である」

上記の主張は正しいが、それに続けて「この性質を持つような、もっと大きな列は知られていない」とも主張していて、それが問題。もっと大きな例を直接構成することは、易しい。

この列は、初項が素数 p で、階差 2, 4, 6, ... を持つのだから、一般項は p に三角数の2倍を足したもの。すなわち k^2 + k + p の形のオイラー風・素数生成式になる。これは熱心に研究されている分野で、同じ性質を持つ P(0) は何千個も記述されている（A191456）。リベラは「より大きなものは知られていない」と言うけれど、2426256797 は何千種類も知られている中で「13番目に小さい雑魚」にすぎない（追記）。

だからといって、速攻 Curios の編集者に連絡するのは「先生に言い付ける」みたいで嫌な感じなので、またメールを書いて直接リベラに連絡したのだが…。「数学的に明らかに間違っている主張」に気付いた場合、瞬時に撤回するのが普通だろう。けれどリベラは「どうするか未定。もう少し考える」という反応。仕方がないので、2426256797 を「無理やり面白く」する方法を考えてみた。

「k=1..9 に対して P(k) = P(k-1) + 2k も素数になるような最小の素数（ただし桁の6個以上が偶数であるとする）は、P(0) = 2426256797」

一応、数学的に正しい主張になった（笑）。リベラには「スペイン語バージョン」の C2 も取り消すようにそれとなく言ってみたのだが（C1 が本当のチャンピオンで C2 はあまり面白くない）、思い出が詰まっているのか、これを消すつもりはないらしい。「フランス語バージョン」については、残念ながら明らかな反例が見つかったので、議論の余地はなく、直接編集者に連絡、リベラの数（実際には第三者の主張を引用したもの）は削除された。フランス語ネタ・スペイン語ネタ・オイラー風ネタと1週間のうちに3回もリベラに文句を言っているが、たまたまタイミング的にそうなっただけ。あくまで数学上の話であり、私生活でけんかをしているわけではない（メール上ではフレンドリーな関係）。それどころか、リベラは「クォーターパウンダー素数113」がたいそう気に入ったらしく、この週末もっと大きい例を探し求め、既に100桁以上の例を得たという。

リベラのサイトは、専門家向けの記事で参照されることもある。そんな大物なのだが、天才にありがちな「大ざっぱな面」も持っているようだ。水も漏らさぬ全数検索や考察・証明をしないで「これが最良」と主張しているケースがある。「昔のパズルの結果を投稿しているので、当時見つかった最善に当たる」とのこと。結果的にそれが刺激になるのだから、悪いことともいえないが…

2019-11-04　2426256797について

リベラの P(0) = 2426256797 についての観察は、単に「k=1..9 に対して P(k) = P(k-1) + 2k も素数である」というだけでなく、「P(k) の直前の素数は P(k-1) である」すなわち「その部分の素数ギャップが 2k である」ということを言いたかったらしい。それにしても「より大きな例は知られていない」という主張はおかしい（1秒で検索できるのだから）。「最も小さな例である」と言うのなら、これはこれで興味深い。

Prime Gaps     2 4  6   8   10    12     14      16       18
Primes       +02 6 12  20   30    42     56      72       90
 2426256797 : PP-P--P---P----P-----P------P-------P--------P
 6430890287 : PP-P--P---P----P-----P------P-------P--------P
 8518049207 : PP-P--P---P----P-----P------P-------P--------P
55065405671 : PP-P--P---P----P-----P------P-------P--------P
55373581421 : PP-P--P---P----P-----P------P-------P--------P
60590486081 : PP-P--P---P----P-----P------P-------P--------P
66945470477 : PP-P--P---P----P-----P------P-------P--------P
76566117071 : PP-P--P---P----P-----P------P-------P--------P
78067026071 : PP-P--P---P----P-----P------P-------P--------P
95748657617 : PP-P--P---P----P-----P------P-------P--------P

#22 と #23 の間で、リベラからメールが来た。#6 を一般化した問題の大きい解だった。お礼として、2426256797 の「間違い」を簡単に検出できる3行の Pari コードを返信。ホイールシーブの一種（追記）。

f(P) = q=P; for( k=0, 9, q+=2*k; if(!isprime(q),return) ); print(P);
g(n) = f(n+11); f(n+17); f(n+41); f(n+101); f(n+137); f(n+167);
forstep( n=0, 3e9, 210, g(n) );

(#24) 2019年12月。1062881: For n=9, 10, 11, 12, 13, the largest prime you can express using only n copies of n without concatenation, with any of the four operations and parentheses, is 2*n^(n-3)-1: (9+9)*9*9*9*9*9-9/9=1062881, (10+10)*10*10*10*10*10*10-10/10=19999999....

2019-12-11　わたしがタイタンになった日

ある理由から f(n) = 2nⁿ⁻³ − 1 の形の素数を考えた。n ≤ 1000 の範囲では n = 4, 6, 7, 9, 10, 11, 12, 13, 129, 169, 645 の11個に限られる。f(7) = 4801, f(9) = 1062881, f(10) = 19999999 は、curios でもある。

副産物として 2 × 645⁶⁴² − 1 が素数と判明。1805桁の素数に行き当たった。「1000桁以上の素数の発見者」は、かつてタイタンと呼ばれた（今どきのコンピューターなら大したことではないが）。f(n) 型の次の素数があるとすれば、恐らく1万桁以上（n > 3000）。

2019-12-12　わたしがタイタンになった日（２）

以下 n ≤ 1000 とする。g(n) = 2nⁿ⁻³ + 1 の形の素数は n = 3, 6, 8, 20, 75, 168 の6個の場合に限られる。

h(n) = nⁿ⁻² + 1 の形の素数は、n = 4, 6 の場合の 17, 1297 に限られる。

h′(n) = nⁿ⁻² − 1 型の素数は、n = 3 の場合の 2 しかない。

そうすると [6, 13] の範囲の8個の自然数 n のそれぞれについて、f(n), g(n), h(n) の少なくとも一つは素数である。特に、[9, 13] の範囲の連続5整数に対して、f(n) は素数。

これらの観察の意味は後回しにして、次に (n + 1)nⁿ⁻⁴ − 1 型の素数を考えると、n = 5, 6, 9, 13, 22, 44, 49, 228, 344, 851 が該当する。最後のケース 852 × 851⁸⁴⁷ − 1 は、2485桁の素数。

2019-12-13　9個の9で書ける最大素数

昨日までの考察は「n 個の n と四則演算・丸括弧だけを使って書き表せる最大の素数は何か」という問題に関連している。例えば「9個の9」の場合、全数検索によれば p = (9+9)*9*9*9*9*9-9/9 = 1062881 が最大の素数。文字列としての連結（例えば2個の9を「99」として使うこと）は許されないものとする。

この「○個の△を使って何々を作れ」は古典的な数学パズルの一つであり、全数検索の原理は、「西暦・平成パズル」を解くアルゴリズムとして紹介されている。それだけなら「ありふれた問題」だろう。

ところで、上記の p は下記の f(9) に他ならない。n = 10 に対する解も、同様に
f(10) = (10+10)*10*10*10*10*10*10-10/10 = 19999999
であり、この性質は n = 11, 12, 13 に対しても成り立つ（しかし n = 14 に対しては不成立）。

簡単な考察によると、n ≥ 4 に対して、もし h(n) が素数ならそれが解であり、そうでなくもし h′(n) が素数ならそれが解であり、そうでなくもし g(n) が素数ならそれが解であり、そうでなくもし f(n) が素数ならそれが解。n = 9, 10, 11, 12, 13 の連続5整数が、この最後の「そうでなくもし」に該当し、その結果「n 個の n で書ける最大素数」が f(n) の形となる。これら4関数の値が一つでも素数になれば、面倒な全数検索を省いて、いきなり解が得られるわけである。この発見（分かってみれば当たり前なのだが）は、小気味よい。

「9個の9」「10個の10」…「13個の13」の最大素数が全部同じパターンになる、という点が、注目に値する。「7個の7」も同じパターンだが、「8個の8」については
g(8) = (8+8)*8*8*8*8+8/8 = 65537
となる。

追記　2020年のCurios

24049

(#25) 2020-07-17: Consider Pell's equation x^2-Dy^2=1. If the least solution is greater than the least solution for any smaller D, then D is usually a prime congruent to 5 mod 8 (e.g., D=61); D=24049 is the largest exception to this, below 10^8.

2020-07-27　ところで、ペル・チャンピオンたち（「その30」参照）のほとんどは 8k+5 型の素数だが、値が小さい場合、例外がある。409 と 24049 は 8k+1 型で、この 24049 が「最後の例外」ではないかと思われる。他のチャンピオンたちの記録は「ケイリーの6乗アシスト」による「ドーピング」（？）の結果だが、24049 は「6乗アシスト」なしで（負のペル経由の「2乗アシスト」だけで、つまり3乗の違いという巨大なハンディーを物ともせず）チャンピオンになってしまうわけで、非常にユニークな素数（言い換えると、その平方根の連分数展開は極めて珍しいパターン）といえるだろう。1週間以上前の話だが、このネタで 24049 の方は、素数キュリオに収録された。

2020-07-29　風変わりなチャンピオン24049

24049のトリビアは次の通り。

「ペル方程式 x^2-Dy^2=1 を考える。その最小解が、より小さなDに対するどの最小解よりも大きいとき、普通、Dは8を法として5と合同の素数である。10^8未満において、D=24049はこれに対する最大の例外である。」

ペルでちょっと遊んだことがある方ならおなじみでしょうが、D=61やD=109の場合、最小の自然数解が記録的にでかくなる。この「記録的」現象が起きるDを調べると、ごく一部の例外を除き、8で割って5余る素数になっている。この観察は「野蛮な解法（その11）」でメモした。大筋では「その29」のケイリーのNoteと結び付いていた。もし x^2-Dy^2=−4 に意味のある解があるなら、その解のざっと6乗がペルの最小解であり、この現象が起きる必要条件（十分条件ではない）が、8で割って5余るということ。細部の説明は少し間違ってるかもしれないが、そんな感じで「8を法として5と合同の素数」が大暴れ。

けれど、それ以外の素数が新記録を作るケースがある。D=24049 は、その最後の例！　…と言い切りたいところだが、残念ながら、野蛮な方法では「10^8未満において」というただし書きを消せない。実際にはDが素数のケースについて、もう一桁上の10^9まで全数検索していて、頑張れば「10^9未満において」まで風呂敷を広げられるが、いずれにしても、単なるブルートフォースでは有限の範囲でしか議論できない…。直観的には、上記「6乗オーダー」をそれ以外の方法で乗り越えるのは無理なので、D=24049 が最後だと予想されるが、本当にそうなのだろうか。議論の便宜上、そうだと仮定するとして、24049 という数は何がそんなに特別なのだろう？

学校の数学には必ず正解があるんだろうけど、こういう「遊びの数学」「野性の数学」は分からない。頑張ればできることなのか、専門知識が必要な難解なことなのか、はたまた専門家にも未解決の問題なのか…

数論は未知がいっぱいの不思議な世界なのに、散策を楽しむ暇もなく、規格化された教科書を与えられ、天下り的に理論を押し付けられ、競争材料にされ、一般的なイメージとしては「数学＝ストレス」になっているのは、切ない。教科書はゲームでいえば攻略本のようなもの。ゲームをしないで攻略本だけ読んでも、楽しいわけない。ひたすら高い山に登りたければ「200年前の数学者が10年かけて見つけたこと」を3行にまとめてもらった方が効率がいいのだが、幸せを味わいたいのなら「200年前の数学者が発見したこと」を独立に再発見する喜びの方が大きい。

63158401

(#26) 2020-08-01: The exact number of seconds in two years can be prime: 63158401. This happened for example in 2011-12, 2012-13, 2016-17, when we had 1 leap day and 1 leap second. The number of days in the last two millennia was also prime (see PC730487).

#25、#26は両方収録されたが、#26に関してはネタがかぶったため（別の方が実質的に同じことを投降し、未編集のままになっていた）、その方の名義で掲載された。わざわざ G. L. Honaker, Jr. からメールで連絡があって、何事かと驚いた（これまで採用になっても、不採用になっても、一度もメール連絡などなかったので）。

199081

(#27) 2020-08-09: The smallest prime of the form p^3+q^3+r^3, where (p,q,r) is a prime triplet.

シンプルでたわいないが、面白いと思ったので…。

追記　2021年・2022年のCurios

2022年2月現在、このうち実際に収録されているのは、#29 と #30。

18379

(#28) 2021-11-25: The smallest integer d such that any of the first 10 primes can be written as |x^2-d*y^2| with some integers x, y.

5231

(#29) 2021-12-22: The smallest p > 5 such that both (p, p+2, p+6) and (prime(p), prime(p+1), prime(p+2)) are prime triples: both (5231, 5233, 5237) and (51131, 51133, 51137) are prime triples, where prime(5231) = 51131, prime(5233) = 51137.

1331

(#30) 2022-02-20: The start of the smallest "anti prime quadruple," where none of {30n+11, 30n+13, 30n+17, 30n+19} is prime.

84

(#31) 2022-06-04: There are exactly 84 Gaussian primes with a norm <= 84. That is, 8*f(84) + 4*g(84) = 84, where f(n) = #{primes 4k+1 <= n} and g(n) = #{primes not 4k+1 <= sqrt(n)}. Although 8*f(84) can't be too different from 4*pi(84) = 92 plus we have 4*g(84) = 12, this equality holds, as Chebyshev's bias is abnormally high at 83.

「数えてびっくりガウス素数」の関連ネタ。カウントでは、単数倍の違いは無視されていない。100、128なども同じ性質を持つが、83におけるバイアスは極端で、この数以下の素数の中にはバニラ素数が40%未満しかない。

---