Ten digit numbers

mathe

那我算41－60吧。呵呵，我就用一台奔腾四

mathe

批量计算还是用我的那个方法比较好。比如你要计算31-40，那么对每个数先用幂运算算出31，然后每次再乘一个。
而不是算完所有数的31次在算32次

无心人

能节约多少时间？

无心人

不过对每个幂， 其起始数字都不一样啊

mathe

看一下我的程序，其实方法很简单，就是对每个整数，计算它的最后一个幂就可以了。（如果某个整数，最后一个幂是h,那么所有比这个数小的，这个幂也不需要算了）

无心人

知道了

你们都倒着来

mathe

原帖由 无心人 于 2009-1-13 08:51 发表
能节约多少时间？

大部分幂运算都可以仅仅用一次乘法运算代替了，这部分时间可以节省。但是数字转化为10进制统计各位数字部分就不能节省了。具体节省多少时间要看两部分计算的比例。如果幂运算占大头，那么节省的运算量会非常客观。

无心人

乘以10位数
如果是64位机器
就简单了

可惜我的服务器内存低
不敢装64位
如果是2G或者4G我就装64位

无心人

你连续算20个幂
每亿数字消耗的时间是？

我这里似乎一个数字的时间在31后
大大增加了
原来估计20分，现在估计60分打不住

gxqcn

改自 33# mathe 的代码：

1. #include < stdio.h >
   
2. #include < time.h >
   
3. #include < HugeCalc.h >
   
4. #include < HugeInt.h >
   
5. 
   
6. #pragma comment( lib, "HugeCalc.lib" )
   
7. 
   
8. #define integer CHugeInt
   
9. 
   
10. 
    
11. #define MIN_EXP    9
    
12. #define MAX_EXP    10
    
13. 
    
14. int cc[MAX_EXP+1] = { 0 };
    
15. 
    
16. void test()
    
17. {
    
18.     time_t t=time(NULL);
    
19.     integer x("10000000002"), z;
    
20.     int i, j, d;
    
21.     int mt=MAX_EXP, c=33334;
    
22.     int count[256];
    
23.     int digits;
    
24.     LPCTSTR p;
    
25. 
    
26.     while( !(!( x-=3 )) )
    
27.     {
    
28.         z.Pow( x, MIN_EXP-1 );
    
29.         for(i=MIN_EXP,d=i*10; i<=mt; i++,d+=10)
    
30.         {
    
31.             z*=x;
    
32. 
    
33.             digits = z.GetDigits();
    
34.             if ( digits == d )
    
35.             {
    
36.                 count['\0']=1;
    
37.                 for(j='0';j<='9';j++)count[j]=i+1;
    
38. 
    
39.                 for( p=z; 0!=--count[*p]; ++p)
    
40.                     ;
    
41. 
    
42.                 if ( '\0' == *p )
    
43.                 {
    
44.                     printf("\nNo.%u\t%s^%d=%s\n\n", ++cc[i], (LPCTSTR)x, i, p-d);
    
45.                     fflush(stdout);
    
46.                 }
    
47.             }
    
48.             else if( digits < d )
    
49.             {
    
50.                 mt=i-1;
    
51.                 if(mt<MIN_EXP)
    
52.                 {
    
53.                     printf("Total cost %dseconds\n", time(NULL)-t);
    
54.                     return;
    
55.                 }
    
56. 
    
57.                 break;
    
58.             }
    
59.         }
    
60. 
    
61.         if( 0==--c )
    
62.         {
    
63.             fprintf(stderr,"\t%s, times %ds\n", (LPCTSTR)x, time(NULL)-t);
    
64.             c = 300000;
    
65.         }
    
66.     }
    
67. }
    
68. 
    
69. int main(int argc, char* argv[])
    
70. {
    
71.     int i;
    
72.     test();
    
73.     for(i=MIN_EXP;i<=MAX_EXP;i++){
    
74.         printf("c[%d]=%d\n",i,cc[i]);
    
75.     }
    
76.     return 0;
    
77. }

复制代码