数学建模社区-数学中国

标题: 3X+1问题(角谷猜想)程序测试数据 [打印本页]

---

作者: hmz20128    时间: 2012-9-15 18:24
标题: 3X+1问题(角谷猜想)程序测试数据
数海聚珠网【程序体验】中的26号程序是专门用来测试3X+1问题的，在没有测试之前，对于一些概念问题先做一个介绍，下文选自：《三思科学》电子杂志，创刊号，2001.07.01文章作者 异调 ！
3x+1问题
一、一个简单的问题
5 N/ D, A" M$ F; p7 S7 o9 h; u, I, ^( C当我们阅读数学史时，会有这样一种印象，数学家们首先研究简单的问题，然后研究越来越复杂的问题。经常性地，高深的数学问题是非常复杂的。只是为了理解问题，我们就得学习非常多的数学知识；而为了解决它，那就得用更复杂的数学知识了。就算我们在学校里的数学考试也是如此，最后一题经常被叫做“最后一大题”，“一大题”是说它表达复杂，里面还有一二三四的小题，要理解题意就得几分钟的时间。弄不好还理解错了，搞得整道题都白白做，被扣去许多分。

0 z$ H5 r, C; V, Q6 _可是数学里不只有这些吓人的“大题”——我是说，数学里还有吓人的“小题”。这样的“小题”理解起来非常容易，却让无数数学家大跌眼镜，怎么冥思苦想也不得其解。3x+1问题大概就是其中最著名而又最简单的一个。它简单到大概任何一个会除2和会乘3的人（比如说，没文化但是经常买菜的老奶奶）都能理解它的意思，但是困难得让数学家至今也没有找到好好对付它的方法。

任取一个自然数，如果它是偶数，我们就把它除以2，如果它是奇数，我们就把它乘3再加上1。在这样一个变换下，我们就得到了一个新的自然数。如果反复使用这个变换，我们就会得到一串自然数。
8 t0 U( q+ A/ c2 I
% J* S8 I* p; n9 W7 K比如说我们先取5，首先我们得到3*5+1=16，然后是16/2=8，接下去是4，2和1，由1我们又得到4，于是我们就陷在4→2→1这个循环中了。
0 I+ }8 Z! u4 m
再举个例子，最开始的数取7，我们得到下面的序列：7→22→11→34→17→52→26→13→40→20→10→5→16→8→4→2→1这次复杂了一点，但是我们最终还是陷在4→2→1这个循环中。
) Q  b( s: T* m% A3 ?9 w7 N' X1 ]3 o
随便取一个其他的自然数，对它进行这一系列的变换，或迟或早，你总会掉到4→2→1这个循环中，或者说，你总会得到1。已经有人对所有小于100*2^50=112589990684262400的自然数进行验算，无一例外。
6 }2 l- R  b/ P0 j
& k& _6 k: p2 p. s8 d那么，是否对于所有的自然数都是如此呢？

这看起来是个多么简单的问题啊！
/ D5 C* @2 D; _# `( g: A& j# b- x, `二、克格勃的阴谋
! m% U1 u0 m  L1 ^7 y2 U这个问题大约是在二十世纪五十年代被提出来的。在西方它常被称为西拉古斯(Syracuse)猜想，因为据说这个问题首先是在美国的西拉古斯大学被研究的；而在东方，这个问题由将它带到日本的日本数学家角谷静夫的名字命名，被称作角谷猜想。除此之外它还有着一大堆其他各种各样的名字，大概都和研究和传播它的数学家或者地点有关的：克拉兹(Collatz)问题，哈斯(Hasse)算法问题，乌拉姆(Ulam)问题等等。今天在数学文献里，大家就简单地把它称作“3x+1问题”。

角谷静夫在谈到这个猜想的历史时讲：“一个月里，耶鲁大学的所有人都着力于解决这个问题，毫无结果。同样的事情好象也在芝加哥大学发生了。有人猜想，这个问题是苏联克格勃的阴谋，目的是要阻碍美国数学的发展。”不过我对克格勃有如此远大的数学眼光表示怀疑。这种形式如此简单，解决起来却又如此困难的问题，实在是可遇而不可求。
( ?0 B) a2 w& }: {
数学家们已经发表了不少篇严肃的关于3x+1问题的数论论文，对这个问题进行了各方面的探讨，在后面我会对这些进展作一些介绍。可是这个问题的本身始终没有被解决，我们还是不知道，“到底是不是总会得到1？”
9 {8 r2 c- p4 s: d3 R# t7 Y2 E7 ?+ s
2 H% [% d  G% J+ o* c在1996年B. Thwaites悬赏1100英镑来解决这个问题。我写一下这个悬赏的文献：Thwaites, B. “Two Conjectures, or How to win ￡1100.”Math.Gaz. 80, 35-36, 1996，好在大家万一证出来时知道跑哪里去领奖。看在钱大爷的份上，3x+1问题于是又多了个名字，叫Thwaites猜想。
% [  G( x, M" Z( C( g, G- w
要是真的有这么一个自然数，对它反复作上面所说的变换，而我们永远也得不到1，那只可能有两种情况。

1)它掉到另一个有别于4→2→1的循环中去了。我们在后面可以看到，要是真存在这种情况，这样一个循环中的数字，和这个循环的长度，都会是非常巨大的；
2)不存在循环。也就是说，每次变换的结果都和以前所得到的所有结果不同。这样我们得到的结果就会越来越大（当然其中也有可能有暂时减小的现象，但是总趋势是所得的结果趋向无穷大）。

因为这是个形式上很简单的问题，要理解这个问题所需要的知识不超过小学三年级的水平，所以每一个数学爱好者都可以来碰碰运气，试试是不是能证明它。不过在这里我要提醒大家的是，已经有无数数学家和数学爱好者尝试过，其中不乏天才和世界上第一流的数学家，他们都没有成功。如果你在几小时内就找到了一个“证明”，那么把它一步一步地严格地写下来，看看是不是严密正确（我可以肯定它是错的，我这样的肯定要冒的危险绝不超过连续中十次彩票头奖的概率，既然我不买彩票，我就没道理不这么肯定:-)）。事实上，在互联网上已经有一些错误的“证明”。据说还有个数学爱好者跑到公证处去公证他的“证明”，生怕别人把他的好主意偷跑了。

* Z. I; m) C% c二十多年前，有人向伟大的数论学家保尔·厄尔多斯(Paul Erdos)介绍了这个问题，并且问他怎么看待现代数学对这问题无能为力的现象，厄尔多斯回答说：“数学还没有准备好来回答这样的问题。”
' e  m$ E1 X0 H0 |/ _' |5 s三、一些概念，一些纪录
& A. [2 h+ x2 h& D虽然证不出猜想，但是数学家们还是得到了许多很可能很有用的结论。让我们先来定义几个概念，然后再来介绍这些结论。
' \1 W- t% L3 i9 d
从一个自然数开始，用上面这个变换，我们可以计算出一串自然数的序列。为了形象起见，我们把这串数列叫做以最初用来开始计算的那个自然数命名的“航班”。
7 `5 |- [5 w0 `" U4 b; l. s8 {5 ~比如说，第6次航班就是 6→3→10→5→16→8→4→2→1
我们把一个航班里的最大数字，叫做这个航班的“最大飞行高度”。比如说，第6次航班的最大飞行高度就是16。我们把航班在数字1“着陆”之前的数字个数（最初的数字包含在内，但1不包含在内），叫做这个航班的“航程”（特别定义第1次航班的航程为0）。第6次航班的航程就是8。如果真有自然数在此变换下永远达不到1，那么这个航班的航程就是无穷了。
' E1 M% Y, k3 p, |: c* z' D: w
接下去的概念稍微有点复杂。我们把从起点开始（但不包括起点）连续的不小于起点的数字的个数，叫作“保持高度航程”。举一个例子来说明这个概念比较方便：
- m3 ?6 C" L9 D' R第11次航班是 11→34→17→52→26→13→40→20→10→5→16→8→4→2→1
我们看到从起点开始，34，17，52，26，13，40，20都不小于起点11，共有7个数字，所以第11次航班的保持高度航程为7。后面的航程中虽然还有数字16大于起始点11，但是它不被算在保持高度航程里了。一个最简单的推论就是，偶数次航班的保持高度航程总是0，因为开始就除以2，跌到较低的高度去了。
. D: h5 B; n: T6 s- Q1 g1 r9 g7 O% f  D
3 Z+ Y$ D: C2 s6 x9 ]  V为什么我们对一个航班的保持高度航程感兴趣？因为如果所有航班的保持高度航程都是有限的话，3x+1问题就成立了。让我们假设已知所有航班的保持高度航程都是有限的，用数学归纳法来证明3x+1问题，也就是所有的航班都在1上“着陆”。我们已经知道第1到第5航班都是在1上着陆的，现在假设对于所有小于n的数字k，第k次航班都在1上着陆，我们来看看第n次航班的情况：由于按假设它的保持高度航程是有限的，所以它迟早会降落在一个比n小的数字上——于是按归纳假设它就会降落在1上！

4 \0 `6 l# E2 M1 w& ^) n8 m我们可以对开始的30班航班列出一个相关数据表来：

3 O' M. N2 o5 j航班  航程  保持高度航程  最大飞行高度
! t/ K4 O6 r- \( v3 P" c1       0           0               1
1 ]& D+ g$ l6 Z, v+ o" m( E. D4 ], U2       1           0               2
1 i2 [' ^# S: N# U" q( b3       7           5              16
6 |6 }: u  Q7 j- V4       2           0               4
: Y$ J! K- a0 V3 d% Z% V5       5           2              16
$ d* [2 S0 u! D+ N6       8           0              16
0 I5 E$ O# F: p7      16          10              52
8      3           0               8
9      19           2              52
5 f$ E+ q" d5 r7 w6 L" `10      6           0              16
11     14           7              52
* n& W5 {+ J+ O3 m# k) }12      9           0              16
13      9           2              40
# O- _" E& b, K- c2 b' x6 \2 s14     17           0              52
15     17          10             160
. Y, Q% Q% t% f8 D16      4           0              16
" e3 P7 {# _9 C/ e" G17     12           2              52
: q4 H" s- B, ]' j; L: e18     20           0              52
19     20           5              88
6 \! `4 _# ~; l2 x, V9 v# s20      7           0              20
- g; c) b0 R; J# `& u/ j21      7           2              64
0 j3 X& D/ [- `5 u. e* V" x22     15           0              52
9 m" J- d6 z: ^2 ]; U% Y9 n* ?( T. [23     15           7             160
24     10           0              24
( y9 O  i7 L% R) }1 b5 M$ E25     23           2              88
  K& k7 h' w( L& h1 D26     10           0              40
27     111         95            9232
6 E* s: \: l; n+ x+ h' V28     18           0              52
29     18           2              88
" O/ n: i( U# V( A7 [- n; J1 a30     18           0             160
: [6 P$ ]9 o! H9 y
下面要说说几个记录。在上面我们已经说过，目前3x+1问题已经被检验到100*2^50=112589990684262400，都没有发现反例。这是葡萄牙阿弗罗(Aveiro)大学的Tomas Oliveira e Silva的工作，用了很巧妙的编程方法。他的主页在http://www.ieeta.pt/~tos/3x+1.html
! j* H, j4 ^1 l
$ i4 O, p# j/ Y# R9 y* ]& l如果一个航班的航程大于所有它前面的航班的航程，我们就把它叫作“航程纪录航班”，比方说第7航班，它的航程是16，比第1到6次航班的航程都长，所以第7航班是个航程纪录航班。今天我们已经知道的航程纪录航班有118个，航程最长的是2234047405400065次航班，它的航程是1871，这是Eric Roosendaal发现的，他有个个人网站http://personal.computrain.nl/eric/wondrous/，里面有各种各样关于3x+1问题的信息，下面的记录也都来自这个网站。
. p* |3 U: I2 r' i
; y- i" u2 R8 F- P同样的，如果一个航班的保持高度航程大于所有它前面的航班的保持高度航程，我们就把它叫作“保持高度航程纪录航班”，比方说从上面的表中我们看到第7航班也是个保持高度航程纪录航班。今天已知的保持高度航程纪录航班有30个，航程最长是1008932249296231次航班，它的保持高度航程是1445（纠正一下：1445是滑翔步数，保持高度航程应是1444）。
9 E0 |9 D! t3 P) o
最大飞行高度记录航班就是那些最大飞行高度记录大于所有它前面的航班的那些航班，现在已知的有76个，最大的是10709980568908647次航班，到达了350589187937078188831873920282244的高度。

对于一个固定航班N，考虑它在1着陆之前所作的变换，如果把其中除以2的变换称为“偶变换”并记为E(N)，而把乘以3再加1的变换称为“奇变换”并记为O(N)。数学家已经证明，O(N)/E(N)<log2/log3。我们注意到，对有些航班来说，O(N)/E(N)非常接近于log2/log3≈0.63092975……。有猜想认为它会越来越接近这个数字（也有相反的猜想，认为不会无限接近），所以大家为此设立了另一个纪录，就是这个比值比所有以前的航班更接近log2/log3的航班。这样的纪录不多，
+ W  ]: {9 }- J* A现在已知的有15个，其中最后一个是N=100759293214567，I(N)/P(N)≈0.604938。值得一提的是N=104899295810901231，它的这个比值还要更靠近，达到0.605413，但是我们不知道它是否是一个纪录，也就是说，我们不知道所有比它小的航班里，是否还有比这个比值更靠近log2/log3的。
/ J& i$ F, M3 D+ ]# h8 r
我们知道，对于任何p，总有至少一个航班，它的航程是p：2p→2p-1→2p-2→……→4→2→1
5 F) l2 u+ m5 ]- n( }但是一般并不需要这么大的航班，就可以达到航程p。在2000年有人提出要找到最小的航班号，使得它的航程恰好是2000。现在最好的纪录是第67457283406188652次航班，但谁都不知道这是不是最小的航程为2000的航班。

计算一个航班的算法是非常简单的——只要除2或乘3加1。但是为了检验大量的和航次巨大的航班，巧妙的编程方法是非常重要的。上面的那些纪录都是由几台类似于我们平时使用的那样的计算机得到的结果。但是如果没有好好地思考和编程，光是硬算，那么使用最先进的计算机恐怕也得不到这样的结果。

为了验证一个航班的确在1上着陆，并不一定需要把结果计算到1。如果你已经验证了所有航次小于n的航班都在1上着陆，那么对于第n次航班，你只要把结果计算到一个小于n的数m就可以了——我们已经验证过第m次航班在1上着陆。事实上，如果我们只要计算到一个以前的航班飞行时到达过的数值就可以了，当然这需要记住以前已经到达过的比较高的高度，这里也必须巧妙地编程使得这样的记忆所使用的内存比较少。

更重要的是使用数学方法去减少计算量。比如说，任何n=4k+1的航班最终都会飞到一个比n更小的高度。首先这是奇数，我们乘3加1得到12k+4，然后连除两次2，就有3k+1<n。所以我们没有必要费功夫去验证4k+1型的航班。另外偶数次航班第一次变换就被除以2，降低了高度，所以同样也不需专门验证。只用这样一个小技巧，我们就使计算量减少到原来的25%。
0 Z8 l4 ^0 D3 q" M0 T% W+ @/ }
1 g8 h" c- p9 C- B如果按照这样的思路下去，我们同样不需要考虑16k+3型的航班，只要考虑到前面的飞行记录：
( a  V7 Z5 e( y1 m- V9 z- Z! ~16k+3→48k+10→24k+5→72k+16→36k+8→18k+4→9k+2→……
而9k+2<16k+3。
8 K& s& ^' V( @
我们可以这样追踪下去，考虑256k+i型的航班，其中i取0到255，那么我们会发现我们需要考虑的类型只有i=27、31、47、63、71、91、103、111、127、155、159、167、191、207、223、231、239、251、255。这样我们要作的计算只有最初的8%不到。

% I5 t' w0 r: y$ K% o. `; {1 Y而Eric Roosendaal得到上面那些纪录的程序，是建立在对65536k+i型航班分析的基础上的，其中只有1729种航班需要真正的检验（只有原来计算量的2.6%）。他的程序还使用了其它的算术技巧，以及可以同时计算好几个航班。Tomas Oliveira e Silva进一步改进了这些技巧，从而使得他成为现在3x+1问题验证的世界纪录保持者（他的计算从1996年8月开始，到2000年4月结束，其间使用了两台133MHz和两台266MHz的DEC Alpha计算机）。Eric Roosendaal还在和其他人一起合作进行计算（包括再次验证以前的结果），如果你愿意加入这个研究项目的话，可以去访问上面给出的他的主页。
9 h0 B9 m, F' S/ l% u8 p6 O" D! X四、理论结果
只要稍微动一下脑筋，我们就知道3x+1问题和下面几个命题都是等价的：
5 p  }% K: n( l8 r3 Z5 i) f1)所有的航班的航程都有限；
2)所有的航班的保持高度航程都有限；
- q/ _+ `& ^7 O: i3)所有的航班中的偶变换的次数都有限；
  W# h' `' n& B, A6 E7 o( ~4)所有的航班中的奇变换的次数都有限；
) r# H$ i' S! ~, A) m  k5)所有的航班的保持高度航程中偶变换的次数都有限；
5)所有的航班的保持高度航程中奇变换的次数都有限。
# A: ^) g1 x9 s9 ^2 Z* P- }
" x1 E7 }% Y7 u& R) v) yR. Terra和C. Everett证明了，“几乎所有的航班都会下降到它的起始点以下”，也就是说“几乎所有的航班的保持高度航程都有限”。
1 F  }" }  w; a2 g这里的“几乎所有”是有确定的数学意义的，它是指：

——存在一个自然数n1，在所有小于n1的航班里，最多只可能有1/10的航班，它们的保持高度航程无限；
9 _9 k. j* [4 @, i7 ]4 @——存在一个自然数n2，它比上面的n1要大，在所有小于n2的航班里，最多只可能有1/100的航班，它们的保持高度航程无限；
% P  c3 I- v% i$ Y/ N( ~' }——存在一个自然数n3，它比上面的n2要大，在所有小于n3的航班里，最多只可能有1/1000的航班，它们的保持高度航程无限；
" i6 g' {* d# ~——等等等等……
7 G) C8 a: l' H- A+ Q& q* A7 r
% J; t' G# n8 o8 N1 z: O0 d) i% I( o这好象很接近证明“所有的航班的保持高度航程都有限”了，于是很接近证明猜想本身了。但是好好想想，这个结论只不过是说明保持高度航程无限的航班会越来越稀少罢了，它们还是有可能存在的……更糟糕的是，这个结论一点也没有排除有其它循环存在的可能。
( B. h* [, u+ T- u, I
) q, L* M6 ~. X% \& v) g6 m对于在1上着陆的航班，数学家们也得到了一些结果。他们证明了，存在一个常数c，当n足够大的时候，在比n小的航班中，能够在1上着陆的航班的个数大于等于nc。在1978年R. Crandal首先给出c=0.05，虽然小了点，但毕竟是开头一步；然后J. Sander给出c=0.3；在1989年I. Krasikov得到c=0.43；1993年G. Wirsching得到c=0.48；最后在1995年D. Applegate和J. Lagarias得到c=0.81。看起来我们越来越接近c=1这个最终目标了。可是我们不知道现在用来得到c的方法是否还可以再用下去，就好象在试图征服哥德巴赫猜想的过程中，陈景润用来证明1+2的方法，似乎不能用来证明1+1了。

2 {5 d$ W; f; t8 N1 R: w1995年的这个证明相当特殊。它使用了计算机程序来解一个十分巨大的方程组，所以这个证明不能用手工来验证。在论文中，我们看见的不是一个关于c=0.81的定理的证明，而是一个关于如何写出这个巨大方程组的说明，和由程序计算出来的结果，以及如何使用这些结果来解释c=0.81。其他的数学家如果想验证这个结果，必须首先看懂关于方程组的证明和那些解释，再按照里面的说明来写一个程序（很复杂的！），运行它，再看看结果是否和文章中的相同。目前四色定理的
* B9 I- L# Q0 m7 N证明也是如此，所以数学家对此很不满意。

& y7 b- A. c7 H' I6 }还有一些结果是关于如果有其他不同于4→2→1的循环存在时，对这样的循环的性质的研究。R. Crandal和N. Yoneda在1978年证明，如果这样一个另外的循环存在的话，那么它的长度（就是在这个循环中数字的个数，比如说循环4→2→1的长度就是3）一定要大于275000。1993年这个体积增大到17087915，最近的结果是102225496。这些结果是通过分析包括我们前面提到的各种纪录得到的，所以这些结果我们还是不能完全通过手工来验证。我们看到，如果真有另外的循环存在的话，那一定是非常非常巨大的！
/ x, R6 G7 i4 o# ^8 L  \五、启发式论证
数学中有一种叫“启发式”的论证方法，建立在估计和概率的手段上。比如说底下的论证方法就是这个类型的：

7 e/ ^( v! g( E- {* `" u4 `" J# X“每个数字要么是奇数要么是偶数，如果随便取一个自然数，碰到奇数和偶数的可能性是一样的。如果我们把一次航班中这一系列数值看作是随机的话，那么使用奇变换和偶变换的可能性也是一样的，所以平均在每两次变换中我们有一次是n→3n+1，有一次是n→n/2。所以平均起来，每次飞行高度的变化就是乘以3/2，于是……就会越飞越高。”
/ j" ~0 t2 H6 z( S$ @
这样的启发式论证就推翻了原来的猜想！但是这个论证显然比较幼稚，因为它没有考虑到，每一次奇变换后随即而来的一定是一次偶变换，因为如果n是奇数的话，3n+1一定是偶数；而每一次偶变换后随即而来的却不一定是一次奇变换。J. Lagarias改进了这个启发式论证。他指出，如果我们把奇变换后再作偶变换考虑在一起，那么这样得到的结果可以看作是真的“很随机”。于是有1/2的可能性它是奇数，有1/4的可能性是一个奇数的2倍，有1/8的可能性是一个奇数的4倍，等等。于是飞行高度的变化就是以下变换的“平均效应”；

* h' u, v% }3 R& ^1 C——n乘以3/2，这有1/2的可能（奇变换后再作偶变换的结果为奇数）；
: m! y; G/ l; Z% z7 `9 L——n乘以3/4，这有1/4的可能（奇变换后再作两次偶变换）；
1 t: Q3 [' ^7 P——n乘以3/8，这有1/8的可能（奇变换后再作三次偶变换）；
…………
0 _* ?- `+ C6 J3 j' p+ `
于是平均来讲，每次变换后高度的变化就是
7 w- `4 N3 Z, a2 |+ C0 F' Oc=(3/2)1/2(3/4)1/4(3/8)1/8(3/16)1/16……=3/4
- W3 r2 j( m0 m- ]. p% x* c# S4 K所以高度在总体上来说应该是越来越低，每次大约低25%，最终降到一个循环上（不过这个论证没有排除有除了4→2→1以外的其他循环）。这个论证可以使我们使用论证中的模型来计算出，从一个自然数开始，平均要多少步的这样的飞行（就是保持高度航程中奇变换的次数），可以使飞行高度降到起始点以下。理论上的数值是3.49265……。如果我们对3到2000000000（二十亿）之间的航班的保持高度航程中奇变换的次数取平均值，我们得到3.4926……。这两个结果惊人的一致
5 t, r$ `* H' B! ~3 e性使我们相信上面的启发性模型是正确的。如果它是正确的，那么就意味着没有保持高度航程无限的航班，于是3x+1猜想就是正确的，至少可以得出没有飞得越来越高的航班的结论。

可是一个启发性论证，就算再有实验证据来表明它是对的，也只不过是个论证，只能使我们对猜想的正确性更充满信心。它不能代替真正的数学证明。比如说，数学家猜想在π的十进位小数表示当中，出现0到9各个数字的可能性是一样的，对π的数值计算也强烈支持这个猜想，可是如果没有数学证明，它还是得被叫做一个猜想，而不是定理。
% |# n' j  A' l& @( Q
用上面这个启发式的概率模型，我们还可以预言，对于第n次航班，它的最大飞行高度不会超过Kn2（对于某个常数K）。数值计算表明对于K=8，这个公式是正确的（同样地，这可以让我们提出猜想，而不是证
8 T* d% X. l8 L5 P. N6 H3 [) z明定理）。
六、会不会永远证不出来？
* q' ?4 v8 S3 D自从哥德尔发表了他的著名的不完备性定理以来，每次碰到一个十分困难的问题时，数学家们就免不了疑神疑鬼——这会不会证不出来？
$ L  K" X6 N/ a5 N' _. n3 E+ `
哥德尔的不完备性定理说，在包含皮亚诺的自然数公理的数学公理系统中，总有不可证明的命题存在。公理系统的这种性质叫不完备性。比如说，如果我们只取欧氏几何的前四条公理，那么平行公理是不能用这前四条公理证明出来的，也就是说只有前四条公理的平面几何是不完备的（这个例子不是很严格，因为欧几里德的公理系统在现代观点下是不严密的，但是我举这个例子只是为了说明不完备性这个概念，所以关系不大）。

3 B5 u: |: h. C( G# ^所以说，如果我们只用皮亚诺的自然数公理，甚至再加上现代的集合论公理系统，也有可能不能证明3x+1问题。甚至即使3x+1猜想其实是错误的，我们也有可能不能证明这一点。比如说，我们可能发现一个航班，我们对它进行计算，发现它飞得越来越高，但是无论如何不能证明它永远也不会回到1上来。
/ q: D: Y; y# a
0 a* E5 o, Y$ i' [当然无论什么数论问题都有可能搞得数学家们这样疑神疑鬼，虽然其实是他们还没有发现证明。不过有一些蛛丝马迹表明我们有必要稍微严肃点看待此问题，因为3x+1问题离不可证明的问题并不太远。
5 O9 u; Z/ B# M1 p
' M3 t6 v( g9 N: E' X% H/ C) VJ. Conway（喜欢数学游戏的朋友可能会记起这个名字来，著名的生命游戏就是他发明的）在1972年考虑了3x+1问题的推广形式。在3x+1问题里，我们把数字除以2，然后得到了2种可能的余数（0或者1），按照余数我们使用2个公式（除以2或者乘3加1）。Conway考虑了除以一个固定的p，按照余数的不同（这时就有p种不同的余数）分别使用p个公式的情况。然后他提出了一个类似“在1着陆”的猜想。他在论文中证明了，这个猜想在集合论公理系统中是不可证的。
, D, _4 g* j, ^8 R3 @
' R# i$ }5 c. v事实上，在任何一个包含了皮亚诺的自然数公理的数学公理系统中，Conway的方法都可以定义一个类似于3x+1问题的不可证命题。当然这不是说有一个在所有公理系统中都不可证的命题。“不可证”总是相对于某公理系统而言的。当然，Conway的方法并没有说明3x+1问题本身是不可证的，也没有说它一定是很困难的（事实上有些3x+1问题的变种是很容易解决的），但是这毕竟说明，有些很象3x+1问题的命题是不可证的，这把事情搞得很可疑。

. v0 I5 G. [1 ?& P& r. W4 ]. j1993年，法国里尔(Lille)的基础信息实验室使用了Conway的方法来演示一套基于逻辑规则的编程形式的威力。同许多数学中的例子一样，先头看上去最没用的课题，会有很具体的用处。
七、各种变种
数学家总喜欢把问题推广，使它更抽象化和一般化，因为这样可以把一种在具体某个问题上使用的方法的威力应用到一般的情况上去，从而得到很有可能是出乎意料的结论。
4 a( ]+ i% i1 H/ D$ _
& W! `5 B/ w) x1 h# B3 D数学家们首先考虑，如果把3x+1问题的规则运用到负整数上去，会产生什么现象。他们发现了三个不同的循环：
1)-1→-2
2)-5→-14→-7→-20→-10
/ W, @- h. Y+ T2 m- G$ B# y3)-17→-50→-25→-74→-37→-110→-55→-164→-82→-41→-122
$ c0 i$ _) y9 q9 t  C→-61→-182→-91→-272→-136→-68→-34
9 h+ s" J2 |5 I% E他们猜想，这就是所有的循环，而所有的负整数都会掉进其中一个里。

他们还提出了5x+1问题，也就是在奇数的情况下用5x+1来取代3x+1。
! v! ], I6 F- K# s6 v4 i1 r这下又有好几个循环：
0 P6 b3 c7 W7 }1 R1)6→3→16→8→4→2→1
! F6 q' `8 i4 r; I  {- Y2)13→66→33→166→83→416→208→104→52→26
7 m' U! r- t" m# E. ?9 m3)17→86→43→216→108→54→27→136→68→34
但是5x+1问题中的第7次航班好象老在那里飞啊飞，怎么也跑不到一个循环里去，但是谁都不能证明的确如此。
# y. @* o# c6 t& r
上面Lagarias的那个启发式论证使得数学家猜想，如果q是大于3的奇数的话，对于qx+1问题，总存在至少一个航程无穷的航班，这看起来很象是一个“反3x+1问题”。

还有许多其他的3x+1问题的推广，一些结果把它们和其它数学领域联系起来，比如说素数理论，某些丢番图方程（求解系数为整数的方程的整数根，比如著名的费尔马大定理就是一个丢番图问题），马尔可夫链（概率论中的递归理论），遍历理论（一种关于函数混合递归的理论）。
$ @, |  c% G' S* n
: H: s  X& Q: a: y7 C# Q# j就算3x+1问题终于被解决了，看看所有这些变种，也够数学家们自娱自乐上几百年的了。

---

作者: hmz20128    时间: 2012-9-15 18:26


: ?+ {: r5 d9 G数海聚珠网【程序体验】中的26号程序运行方式分两种类型：
" ~' l; G6 x" U9 @) _' i; g, W类型1，运行某一个数（某个航班）归一的过程及最后输出这个航班的航程（归一步数）i、保持高度航程b（滑翔步数=b+1）、最大飞行高度(最高点)z、奇变换j、偶变换u。
类型2，运行某一个数列（一列航班）的航程（归一步数）i、保持高度航程b（滑翔步数=b+1）、最大飞行高度(最高点)z、奇变换j、偶变换u。
在类型2中还可以用i、b、z、j、u中的某一项排序输出。
例：选择类型1的选项，输入k=3，o=1，v=1（最后停止的数），C=2234047405400065，D=1，b=1，r=0（航班是自然数型，C输入航班，则D=1、b=1，r=0）
" F* g* F3 @6 V输出：
: t" {# V) v  J" ]T=1017299500769
$ ?; p7 V: x) u% J7 e% S) h" K3 KT=3051898502308
, `/ Y% B" l- RT=1525949251154
T=762974625577
+ D0 W. L$ s" x) \& Q! \1 rT=2288923876732
8 c: J: ~8 q- t  t$ sT=1144461938366
3 O' h" U5 k/ \$ eT=572230969183
T=1716692907550
T=858346453775
……
T=155413970848165054
T=77706985424082527
T=233120956272247582
7 ?# G- E, Z4 a+ y8 @' TT=116560478136123791
+ P; W8 [; q5 Q- A" QT=349681434408371374
- l0 L% e0 P' r$ H) V/ J- i9 IT=174840717204185687
T=524522151612557062
T=262261075806278531
T=786783227418835594
T=393391613709417797
7 w; i! _) X# AT=1180174841128253392
……
) K5 Q5 S# d  z2 J. [T=160
T=80
T=40
T=20
3 L% ]) I% k+ O+ vT=10
; W4 G8 b" `* b  k( g1 d  |T=5
+ e$ Y: p4 Q( n& K8 m; e/ |* [) KT=16
T=8
T=4
T=2
8 J+ a% I! `) yT=1
, }5 _: [' y5 V  i=1871
3 t( t  H& g, b  b=2
$ O% Y) K  g. [9 Y9 z6 r* A+ d6 h3 K' O  z=1180174841128253392
  u=1167，j=704
, V! q- E0 A9 R* X! g输出的结果说明，加上省略的部分到T=1是航班2234047405400065的归一过程，下面的 i=1871表示航程，b=2表示保持高度航程，z=1180174841128253392表示的是最大飞行高度，u=1167表示偶变换次数，j=704表示奇变换次数。

---

作者: hmz20128    时间: 2012-9-21 19:28
一楼中的介绍：如果一个航班的保持高度航程大于所有它前面的航班的保持高度航程，我们就把它叫作“保持高度航程纪录航班”，比方说从上面的表中我们看到第7航班也是个保持高度航程纪录航班。今天已知的保持高度航程纪录航班有30个，航程最长是1008932249296231次航班，它的保持高度航程是1444（滑翔步数是1445）
7 F, y9 R  }/ ~9 f5 V用26号程序。输入：k=3，o=1，v=1（最后停止的数），C=1008932249296231，D=1，b=1，r=0（航班是自然数型，C输入航班，则D=1、b=1，r=0）
; B) X, [- @4 h输出：
6 R/ ^1 C) q; m/ z( ]T=1008932249296231
T=3026796747888694
& d5 }8 Y) g3 \) G' ?T=1513398373944347
, F# `  ?; m% ]6 V0 eT=4540195121833042
5 n- l0 X' k; J* t, ]T=2270097560916521
T=6810292682749564
T=3405146341374782
+ n# V4 f, D- }4 GT=1702573170687391
! S; {2 l7 J- j$ O* r7 h8 aT=5107719512062174
T=2553859756031087
1 [% H! c- s0 V5 m# RT=7661579268093262
' U" m: G1 s* Z: @" w……
T=3630702040938956164751
3 }/ \, }' m. u( I/ {$ E. Y- h4 i* QT=10892106122816868494254
9 ?* I6 F9 R: L7 |3 QT=5446053061408434247127
T=16338159184225302741382
  I# p; ?) a( o& z5 NT=8169079592112651370691
7 X$ K( B  Z. M  y$ U$ E' f2 B9 cT=24507238776337954112074
T=12253619388168977056037
. k9 L! B; J* c' [- _T=36760858164506931168112
- p/ n$ |& c5 i; i……
, H# c5 _' N& {5 ?T=106
T=53
; [: B3 y4 K6 a# P8 \  R; p: G% OT=160
, F4 ]- X$ D! a! eT=80
" p- S. t/ R; j$ _. A; g* G) AT=40
1 b' k4 x7 W8 YT=20
# y4 |; D& S8 p& c: K) s6 hT=10
$ D& D1 m6 L5 w0 OT=5
T=16
. Y% i  p% r5 m$ C) h- a; YT=8
T=4
- x$ |% b- d' c5 AT=2
T=1
  i=1831
  b=1444
  z=36760858164506931168112
  u=1142,j=689
1 v1 @( u1 l& ?: V! y" d( F. }输出的结果说明，加上省略的部分到T=1是航班1008932249296231的归一过程，下面的 i=1831表示航程，b=1444表示保持高度航程，z=36760858164506931168112表示的是最大飞行高度，u=1142表示偶变换次数，j=689表示奇变换次数。

---

作者: hmz20128    时间: 2012-9-21 19:29
一楼中的介绍：最大飞行高度记录航班就是那些最大飞行高度记录大于所有它前面的航班的那些航班，现在已知的有76个，最大的是10709980568908647次航班，到达了350589187937078188831873920282244的高度。
用26号程序。输入：k=3，o=1，v=1（最后停止的数），C=10709980568908647，D=1，b=1，r=0（航班是自然数型，C输入航班，则D=1、b=1，r=0）
输出：
8 y- ?  Q+ J; t$ z- dT=10709980568908647
& T" t/ O  o  k. l0 G8 LT=32129941706725942
T=16064970853362971
8 B* t/ M3 h0 m3 j2 rT=48194912560088914
T=24097456280044457
9 u0 N9 Q8 u; S7 `6 w" vT=72292368840133372
T=36146184420066686
; G2 P; i  V% j$ S% V# sT=18073092210033343
……
% d2 l& I' h* X+ }& r# NT=69252185271521617547036823759454
  `6 n- q; v  v4 {T=34626092635760808773518411879727
T=103878277907282426320555235639182
T=51939138953641213160277617819591
2 b8 t- e3 V$ |, E7 J( YT=155817416860923639480832853458774
0 ?* N/ d! t' d8 e* \. O1 T* YT=77908708430461819740416426729387
: b0 `6 v9 k. U8 U/ g6 K8 J$ Q4 jT=233726125291385459221249280188162
T=116863062645692729610624640094081
T=350589187937078188831873920282244
……
* C% v0 @' ~0 B4 J6 w5 ZT=106
5 \( c9 M) B/ S" b: y2 s7 qT=53
T=160
T=80
T=40
# b/ l1 v* D* l# u8 p; qT=20
T=10
" p# j' W2 C  k5 E3 r  ?4 v/ mT=5
& e  t) k, C, a1 BT=16
T=8
) m- ^! m2 n# |6 }* FT=4
$ ]9 H, D2 `4 E: D6 z% N$ Y, LT=2
T=1
  i=1165
* f% L) J) b% m: q6 j* n( J  b=790
  z=350589187937078188831873920282244
. g  X3 J, g- X  u=735,j=430
输出的结果说明，所有的T加上省略的部分到T=1是航班10709980568908647的归一过程，下面的 i=1165表示航程，b=790表示保持高度航程，z=350589187937078188831873920282244表示的是最大飞行高度，u=735表示偶变换次数，j=430表示奇变换次数。

---

作者: hmz20128    时间: 2012-9-25 20:38
关于3X+1问题的滑翔步数、归一步数、与奇变换和偶变换有关的剩余因子、完整度、最高点等数据记录，国外的一些爱好者都有深入的研究，目前介绍最详细并且被翻译成中文的网址是：http://www.equn.com/3x+1/;WkkY*，里面归纳的数据非常完整，介绍的非常详细，值得还在继续研究3X+1的爱好者一看。

数海聚珠网【程序体验】是专门为数论爱好者提供的免费测试工具，其中26程序是专门为3X+1爱好者提供的运行测试某一个航班的归一过程和一列航班的航程（归一步数）i、保持高度航程b（滑翔步数=b+1）、最大飞行高度(最高点)z、奇变换j、偶变换u的数据的，有兴趣的也可用26号程序寻找目前还没有填充的空白。

---

作者: hmz20128    时间: 2012-10-17 23:24
用26号程序运行某一范围等差级数110X+1航程（归一步数）i、保持高度航程b（滑翔步数=b+1）、最大飞行高度(最高点)z、奇变换j、偶变换u
6 G" V: g$ {! C( `2 \输入：K=3，o=1，v=1，h=1，L=111（后面两空都填1），K1=110，r=1
& f: b+ ^& R# b" v8 x6 u" ^选择第二选项T=K1*b+r
输出：
T=111
  i=69
  b=49
  z=9232
  u=45,j=24
9 }( y# O* c2 j. Y" T4 FT=221
6 ?/ V3 s/ U* r( G! k  i=114
8 L3 p& Z; K7 Z' C! V. j  b=2
3 N! c2 i4 k1 i/ a: O  z=9232
  u=73,j=41
" U3 s; Q2 S. L! L/ g. Q& XT=331
  i=24
% l5 a# F. h" c6 E  b=7
  z=1492
  u=18,j=6
9 h, ]" x2 B( x# ~+ @/ vT=441
) I5 i8 \+ d- _. w" M  i=27
  b=2
: ]4 x% r! i2 B- I" o  z=1492
  u=20,j=7
6 A# z- z8 S6 X$ S3 ]$ z# w2 kT=551
" X  G3 K2 I& `9 Y0 m# c4 N) C  i=43
( z! `) \+ H* v6 W4 }  e  b=12
, {( ?" b8 v+ V. E# u- N. A  z=4192
7 t2 X7 }, |  _" Q4 W, Y  u=30,j=13
$ Z, {7 N* q5 `: ?, ~* W4 A! UT=661
" E, w) @1 _+ N0 [" _  i=113
  b=2
  z=9232
" S; ]- I9 T' ~7 S" h$ a. F  u=73,j=40
. p6 [. j  C  P6 W2 ]T=771
: _! Z4 u% w& ]6 Z3 v8 a3 I  i=33
, ?2 m1 @$ S) F' ^, |- T9 \  b=5
+ @9 _$ l( V+ B  {$ [  z=3472
  u=24,j=9
T=881
  i=116
0 Y' |! j5 o! }7 e9 H. I/ k; f( h  b=2
  z=9232
  u=75,j=41
T=991
  i=98
  b=41
! G# k8 w" W! J( x! J  z=15064
  u=64,j=34
/ J2 x5 L7 s  a6 \) ^+ D7 J* d$ ?T=1101
  e6 H: C4 k; l/ U+ p  i=93
  b=2
  z=9232
/ I' w  R% _$ i& ~* P' x  u=61,j=32
4 r9 a# Z9 ~+ k: I: }/ b……
  E0 S9 M2 ]! M2 ?T=11881
  i=99
4 z8 ]# A7 f5 K# U0 a1 T2 ?  b=2
  z=114208
" h5 A2 a; U. ?0 M0 ^/ D( z! b  u=66,j=33
T=11991
8 W+ ~8 V- R$ E0 f  i=81
  b=7
  z=80944
4 H% d3 q( p* ~8 [& v; L  u=55,j=26
7 ]$ M! R$ p2 H0 W0 \% }7 u4 x4 `3 x1 ST=12101
; R! |7 T" H% C8 t+ Y  i=68
  b=2
3 _" z) Y$ m2 r! m  z=36304
4 `% W( r. V9 }9 f! d0 w  u=47,j=21
T=12211
  i=112
  b=5
  z=54952
0 D- P* L* c4 ]+ I( V. Y3 U% f  u=74,j=38
T=12321
1 s( m! R  L; x! E1 _+ c/ a+ B  i=156
) ]  h; n- j. i  b=2
  z=36964
. O0 k* [5 G  k) e3 H; ^  u=101,j=55
( v1 ~9 W- o, O$ l: d4 U. G从保持高度航程的b来看，有一半数据都是b=2，若去掉这一半的数据，可输入K1=220，r=111
输出：
4 u4 e( v7 j1 z% ^T=331
  i=24
  f0 m8 A9 B. \+ v  b=7
  z=1492
( v- v; {% {. m0 M& t/ l# B; M5 M5 s  u=18,j=6
T=551
  i=43
& P2 W$ r- W0 b& T' v3 Q  b=12
6 k; B  I# n6 W8 W  z=4192
9 F( n% s& d9 _+ Z* ^1 V( l0 N' `' n  u=30,j=13
! w) L2 |1 q' p- o# S& DT=771
  i=33
! K6 e6 Z' D; w5 K3 c  b=5
) u$ s& ~" i( @  z=3472
  u=24,j=9
T=991
* D, X& k- ~* t3 b1 [$ a  i=98
  b=41
9 \$ l& ^, y! r3 k' p  z=15064
$ G5 f  C  y, y. H. {0 o  u=64,j=34
T=1211
8 G. {: E8 ~; k& {  i=70
  b=10
  z=39364
, }& U9 h6 w; l  u=47,j=23
: T' G9 {; T8 ?T=1431
  i=26
, W/ @% H& B' |8 Z  b=7
2 c4 p% i* g$ G  z=9664
7 W$ F" L8 M* f. [9 K  u=20,j=6
T=1651
  i=135
  b=5
  z=9232
  u=87,j=48
T=1871
  i=130
  b=12
2 w% P, c& h! \. V/ x  l9 A  z=18952
# t5 i4 z# w( G  S2 J  u=84,j=46
+ U' `3 k2 M4 R( X5 p( }T=2091
% B" ~1 q4 ]! [  [  i=37
8 a0 a% o% t$ [' x  b=7
9 I$ F+ r/ S. ?% ?- {5 V- ?  z=9412
  u=27,j=10
) O# u" w) h) g! ?5 M3 q- I( p……
" n- `8 o+ u( k0 R4 `. cT=24091
1 a, h3 \8 h; [9 I: [1 P+ x/ y  i=250
  b=113
  z=9635920
! S- M) N' Z" c$ J' U" a  u=159,j=91
0 _7 n& p9 ?6 O! zT=24311
. N1 y. J) n# O1 J) k  i=188
  b=7
9 j( G2 G$ I! H# n  ]: C1 b  z=164104
  u=121,j=67
3 [, C' S8 l, ~9 g; n2 V1 zT=24531
  i=69
( v% y( |. k& ^0 r0 ?/ K# X- `! v/ x  b=5
  z=157192
  u=48,j=21
T=24751
  i=188
  b=12
4 b6 f# U5 F7 {2 O6 o% ~& l' f  z=338740
2 S$ ?+ i) F; j& M$ o  M  u=121,j=67
若要找b在220X+111中X在1-111的范围内保持最高纪录(后面的b大于前面的b)
* M' Z4 w! y, z在前面输入的基础上再选择附选框,并选择下拉选项b>G=7(7是第一个b值)
1 n! W! S  V6 j. B  g5 n输出：
! U0 E4 U  w+ Q; D+ c6 W1 tT=551
- ?3 g- L! T) J  i=43
  b=12
  z=4192
) s" |$ \+ `- z+ @& ]  w; i+ r  u=30,j=13
: e% v, M: T9 D' v: hT=991
  i=98
  b=41
  z=15064
  u=64,j=34
& \8 R+ R" c+ N9 S* Y  PT=2751
  i=172
$ Q. q/ G- \; @# f# @  b=69
  z=225988
  u=110,j=62
, v3 [; J' L8 I) d" p5 {# G7 W+ T  _T=10011
  i=210
  b=82
5 j- t! g6 r, M' Q, U  z=2633632
  u=134,j=76
$ d  P8 S* K9 ~9 ST=18591
  i=154
/ l. ]6 ^/ a( H7 Z: q, @  b=100
  z=2035348
  u=100,j=54
T=21231
  i=193
  b=113
  z=11927896
  u=124,j=69
还要可以选择其他内容做为排序对象。
6 q9 h4 V6 k4 l' b, V4 y# o9 q1 V- E例：以航程（归一步数）i为标准排序，在上面的基础上，选择下拉选项i>G=69(69是第一个i值)
- D; E3 k: Q) `2 @* F# ?( z输出：
% j5 m/ [  h8 xT=991
; |5 A. z" ?/ l& n0 C; E  i=98
  b=41
  z=15064
  u=64,j=34
T=1651
' s0 C! G  t# `6 ]2 e) n7 G- `% ^- x  i=135
  b=5
  z=9232
  u=87,j=48
T=2311
' O) l! G! I* |4 }* `' H. C  i=151
  b=10
# {* c1 o% {3 |  m5 Y* L1 ?  z=15604
( Y/ `4 @& `0 {9 n& C+ B  u=97,j=54
* Y7 s. x8 R& h2 w& aT=2751
  i=172
  b=69
( W5 q+ v5 U; N- `% X  z=225988
) p# g# @; [8 s+ \; ^* n  u=110,j=62
1 A) U$ Z, u$ qT=7591
  i=176
  b=15
  z=57652
  u=113,j=63
) C/ v- K6 H1 r1 jT=10011
% [. I, K" t2 K! M% B8 M. O  i=210
  b=82
  z=2633632
  u=134,j=76
$ E, z: h( D0 VT=24091
) E  D. u4 a! h' L# Q( ~# N  i=250
$ M+ ?, @( o6 @  H8 X; g. S  b=113
  z=9635920
3 U) F) m+ R9 Y: H- P  u=159,j=91
& n( A0 X8 t/ `! T3 C$ o+ S

---

作者: 1097908652    时间: 2012-12-23 13:04
东西太多了，感到飘过。还不如给个附件哟

---

作者: 1097908652    时间: 2012-12-23 13:05
发完之后意识到原来我抢到了一个沙发呀。。。纪念一下咯

---

作者: 弘道    时间: 2014-7-28 22:50
谢谢楼主……辛苦啦！………………

---

作者: 弘道    时间: 2014-7-28 22:50
谢谢楼主……辛苦啦！………………

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5