主题:[讨论]冰雹数字,有兴趣的来看看
BlackBlizzard
[专家分:50] 发布于 2006-07-18 00:55:00
[em14]对于任意的一个正整数
1.若n是偶数,n=n/2;
2.若n是奇数,n=3*n+1;
3.若n是1,则算法结束;
我的程序:
#include <dos.h>
#include <stdio.h>
int main()
{
int n,i=1;
system("cls");
printf("please input the hailstone generate number\n");
scanf("%d",&n);
printf("Hailstone generated by :%d\n",n);
printf("%4d ",n);
while(1)
{
if(n%2==0)
{n=n/2;
i++;
printf("%4d ",n);
}
if(n%2==1&&n!=1)
{ n=3*n+1;
printf("%4d ",n);
i++;
}
if(n==1)
break;
}
printf("\n");
printf("Number of hailstone generated :%d",i);
return 0;
}
结果:
please input the hailstone generate number
3215
Hailstone generated by :3215
3215 9646 4823 14470 7235 21706 10853 32560 16280 8140 4070 2035 6106 3053 9160
4580 2290 1145 3436 1718 859 2578 1289 3868 1934 967 2902 1451 4354 2177 6532
3266 1633 4900 2450 1225 3676 1838 919 2758 1379 4138 2069 6208 3104 1552 776
388 194 97 292 146 73 220 110 55 166 83 250 125 376 188 94
47 142 71 214 107 322 161 484 242 121 364 182 91 274 137 412
206 103 310 155 466 233 700 350 175 526 263 790 395 1186 593 1780
890 445 1336 668 334 167 502 251 754 377 1132 566 283 850 425 1276
638 319 958 479 1438 719 2158 1079 3238 1619 4858 2429 7288 3644 1822 911
2734 1367 4102 2051 6154 3077 9232 4616 2308 1154 577 1732 866 433 1300 650
325 976 488 244 122 61 184 92 46 23 70 35 106 53 160 80
40 20 10 5 16 8 4 2 1
Number of hailstone generated :168
你们可以试一下,都是收敛的,后面的都是16 8 4 2 1
回复列表 (共23个回复)
沙发
punkrocker [专家分:720] 发布于 2006-07-18 09:08:00
好象是日本某个伙计研究的OOXX猜想?
板凳
euc [专家分:4310] 发布于 2006-07-18 13:55:00
有点像随机数哈
ps: [url=http://www.hnqx.gov.cn/article/Articledisp.asp?articleID=589]冰雹是怎样形成的[/url]
3 楼
BlackBlizzard [专家分:50] 发布于 2006-07-18 19:59:00
对于任意的一个正整数
1.若n是偶数,n=n/2;
2.若n是奇数,n=3*n+1;
3.若n是1,则算法结束;
随机输入一个在正整数按照上面的规则
即产生hailstone
4 楼
fwjmath [专家分:80] 发布于 2006-07-19 12:06:00
这个还是可以吧~~~不过就是可能有点太慢了~~~
6 楼
fwjmath [专家分:80] 发布于 2006-07-19 17:05:00
我就写一下我的计算里边的主要代码吧~~~不过我不是生成整个序列的数字的~~~只是整个序列的长度~~~
asm mov ax, n
asm xor si, si
start:
asm cmp eax, 4
asm jb ending
asm test ax, 1
asm jnz evenpre
asm test ax, 2
asm jnz res12
asm add si, 4
asm mov ebx, eax
asm shr ebx, 1
asm shl eax, 1
asm jc overflow
asm inc ax
asm stc
asm adc eax, ebx
asm jc overflow
asm jnc start
evenpre:
asm test ax, 2
asm jnz res12
asm add si, 2
asm shr eax, 2
asm jnc start
res12:
asm mov ebx, eax
asm shr ebx
asm sub eax, ebx
asm jmp start
overflow:
printf("overflow");
ending:
asm mov n, ax
asm mov i, si
switch(n){
case 2: i+=1;
case 3: i+=7;
case 4: i+=2;
}
7 楼
lifnmm [专家分:90] 发布于 2006-07-19 17:12:00
main()
{
int i,sum;/*sum is for the count*/
int n;
printf("put in the init num n:\n");
scanf("%d",&n);
sum=1;
printf("%5d ",n);
while(n!=1)
{
if(n%2==0)
{
n=n/2;
printf("%5d ",n);
}
else
{
n=3*n+1;
printf("%5d ",n);
}
sum++;
if(sum%10==0)
printf("\n");
}
printf("\n");
printf("end\n");
}
8 楼
BlackBlizzard [专家分:50] 发布于 2006-07-19 18:23:00
6楼的兄弟,用的汇编啊?我不懂那个东西也
7楼这个程序~时间上有提高?我怎么看不出来,还有运行成功了?
9 楼
fwjmath [专家分:80] 发布于 2006-07-19 18:35:00
的确是汇编~~~
不过我猜想把n=3*n+1换成n*=3;n++,把n/=2换成n<<=1,把n%2==0换成n&1,然后像我的程序那样分开讨论被四除的余数的话,应该速度也会有所提高的~~~
10 楼
rickone [专家分:15390] 发布于 2006-07-19 18:46:00
建议,不要讨论这里的速度(或者效率)问题,你们编的都很棒。
我来回复