求大数40000阶乘

近来在看算法导论，编译原理方面的书籍，想起写些小玩意就当休闲了，以下是截图

 

C

 

代码如下：

2. #include <stdio.h>
3. #include <malloc.h>
4. #include <time.h>
5. void main(void)
6. {
7.     
8.     int* Data = malloc(50000*sizeof(int));  //大数的阶乘位数数组如果用数组来存已超过1M堆栈的空间，所以这里用动态存储,预计计算结果不超过100000个数据元素，一个数据元素中最大存储值不大于46340，即最终最大的阶乘值最多可以达到231700000位(天文数字)
9.     int Digit;  //数据位数变量
10.     int i,j,r,k;    //循环计数器变量
11.     int N;  //用户输入要求阶乘的值
13.     time_t starttime;
14.     time_t endtime;
15.     time_t elapsetime;
16.     
17.     for (i = 1; i < 50000+1; i++)
18.     {
19.         Data[i] = 0;
20.     }
22.     Data[0] = 1;    //0位为标志位1
23.     Data[1] = 1;    //1位为标志位1
24.     Digit = 1;  //数据位为1，代表大数组中已有几个数组元素作为数据位了
26.     printf("Enter a number what you want to calculus : ");
27.     scanf("%d", &N);    //读取欲求N值
29.     time(&starttime);   //time start
31.     for (i = 1; i < N+1; i++)
32.     {
33.         for (j = 1; j < Digit+1; j++)
34.         {
35.             Data[j] *= i;   //数组元素中已有位（数据位）的内容运算
36.         }
38.         //转换为46340进制
39.         for (j = 1; j < Digit+1; j++)
40.         {
41.             if (Data[j] > 46340)    //每一位达到46340时，就考虑进位问题（即46340进制），因为在x86 32位os中int型的最大正值是2的31次方减一即2147483647，考虑到最大一个数据位乘46340以上时将超过int的最大范围（这里也说明可以求最大46340左右的数的阶乘），在这一块上可以进行优化（这里一个元素存储46340是否为了尽量使用一个int型所能存储的信息量，即节约了空间），例如修改成long型元素。
42.             {
43.                 for (r = j; r < Digit+1; r++)
44.                 {
45.                     if (Data[Digit] > 46340)    //动态增加一位，因为可能遇到最高位数据位是46339的情况，最多最高位位数计数器加一
46.                     {
47.                         Digit++;    
48.                     }
49.                     Data[r+1] += Data[r]/46340; //进位
50.                     Data[r] = Data[r]%46340;    //留余
51.                 }
52.             }
53.         }
54.     }
55.     printf("%d! = ", N);
56.     for (k = Digit; k > 0; k--)
57.     {
58.         printf("%d", Data[k]);
59.     }
60.     printf("/n");
61.     time(&endtime); //time end
62.     printf("start time : %d/n", starttime);
63.     printf("end time : %d/n", endtime);
64.     elapsetime = endtime - starttime;
65.     printf("elapsetime : %d/n", elapsetime);
67. }

大数阶乘，大进制，数组，分治策略。思路可以借鉴到其他大数预算操作上。

 

http://blog.csdn.net/jrckkyy
http://hi.baidu.com/jrckkyy

[e-mail]