nn:=300:T:=阵列（1..nn）：
对于从1到15的n，do:jj:=0:
对于从2到10^5的m，当（jj=0）do:p:=ithprime（m）：
对于i从1到nn，当（jj=0）do：
T[i]：=0:od:a:=1:T[1]：=p:x:=p:
当（x>1）从1到nn时，执行以下操作：
如果irem（x，2）=0，则
x：=x/2：
其他的
a： =a+1：温度[a]：=x：
x:=3*x+1:fi:
日期：
jj:=0:aa:=a:itr:=0:
对于从2到n+1的j，执行以下操作：
如果type（T[j]，prime）=true，则
itr:=itr+1：
其他fi：
日期：
如果itr=n，则
jj:=1:printf（“%d%d\n”，n，p）：
其他的
图1：
日期：
日期：
