|
楼主 |
发表于 2013-9-10 12:45
|
显示全部楼层
【讨论】将 3x+1 问题推广成 qx+1 问题,结果如何?
[这个贴子最后由天山草在 2013/09/10 00:47pm 第 1 次编辑]
写一个程序,对于 q = 3,5,7,9,……,999 这些奇数分别进行试验,给定的自然数取 1000 以内,看看每一个奇数能得到多少“非平凡解”。程序如下:
For[q = 3, q <= 999, q = q + 2,
ss = 0;
For[k = 1, k <= 1000, k++,
 x = k;
 s = 0;
 While[x != 1 && s < 500, s = s + 1;
  If[EvenQ[x], x = x/2, x = q*x + 1]]; If[x == 1, ss = ss + 1]];
If[ss - 10 != 0,
 rint["q=", q, ",  ", "1000 以内能收敛到 1 的非平凡解个数 = ", ss - 10]]]
运行结果:
q=3,  1000 以内能收敛到 1 的非平凡解个数 = 990
q=5,  1000 以内能收敛到 1 的非平凡解个数 = 71
q=7,  1000 以内能收敛到 1 的非平凡解个数 = 30
q=9,  1000 以内能收敛到 1 的非平凡解个数 = 31
q=11,  1000 以内能收敛到 1 的非平凡解个数 = 5
q=13,  1000 以内能收敛到 1 的非平凡解个数 = 2
q=15,  1000 以内能收敛到 1 的非平凡解个数 = 18
q=17,  1000 以内能收敛到 1 的非平凡解个数 = 16
q=21,  1000 以内能收敛到 1 的非平凡解个数 = 12
q=23,  1000 以内能收敛到 1 的非平凡解个数 = 4
q=31,  1000 以内能收敛到 1 的非平凡解个数 = 12
q=33,  1000 以内能收敛到 1 的非平凡解个数 = 13
q=35,  1000 以内能收敛到 1 的非平凡解个数 = 4
q=39,  1000 以内能收敛到 1 的非平凡解个数 = 4
q=43,  1000 以内能收敛到 1 的非平凡解个数 = 2
q=45,  1000 以内能收敛到 1 的非平凡解个数 = 4
q=51,  1000 以内能收敛到 1 的非平凡解个数 = 8
q=63,  1000 以内能收敛到 1 的非平凡解个数 = 9
q=65,  1000 以内能收敛到 1 的非平凡解个数 = 10
q=73,  1000 以内能收敛到 1 的非平凡解个数 = 8
q=85,  1000 以内能收敛到 1 的非平凡解个数 = 10
q=89,  1000 以内能收敛到 1 的非平凡解个数 = 6
q=91,  1000 以内能收敛到 1 的非平凡解个数 = 5
q=93,  1000 以内能收敛到 1 的非平凡解个数 = 7
q=105,  1000 以内能收敛到 1 的非平凡解个数 = 5
q=117,  1000 以内能收敛到 1 的非平凡解个数 = 5
q=127,  1000 以内能收敛到 1 的非平凡解个数 = 7
q=129,  1000 以内能收敛到 1 的非平凡解个数 = 7
q=151,  1000 以内能收敛到 1 的非平凡解个数 = 3
q=195,  1000 以内能收敛到 1 的非平凡解个数 = 6
q=217,  1000 以内能收敛到 1 的非平凡解个数 = 3
q=255,  1000 以内能收敛到 1 的非平凡解个数 = 5
q=257,  1000 以内能收敛到 1 的非平凡解个数 = 5
q=273,  1000 以内能收敛到 1 的非平凡解个数 = 7
q=315,  1000 以内能收敛到 1 的非平凡解个数 = 7
q=341,  1000 以内能收敛到 1 的非平凡解个数 = 9
q=381,  1000 以内能收敛到 1 的非平凡解个数 = 5
q=399,  1000 以内能收敛到 1 的非平凡解个数 = 1
q=455,  1000 以内能收敛到 1 的非平凡解个数 = 7
q=511,  1000 以内能收敛到 1 的非平凡解个数 = 3
q=513,  1000 以内能收敛到 1 的非平凡解个数 = 3
q=585,  1000 以内能收敛到 1 的非平凡解个数 = 8
q=657,  1000 以内能收敛到 1 的非平凡解个数 = 2
q=771,  1000 以内能收敛到 1 的非平凡解个数 = 4
q=819,  1000 以内能收敛到 1 的非平凡解个数 = 8
|
|