コンテンツにスキップ

コラッツ予想:3行で書けるのに誰も解けない問題

前提:ゼロで割ってはいけない理由:0 に逆数を与えると数の世界が 1 点に潰れる

生 Markdown
  • ルールはたった 2 行です。自然数 nn が偶数なら n/2n/2、奇数なら 3n+13n+1 に置き換える。これを繰り返すと必ず 11 にたどり着く——これがコラッツ予想で、1930 年代に提出されて以来未解決です。
  • 小さい数で試すとすぐ 11 に落ちます。ところが 2727 から始めると 111111 手かかり、途中で 92329232 まで跳ね上がります。「単純=易しい」ではないことを、この 1 例が教えてくれます。
  • 何も分かっていないわけではありません。「奇数が 1 個または 2 個だけの巡回は自明なもの以外に存在しない」「少なくとも 4 分の 3 の自然数は 3 手以内に元より小さくなる」といった主張は、この記事の中で完全に証明できます。
  • 「たぶん正しい」と信じられている理由は、1 手あたり平均して 3/40.866\sqrt{3/4} \approx 0.866 倍に縮むという確率的な見積もりです。ただしこれは証明ではありません。
  • 3n+13n+13n13n-1 に変えるだけで反例(5147201055 \to 14 \to 7 \to 20 \to 10 \to 5)が現れます。つまり正しい証明は「+1+11-1 か」を区別できるほど繊細でなければならず、そこが最大の壁です。

1. 動機:ルールを 3 行で書ける未解決問題

Section titled “1. 動機:ルールを 3 行で書ける未解決問題”

数学の未解決問題というと、ふつうは問題文を理解するだけでひと苦労します。リーマン予想を人に説明しようとすれば、まず複素数とゼータ関数の話から始めなければなりません。

コラッツ予想はそうではありません。ルールはこれだけです。

  1. 自然数を 1 つ選ぶ。
  2. それが偶数なら 2 で割る。奇数なら 3 倍して 1 を足す。
  3. 2 に戻る。

そして予想はこうです。どんな自然数から始めても、いつかは 11 に到達する。

小学生でも遊べます。実際に 66 から始めてみましょう。

631051684216 \to 3 \to 10 \to 5 \to 16 \to 8 \to 4 \to 2 \to 1

88 手で着地しました。11 に着いたあとは 14211 \to 4 \to 2 \to 1 とぐるぐる回るので、そこで打ち切ります。

この問題は、1937 年頃にドイツの数学者ローター・コラッツが提出したとされます。以来、角谷静夫にちなむ「角谷の問題」、ウラム、ハッセ、スウェイツなど多くの名前が付き、「シラキュース問題」「3n+13n+1 問題」とも呼ばれてきました。名前が多いというのは、それだけ多くの人が独立に取り憑かれたということです。

エルデシュ・パールはこの問題に 500 ドルの懸賞をかけたうえで、現代の数学はまだこの種の問題を扱えるほど成熟していない、という趣旨のことを述べたと伝えられています。2021 年には日本の企業が 1 億 2000 万円の懸賞金を出したことも話題になりました。それでもまだ解けていません。

2. 準備:コラッツ写像・軌道・停止時間

Section titled “2. 準備:コラッツ写像・軌道・停止時間”

遊びを数学にするために、言葉を決めます。以下、N={1,2,3,}\mathbb{N} = \{1, 2, 3, \ldots\} とし、00 は含めません。

定義 2.1コラッツ写像

写像 C:NNC : \mathbb{N} \to \mathbb{N}

C(n)={n/2(n が偶数)3n+1(n が奇数)C(n) = \begin{cases} n/2 & (n \text{ が偶数}) \\ 3n+1 & (n \text{ が奇数}) \end{cases}

で定める。これをコラッツ写像と呼ぶ。

nn が奇数のとき 3n+13n+1 は必ず偶数なので、奇数の次には必ず偶数が来ます。そこで「奇数のときは 3n+13n+1 してすぐ 2 で割る」とまとめた写像も便利です。

定義 2.2短縮コラッツ写像

写像 T:NNT : \mathbb{N} \to \mathbb{N}

T(n)={n/2(n が偶数)(3n+1)/2(n が奇数)T(n) = \begin{cases} n/2 & (n \text{ が偶数}) \\ (3n+1)/2 & (n \text{ が奇数}) \end{cases}

で定める。これを短縮コラッツ写像と呼ぶ。

TTCC の 1 手または 2 手をまとめたものなので、CC11 に到達することと TT11 に到達することは同値です。証明では計算が短くなる方を使います。

定義 2.3軌道・停止時間・巡回

nNn \in \mathbb{N} に対し、列 n,C(n),C2(n),n,\, C(n),\, C^2(n),\, \ldotsnn軌道と呼ぶ。

σ(n)=min{k1:Ck(n)<n}\sigma(n) = \min\{\, k \ge 1 : C^k(n) < n \,\}

nn停止時間と呼ぶ(そのような kk が存在しないときは σ(n)=\sigma(n) = \infty と定める)。

また、Ck(m)=mC^k(m) = m となる k1k \ge 1 が存在するとき、mm を含む有限集合 {m,C(m),,Ck1(m)}\{m, C(m), \ldots, C^{k-1}(m)\}巡回と呼ぶ。

停止時間は「元の数より小さくなるまでに何手かかるか」です。11 に着くまでの手数(総停止時間)より扱いやすく、後で見るようにこちらだけで予想を言い換えられます。

flowchart LR
A["自然数 n"] --> B&#123;"n は偶数か"&#125;
B -- "はい" --> C["n / 2 に置き換える"]
B -- "いいえ" --> D["3n + 1 に置き換える"]
C --> E&#123;"n = 1 か"&#125;
D --> E
E -- "いいえ" --> A
E -- "はい" --> F["終了"]
コラッツ写像の 1 手

理屈より先に、数を眺めます。

例 3.11 から 12 までの成績表

nn について、11 に到達するまでの手数(総停止時間)と、軌道の最大値を並べます。

nn手数軌道の最大値
101
212
3716
424
5516
6816
71652
838
91952
10616
111452
12916

たとえば 77 の軌道は

72211341752261340201051684217 \to 22 \to 11 \to 34 \to 17 \to 52 \to 26 \to 13 \to 40 \to 20 \to 10 \to 5 \to 16 \to 8 \to 4 \to 2 \to 1

で、確かに 16 手、最大値 5252 です。隣り合う 6677 で手数が 8816168899331919 と、まったく揃いません。この不規則さがこの問題の本質です。

例 3.227 という暴れ者

2727 から始めると、軌道はこう動き出します。

278241124623194471427121410727 \to 82 \to 41 \to 124 \to 62 \to 31 \to 94 \to 47 \to 142 \to 71 \to 214 \to 107 \to \cdots

下がったかと思うとまた上がる、を延々と繰り返し、途中で最大値 92329232 に達し、111111 手かけてようやく 11 に着きます。出発点は 2727、つまり 22 桁です。それが 44 桁まで登るのです。

2727 の隣の 26261010 手、28281818 手で終わります。2727 だけが突出しているわけで、「小さい数だから短いはず」という直感はここで完全に壊れます。

手で追うのは大変なので、プログラムに任せます。次のコードは標準ライブラリだけで動きます。

def total_stopping_time(n):
"""n が 1 に到達するまでの手数と軌道の最大値を返す。"""
steps, peak = 0, n
while n != 1:
n = n // 2 if n % 2 == 0 else 3 * n + 1
peak = max(peak, n)
steps += 1
return steps, peak
print(total_stopping_time(27)) # (111, 9232)
print(max(range(1, 10**6), key=lambda n: total_stopping_time(n)[0]))
# 837799 (100 万未満で最も手数が多い数。524 手)

最後の行を実行すると 837799837799 が返ります。手数は 524524 です。100100 万未満の数で最悪でも 524524 手、というのは「意外と小さい」と感じるでしょうか。それとも「11 に着く保証がないのに 524524 手も彷徨うのか」と感じるでしょうか。どちらの感想も正しいのが、この問題の面白いところです。

4. 定義からすぐ証明できること

Section titled “4. 定義からすぐ証明できること”

「何も分かっていない」わけではありません。ここでは紙と鉛筆だけで証明できる事実を 3 つ挙げます。

定理 4.12 の冪は素直

k0k \ge 0 を整数とすると、Ck(2k)=1C^k(2^k) = 1 である。すなわち 2k2^k はちょうど kk 手で 11 に到達する。

証明(定理 4.1)

kk についての帰納法で示します。k=0k = 0 のとき 20=12^0 = 1 で、C0(1)=1C^0(1) = 1 ですから成立します。

kk で成立するとします。2k+12^{k+1} は偶数なので、定義 2.1 より C(2k+1)=2k+1/2=2kC(2^{k+1}) = 2^{k+1}/2 = 2^k です。よって

Ck+1(2k+1)=Ck(C(2k+1))=Ck(2k)=1C^{k+1}(2^{k+1}) = C^{k}\bigl(C(2^{k+1})\bigr) = C^k(2^k) = 1

となり、最後の等号で帰納法の仮定を使いました。以上で k+1k+1 でも成立します。

つまり 2,4,8,16,32,2, 4, 8, 16, 32, \ldots は一直線に落ちます。予想が難しいのは、こういう素直な数のせいではありません。

命題 4.24 で割って 1 余る数は 3 手で小さくなる

nn を偶数とすると σ(n)=1\sigma(n) = 1 である。また n5n \ge 5n1(mod4)n \equiv 1 \pmod 4 を満たすならば σ(n)=3\sigma(n) = 3 であり、しかも

C3(n)=3n+14C^3(n) = \frac{3n+1}{4}

が成り立つ。したがって、σ(n)3\sigma(n) \le 3 を満たす自然数の(自然密度の意味での)割合は少なくとも 3/43/4 である。

証明(命題 4.2)

nn が偶数なら 定義 2.1 より C(n)=n/2<nC(n) = n/2 < n なので σ(n)=1\sigma(n) = 1 です。

次に n1(mod4)n \equiv 1 \pmod 4n5n \ge 5 とし、n=4k+1n = 4k+1k1k \ge 1)と書きます。nn は奇数なので

C(n)=3(4k+1)+1=12k+4.C(n) = 3(4k+1)+1 = 12k+4 .

これは偶数なので C2(n)=6k+2C^2(n) = 6k+2、これも偶数なので C3(n)=3k+1C^3(n) = 3k+1 です。ここで

3n+14=12k+3+14=3k+1\frac{3n+1}{4} = \frac{12k+3+1}{4} = 3k+1

なので、主張の等式が確かめられました。

σ(n)=3\sigma(n) = 3 であることを見ます。まず C3(n)=3k+1<4k+1=nC^3(n) = 3k+1 < 4k+1 = nk1k \ge 1 から従います。一方、途中の 2 手では小さくなりません。実際 C(n)=12k+4>4k+1=nC(n) = 12k+4 > 4k+1 = n8k+3>08k+3 > 0 より)、C2(n)=6k+2>4k+1=nC^2(n) = 6k+2 > 4k+1 = n2k+1>02k+1 > 0 より)です。よって nn より小さくなる最初の時刻はちょうど 33 です。

最後に密度を数えます。11 以上 NN 以下の自然数のうち、偶数はおよそ N/2N/2 個、44 で割って 11 余るものはおよそ N/4N/4 個あり、両者は重なりません。したがって σ(n)3\sigma(n) \le 3 を満たすものは少なくともおよそ N/2+N/4=3N/4N/2 + N/4 = 3N/4 個あり、NN \to \infty として割合 3/43/4 以上を得ます。

この計算は「もっと細かい剰余で調べれば、もっと多くの nn を捕まえられるのでは」という発想につながります。実際そのとおりです。なお、以下でも使う ab(modm)a \equiv b \pmod m という合同の記法については 合同の定義(定義 3.2)[数学はなぜ難しいのか] を参照してください。

例 4.316 で割って 3 余る数も落ちる

n3(mod16)n \equiv 3 \pmod{16}、すなわち n=16k+3n = 16k+3k0k \ge 0)とします。CC を 6 回適用します。

16k+348k+1024k+572k+1636k+818k+49k+216k+3 \to 48k+10 \to 24k+5 \to 72k+16 \to 36k+8 \to 18k+4 \to 9k+2

最後の 9k+29k+216k+316k+3 より小さい(差は 7k+1>07k+1 > 0)ので、σ(n)6\sigma(n) \le 6 です。n3(mod16)n \equiv 3 \pmod{16} の数は全体の 1/161/16 を占め、これは 命題 4.2 で捕まえた偶数とも 1mod41 \bmod 4 とも重なりません(3mod163 \bmod 16 の数は 3mod43 \bmod 4 の奇数だからです)。合わせて割合は 1/2+1/4+1/16=13/161/2 + 1/4 + 1/16 = 13/16 以上になります。

この手続きを mod2k\bmod 2^k でどこまでも続けると、次の定理が得られます。

注意 4.4

テラス(R. Terras, 1976)は、σ(n)<\sigma(n) < \infty を満たす nn の自然密度が 11 であることを証明しました。つまり「ほとんどすべての自然数は、いつか自分より小さくなる」までは分かっています。これは 命題 4.2例 4.3 の手続きを mod2k\bmod 2^k で押し切り、kk \to \infty としたものです。ただし密度 11 は「例外が有限個」を意味しません。例外の集合が無限にあっても密度は 00 になり得ます。

停止時間だけを見ればよい、というのは次の言い換えから正当化されます。

定理 4.5停止時間による言い換え

次の 2 つは同値である。

(a) すべての nNn \in \mathbb{N} の軌道は 11 を含む(コラッツ予想)。

(b) すべての n2n \ge 2 に対して σ(n)<\sigma(n) < \infty である。

証明(定理 4.5)

(a) \Rightarrow (b)。 n2n \ge 2 とします。(a) より、ある k0k \ge 0Ck(n)=1C^k(n) = 1 です。n2>1n \ge 2 > 1 なので k1k \ge 1 であり、Ck(n)=1<nC^k(n) = 1 < n です。よって nn より小さくなる時刻が少なくとも 1 つ存在し、その最小値 σ(n)\sigma(n) は有限です。

(b) \Rightarrow (a)。 nn についての強い帰納法で「nn の軌道は 11 を含む」を示します。

n=1n = 1 のときは軌道の最初の項が 11 なので成立します。

n2n \ge 2 とし、nn より小さいすべての自然数について主張が成り立つと仮定します。(b) より k=σ(n)k = \sigma(n) は有限で、m=Ck(n)m = C^k(n) とおくと m<nm < n、かつ m1m \ge 1 です。帰納法の仮定より mm の軌道は 11 を含み、ある j0j \ge 0Cj(m)=1C^j(m) = 1 となります。すると

Ck+j(n)=Cj(Ck(n))=Cj(m)=1C^{k+j}(n) = C^j\bigl(C^k(n)\bigr) = C^j(m) = 1

なので、nn の軌道も 11 を含みます。

5. ぐるぐる回る危険:巡回は存在するか

Section titled “5. ぐるぐる回る危険:巡回は存在するか”

予想が破れるとしたら、破れ方は 2 通りしかありません。

  1. どこかの nn の軌道が 11 を含まない巡回に入る。
  2. どこかの nn の軌道が無限に大きくなり続ける(正の無限大に発散する(定義 5.2)[ゼロで割ってはいけない理由])。

このうち 1 については、かなりのことが証明できます。

補題 5.1巡回は奇数を含む

CC の任意の巡回は、少なくとも 1 つの奇数を含む。

証明(補題 5.1)

巡回 {m,C(m),,Ck1(m)}\{m, C(m), \ldots, C^{k-1}(m)\} がすべて偶数からなると仮定します。すると 定義 2.1 より各項は前の項の半分なので、Ck(m)=m/2kC^k(m) = m/2^k です。巡回の定義から Ck(m)=mC^k(m) = m なので m=m/2km = m/2^k、すなわち m(2k1)=0m(2^k - 1) = 0 となります。k1k \ge 1 より 2k112^k - 1 \ge 1 なので m=0m = 0 ですが、これは mNm \in \mathbb{N} に反します。

巡回に含まれる奇数を順に n1,n2,,nrn_1, n_2, \ldots, n_r(互いに相異なる)とします。nin_i が奇数なら 3ni+13n_i + 1 は偶数で、そこから偶数が続く限り 2 で割られ、次の奇数 ni+1n_{i+1}(添字は rr の次を 11 と読む)に着きます。割った回数を ai1a_i \ge 1 とすると

3ni+1=2aini+1(i=1,,r)3 n_i + 1 = 2^{a_i} n_{i+1} \qquad (i = 1, \ldots, r)

が成り立ちます。この関係式が巡回を調べる出発点です。

定理 5.2奇数が 2 個以下の巡回

CC の巡回で、含まれる奇数の個数が 11 個または 22 個であるものは、{1,4,2}\{1, 4, 2\} に限る。

証明(定理 5.2)

補題 5.1 より奇数は少なくとも 1 個あります。上で導いた関係式を使います。

奇数が 1 個の場合。 r=1r = 1 なら n2=n1=nn_2 = n_1 = n と読み替えて 3n+1=2an3n + 1 = 2^a n、すなわち

n(2a3)=1.n (2^a - 3) = 1 .

nn2a32^a - 3 はともに整数で積が 11、かつ n1n \ge 1 なので n=1n = 1 かつ 2a3=12^a - 3 = 1、つまり 2a=42^a = 4a=2a = 2 です。このとき巡回は 14211 \to 4 \to 2 \to 1 で、集合として {1,4,2}\{1, 4, 2\} です。

奇数が 2 個の場合。 相異なる奇数 n1n2n_1 \ne n_2 について

3n1+1=2a1n2,3n2+1=2a2n13n_1 + 1 = 2^{a_1} n_2, \qquad 3n_2 + 1 = 2^{a_2} n_1

が成り立つとします。辺々掛けて A=a1+a22A = a_1 + a_2 \ge 2 とおくと

(3n1+1)(3n2+1)=2An1n2,(3n_1+1)(3n_2+1) = 2^{A} n_1 n_2 ,

左辺を展開して

9n1n2+3(n1+n2)+1=2An1n2,9 n_1 n_2 + 3(n_1 + n_2) + 1 = 2^A n_1 n_2 ,

すなわち

(2A9)n1n2=3(n1+n2)+1.(2^A - 9)\, n_1 n_2 = 3(n_1 + n_2) + 1 .

右辺は正なので左辺も正で、2A>92^A > 9 です。22 の冪で 99 を超える最小のものは 1616 なので 2A162^A \ge 16、よって 2A972^A - 9 \ge 7 です。n1n2>0n_1 n_2 > 0 なので

7n1n2(2A9)n1n2=3(n1+n2)+1.7 n_1 n_2 \le (2^A - 9)\, n_1 n_2 = 3(n_1 + n_2) + 1 .

一般性を失わず n1<n2n_1 < n_2 とします(相異なるので等号は起きません)。

n1=1n_1 = 1 のとき、この不等式は 7n23+3n2+17 n_2 \le 3 + 3 n_2 + 1、すなわち 4n244 n_2 \le 4 となり n21n_2 \le 1 です。これは n2>n1=1n_2 > n_1 = 1 に矛盾します。

n13n_1 \ge 3 のとき(n1n_1 は奇数なので 11 の次は 33 です)、左辺は 7n1n221n2=3n2+18n27 n_1 n_2 \ge 21 n_2 = 3 n_2 + 18 n_2 と下から評価できます。一方 n2>n13n_2 > n_1 \ge 3 より 18n2>18n13n1+118 n_2 > 18 n_1 \ge 3 n_1 + 115n145>115 n_1 \ge 45 > 1 だから)です。したがって

7n1n23n2+18n2>3n2+3n1+1=3(n1+n2)+17 n_1 n_2 \ge 3 n_2 + 18 n_2 > 3 n_2 + 3 n_1 + 1 = 3(n_1+n_2) + 1

となり、先ほどの不等式に矛盾します。

以上より奇数が 2 個の巡回は存在せず、巡回は {1,4,2}\{1,4,2\} に限られます。

この議論は奇数の個数 rr を増やすと急激に難しくなりますが、同じ関係式を精密化する方向で研究が進んでいます。現在では、11 を含まない巡回が存在するとすれば、その長さは 1010 億手を優に超えることが示されています。

flowchart LR
n12["12"] --> n6["6"] --> n3["3"] --> n10["10"] --> n5["5"] --> n16["16"] --> n8["8"] --> n4["4"] --> n2["2"] --> n1["1"]
n80["80"] --> n40["40"] --> n20["20"] --> n10
n13["13"] --> n40
n21["21"] --> n64["64"] --> n32["32"] --> n16
n128["128"] --> n64
1 に流れ込む数たち(矢印は 1 手の適用)

この図の枝をどこまでも伸ばしていったとき、すべての自然数がこの 1 本の木の中に現れるか——それがコラッツ予想です。図を見ると 1616 には 323255 の 2 本、4040 には 80801313 の 2 本が流れ込んでいます。偶数 mm には必ず 2m2m が流れ込み、さらに m4(mod6)m \equiv 4 \pmod 6 のときは奇数 (m1)/3(m-1)/3 も流れ込みます。木は上に向かって指数的に広がるので(指数的な増大がどれほど速いかは 紙を 42 回折る話(例 6.1)[数学はなぜ難しいのか] が分かりやすい例です)、「全部を尽くしているか」を確かめるのは目で見るほど簡単ではありません。

6. なぜ「たぶん正しい」と思われているのか

Section titled “6. なぜ「たぶん正しい」と思われているのか”

証明はないのに、多くの数学者はコラッツ予想が正しいと考えています。根拠は次の見積もりです。

奇数 nn に短縮写像 定義 2.2 を当てると (3n+1)/21.5n(3n+1)/2 \approx 1.5 n、偶数に当てると 0.5n0.5 n です。ここで「軌道に現れる数の偶奇はコイン投げのようにランダムだ」と仮定してみます。すると 1 手あたりの倍率の幾何平均

1.5×0.5=0.750.866\sqrt{1.5 \times 0.5} = \sqrt{0.75} \approx 0.866

です。11 より小さい。つまり典型的な軌道は、増えたり減ったりしながら平均としては 1 手ごとに約 13% ずつ縮んでいくはずだ、というわけです。縮み続けるなら、いつかは小さな数に落ちる。落ちればあとは 定理 4.5 の帰納法が効きます。

この見積もりは手数の予測まで与えます。CC で考えると、奇数 1 個につき平均 2 回の halving が続く(3n+13n+12j2^j でちょうど割り切れる確率が 2j2^{-j} なので期待値が 22)ので、奇数 kk 個・偶数 2k2k 個で全体の倍率は 3k/22k=(3/4)k3^k / 2^{2k} = (3/4)^k です。これが 1/n1/n になるのは k=lnn/ln(4/3)k = \ln n / \ln(4/3) のときで、総手数は 3k10.4lnn3k \approx 10.4 \ln n と予測されます。nn100100 万程度なら 10.4×13.814410.4 \times 13.8 \approx 144 手。実測でもその近辺に集中します。悪くない予測です。

証明できている最良の結果は、この確率的直観を厳密化する方向にあります。テレンス・タオは 2019 年に、任意の発散関数 fff(n)f(n) \to \infty)について、ほとんどすべての nn の軌道が f(n)f(n) より小さい値を取ることを示しました。「ほとんどすべての nn はいったんは非常に小さくなる」ところまでは来ているわけです。それでも「すべての nn11 に着く」との間には、依然として深い溝があります。

7. なぜ難しいのか:3n−1 という残酷な反例

Section titled “7. なぜ難しいのか:3n−1 という残酷な反例”

素朴な証明方針が失敗する理由を、いちばん鮮明に示す例を挙げます。ルールの +1+11-1 に変えてみます。

命題 7.13n−1 版には自明でない巡回がある

写像 D:NND : \mathbb{N} \to \mathbb{N} を、nn が偶数なら D(n)=n/2D(n) = n/2、奇数なら D(n)=3n1D(n) = 3n - 1 で定める。このとき DD{1,2}\{1, 2\} 以外に巡回

5147201055 \to 14 \to 7 \to 20 \to 10 \to 5

をもつ。

証明(命題 7.1)

順に計算します。55 は奇数なので D(5)=351=14D(5) = 3 \cdot 5 - 1 = 141414 は偶数なので D(14)=7D(14) = 777 は奇数なので D(7)=371=20D(7) = 3 \cdot 7 - 1 = 202020 は偶数なので D(20)=10D(20) = 101010 は偶数なので D(10)=5D(10) = 5。よって D5(5)=5D^5(5) = 5 であり、{5,14,7,20,10}\{5, 14, 7, 20, 10\} は巡回です。11 を含まないので {1,2}\{1,2\} とは異なります。

なお D(1)=2D(1) = 2D(2)=1D(2) = 1 なので {1,2}\{1,2\} も巡回です。

これが効いてきます。第 6 節の確率的な見積もりは、+1+11-1 に変えてもまったく同じです(3n13n3n-1 \approx 3n ですから)。つまり「平均して縮むから 11 に落ちる」という議論は、反例のある問題に対しても同じ結論を出してしまうのです。したがってこの方針は、そのままでは決して証明になりません。

例 7.2負の数まで含めるともっとひどい

DD の巡回は上のものだけではありません。

1750257437110551648241122611829127213668341717 \to 50 \to 25 \to 74 \to 37 \to 110 \to 55 \to 164 \to 82 \to 41 \to 122 \to 61 \to 182 \to 91 \to 272 \to 136 \to 68 \to 34 \to 17

も巡回です(各手を 命題 7.1 と同じ要領で確かめてください)。

さらに、nnn \mapsto -n と置き換えると DD は元のコラッツ写像 CC そのものになります。したがって CC を負の整数まで拡張すると、121-1 \to -2 \to -1514720105-5 \to -14 \to -7 \to -20 \to -10 \to -5、そして 17-17 から始まる長さ 1818 の巡回、と 3 つの巡回が現れます。CC の式は正負で何も変わらないのに、正の側だけ巡回が 11 つしかない(と予想される)——この非対称性を説明できる理屈が、今のところ誰にも見つかっていません。

もう 1 つ、深いところからの警告があります。ジョン・コンウェイは 1972 年、コラッツ写像を「nnmm で割った余りごとに別々の一次式を当てる」形に一般化した写像の族について、与えられた出発点が 11 に到達するかどうかを判定する一般的アルゴリズムが存在しない(アルゴリズム的に決定不能である)ことを示しました。コラッツ予想そのものが決定不能だと言っているわけではありませんが、「この種の問題に万能の解法はない」ことが証明されている、というのは重い事実です。

注意 7.3計算による検証の現状

コンピュータによる検証は、2682.95×10202^{68} \approx 2.95 \times 10^{20} までのすべての自然数が 11 に到達することを確かめています(Barina, 2021)。これは膨大ですが、無限に比べれば何でもありません。実際、最小の反例が 1010010^{100} のあたりに潜んでいたとしても、現在の検証はそれを一切排除できません。有限個の確認をいくら積み上げても証明にならないことは、40 回当たって 41 回目に外れる公式(例 4.2)[数学はなぜ難しいのか] が端的に示すとおりです。

一方、理論的な側からは、xx 以下の自然数のうち 11 に到達するものの個数が少なくとも x0.84x^{0.84} 個であること(Krasikov–Lagarias, 2003)などが示されています。x0.84x^{0.84}xx に比べればずっと少ない——「ほとんど全部」を数え上げる作業が、いかに難しいかが分かります。

演習 8.1

CC定義 2.1 の写像とする。n=9n = 9 の軌道を書き下し、11 に到達するまでの手数と軌道の最大値を求めよ。

解答

順に計算します。

9281472211341752261340201051684219 \to 28 \to 14 \to 7 \to 22 \to 11 \to 34 \to 17 \to 52 \to 26 \to 13 \to 40 \to 20 \to 10 \to 5 \to 16 \to 8 \to 4 \to 2 \to 1

99 を第 00 項として数えると、11 は第 1919 項です。よって手数は 1919 です。最大値は 5252 です。

例 3.1 の表で 77 の手数が 1616 でしたが、上の軌道は 3 手目で 77 に落ちています。3+16=193 + 16 = 19 で一致します。このように「途中で既知の数に合流したら、そこから先は計算し直さなくてよい」というのが、コラッツ数列を高速に計算するときの基本テクニックです。

演習 8.2標準

nn を奇数とする。4n+14n+1 もまた奇数であることを確かめたうえで、

C3(4n+1)=C(n)C^3(4n+1) = C(n)

を示せ。これは何を意味するか、n=5n = 5 の場合で確かめよ。

解答

nn が奇数なら 4n4n は偶数なので 4n+14n+1 は奇数です。よって

C(4n+1)=3(4n+1)+1=12n+4.C(4n+1) = 3(4n+1) + 1 = 12n + 4 .

これは偶数なので C2(4n+1)=6n+2C^2(4n+1) = 6n + 2、これも偶数なので

C3(4n+1)=3n+1.C^3(4n+1) = 3n + 1 .

一方 nn は奇数なので C(n)=3n+1C(n) = 3n+1 です。両者は一致します。

つまり 4n+14n+1 の軌道は 3 手で nn の軌道に合流します。

n=5n = 5 のとき 4n+1=214n+1 = 21 です。2164321621 \to 64 \to 32 \to 16 で、C(5)=16C(5) = 16 に確かに合流しています。この事実は、定理 5.2 の後の図で 21216464 を通って 1616 に流れ込んでいたことに対応します。奇数 nn に対して n,4n+1,4(4n+1)+1=16n+5,n, 4n+1, 4(4n+1)+1 = 16n+5, \ldots が全部同じところに合流するので、コラッツの木には無限に続く「そっくりな枝」がたくさんあります。

演習 8.3標準

n11(mod16)n \equiv 11 \pmod{16} を満たす nn は、例 4.3 と同じ 6 手では nn より小さくならないことを示せ。

解答

n=16k+11n = 16k + 11k0k \ge 0)とおいて計算します。nn は奇数なので

C(n)=48k+34.C(n) = 48k + 34 .

偶数なので C2(n)=24k+17C^2(n) = 24k + 17、これは奇数なので C3(n)=72k+52C^3(n) = 72k + 52、偶数なので C4(n)=36k+26C^4(n) = 36k + 26、偶数なので C5(n)=18k+13C^5(n) = 18k + 13、これは奇数なので

C6(n)=54k+40.C^6(n) = 54k + 40 .

54k+40>16k+1154k + 40 > 16k + 11(差は 38k+29>038k + 29 > 0)なので、6 手では小さくなっていません。途中の値もすべて 16k+1116k+11 より大きいことは、たとえば C5(n)=18k+13>16k+11C^5(n) = 18k+13 > 16k+11(差 2k+2>02k+2 > 0)のように順に確かめられます。

つまり 例 4.3 の手続きは剰余類ごとに「当たり外れ」があり、外れたものは mod32\bmod 32mod64\bmod 64 とさらに細かく分けて調べ直すことになります。注意 4.4 のテラスの定理は、この細分をどこまでも続けると取りこぼしの割合が 00 に近づく、という主張です。

演習 8.4

TT定義 2.2 の短縮コラッツ写像とする。k0k \ge 0m1m \ge 1 に対して

Tk(2km1)=3km1T^k(2^k m - 1) = 3^k m - 1

が成り立つことを kk についての帰納法で示せ。これを使って、n=2k1n = 2^k - 1 の形の数が最初のうち大きく増えることを説明せよ。

解答

k=0k = 0 のとき。 左辺は T0(m1)=m1T^0(m - 1) = m - 1、右辺は m1m - 1 で一致します。

kk で成立すると仮定して k+1k+1 を示す。 N=2k+1m1N = 2^{k+1} m - 1 とおきます。2k+1m2^{k+1} m は(k+11k + 1 \ge 1 なので)偶数なので NN は奇数です。よって 定義 2.2 より

T(N)=3(2k+1m1)+12=32k+1m22=32km1=2k(3m)1.T(N) = \frac{3(2^{k+1} m - 1) + 1}{2} = \frac{3 \cdot 2^{k+1} m - 2}{2} = 3 \cdot 2^{k} m - 1 = 2^k (3m) - 1 .

これは 2km12^k m' - 1 の形(m=3m1m' = 3m \ge 1)なので、帰納法の仮定が使えて

Tk+1(N)=Tk(T(N))=Tk(2km1)=3km1=3k3m1=3k+1m1.T^{k+1}(N) = T^k\bigl(T(N)\bigr) = T^k(2^k m' - 1) = 3^k m' - 1 = 3^k \cdot 3m - 1 = 3^{k+1} m - 1 .

以上で帰納法が完成します。

意味。 m=1m = 1 とすると、n=2k1n = 2^k - 1kk 手の短縮ステップで 3k13^k - 1 になります。倍率はおよそ (3/2)k(3/2)^k で、kk が大きいほど激しく増えます。第 6 節の確率的な見積もりでは「1 手あたり平均 0.8660.866 倍」でしたが、この形の数は最初の kk 手すべてが奇数ステップで、平均から大きく外れるのです。

たとえば k=5k = 5 なら n=31n = 31 で、

31477110716124231 \to 47 \to 71 \to 107 \to 161 \to 242

と 5 手で 242242 まで登ります(351=2423^5 - 1 = 242 です)。例 3.22727 が暴れたのも同じ理由で、2727 の軌道は途中で 3131 を通ります。「たまたま奇数が連続する数」はいくらでも作れるので、確率的な議論だけで予想を証明できない理由がここにも見えます。

  • J. C. Lagarias, “The 3x+1 problem and its generalizations”, American Mathematical Monthly 92 (1985), 3–23. — この問題の標準的な入門サーベイ。歴史的経緯と主要結果がまとまっています。
  • J. C. Lagarias (ed.), The Ultimate Challenge: The 3x+1 Problem, American Mathematical Society, 2010. — 上のサーベイを含む論文集。現在の到達点を知るならまずこれです。
  • R. Terras, “A stopping time problem on the positive integers”, Acta Arithmetica 30 (1976), 241–252. — 停止時間が有限な自然数の密度が 11 であることの原論文。
  • I. Krasikov and J. C. Lagarias, “Bounds for the 3x+1 problem using difference inequalities”, Acta Arithmetica 109 (2003). — xx 以下で 11 に到達する数の個数の下界。
  • T. Tao, “Almost all orbits of the Collatz map attain almost bounded values”, arXiv:1909.03562 (2019). — 確率的直観を厳密化した現時点で最強の結果。
  • D. Barina, “Convergence verification of the Collatz problem”, The Journal of Supercomputing 77 (2021). — 2682^{68} までの計算機検証。
  • J. H. Conway, “Unpredictable iterations”, Proceedings of the 1972 Number Theory Conference, University of Colorado, 1972. — 一般化コラッツ写像の決定不能性。

Appendix: 手を動かすためのヒント

Section titled “Appendix: 手を動かすためのヒント”

自分で実験するときの注意。 第 3 節のコードは素朴な実装なので、10710^7 を超えるあたりから遅くなります。速くするには、演習 8.1 の解答で触れたメモ化(一度計算した数の手数を辞書に保存する)を入れるのが定石です。ただし辞書が巨大になるので、「nn 未満の値に落ちたら打ち切る」方式(つまり停止時間 σ(n)\sigma(n) だけを計算し、定理 4.5 に頼る)のほうがメモリ効率は良くなります。

見つけても喜びすぎないこと。 「反例を見つけた」と思ったときは、まず桁あふれを疑ってください。C 言語や Java の 64 ビット整数は 9.2×10189.2 \times 10^{18} 程度で溢れます。272792329232 まで登ったように、コラッツ軌道は出発点の何桁も上まで行くので、101810^{18} 台の数を試すと途中で簡単に溢れます。Python の整数は自動的に多倍長になるので、この点では安心です。

それでも挑戦したい人へ。 この記事で見たとおり、素朴な確率的議論は 3n13n-1 版(命題 7.1)を排除できず、剰余類による議論は 注意 4.4 の壁を越えられません。新しい証明は、+1+1 という項が正の整数の世界で果たしている特別な役割を掴まえるものでなければならないはずです。それが何なのかは、まだ誰も知りません。数学の問題がなぜ難しくなるのかについては 数学はなぜ難しいのか も参考にしてください。計算機に大きく頼って解決した例としては 四色定理定理 5.1[四色定理])が対照的で、こちらは「有限個の場合に帰着できた」こと(不可避集合と可約配置(定義 5.3)[四色定理])が決め手でした。コラッツ予想には、その有限化の手段がまだ見つかっていないのです。

この記事の誤りを報告する ・運営: 夢現技研合同会社料金プラン利用条件特定商取引法に基づく表記

© 2026 夢現技研合同会社 ・本文の LLM への入力は自由です。コード例は MIT ライセンスです。