コラッツの問題

コラッツの問題は、数論の未解決問題のひとつである。問題の結論の予想を指してコラッツ予想と言う。伝統的にローター・コラッツの名を冠されて呼ばれるが、固有名詞に依拠しない表現としては3n+1問題とも言われ、また初期にこの問題に取り組んだ研究者の名を冠して、角谷の問題、米田の予想、ウラムの予想、シラキュース問題などとも呼ばれる。

数学者ポール・エルデシュは「数学はまだこの種の問題に対する用意ができていない」と述べた。また、ジェフリー・ラガリアスは2010年に、コラッツの予想は「非常に難しい問題であり、現代の数学では完全に手が届かない」と述べたLagarias。

2019年9月、テレンス・タオはコラッツの問題がほとんどすべての正の整数においてほとんど正しいとするプレプリントを発表し、論文は2022年5月に出版された。

1. 概要
-----

任意の正の整数 n}} に対して、以下で定められる操作について考える。

//n// が偶数の場合、//n// を 2 で割る
//n// が奇数の場合、//n// に 3 をかけて 1 を足すこのとき、「どんな初期値から始めても、有限回の操作のうちに必ず 1 に到達する」という主張が、コラッツの予想である。

モジュラ計算の表記を用いて、上記の操作を関数 f}} として定義する。

任意の正の整数 n}} に対して、この関数 f}}  を繰り返し適用することで整数の無限列 //a////i//}//i// ≧ 0}}を得る。数式としての定義は次のようになる。

このとき、コラッツ予想を形式的に書くと、以下のようになる。

すなわち、任意の正整数 n}} に対して f}} を繰り返し適用していると、1=//f//…)) = 1}} がいつかは成り立つという主張がコラッツ予想である。1に到達して以降は、この数列は 1, 4, 2 を繰り返すのみとなるため、以降では最初に1に到達して以降の数は表記しないことにする。

コラッツ予想が誤りであるなら、1 を含まない数列を生成する初期値が存在するということになる。そのような数列は、1 を含まない繰り返し数列に突入するか、もしくは際限なく増大していくかのいずれかである。そのような数列はいまだ見つかっていない。

コンピュータを用いた計算により、268 までの初期値には反例がないことが確かめられている。

2. 例
----

初期値を 6 にしたとき、6, 3, 10, 5, 16, 8, 4, 2, 1 という1に到達する数列を得る。
初期値を 11 とすると、11, 34, 17, 52, 26, 13, 40, 20, 10, 5, 16, 8, 4, 2, 1 となり、1 に到達するまでにより長くかかる。
初期値として 27 を選ぶと、以下のように数列の項数は112にまで及び、その値は最終的に 1 に到達する前に最大 9232 まで増大する。

各整数において1になるまでの回数はA006577}}を，各整数に対する最大値はA025586}}を参照。
「3を掛けて1を加えた後、奇数になるまで２で割る」という1回の作業で1となる数値は、1 と無限に存在する。

3. 歴史
-----

1930年代ごろ、学生であったコラッツは数論的関数を有向グラフとして表現して、そのグラフの性質が関数のふるまいとどのように関連付けられるかに興味を持った。この時の彼のノートには、現在知られるものとは異なる関数について調べられている。コラッツはこれらの関数についての問題を出版することはなかったが、1950年の国際数学者会議をきっかけに広くこの問題が知られるようになり、1963年に初めて  紙面に取りあげられた。その後、この問題に興味を持った数学者によって研究される間に、様々な別名を得ることとなった。「ハッセのアルゴリズム」はヘルムート・ハッセに、「シラキュース問題」はシラキュース大学に、「角谷の問題」は角谷静夫に、「ウラムの問題」はスタニスワフ・ウラムに由来する。

コラッツ予想の解決にはいくつかの懸賞金がかけられており、コクセターが50ドル、エルデシュが500ドル、ブライアン・スウェイツBryan Thwaites}}は1000ポンドを申し出ている。

2011年度大学入試センター試験数学IIB第6問に題材として取り上げられた。

2019年、テレンス・タオは「1=lim//N//→∞ //f// = +∞}} を満たす任意の関数 //f// : ℕ+ → ℝ}} について、N}} から始まるコラッツ数列の最小値は  ほとんど全ての自然数 N}} において //f//}} より小さい」ことを証明し、コラッツ予想は「『ほとんど』全ての自然数で『ほとんど』正しい」ことが示された。帰結として例えば、N}} から始まるコラッツ数列の最小値はほとんど全ての自然数 N}} において log log log log //N//}} より小さくなることが従う。

2021年7月7日、株式会社音圧爆上げくんが、「コラッツ予想の真偽を明らかにした方に懸賞金1億2000万円を支払います」と発表した。

4. この予想を支持する論拠
--------------

コラッツの予想は未解決だが、この問題を研究した多くの数学者は正しいだろうと考えている。そう予想される理由を以下に挙げる。

4.1. 経験的証拠
----------

この予想は、初期値が268 ≈ 2.95×1020までは成り立つことがコンピュータで確認されている。この数値は、コンピュータのスピードアップとテスト技術の向上に伴って更新され続けている。

しかしこういったコンピューターでの検証は、予想が真であるという証拠にはならない。ポリア予想、メルテンス予想、スキューズ数の場合に示されているように、非常に大きな数において予想の反例が見つかることがある。

4.2. ヒューリスティクス
--------------

厳密な議論ではないヒューリスティクスであるが、ステップを経るごとに数の大きさがどのようになるかを考察する。今、//n// が偶数ならば、次のステップで大きさは半分になる。また、//n// が奇数ならば、次のステップで 3//n// + 1 になるが、これは必ず偶数であるから、その次に /2 になることまでは確定している。ここまでを一つのステップと解釈すれば、このステップで大きさは約 3/2 倍になる。1回のステップを経た後に偶数になるか、奇数になるかが半々であると考えると、1/2 の確率で数の大きさは 1/2 倍になり、残る 1/2 の確率で数の大きさは 3/2 倍になる。よって、1回のステップで、数の大きさは 1/2 × 1/2 倍になると期待される。これは 1 より小さな値であるから、ステップを経るごとに「確率的に」小さくなると考えられる。この意味で、いつかは 1 に到達するとの予想は確からしい。

確率論の言葉を用いるとこれは無限のステップ数を取る極限で1に平均収束するということであり、厳密には予想の確からしさとは無関係である。

5. コラッツの数列を計算するプログラミング例
-----------------------

以下は擬似コードによるプログラミング例である。与えられた初期値に対するコラッツ数列を計算できる。def collatzshow nif n.odd? and n > 1 collatzelse if n.even? collatz

このプログラムは 1, 4, 2, 1 の無限サイクルを省くために 1 に到達すると停止するようになっている。もしコラッツ予想が正しければ、いかなる正の整数の初期値を与えても停止する。

この問題を、初期のコンピュータで計算することで研究したのがスタニスワフ・ウラムであった。

6. サイクル
-------

n}} が奇数であるときに 3//n// + 1}} は必ず偶数になることから、コラッツ写像の「ショートカット」形式が次の定義で与えられる：

コラッツ数列のサイクルとは相異なる正整数からなる数列 }} であって、, , ..., を満たすものである。唯一知られているサイクルは長さ2の で、これは自明なサイクルと呼ばれる。

6.1. サイクルの長さ
------------

自明ではないサイクルの長さは少なくとも 17087915 である。Eliahouは1993年の論文で、サイクルの最小値が 240}} を超えるならば、周期の長さ p}} が

となることを示した。ここでは非負整数で、 かつ である。この結果は、の 連分数展開と関連している。最新の検証では、268}}までコラッツ予想が正しいと判明しているので、同様の論法によりサイクルの長さの下限は114208327604}}  であると言える。

6.2. //k//-サイクル
---------------

k}}-サイクルとは、奇数が連続する増加列 k}} 個と偶数が連続する減少列 k}} 個の 2//k//}} 個の部分で全体が構成されたサイクルである。例えばサイクルが、奇数による増加列1個と、偶数による減少列1個で構成されたなら、1-サイクルと呼ぶ。Steiner は、1-サイクルが、自明なサイクルのみであることを証明した。ジョン・サイモンズ      }}は、Steinerの方法を使って、2-サイクルが存在しないことを証明した。サイモンズとde Wegerはこの証明を拡張して、1=//k// = 68}} までの場合において、非自明な//k}}//-サイクルが存在しないことを示した。また、k}} が69以上の場合においては、存在した場合のサイクルの最小値の範囲を与えた。例えば 1=//k// = 69}} の場合は、サイクルの最小値は 3.3389より大きく 6.4877より小さいことが示されている。

7. コラッツ予想の他の形式
--------------

7.1. ボトムアップ方式
-------------

別のアプローチをとってみる：ボトムアップ方式でいわゆる//コラッツグラフ//を作成してみる。 コラッツグラフは、以下のように操作を逆にして定義する。

したがって、すべての正の整数が最終的に1になることを証明する代わりに、1がすべての正の整数に逆方向につながることを証明すればよい。 任意の整数n}}, //n// ≡ 1 }} ⇔ 3//n// + 1 ≡ 4 }} である。ゆえに、 sfrac ⇔ //n// ≡ 4 }}である。 推測的に、この逆関係は、1–2–4ループを除いてツリーを形成する。

操作3//n// + 1}}を、ショートカット操作sfracに置き換えた場合は、 コラッツグラフの定義は以下のようになる。

任意の整数n}}, //n// ≡ 1 }} ⇔ sfrac。ゆえに、 sfrac ⇔ //n// ≡ 2 }}である 。推測的に、この逆関係は、1–2ループを除いてツリーを形成する。

7.2. パリティシーケンス
--------------

本節では、コラッツ関数を少し変形したものを考える：

//n//が奇数の場合には3//n//   +   1が必ず偶数になるので上記のようにできる。

Pをパリティ数とする。P   =   0 で、P =   1 である。整数//n//のパリティシーケンスを、//pi// = P, ただし//a//0 = //n//, and //a////i//+1 = //f// と定義する。 /2 または //n///2、どちらの操作が適用されるかは、パリティに依存する。パリティシーケンスは//fによる//操作の場合分けに等しい。//f//に対してこの形式を適用すると、２つの整数//m// と//n//のパリティシーケンスは、//m// と//n//が2//k//を法として合同の場合のみ、最初の//k//項で一致するこが示される。これは、すべての整数がパリティシーケンスにより一意に識別されること意味し、さらに複数のコラッツ数列がある場合、対応するパリティシーケンスが異なる必要があることを意味する。

//n//=//a//·2//k//   +   //b// に関数 //f// を //k// 回適用すると、//a//·3//c//   +   //d//となる。ここで//d//は//b//に関数//f//を//k//回適用した結果で、//c//はその過程で3倍の演算を行った回数である。 //b//が2//k//   -   1の場合には、 //k//回の増加があり、結果は2·a·3//k//   -   1となる。//a//に掛かる係数は、//a//には無関係で、//b//にのみ依存する。これにより、特定の形式の数値が特定の反復回数の後、常により小さい数値になることを予測できます。例えば、4//a//   +   1 は、2回の//f//操作により3//a//   +   1となり、16//a//   +   3は4回の//f//操作により9·//a//   +   2となる。これらの小さくなった数が1へとつながるかどうかは、 //a//の値に依存する。

8. 問題の拡張
--------

8.1. 整数全体への拡張
-------------

正整数だけではなく、負数も含む任意の整数から開始するコラッツ数列について考えることができる。このとき、自明なサイクル0→0を除いて、合計4つの既知のサイクルがあり、知られている限りでは任意の0でない整数は以下に示された4つのサイクルのいずれかに到達する。

以下の表において、奇数の値は太字で示されている。また各サイクルは、最小絶対値の要素を先頭にしている。

一般化されたコラッツの予想は、f}} による反復の下で、すべての整数が最終的に上記の4つのサイクルの1つ  に入るという主張である。

8.2. 分母が奇数である有理数への拡張
--------------------

コラッツ写像は、分母が奇数である既約分数に拡張できる。ここで偶奇性の判定は分子の偶奇性によって判定し、それ以外は定義域が整数の場合と同様に行う。すなわち、分子が偶数であれば2で割り、奇数であればその有理数に3を掛けて1を足すこととする。「ショートカット」した定義を使用する場合、任意の周期的かつ既約なパリティ列は一意な有理数によって生成されることが知られている。逆に、分母が奇数であるすべての有理数は、周期的なパリティ列を持っていると予想されている 。

長さ n}} のパリティ列が奇数を正確に m}} 個含むとして、そのインデックスを //k//0 < … < //k////m//−1}} で表すとすると、パリティサイクルを生成する有理数が次のように決定される：

例えば、長さ7のパリティサイクル }} は、0, 2, 3, 6 の位置にパリティ1  が存在している。従って 1=//m// = 4}}, 1=//k////i// = 0, 2, 3, 6 }} と置くと、周期列を生成する有理数が次のように求まる。

これは、下記の有理数のサイクルとなり、パリティは確かに }} の周期となっていることがわかる。

}} の巡回置換は、上記の分数の1つに関連付けられる。例えば、サイクル }} から同様の手順で有理数を得ると

となり、先の例で得たサイクルの2番目に対応している。

コラッツ予想が正しいことを仮定すると、}} と }} が正の整数によって生成される唯一のパリティサイクルであることが示される。

有理数の分母が3の倍数でない場合、この拡張されたコラッツ数列はすべて同じ分母を持ち、このとき分子の列はコラッツ関数の「3//n// + //d//}}」への一般化を適用することによって得ることができる。

8.3. 2-進整数環への拡張
---------------

有理数と類似の拡張として、コラッツ写像の2-進整数環 への拡張が挙げられる。この環は部分環として奇数分母の有理数全体からなる環を含んでいる。このとき写像は

で与えられ、これは 上でwell-definedな連続写像であり、かつハール測度を保存する。加えて、その力学系はエルゴードであることが知られている。

8.4. 実数・複素数への拡張
---------------

コラッツ写像は、次の式で定義される滑らかな関数を整数に制限したものと見なせる：

実直線上でのこの写像の反復については、マーク・チェンバーランド      }}の研究によると、この関数は実直線上に無限に多くの不動点を持ち、また反復合成について単調に無限大に発散する軌道も持つ。

//f//}} の正の不動点を小さい方から 0 < μ1 < μ2 < μ3 < …}} とおく。このとき、[μ1, μ3]}} の範囲には周期アトラクターがちょうど2つ存在して、1つは 1, 2}}}、もう1つは 1.192531907..., 2.138656335...}   }}である。 Letherman, Schleicher および Wood はこの関数の研究を複素平面へと拡張した。ほとんどの複素平面上の点は無限遠へと発散するが、境界となるジュリア集合はフラクタル図形を描く。

9. 検証計算の最適化
-----------

9.1. 時間と空間のトレードオフ
-----------------

上記の //パリティシーケンス// により、コラッツ数列の計算の高速化が可能である。

k}} ステップ先にジャンプするために、まず現在の数値を2つに分割する： b}}と、a}}。k}} ステップをジャンプした結果は以下になる：

.//c//}} と //d//}} 配列は、k}} ビットの数すべてについて事前計算しておく。//d//}} は b}} に //f//}} 関数を k}} 回行った数で、//c//}} はその間に登場した奇数の数である。例えば k =}} 5}} なら5ステップのジャンプが可能で、そのために下位5ビットを分割して、下記配列を使う：

これは、2k}} 個の事前計算とストレージが要求される。これにより計算速度が k}} 倍高速化出来るが、時間と空間のトレードオフである。

9.2. Modular restrictions
-------------------------

コラッツ予想の反例を探すにあたって、上記の事前計算を使うと加速させることができる。もし、与えられた b}} と k}} で、不等式

が成り立ったとすると、すべての a}} についても成立する。そうすると、最初の反例は、 2k}} を法として b}} と合同ではないと言える。

例えば、最初の反例は必ず奇数である。なぜならば //f//  =}} n}} で、2n}} よりも小さいからである。同様に、最初の反例は mod 4}} で3と合同である。なぜならば、//f// 2 =}} 3n + 1}} で、4n + 1}} より小さいからである。コラッツ予想の反例でないすべての a}} には、上記不等式が成立する k}} が存在する。なので、ある数のコラッツ予想の検証するとき、合同なクラスを検証するとよい。k}} が増加したときの検証は、より小さい k}} で除外されなかった b}} についてのみ検証すればよい。指数関数的に小さい割合だけが残ることになるname="Lagarias  Th. D"。例えば mod 32}} では、b =}} 7, 15, 27, 31}} だけが残る。

10. 類似の問題
---------

関数 //f// の定義を少し変えることにより、コラッツの問題の類似を考えることができる。

10.1. 変数//n//が奇数の時の乗数の奇数一般への拡張による類似問題
-------------------------------------

例えば「任意の正の整数 //n// に対して

//n// が偶数の場合、//n// を 2 で割る
//n// が奇数の場合、//n// に 2//m// – 1 をかけて 1 を足すという操作を繰り返すと、有限回で 1 に到達する」という命題を考える。//m// = 1 のときは、この命題が正しいことを簡単に証明できる。//m// = 2 の場合が上述のコラッツの問題である。//m// ≥ 3 の場合は、//m//の値と//n//の初期値によっては、1を含まない繰り返し数列、もしくは際限なく増大していく数列が得られるため、この命題は一般に成り立たない。たとえば //m// = 3 の場合、nの初期値を13に設定すると、13, 66, 33, 166, 83, 416, 208, 104, 52, 26, 13 という1を含まない数列のサイクルが得られる。これは上記のヒューリスティクスの観点からして、//m//が大きくなるほど1に到達する可能性は低くなると予想されることとも符合する。

10.2. 変数//n//が奇数の時の加算数の奇数一般への拡張による類似問題
--------------------------------------

また、もう一つの類似として、「任意の正の整数 //n// に対して

//n// が偶数の場合、//n// を 2 で割る
//n// が奇数の場合、//n// に 3をかけて 2//l// - 1  を足すという操作を繰り返すと、有限回で 1 に到達する」という命題を考える。 ここで、//l// = 1 のときが上述のコラッツの問題である。しかし、//l// ≥ 2の場合、1を含まない繰り返し数列が得られる場合があるので、この命題は一般に成り立たない。

たとえば、//l// = 2として、初期値//n// = 43を与えた場合、43, 132, 66, 33, 102, 51, 156, 78, 39, 120, 60, 30, 15, 48, 24, 12, 6, 3, 12, 6, 3という数列が得られ、この命題は成り立たない。初期値//n//が1, 2などなら有限回で1に到達するが、他の初期値に対しては3, 12, 6, 3と、3を繰り返すサイクルになると思われる。そこで//l// = 2に対してコラッツの予想を応用し、「任意の正の整数 //n// に対して、上記の操作を行えば、有限回で1または3に到達する」という命題を代わりに立てれば、これが成り立つと予想される。

この二つの予想を一般化して、「任意の正の整数 //n// に対して

//n// が偶数の場合、//n// を 2 で割る
//n// が奇数の場合、//n// に 3をかけて 2//l// – 1  を足すという操作を繰り返すと、有限回で1または2//l// – 1  に到達する」という命題を立てたとしても、//l// ≥ 3以上の場合には、この命題は一般に成り立たない。たとえば//l// = 3の場合、任意の自然数//n//が1または5に到達するという命題になるが、//n//＝13の時、13, 44, 22, 11, 38, 19, 62, 31, 98, 49, 152, 76, 38, 19と、19を繰り返す無限ループになり、1にも5にも到達しない。

ただし、上の、2//l// – 1  が、0以上の整数//a//を用いて3//a//–1  で表されるときには、上記のプロセスを繰り返せば、有限回数で1または3//a//–1  に到達することは予想される。//a// = 1の場合がコラッツの問題である。//a// = 2の場合は、上記で//l// = 2のケースである。

10.3. 変数//n//が奇数の時の乗数と加算数双方の、奇数一般への拡張による類似問題
--------------------------------------------

以上のことから、一般化は困難ではあるが、個別に考えるなら、さらに進んで、「任意の正の整数 //n// に対して

//n// が偶数の場合、//n// を 2 で割る
//n// が奇数の場合、//n// に 2//m// – 1  をかけて、2//l// – 1  を足すという操作を繰り返すとき、//n//、//m//、//l//の値に応じてどのような数列が展開されるか」

という問題にも拡張できるなど、コラッツの問題の類似問題の幅は広い。

11. ギャラリー
---------

Collatz_orbits_of_the_all_integers_up_to_1000.svg|最初の1000個の数字の数列を示す有向グラフ。 CollatzConjectureGraphMaxValues.jpg|x}}軸は初期値を表し、y}}軸は、1に至るまでの間での最大値を示す。このプロットでは、y}}軸の範囲を制限している：いくつかのx}}では、途中の最大値が2.7e7}}にもなる。 All Collatz sequences of a length inferior to 20.svg|20ステップ未満のすべてのツリー。

12. 参考文献
--------

13. 関連文献
--------

『数学100の問題』 日本評論社 978-4535606142}} 　「角谷の問題」として取り上げられている。
『数学の魔法の宝箱』 イアン・スチュアート  ソフトバンク クリエイティブ 978-4-7973-5982-4}}：「コラッツ＝シラキュース＝ウラム問題」として紹介されている。

14. 関連項目
--------

タグシステム
FRACTRAN

15. 外部リンク
---------

url=http://www.geocities.co.jp/HeartLand/2327/rada/korattsu.html コラッツ予想計算
[1] 　コラッツ予想計算、拡張問題にも対応可能。
[2]}}
[3] コラッツの問題を検証する為の分散コンピューティング

16. 脚注
------