|
c： =0:
对于n从3×2到800 do：
x:=2:lst:={n}:lst1:={}:x:=n:
对于k从1到120，而（x>1），则：
a:=0:
如果类型（x，'even'），则
x:=x/2:lst:=lst联合{x}:a:=a+1:
其他
x:=3*x+1:lst:=lst并集{x}:a:=a+1:
金融机构：
外径：
n1:=无（lst）：
对于从1到n1的u，do：
如果irem（lst[u]，2）=1，则
lst1:=lst1联合{lst[u]}：
其他
金融机构：
外径：
m1:=最大值（op（lst1））：n1:=无（lst1）：
对于从2乘2到m1+1的p，do：
lst2:={}：
对于从1到n1的q，do：
lst2:=lst2联合{irem（lst1[q]，p）}：
外径：
lst3:={}：n2:=nops（lst2）：kkk:=0:
对于从1到n2的i，do:jjj:=0:
对于从1到n2的j，do：
z:=irem（lst2[i]*lst2[j]，p）：lst3:=lst3联合{z}：
如果z=1，则jjj:=1：否则fi
外径：
如果jjj=0，则kkk:=1:否则fi:
外径：
n3:=nops（lst3）：iii:=0:
对于b，从1到n3，而（iii=0，n2=n3，kkk=0）
执行：
如果lst2[b]<>lst3[b]，则
iii:=1:其他
金融机构：
外径：
如果iii=0，n2=n3，kkk=0，则c:=c+1：
printf（“%d%d\n”，c，p）：
其他
金融机构：
外径：
x:=2:
外径：
| 