コラッツ予想がとけたらいいな

**１３２人目の素数さん** · 2012/10/14(日) 10:32:39.71

525 名前：１３２人目の素数さん[sage] 投稿日：2012/09/03(月) 18:24:27.22
http://d.hatena.ne.jp/righ1113/
コラッツ予想について、証明を考えてみました。
ご指摘ご意見ご感想など、ぜひよろしくお願いします。

**righ1113** · 2013/01/18(金) 15:40:11.54

>>45
初期値1のコラッツ・パターンはそれで合っています。
これで良いのです。

「元のパターン」の右端は、『コラッツ操作が1に収束するまでは』傾きlog[2](3/2)の直線をとります。
でした。言葉足らずでした。

そしてこれが背理法の要です。
>>42の
> 二つの傾きを重ねると、あるところで交差・逆転します。
> すると、大小関係が逆転するので(ｲ)式と矛盾します。
> これを回避するには、二つの傾きが交差する前に、コラッツ操作が1に収束するしかありません。

は1行追加して、

二つの傾きを重ねると、あるところで交差・逆転します。
すると、大小関係が逆転するので(ｲ)式と矛盾します。
これを回避するには、二つの傾きが交差する前に、コラッツ操作が1に収束して、
『「元のパターン」の右端傾きが1になるしかありません。』
です。

**righ1113** · 2013/01/18(金) 15:43:42.91

こうなると大小関係が逆転して矛盾するから、
o　　　x
　o　　x
　　o　　x
　　　o　x
　　　　o　x
　　　　　ox
　　　　　　xo
こうなるしかありません。
o　　　x
　o　　x
　　o　　x
　　　o　x
　　　　o　xここで1に収束
　　　　　o　　x
　　　　　　o　　　x
　　　　　　　o　　　　x

**１３２人目の素数さん** · 2013/01/18(金) 17:34:29.63

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　このスレは馬と鹿と豚さんばかりね。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　　　
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

**１３２人目の素数さん** · 2013/01/19(土) 21:22:38.87

1だけが特別ってのも変だと思う。
他の数でも+1が上手く作用して、部分的に傾きが大きくなることもある。
27なんかはそれが多く起こって長く続くんだろう。

**righ1113** · 2013/01/20(日) 16:01:48.02

うーん変ですか……
1だけが特別というのは
1だけが4-2-1ループする唯一の奇数という
コラッツ予想の特徴をよく表していると思うのですが……

**１３２人目の素数さん** · 2013/01/29(火) 20:15:10.69

あ

**１３２人目の素数さん** · 2013/01/29(火) 22:53:37.56

傾きの大小だけで交わるかどうかを判断することはできない。
例えば、y=(x-1)/xはx>0の範囲で増加するが、xをどれだけ大きくしても1を超えることはない。
実際、「差のパターン」の差分を計算してみると、確かにlog[2](3/2)より大きくはなるがlog[2](3/2)に収束しそうな感じがした。

**righ1113** · 2013/01/31(木) 18:42:04.22

がーん　そんなあ。
確かにy=(x-1)/xの傾き1/x^2は、y=1の傾き0よりも大きいけど、
二つのグラフは交わりませんね。

ということは
「差のパターン」は「元のコラッツパターン」に漸近するんですかね？

**１３２人目の素数さん** · 2013/02/01(金) 12:18:46.16

だんだん平行になっていく、ぐらいじゃないかなあ
直感的には「差のパターン」は「元のコラッツパターン」よりずいぶん小さい気がする

**１３２人目の素数さん** · 2013/02/01(金) 21:31:47.57

数学板でちゃんとした投稿が大半を占めている貴重なスレなので
お礼にこのスレに関連する情報を１つ書いておくね。知ってたらゴメン。

コラッツ予想に関して、最近までの主要な結果（一般化も含め）や昔の論文の復刻を纏めた次の本が
２年ちょい前にアメリカ数学会(American Mathematical Society)から出版された。

Jeffery C. Lagarias "The Untimate Challenge: The 3x+1 Problem", ISBN: 978-0-8218-4940-8

AMSの会員ならAMSのホムペから少し安く買えるけど、会員以外でもアメリカのAmazonとかで買えるはず。
CoxeterやConwayやRichard Guyら錚々たる連中の昔の論文もリプリントされていて読めるので
コラッツ大好きな人はこの本を持ってるとちょっとハッピーだと思う。

**righ1113** · 2013/02/04(月) 13:42:19.67

>>54
うわああああああああ

>>55
情報ありがとうございます！
さっそく見てみます！

**１３２人目の素数さん** · 2013/02/10(日) 19:46:59.49

有理数への拡張を詳しく考えてみた

分母が奇数のものだけを考え、分子が奇数なら「3倍して1足して2で割る」、分子が偶数なら「2で割る」という操作をする。
例えば1/5から始めると、
1/5→4/5→2/5→1/5→…
とループする。

ここで、分子に奇数が現れたら○、偶数が現れたら×と書くことにすると、上記のループは「○××」が繰り返してることになる。
そこで、逆に「○××」で元に戻る数が他にあるかを考える。
初期値をxとすると、
x→(3x+1)/2→(3x+1)/4→(3x+1)/8
と変化するはずだから、方程式
(3x+1)/8=x
が得られ、これを解くと解はx=1/5のみ。
よって、「○××」が繰り返されるループは1/5を初期値とするもののみであることがわかる。

同様にして、○×からなる有限列を一つ指定すれば、それに応じて有理数が一つ定まる。
得られた有理数が実際に指定してループを辿るかは非自明だが、おそらく辿る。

**１３２人目の素数さん** · 2013/02/10(日) 19:50:35.62

「×○×」や「××○」を指定しても、初期値が変わるだけで「○××」と同じループが得られる。
そこで、これは「○××」が円形に並んでいるものと考えられる。
この考えから、次のことがわかる。

定理
n個の奇数とm個の偶数からなるループの個数は、n個の○とm個の×からなる円順列の個数に等しい。

例
奇数2個、偶数4個からなるループを考える。
○2個、×4個からなる円順列は、
○×　○×　○×
○×　××　××
××　○×　×○
(半時計周りに並んでいるとする)
の3通り。左上の○に対応する数は、それぞれ1/11,7/55,1/5となる。
これらは、実際に指定した偶奇を辿ってループすることが確かめられる。

奇数n個、偶数m個からなるループについて方程式を作ると、
(3^n*x+a)/2^(n+m)=x　(aは正の整数)
という形になる。これを解くと、
x=a/(2^(n+m)-3^n)
となる。したがって、ループを構成する数は、分母が2^(n+m)-3^n(の約数)であるような分数として書ける。

とりあえずここまで

**righ1113** · 2013/02/14(木) 18:59:31.92

http://cdn-ak.f.st-hatena.com/images/fotolife/r/righ1113/20130214/20130214184540.jpg

7から始まるコラッツ数列を強制的にループさせてみると、
3周目で「左端を伸ばすパターン」と「差のパターン」が交差した。
ここから何か言えないだろうか。

**１３２人目の素数さん** · 2013/02/14(木) 20:32:47.98

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　このスレには馬と鹿と豚さんしかいないのね。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　　　
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

**１３２人目の素数さん** · 2013/02/27(水) 21:26:07.29

http://science2.2ch.net/test/read.cgi/math/1072595845/
http://science6.2ch.net/test/read.cgi/math/1240289175/

**righ1113** · 2013/02/28(木) 22:02:47.37

コラッツ操作を式であらわすと
初期値をx、ステップ数をn、nステップ後の奇数をx_nとして
x(3/2)^n + 3^(n-1)/2^n + p_1*3^(n-2)/2^n-1 +...+ p_1*p_2*...p_(n-2)*3^1/2^2 + p_1*p_2*...p_(n-1)/2
　= p_1*p_2*...p_n*x_n
となる。
左辺第一項x(3/2)^nが「左端を伸ばすパターン」、
左辺残りが「差のパターン」、
右辺p_1*p_2*...p_n*x_nが「元のコラッツパターン」に対応している。

p_nはコラッツ操作で偶数が2回以上続いた時の2のべき。
「差のパターン」のp_nの積の部分は例えば
1,1,2,2,2,4,4,8,8,8,8,16...
のように増加していく。

ここから先に進まない……

**righ1113** · 2013/03/08(金) 21:56:02.23

「差のパターン」3^(n-1)/2^n + p_1*3^(n-2)/2^n-1 +...+ p_1*p_2*...p_(n-2)*3^1/2^2 + p_1*p_2*...p_(n-1)/2
をf(n)とおく。
「差のパターン」の右端グラフは2進数目盛上にあるので
log{f(n)}になる。底は2。

これの2階差分u_(n+1) - 2u_n + u_(n-1)を取ると、
log{f(n+1)} - 2log{f(n)} + log{f(n-1)}
= log{ f(n+1)*f(n-1) / f(n)^2 }
f<0だから、この式は正。
2階差分が正だから、「差のパターン」の右端グラフは下(左)に凸。

この結果と先にあげた
「左端を伸ばすパターン」右端傾きより「差のパターン」右端傾きが大きい　を合わせて
二つの右端グラフが交差する、と言えるのではないだろうか。

**righ1113** · 2013/03/08(金) 21:59:18.99

×　f<0
○　f>0

**righ1113** · 2013/03/09(土) 00:18:05.28

あ、やべ。
logの中が1より小さかったらlogは負になるんだ。
考え直します。

**righ1113** · 2013/03/10(日) 13:24:20.36

log{ f(n+1)*f(n-1) / f(n)^2 }の中が1より大きいか考える。
n=5でf(n)が5項ある場合を考える。nが増えても同様。

(a+b+c+d+e+f)*(a+b+c+d) / (a+b+c+d+e)^2
a+b+c+d+e=Aと置くと
(A+f)*(A-e) / A^2 = {A^2+A(f-e)-fe} / A^2

e = (1/2)*p_1*p_2*...p_(n-1)
f = (1/2)*p_1*p_2*...p_(n-1)*p_n
なので、p_n=1の場合、分子がA^2-feとなって分数式 < 1。
p_n=2の場合、f=2eだから分子がA^2+Ae-fe、A>fだから分数式 > 1。
p_n=4,8,16,...も同様に分数式 > 1。

よってlog{ f(n+1)*f(n-1) / f(n)^2 }は
p_n=1の場合、負。
p_n=2以上の場合、正。
「差のパターン」の右端グラフは
偶数が1回しかない時、上に凸で、偶数が2回以上の時、下に凸という、
なにも進展しない結果になった。

コラッツ操作で、偶数が2回以上続くほうが多い、とか言えないかなあ。。。

**righ1113** · 2013/03/14(木) 17:54:50.43

ネットを見てると、よく
「4x+1となる数については調べなくてよい」
とありますが、これは数学的帰納法の一部分であって、
「4x+1はコラッツ操作で1に収束する」が証明された訳ではないですよね。

何故こんな事を聞くかというと、4x+3もコラッツ操作で4x'+1に接続されるので、
予想が解けちゃうなー、と思ったので。

**１３２人目の素数さん** · 2013/03/18(月) 21:39:37.85

科学技術フェスタで、奇数の時3n+1にする代わりに3n-1にすると、ループが3種類できて
その3種類の頻度はほぼ同じになるっぽい（証明はできてない）、という発表を高校生がやっとった。

**１３２人目の素数さん** · 2013/04/05(金) 18:43:02.40

い

**righ1113** · 2013/04/15(月) 20:42:39.70

全ての数がコラッツ操作で小さくなれば、
コラッツ予想は証明されます。
値が2xの場合は2で割って小さくなります。
同様に値が4x+1や16x+3でもコラッツ操作を続けると小さくなります。
これを全ての数で言えないかを今考えています。

**righ1113** · 2013/04/19(金) 20:36:36.42

こんな数学的帰納法を考えています。

・x = 1で成り立つ
・x < kで成り立つと仮定したとき、
　x = kでも成り立つ
値が4x+1や16x+3の場合はコラッツ操作をおこなうと小さくなるので、
上記の方法が使えます。
すぐ小さくならない値k1については、
・x=k1とそれに連結される数を除くx<k2で成り立つと仮定したとき、
　x=k2でも成り立つ
　するとk1はk2に連結するのでx=k1でも成り立つ
という方法が考えられます。

図にするとこんな感じです。
　　k2
　／＼
／　　k3
k1

k2>k3となるk2は、全ての数が4x+1を通過することから
すぐ見つけられます。
あとはk1とk3がつながってループする事がなければ良いわけです。

**righ1113** · 2013/04/30(火) 18:13:18.50

あと、すぐ小さくならない値が複数個連結される可能性もあります。
図にするとこんな感じです。
　　　　　　　　　　　／＼
　　　　　　／＼...／　　　↓
　／＼...／　　　　　　　　｜
／　　　　　　　　　　　　　｜
↑................................................┘最後はループする

**righ1113** · 2013/04/30(火) 18:47:22.79

式であらわします。
すぐ小さくならないk1がすぐ小さくならないk2に連結されるとします。
k1の側は、例えばk1=27+2^5*xとおくと、2ステップ後にちょびっと小さくなるので、その値は
　((3*k1+1)/2*3+1)/4 = (9*k1+5)/8 = 31 +3^2*2^2*x1　　　－－－①
となります。

k2の側は、k2から後ろ向きにnステップ伸びるとして、
コラッツ逆操作「2のべき乗をかけて1引いて3で割る」をおこないます。
n=1の場合は
　(k2*2^p1 -1)/3
n=2の場合は
　((k2*2^p1 -1)/3*2^p2 -1)/3 = k2*2^(p1+p2)/3^2 -2^p2/3^2 -1/3
一般化して
　(1/3^n)( k2*2^(p1~pn) -2^(p2~pn) -3*2^(p3~pn) -...-3^(n-2)*2^pn -3^(n-1) )　　　－－－②
となります。（ p1+p2+...+pn を p1~pn と略記）

連結されるので、①と②をイコールで結びます。
　3^n*( 31 +3^2*2^2*x1 )
　　 = k2*2^(p1~pn) -2^(p2~pn) -3*2^(p3~pn) -...-3^(n-2)*2^pn -3^(n-1)

**righ1113** · 2013/04/30(火) 20:22:52.00

n=1の場合は
　2*47 +3^3*2^2*x1 = k2*2^p1 => p1=1となって
　　47 +3^3*2*x1 = k2
n=2の場合は
　3*(2*47 +3^3*2^2*x1) = k2*2^(p1+p2) -2^p2 => p2=1
　2*71 +3^4*2*x1 = k2*2^p1 => p1=1となって
　　71 +3^4*x1 = k2
第一項は31から始まるコラッツ数列になること、
pnはそのときの2で割る回数になるところがポイントです。
一般化すると、以下になります。　CO31 は31から始まるコラッツ数列です。
　CO31 +3^(n+2)*x1/2^(p1~p(n-2)) = k2

**righ1113** · 2013/04/30(火) 20:37:41.06

k2を固定します。ここではk2=71+2^7*x2とおいてみます。
　CO31 +3^(n+2)*x1/2^(p1~p(n-2)) = 71+2^7*x2
　2^(p1~p(n-2))*(CO31 -71) = 2^(7+p1~p(n-2))*x2 -3^(n+2)*x1

この式の形　c = 2^a * x1 - 3^b * x2　を考えることによって
何か言えるのではと思うわけです。

**righ1113** · 2013/05/07(火) 23:58:33.04

>>74の後半は間違っていました。
第一項が分数になることもあります。
コラッツ数列ぽいことはぽいのですが......

**righ1113** · 2013/05/15(水) 19:03:20.12

ここでちょっと内容を変えて、
コラッツ操作ですぐ小さくなる値とそうでない値をまとめます。
すぐ小さくなる値を「良い値」、そうでない値を「悪い値」と呼びましょう。
全ての数を2進数の下位ビットで場合分けして、良い/悪いを調べます。
（2進数は左が下位）

2進数　良い/悪い　どう小さくなるか
0…　良い　　2x　→　x
10…　良い　　4x+1　→　3x+1
1100…　良い　　16x+3　→　3^2x+2
1101…
　　11010…　良い　　32x+11　→　3^3x+10
　　11011…　★悪い　　27+2^5x

**righ1113** · 2013/05/15(水) 19:04:12.24

1110…
　　11101…　良い　　32x+23　→　3^3x+20
　　11100…
　　　　111000…
　　　　　　1110000…　良い　　2^7x+7　→　3^4x+5
　　　　　　1110001…　★悪い　　71+2^7x
　　　　111001…
　　　　　　1110010…
　　　　　　　　11100101…　★悪い　　167+2^8x
　　　　　　　　11100100…　良い　　2^8x+39　→　3^5x+38
　　　　　　1110011…　★悪い　　103+2^7x
1111…
　　11110…
　　　　111101…　★悪い　　47+2^6x
　　　　111100…
　　　　　　1111000…　良い　　2^7x+15　→　3^4x+10
　　　　　　1111001…
　　　　　　　　11110010…　良い　　2^8x+79　→　3^5x+76
　　　　　　　　11110011…　★悪い　　207+2^8x
　　11111…　★悪い　　31+2^5x

**righ1113** · 2013/05/15(水) 20:12:42.62

「11011…」「1110001…」「11100101…」「1110011…」
「111101…」「11110011…」「11111…」が悪い値ですが、
「1110001…」と「11100101…」は、次ステップで「11011…」になるので
考えなくて良いです。

よって、「11011…」「1110011…」
「111101…」「11110011…」「11111…」
の五つの場合を考えれば良い事になります。

**１３２人目の素数さん** · 2013/05/29(水) 21:16:17.31

う

**righ1113** · 2013/06/03(月) 02:29:13.29

一つ定理が出来たので書きます。

コラッツ操作でxがsステップで小さくなれば、
3^s < 2^lを満たすx+2^l*yもsステップで小さくなる。

コラッツ操作で奇数→奇数までを1ステップと数えます。
証明は、まず、xがsステップで小さくなれば、その時の2で割った合計をl0と置くとき、
3^s < 2^l0が成り立つ事を言います。―――①
sステップ後の数は、それぞれのステップで2で割った回数をpnとすると、
(3^s/2^(p1~ps))*x + 3^(s-1)/2^(p1~ps) + 3^(s-2)/2^(p2~ps) + ... + 1/2^ps < x
(3^s/2^(p1~ps))*x < x
3^s < 2^(p1~ps) = 2^l0　　　　よって①は成り立つ。

**righ1113** · 2013/06/03(月) 02:52:24.88

次に、3^s < 2^lを条件（②）としたx+2^l*yもsステップで小さくなる事を証明する。
①と②を重ねると次の三つの場合がある。
　　2^l0 < 2^l, 2^l0 = 2^l, 2^l0 > 2^l

2^l0 < 2^lの場合、xのsステップ後をx1とおいて、x+2^l*y―(sステップ後)→x1+3^s*2^(l-l0)y
　x > x1, 2^(l-(l-l0))*y > 3^s*yだから、x+2^l*y > x1+3^s*2^(l-l0)yが成り立つ。

2^l0 = 2^lの場合、x+2^l*y―(sステップ後)→x1+3^s*y
　x > x1, 2^l*y > 3^s*yだから、x+2^l*y > x1+3^s*yが成り立つ。

2^l0 > 2^lの場合、x+2^l*y―(sステップ後)→2^(l0-l)*x1+3^s*y
　2^l*y > 3^s*y。x > 2^(l0-l)*x1を言いたいので、次の補題を考える。

3^s < 2^l が成り立てば、xはsステップ後、計l回2で割った時点で小さくなる。―――③

**righ1113** · 2013/06/03(月) 03:09:29.21

③を証明する。コラッツパターンを使う。
http://cdn-ak.f.st-hatena.com/images/fotolife/r/righ1113/20130603/20130603020241.jpg

コラッツパターンより、
　　log[2]x > log[2]x -(l0-s) + s*log[2](3/2)―――④
　　log[2]x > log[2]x+1 -(l0-s) + s*log[2](3/2)―――⑤
が言えればよい。

3^s < 2^l0　->　s*log3 < l0*log2　->　log[2]3 < l0/s　とlog[2]3 -1 = log[2](3/2)から、
l0/s -1 > log[2](3/2)　->　l0-s > s*log[2](3/2)
log[2]x > log[2]x -(l0-s) + s*log[2](3/2)　　　④が言えた。

④が言えたから条件3^s < 2^lより、log[2]x > log[2]x -(l-s) + s*log[2](3/2)が言える。
これとl0-1 ≧ lより、log[2]x > log[2]x+1 -(l0-s) + s*log[2](3/2)　　　⑤が言えた。

④⑤が証明できたので、③も成り立ち、
2^l0 > 2^lの場合も　x+2^l*y > 2^(l0-l)*x1+3^s*yが成り立つ。
以上です。

**１３２人目の素数さん** · 2013/06/13(木) 16:41:12.70

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　このスレは馬と鹿と豚さんばかりね。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　　　
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

**１３２人目の素数さん** · 2013/06/13(木) 19:17:37.40

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　このスレは馬と鹿と豚さんばかりね。
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　　　
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

**righ1113** · 2013/07/07(日) NY:AN:NY.AN

③の証明に不備があったので修正します、③を少し弱めます。
　　3^s < 2^l < 2^l0が成り立てば、xはsステップ後、
　　計l回2で割った時点で同じか小さくなる。―――⑥
これを証明します。
logの底は2です。

sステップ後のコラッツ値をxsとおいて対数を取ると
以下が成り立ちます。
　　log(xs) = log(x) -(l-s) +s*log(3/2) +log(1+1/3x)…(1+1/3x_s-1)

>>83より
　　log(x) > log(x) -(l-s) + s*log(3/2)―――④
　　log(x) > log(x)+1 -(l-s) + s*log(3/2)―――⑤
が成り立ちます。切り上げて
　　[log(x)] ≧ [log(x) -(l-s) + s*log(3/2)]―――⑦
　　[log(x)] ≧ [log(x)+1 -(l-s) + s*log(3/2)]―――⑧

コラッツパターンより
　　[log(x) -(l-s) + s*log(3/2)] = [log(x) -(l-s) +s*log(3/2) +log(1+1/3x)…(1+1/3x_s-1)]
　　[log(x) -(l-s) + s*log(3/2)] +1 = [log(x) -(l-s) +s*log(3/2) +log(1+1/3x)…(1+1/3x_s-1)]

⑦⑧に代入
　　[log(x)] ≧ [log(x) -(l-s) +s*log(3/2) +log(1+1/3x)…(1+1/3x_s-1)]

切り上げを外して
　　log(x) ≧ log(x) -(l-s) +s*log(3/2) +log(1+1/3x)…(1+1/3x_s-1)

以上で⑥が証明できました。

**righ1113** · 2013/07/10(水) NY:AN:NY.AN

話題がコロコロ変わってすみません。おさらいです。
コラッツ予想を2進数で考えます。
コラッツ数列の奇数のみを並べると以下のようなパターンができます。
これをコラッツパターンと名付けます。
（下位ビットが左）
111　　　　7
1101　　　　11
10001　　　　17
01011　　　　13
000101　　　　5
0000001　　　　1

コラッツパターンは以下のルールで下へ伸びていきます。
(1)「1」の塊は、次ステップで両端が離れる
　　「11」は「1001」に、「111」は「10101」になります。
(2)単独の「1」は、次ステップで「11」になる
(3)「11011」のような、次ステップで左「1」と右「1」が
重なる場合は、右(上位)へ繰り上がる
　　「11011」は次ステップで「1000101」になります。
(4)最後に、左端に+1する

sステップ後値の初期値0位置からの距離をLnとおきましょう。
例の5ステップ目はLn=7となります。

**righ1113** · 2013/07/10(水) NY:AN:NY.AN

次に、ルール(4)を削除したパターンを考えてみます。
左端がえんえんと左へ伸びていきます。
00000111
000010101
000111111
0010111101
01110110001
10100101011
左端を伸ばすパターンと名付けます。

sステップ後値の初期値0位置からの距離をLlとおきましょう。
例の5ステップ目はLl=6となります。

**righ1113** · 2013/07/13(土) NY:AN:NY.AN

例からも分かるように、Ln=Llの時とLn=Ll+1の時があります。
これ以外はないことを証明します。
コラッツパターンのルールより、1ステップ後の右端は+0 or +1です。
左端に+1しているから、Ln≧Llです。
なので、Ln=Llから1ステップ後Ln=Ll+1になるところを考えてみます。

http://cdn-ak.f.st-hatena.com/images/fotolife/r/righ1113/20130713/20130713165031.jpg
図より、Ln=Ll→Ln=Ll+1になることはありえるが、その次のステップで
Ln=Ll+1→Ln=Llになるので、2以上ずれることはありません。
Ln=Ll or Ln=Ll+1です。

**righ1113** · 2013/07/13(土) NY:AN:NY.AN

sステップ後値の初期値0位置からの距離Lnは、
コラッツパターンを2進数で書いているのでlog[2]の対数目盛と見なして、
sステップ後コラッツ値のlog[2]を取れば良いことになります。

コラッツ操作27→41を変形すると
　　(27*3+1)/2 = 41 = (27+1/3)*3/2 = 27*(1+1/(3*27))*3/2
logをとって
　　log41 = log27 +log(1+1/(3*27)) +log(3/2)
コラッツ操作41→31を変形すると
　　(41*3+1)/4 = 31 = (41+1/3)*3/4 = 41*(1+1/(3*41))*3/4
logをとって
　　log31 = log41 +log(1+1/(3*41)) +log(3/2) -log2
　　　　　　= log27 + log(1+1/(3*27))(1+1/(3*41)) +2*log(3/2) -log2
よって一般化するとLnは以下になります。
引き算の部分はコラッツパターンでは右によせているので消えます。
　　Ln = [log(x) +s*log(3/2) +log(1+1/3x)…(1+1/3x_s-1)]

左端を伸ばすパターンの式は、初期値に次々と3/2をかければ良いので
　　log( x*(3/2)^n ) = log(x) +s*log(3/2)
となります。切り上げてLl = [log(x) +s*log(3/2)]です。

**righ1113** · 2013/07/13(土) NY:AN:NY.AN

よって、>>89より最大でもLn=Ll+1なので、
　　[log(x) +s*log(3/2) +log(1+1/3x)…(1+1/3x_s-1)] = [log(x) +s*log(3/2)] +1
切り上げを外して
　　log(1+1/3x)…(1+1/3x_s-1) < 2
logをとって
　　(1+1/3x)…(1+1/3x_s-1) < 4
となります。

**righ1113** · 2013/07/13(土) NY:AN:NY.AN

もしコラッツ予想で4-2-1以外のループがあったら
　　(1+1/3x)…(1+1/3x_s-1)
の中のループ1周期の積をXとおいて
　　X*X*X*…
となりますが、X>1なので、いずれ
　　X*X*X*… > 4
となって
　　(1+1/3x)…(1+1/3x_s-1) < 4
と矛盾します。

よって
　　コラッツ予想で4-2-1以外のループは存在しない
ことが証明できました。

**１３２人目の素数さん** · 2013/07/19(金) NY:AN:NY.AN

>>89の画像の上から2番目の図で、
コラッツパターンは11、左端を伸ばすパターンは01
となることはないか？

**righ1113** · 2013/07/20(土) NY:AN:NY.AN

ちょっとまってください

**１３２人目の素数さん** · 2013/07/20(土) NY:AN:NY.AN

一通り検証したけど、そこ以外は間違いはなさそう
本質的に難しいのはここなのかも
直接修正できなくても、Ln-Llが上から抑えられさえすればおｋ

**righ1113** · 2013/07/22(月) NY:AN:NY.AN

>>93
>>95
ありがとうございます。

指摘の部分ですが、すぐにできそうにないです。
>>89の画像の上から2番目の図で、
コラッツパターンが0011、左端を伸ばすパターンは**01の時は、
同じように次々ステップでずれはなくなります。
コラッツパターンが0111、1011、1111の時はどうしよう……

**righ1113** · 2013/08/03(土) NY:AN:NY.AN

修正できました。流れは以下です。

初めてコラッツパターンと左端を伸ばすパターンがずれる所を考える
ずれるステップをsとおく
　　↓
s-1,s-2,s-3のずれはない
　　↓
特定のパターン（2つ）しかあらわれない
　　↓
その特定のパターンはsでずれて、s+1,s+2,s+3ではずれない
　　↓
次にずれる時s2も、s2-1,s2-2,s2-3のずれはない

**righ1113** · 2013/08/03(土) NY:AN:NY.AN

特定のパターン1つ目です。
　　　　コラッツパターン　左端を伸ばすパターン
s-3  　**1　　　　　  　***1　　　　　　　　
s-2  　*11　　　　　  　**11　　　　　　　　
s-1  　0101　　　　　　　　**101　　　　　　　　　　
s　   　0000[1]　　　　　*1111　　　　　　　　　　
s+1  　**011　　　　  　101101　　　　　　　
s+2  　**1001　          　***0001　　　　      
s+3  　*11011　          　*****11　　　　      

特定のパターン2つ目です。s+2でまたずれるのでs'と置きなおしています。
　　　　　　コラッツパターン　左端を伸ばすパターン
s-3　　     　**1　　　　　  　***1　　　　　　　　
s-2　　     　*11　　　　　  　**11　　　　　　　　
s-1　　     　0101　　　　　　　　*1001　　　　　　　　　　
s　　　     　0000[1]　　　　　11011　　　　　　　　　　
s+1　　     　**011　　　　  　***101　　　　　　　
s+2 -> s'　　**100[1]　        　**1111　　　　　　　
s+3 -> s'+1　*11011　          　*101101　　　　      
s+4 -> s'+2　**00101　　　  　*******1　　　　　　
s+5 -> s'+3　***1111　　　  　******11　　　　　　

Ln=Ll or Ln=Ll+1が言えます。

**righ1113** · 2013/08/03(土) NY:AN:NY.AN

ぐじゃりましたね。こうです。
http://cdn-ak.f.st-hatena.com/images/fotolife/r/righ1113/20130803/20130803214250.jpg

**righ1113** · 2013/08/20(火) NY:AN:NY.AN

無限大に発散するほう、いけるか……

**なんとなくな一考察** · 2013/08/22(木) NY:AN:NY.AN

わかったよ
証明できるかも

学校の先生に聞いてみる

**righ1113** · 2013/08/25(日) NY:AN:NY.AN

無限大に発散するほう、むずかしいお......

**righ1113** · 2013/09/16(月) 04:56:49.18

無限大に発散するほう、できました。
コラッツ値xsが無限大に発散するとします。
xs = x0 *3^s/2^l *(1+1/3x0)…(1+1/3x_s-1)
xs < x0 *(3/2)^s s<lなので

かっこの部分を考えます。
(1+1/3x0)…(1+1/3x_s-1)… > (1+1/3x0)…(1+1/(3x0*(3/2)^(s-1)))…
> 1 +1/3x0 +1/3x0(3/2) +…+1/3x0(3/2)^(s-1) +…　　　等比数列の和

1+1/x0 < (1+1/3x0)…(1+1/3x_s-1)…
左辺第二項を大きくして、イコールになるところをα0とおく
1+α0/x0 = (1+1/3x0)…(1+1/3x_s-1)…　　　　①
同様にx1からスタートして
1+α1/x1 = (1+1/3x1)…(1+1/3x_s-1)…　　　　②

**righ1113** · 2013/09/16(月) 04:58:38.08

②を①に代入
(1+1/3x0)(1+α1/x1) = 1+α0/x0
きれいにして
(3α0-1)x1 = α1(3x0+1)

x1=(3x0+1)/2^qを代入して
(3α0-1)(3x0+1) = 2^q * α1(3x0+1)
x0の部分が消えて
α1 = α0 * 3/2^q * (1-1/3α0)

xsが発散する
　→　q=1が多い、αsが大きくなるとかっこの効果も弱まる
　　→　αsも発散する

**righ1113** · 2013/09/16(月) 05:16:00.29

一方、無限大に発散するコラッツ列でx0を最小値にとれば、
x0からx∞まで等比数列で下から押えられる
と仮定する。　xs > x0*a^s　　　　　　　　　③

(1+1/3x0)…(1+1/3x_s-1)… < (1+1/3x0)…(1+1/(3x0*a^(s-1)))…
≒ 1 +1/3x0 +1/3x0a +…+1/3x0a^(s-1) +…
不等号が成り立つようにaを少し小さくする

(1+1/3x0)…(1+1/3x_s-1)… = 1+α0/x0 < 1+1/x0 *a/3(a-1)

コラッツ列をグラフにした時に、下に凸な点を考える。
x0より後でx0の次に小さいxm1で
1+αm1/xm1 < 1+1/xm1 *a/3(a-1)
xm1より後でxm1の次に小さいxm2で
1+αm2/xm2 < 1+1/xm2 *a/3(a-1)
このプロセスはいくらでも続けられるので、
αmは発散するが、a/3(a-1)は一定なので矛盾する。

**righ1113** · 2013/09/16(月) 05:53:15.99

あとは③を証明すれば良い。
傾きaは、x0 *a^s < xsを満たすので、
log(a) < log(xs/x0)/s = log(x0 *3^s/2^l *(1+1/3x0)…(1+1/3x_s-1) /x0)/s
= log3 -l/s +logA/s　　　　　　　　(1+1/3x0)…(1+1/3x_s-1) = Aとおく

コラッツパターンより　l-s < log(xs)なので、
l/s < log(xs)/s +1 = log(x0)/s +log3 -l/s +logA/s +1
l/s < log(x0)/2s +log(3)/2 +logA/2s +1/2

l/sはsが大きくなると小さくなるので、傾きaはsが大きくなると大きくなる。
（直線x0-xm1の傾きより、x0-xm2、x0-xm3…の傾きのほうが大きい。
　　aを直線x0-xm1の傾きにとれば、コラッツ値はそれより上にある。）

以上で
　　コラッツ予想で無限大に発散する数はない
ことが証明できました。

**righ1113** · 2013/09/19(木) 20:24:57.81

>>106の
l-s < log(xs)
は自明じゃなかったです。
修正します。

**righ1113** · 2013/09/29(日) 01:07:15.53

修正できました。証明したい補題は以下です。

無限大に発散するコラッツ列でx0を最小値にとれば、
x0からx∞まで等比数列で下から押えられる
　xs > x0*a^s　　　　　　　　　③

コラッツパターンにおいて、左端の傾きをd1、
右端の傾きをd2(=log1.5)とおきます。
sステップ後の左端までの距離l-sは
　　l-s = s*d1
で、左端から右端までの距離log(xs)は
　　log(xs) = log(x0) +s*d2 -s*d1
です。
s*d1≦log(x0)の区間では
　　s*d1≦log(x0) +sd1 -sd1
　　　　< log(x0) +sd2 -sd1　　∵ d1 < d2
　　l-s < log(xs)
が成り立ちます。

**righ1113** · 2013/09/29(日) 01:09:27.19

s0 < s < sgまでxs > x0*a^sが成り立つ事がわかりました。
あとはこれをs∞まで広げれば良いわけです。

s0からsgの間に傾きaを上回る二点xf,xf+1が存在する
事が言えます。もしなかったらコラッツ値のグラフが傾きa直線とぶつかって
矛盾するからです。

xf,xf+1でも同様に
xs > xf * b^t　　　s < sh　　　が言えます。変形して
　　> x0 * a^f * b^t
　　> x0 * a^(f+t)　　∵ a<b

成り立つ区間がs < sgからs < sg < shにのびました。
このプロセスを繰り返せばs∞までxs > x0*a^sが言えます。
③が証明できました。

**righ1113** · 2013/10/06(日) 03:19:16.82

>>109の
s0からsgの間に傾きaを上回る二点xf,xf+1が存在する
をちゃんと説明すると、

http://cdn-ak.f.st-hatena.com/images/fotolife/r/righ1113/20131006/20131006025445.jpg

図のように
s0からsgの間に、コラッツ値x0,x1,x2,x3をとる
x1,x2,x3とx0との傾きは、sが大きくなるに従って大きくなるので、
　　x0-x1傾き < x0-x2傾き < x0-x3傾き
　　⇒ x0-x1傾き < x2-x3傾き
よって
　　s0 < s < sgまでxs > x0*a^s
　　s < sh　　　　　xs > xf * b^t
において
　　a < b
が言えます。

**１３２人目の素数さん** · 2013/10/22(火) 04:16:03.55

え

**１３２人目の素数さん** · 2013/11/06(水) 01:58:23.54

お

**１３２人目の素数さん** · 2013/11/19(火) 23:24:04.86

か

**１３２人目の素数さん** · 2013/12/03(火) 00:35:32.97

き

**１３２人目の素数さん** · 2013/12/21(土) 05:55:12.66

9232ってどうして頻発するの？

**１３２人目の素数さん** · 2014/01/04(土) 23:57:48.91

く

**righ1113** · 2014/01/05(日) 00:07:19.41

>>115
コラッツ操作で9232を通過する数がどうして多いか、ということですか。
自分は偶数を省いてやってたので気づきませんでした。
調べてみます。

**righ1113** · 2014/01/15(水) 21:11:20.77

うーんわかんないです

**１３２人目の素数さん** · 2014/01/16(木) 00:03:14.86

そうですか
寝てください

**righ1113** · 2014/01/18(土) 00:53:38.62

xを最大値に持つ数の個数をコンピュータで調べました。
500000くらいまで調べました。以下の事が分かりました。
・ほとんど0個
・1/5ぐらいで5個とか
・1/3000ぐらいで50個とか
・225988を最大値に持つ数は386個
・250504を最大値に持つ数は1759個
・560356を最大値に持つ数は500個
・575728を最大値に持つ数は550個
・695464を最大値に持つ数は612個
9232の1579個を超える数も見つかりました。
ごくまれに大きい個数が出てくるみたいです。

なぜこんなに偏っているのかは謎です……

**righ1113** · 2014/01/18(土) 00:54:28.34

ソースコードです。Haskellです。
Prelude> :l tree
*CTree> map colmaxcnt [1..100]
のように使います。

-- tree.hs start
module CTree where
data CTree = Leaf Int | Node CTree Int CTree deriving (Eq,Show)
collatz :: Int -> Int
collatz 1 = 1
collatz x | odd x　　 = x * 3 + 1
　　　　　| otherwise = x `div` 2
-- 木を引数まで成長させる
growm :: Int -> CTree -> CTree
growm _ (Leaf 0) = Leaf 0
growm y (Leaf x) | x > y 　　　 = Leaf 0
　　　　　　　　 | even(x)&&((mod (x-1) 3)==0)
　　　　　　　　　　　　　　　　= (Node (Leaf $ div (x-1) 3) x (Leaf $ x*2))
　　　　　　　　 | otherwise　　= (Node (Leaf 0) x (Leaf $ x*2))
growm y (Node t1 x t2) = (Node (growm y t1) x (growm y t2))

**righ1113** · 2014/01/18(土) 00:55:03.05

flatten :: CTree -> [Int]
flatten (Leaf x) = [x]
flatten (Node t1 x t2) = flatten(t1) ++ [x] ++ flatten(t2)
-- 空でないリストから収束した値を返す
conver :: Eq a => [a] -> a
conver [x] = x
conver (x1:x2:xs) = if x1==x2 then x1 else conver (x2:xs)
-- 最終結果
colmaxcnt :: Int -> Int
colmaxcnt 4 = 2
colmaxcnt x = if (any (x<) col) then 0 else chk
　　where col = takeWhile (1/=) (iterate collatz x)
　　　　　chk = length $ filter (\x->x/=0) $ flatten $ conver $ iterate (growm x) (Leaf x)

**１３２人目の素数さん** · 2014/02/21(金) 13:25:22.02

9232 が目立ってみえるのは、単に小さい初期値でしか調べてないから

**righ1113** · 2014/02/25(火) 00:41:37.58

僕もそれがいいたかったんです。ほんとです。

**１３２人目の素数さん** · 2014/02/25(火) 04:21:25.45

位数1024の群みたいなものか

**１３２人目の素数さん** · 2014/02/28(金) 23:35:58.04

9232未満の6分の1程度が9232に恋をして散ってしまう、というのは特筆すべきことだと思うけど。

**righ1113** · 2014/03/10(月) 18:37:13.99

無限大の証明ですが、間違っていました(>_<)
αは有限値をとるみたいです。
>>103-110は無しでお願いします。
新しい証明を考え中です。

**righ1113** · 2014/04/01(火) 20:46:55.54

あ、でも>>87-99の
4-2-1以外のループは存在しない証明は自信あります。

**righ1113** · 2014/04/22(火) 11:07:09.03

【検証】コラッツの予想(1-1000)
http://r-2ch.com/t/math/1240289175/
（>>12と同じ）にあった割数列というのを調べている。
これで全ての3の倍数の奇数を表わせないかなと。

気になるレスを抜き出してみる。

**righ1113** · 2014/04/22(火) 11:11:01.21

106
4 年前
ざっと計算機を回してみた感じでは、

任意の３の倍数でない奇数xに対して3の倍数yが存在して、
x∈collatz_set(y)
が成り立ちそうに見える。
80 ID: 2009/05/15 20:59

**righ1113** · 2014/04/22(火) 11:15:34.73

108
4 年前
もし>>106が成り立てば、コラッツの予想は３の倍数だけ調べればいいってことになって、
コラッツ素数の概念にもそれなりに意味が出てくる。
個人的には、そうであって欲しいところだ。
80 ID: 2009/05/15 21:25

**righ1113** · 2014/04/22(火) 11:20:26.74

110
4 年前
3で割り切れない数を9で割った余りは、1,2,4,5,7,8のどれかだけど、これは

1*2^6≡1
2*2^5≡1
4*2^4≡1
5*2^1≡1
7*2^2≡1
8*2^3≡1 (mod 9)

のように、どれも2を適当な回数掛けることで、9で割ると1余る偶数にできる。
ここから1を引いて3で割れば3の倍数である奇数になる。
>>106の予想は正しい。
１３２人目の素数さん ID: 2009/05/16 00:38

**righ1113** · 2014/04/22(火) 11:23:45.07

285
4 年前
>>283
ありがとう！ここから何か出てこないかな・・

一つ、完全割数列→完全割数列に関しての操作を見つけた。
長さｎの完全割数列→長さｎ＋１の完全割数列

まず、長さｎの完全割数列を、初項に０をつけたｎ＋１型で表す。
長さｎの完全割数列でできる最終値を９で割ったあまりが・・
３　・・　[4,+2]or[1,-2]をつける
６　・・　[2,+2]or[3,-2]をつける
０　・・　[6,+2]or[5,-2]をつける
分かりづらいと思うので例を。

21≡3(mod 9)　21=[0,6]
このとき、[4,6+2]と[1,6-2]が存在する。
ちなみに、[0,1,…,]みたいに、２項目（本来の初項）が１か２のときは２で引けない。
このときは本来の初項に６を足した[0,7,…,]から考える。本来の初項に６の倍数を加減してもOKなので。
170 ID: 2009/07/22 22:04

**righ1113** · 2014/04/22(火) 11:27:19.74

289
4 年前
全ての完全割数列を列挙できるかはゴメン、証明してないや。
でも、ちょっとやればできる気がする。今度時間ができたとき確信を得てみるよ。

＞どの数から始めても、コラッツ数列は一意に定まるからね。
確かにそうだね。これ書いた時は、一意でないものがあればそいつは１４２１以外のループをもつのかと
なんとなく思っていたけど、割数列を基に１スタートで逆にたどっていくならば
どこかでループしてしまうような値にはたどり着かないもんね。
ループしてしまうならば１にたどり着かないのだから。
「完全割数列で全ての３の倍数の奇数を表せる」だけ分かれば良いか。
170 ID: 2009/07/22 23:40

**１３２人目の素数さん** · 2014/04/22(火) 21:50:24.56

この問題が長年解かれないのは解こうとするのがアマばかりなのも一因だと思うよ

**１３２人目の素数さん** · 2014/04/22(火) 23:23:06.27

フェルマー・ワイルズの定理みたいに、これが解けたら重要な予想が証明できるってことがあるといいんだけど。何かあるんだったっけ？

**１３２人目の素数さん** · 2014/04/23(水) 02:43:17.47

>>135
違うよ。もうアマしか残ってないだけだよ。

かつて簡単に解けるだろうとコラッツ予想に手を出して時間を浪費した、
数多の研究者の屍で山が築かれたから。

**１３２人目の素数さん** · 2014/04/23(水) 17:26:51.64

　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　／:::::; -‐''"　　　　　　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　　　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　　　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　
.　　　　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　
　　　　 | 　　　　l＾,人| 　` `-' 　　ゝ　 |　　　　　　　　
　　　 | 　　　　 ` -'＼　　　　　ｰ' 　人　　　　　　　　　　　私は死なないわよ。
　　　　|　　　　　　　 /(l 　　　＿＿／　ヽ、　　　　　　　　　　　でも最近一寸太ったかしら。
　　　 |　　　　　　　（:::::`‐-､__ 　|::::`､　　ﾋﾆﾆヽ、　　　　　　　　　　Windows ver.10 で　　　　
　　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　元の痩せた姿にしてよね。
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　ﾋﾆ二、　＼　　　　　　　　　　　　　
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼　 |　'､　　＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/

**righ1113** · 2014/04/24(木) 00:14:09.98

まず、>>133の各変換に名前をつけ、式であらわす。

　A:[4,2]　　　B:[1,-2]
　C:[2,2]　　　D:[3,-2]
　E:[6,2]　　　F:[5,-2]
そしてG:[+6]　　　H:[-6]

981[7,1...]⇒15[1,1...]⇒81[2,3,1...]のように、
まずHをおこなってC（A,E）をおこなう変換は頭にHをつけて
HA,HC,HEであらわす。

変換前のコラッツ値をxとおくと変換後は、
　HA:[4,2]　x/3-2　　　B:[1,-2]　x/3/2-1/2
　HC:[2,2]　x/3/4-3/4　　　D:[3,-2]　2x/3-1
　HE:[6,2]　4x/3-7　　　F:[5,-2]　8x/3-3
　G:[+6]　　　64x+21
となる。

◆注意点１
　HC:117[5,1...]⇒ﾅｼ[-1,1...]⇒9[2,1...]のように、
H後にマイナスになる場合がある。

◆注意点２
　HC:981[7,1...]⇒15[1,1...]⇒C:⇒81[2,3,1...]のように、
A,C,EはHA,HC,HEとしてもあらわれる。

**righ1113** · 2014/04/24(木) 00:38:42.01

各変換でどのような数があらわれるか見ていく。

B:[1,-2]　x/3/2-1/2は21+72xを【3+12x】にうつす。
（例　21/3/2-1/2=3）
HC:[2,2]　x/3/4-3/4は117+288xを【9+24x】にうつす。
D:[3,-2]　2x/3-1は69+72xを【45+48x】にうつす。
HA:[4,2]　x/3-2は213+288xを【69+96x】にうつす。
F:[5,-2]　8x/3-3は45+72xを【117+192x】にうつす。
HE:[6,2]　4x/3-7は309+288xを【405+384x】にうつす。
G:[+6]　64x+21は3+6xを【213+384x】にうつす。

【3+12x】【9+24x】【45+48x】【69+96x】【117+192x】【405+384x】【213+384x】
と割数列1項[6]であらわされる【21】を加えて、
全ての3の倍数の奇数は完全割数列で表わされる。

**righ1113** · 2014/05/25(日) 18:04:56.82

変換でマイナス値を経由するとか、あやしいところもあるので
>>133を書き換えます。

長さｎの完全割数列→長さｎ＋１の完全割数列
まず、長さｎの完全割数列を、初項に０をつけたｎ＋１型で表す。
長さｎの完全割数列でできる最終値を９で割ったあまりが・・
３　・・　A:[6,-4]orB:[1,-2]をつける
６　・・　C:[4,-4]orD:[3,-2]をつける
０　・・　E:[2,-4]orF:[5,-2]をつける
元の初項が負になる場合はあらかじめG:[+6]をおこなう。

**righ1113** · 2014/05/26(月) 22:00:52.78

>>141を証明します。
A～FとGの各変換で[3の倍数の奇数]を[3の倍数の奇数]に写す
事を証明します。

A:[6,-4]
3 mod 9かつ奇数から変換前の数xは　18t+3
さらに割数列の初項が4以下は変換できないから
　(18t+3)*3+1= 54t+10　tが奇数の場合を除外、
さらに3=[1,…]だからt=0も除いて、
　x=36t+3、t=1,2,3,…　これが変換前の数。

A[6,-4]の変換関数は
(((3x+1)*2^-4-1)/3*2^6-1)/3 = 4x/3-7。

変換後の数は4x/3-7 にx=36t+3 を代入して
48t-3= 45,93,…　は3の倍数の奇数である。

**righ1113** · 2014/06/15(日) 13:17:22.81

A書き直します。
A:[6,-4]
3 mod 9かつ奇数から変換前の数xは　18t+3
さらに3=[1,4]だからt=0も除いて、
　x=18t+3、t=1,2,3,…　これが変換前の数。

A[6,-4]の変換関数は
(((3x+1)*2^-4-1)/3*2^6-1)/3 = 4x/3-7。

変換後の数は4x/3-7 にx=18t+3 を代入して
24t-3= 21,45,69,…　は3の倍数の奇数である。

**righ1113** · 2014/06/15(日) 13:18:44.73

B:[1,-2]
3 mod 9かつ奇数から変換前の数xは　18t+3
さらに割数列の初項が4以下は変換できないから
　(18t+3)*3+1= 54t+10⇒27t+5　tが偶数の場合を除外、
　x=18(2t+1)+3、t=0,1,2,3,…　これが変換前の数。

B[1,-2]の変換関数は
(((3x+1)*2^-2-1)/3*2^1-1)/3 = x/6-1/2。

変換後の数はx/6-1/2 にx=18(2t+1)+3 を代入して
6t+3= 3,9,15,…　は3の倍数の奇数である。

**righ1113** · 2014/06/18(水) 18:21:46.37

C:[4,-4]
6 mod 9かつ奇数から変換前の数xは　x=9(2t+1)+6、t≧0

C[4,-4]の変換関数は
(((3x+1)*2^-4-1)/3*2^4-1)/3 = x/3-2。

変換後の数はx/3-2 にx=9(2t+1)+6 を代入して
3(2t+1)+0= 3,9,15,…　は3の倍数の奇数である。