回 帖 发 新 帖 刷新版面

主题:[讨论]冰雹数字,有兴趣的来看看

[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个回复)

沙发

好象是日本某个伙计研究的OOXX猜想?

板凳

有点像随机数哈

ps: [url=http://www.hnqx.gov.cn/article/Articledisp.asp?articleID=589]冰雹是怎样形成的[/url]

3 楼

对于任意的一个正整数
1.若n是偶数,n=n/2;
2.若n是奇数,n=3*n+1;
3.若n是1,则算法结束;
随机输入一个在正整数按照上面的规则
即产生hailstone

4 楼

这个还是可以吧~~~不过就是可能有点太慢了~~~

5 楼

那你有什么好的改进?

6 楼

我就写一下我的计算里边的主要代码吧~~~不过我不是生成整个序列的数字的~~~只是整个序列的长度~~~

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 楼

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 楼

6楼的兄弟,用的汇编啊?我不懂那个东西也
7楼这个程序~时间上有提高?我怎么看不出来,还有运行成功了?

9 楼

的确是汇编~~~
不过我猜想把n=3*n+1换成n*=3;n++,把n/=2换成n<<=1,把n%2==0换成n&1,然后像我的程序那样分开讨论被四除的余数的话,应该速度也会有所提高的~~~

10 楼

建议,不要讨论这里的速度(或者效率)问题,你们编的都很棒。

我来回复

您尚未登录,请登录后再回复。点此登录或注册