实现3x+1问题的程序代码

/相

可以使用任何程序写代码
要求有迭代的中间结果,以及迭代的次数

epezzaqegazeu

我用mathematica写的

1. (*这个代码使用vim编写*)(*3x+1问题的mathematica代码*)Clear["Global`*"];(*清除所有变量*)(*定义了两个纯函数*)(*第一个纯函数表示迭代的方式,如果是偶数则除以2;如果是奇数则乘以3再加上1*)(*第二个纯函数表示判断条件,如果不等于1则一直迭代下去*)num=20;(*被考察的整数*)p=NestWhileList[If[Mod[#,2]==0, #/2,3*#+1]&, num, #!=1&](*得到的迭代结果*)Length[p](*迭代结果的长度*)

复制代码

当然出于美观的考虑,我加上代码的图片,使得其有颜色.

ogemama

整数200的结果是
{200, 100, 50, 25, 76, 38, 19, 58, 29, 88, 44, 22, 11, 34, 17, 52, \
26, 13, 40, 20, 10, 5, 16, 8, 4, 2, 1}
和27

ifequwu

2# mathematica

jnetoji

可以仅考察奇数，如果是偶数则一次性除尽所有的2，所以迭代序列中只出现奇数。

1. num=25;(*被考察的整数*)p=NestWhileList[(3*#+1)/2^Integerexponent[#,2]&, num, #!=1&](*得到的迭代结果*)Length[p](*迭代结果的长度*)

复制代码

25的迭代结果是

1. {25, 19, 29, 11, 17, 13, 5, 1}8

复制代码

asukrimukozor

刚搜索了一下3x+1问题，居然简单到连我都可以理解。
Forcal代码：

1. i: SetIntStackMax[1000];i: x3(x:static,max) =   max++, printff{"{1,i}  ",x},   which[x==1, return(max) : x%2, x3(3*x+1) : x3(x/2)];i: printff{"\r\n结果："}, x3[200];

复制代码

结果：200  100  50  25  76  38  19  58  29  88  44  22  11  34  17  52  26  13  40  20  10  5  16  8  4  2  1  
i: 27

ameahopiyoyif

按5楼的算法，Forcal代码：

1. i: SetIntStackMax[1000];i: x3(x:static,max) = which{x%2, {printff("{1,i}  ",x), max++, which[x==1, return(max) : x3(3*x+1)]} : x3(x/2)};i: printff{"\r\n结果："}, x3[25];

复制代码

结果：25  19  29  11  17  13  5  1  
i: 8

ucobiit

7＃代码改成如下更好：

1. i: SetIntStackMax[1000];i: x3(x:static,max) = which{x%2, {printff("{1,i}  ",x), max++, which[x==1, return(max) : x3(3*x+1)]} : {while[!(x%2),x=x/2], x3(x)}};i: printff{"\r\n结果："}, x3[25];

复制代码

结果：25  19  29  11  17  13  5  1  
i: 8

ozjazit

>>Collatz(200)
in> Collatz(200)
200-->100-->50-->25-->76-->38-->19-->58-->29-->88-->44-->22-->11-->34-->17-->52-->26-->13-->40-->20-->10-->5-->16-->8-->4-->2-->1

out> 26

------------------------

By Calculator
使用手册
下载