大整数阶乘的计算

 1 #include <stdio.h>
 2 
 3 #define mo 10000
 4 int ans[1000001];
 5 //这里高精度整数存储格式是：ans[0]存储位数，低位在前，高位在后，每一个ans[i]存储四位整数
 6 
 7 void mul(int *,int);
 8 void print_arr(int *);
 9 void fact(int n);
10 
11 int main()
12 {
13     int n,i;
14     scanf("%d",&n);
15     printf("%d!=",n);
16     fact(n);
17     print_arr(ans);/**/
18 
19     /*ans[0]=1;
20     ans[1]=1000;
21     n=15;
22     mul(ans,n);
23     print_arr(ans);
24     for(i=0;i<=ans[0];i++) printf("%d ",ans[i]);*/
25 
26     return 0;
27 }
28 void mul(int * a,int c)
29 {
30     int i,x=0;
31     for (i=1;i<=a[0];i++)
32     {
33         a[i]=a[i]*c+x;
34         x=a[i]/mo;  //进位
35         a[i]=a[i]%mo;
36     }
37     while(x>0)      //向更高位进位
38     {
39         a[++a[0]]=x%mo;
40         x=x/mo;
41     }
42 }
43 
44 void print_arr(int * a)
45 {
46     int i;
47 
48     //单独输出高位。避免在循环体内用%04d的格式输出从而产生无意义的前缀0。
49     printf("%d",a[a[0]]);
50 
51     for (i=a[0]-1;i>0;i--)
52         printf("%04d",a[i]);
53     putchar('\n');
54 }
55 
56 void fact(int n)
57 {
58     if (n==1)
59     {
60         ans[0]=1;
61         ans[1]=1;
62         return ;
63     }
64     fact(n-1);   //先求（n-1）的阶乘
65     mul(ans,n);  //然后求n的阶乘
66 }