阶乘的计算以及大数的表示

一、精确计算1000!的阶乘

1000!有多大呢？拿微软自带的计算器一算，结果是4.02*10^2567，共有2568位。

在C语言中我们没有能够精确表示这个数字的数据类型。

如果非要计算，那么只能以数组的形式存放每一位数字。

代码不太难，如下：

 1 #include <stdio.h>
 2 #include <string.h>
 3
 4 #define maxn 3000
 5
 6 int f[maxn];
 7
 8 int main()
 9 {
10     int i, j, n;
11     scanf("%d", &n);
12     memset(f, 0, sizeof(f));
13     f[0] = 1;
14     for (i = 2; i <= n; i++) {
15         int c = 0;
16         for (j = 0; j < maxn; j++) {
17             int s = f[j] * i + c;
18             f[j] = s % 10;
19             c = s / 10;
20         }
21     }
22     for (j = maxn-1; j >= 0; j--) if (f[j]) break;
23     printf("total bits：%d\n", j+1);
24     for (i = j; i >= 0; i--) printf("%d", f[i]);
25     printf("\n");
26
27     return 0;
28 }

如果你不想像第16行一样，遍历数组里的每一位来做乘法，那么可以加上一个跟踪参数validBits来设置循环的阈值。这个参数初始化为1，当j达到validBits-1而且进位不为零的时候，validBits自增1。

如果现在要我们计算10000000!的阶乘，怎么办？刚才我们用计算器计算1000的阶乘，可以得到一个科学表示法的结果。当我们计算10000的阶乘时，Windows的计算器提示“溢出”。

那么我们现在简化问题，先计算10000000的阶乘有多少位有效数字。

二、计算阶乘的位数

先给出一个结论：有一个正整数n，则它的位数为log10(n) + 1。

解释如下：

对于正整数n，我们有10^(x-1) <= n < 10^x，其中x为正整数。则我们可以指出n的位数为x。

对不等式取对数：log10(10^(x-1)) <= log10(n) < log10(10^x)，简化得到：x-1 <= log10(n) < x。

则我们有log10(n) + 1 >= x，且log10(n) < x，

则对log10(n)取整之后，我们有(int)log10(n) + 1 = x。

现在我们要求n!的位数，则相当于求(int)log10(n!) + 1。

而log10(n!)=log10(n) + log10(n-1) + log10(n-2) + ... + log10(1)

我们很容易编程解决这个问题。

 1 #include <stdio.h>
 2 #include <math.h>
 3
 4 int
 5 main(int argc, char** argv)
 6 {
 7     double result;
 8     int input;
 9     int j;
10
11     result = 0;
12     scanf("%d", &input);
13
14     for (j = 1; j <= input; j++) {
15         result += log10(j);
16     }
17
18     printf("%d\n", (int)result + 1);
19
20     return 0;
21 }

运行程序，输入10000000，输出结果为65657060。

所以10000000的阶乘的结果有65657060位数字。