阶乘的精确值

例如输入不超过1000的正整数n，输出n!=1X2X3...Xn的精确结果。例如输入30输出 265252859812191058636308480000000

当输入1000时1000的阶乘约4*10^(2567)，大概3000位数字，如果采用普通的阶乘算法很明显溢出：

int fun(int n){
    int i;
    int s=1;
    for(i=1;i<=n;i++)
    s*=i;
    return s;
    } //求n的阶乘的普通算法

故需要引入高精度算法。引入一个大小3000的数字a_int，让a_int[0]保留个位，让a_int[1]保留十位，让a_int[2]保留百位。。。之所以这么做我们可以按照以下分析：

当n=1时，最终结果1，a_int[0]=1，其他位0即可；

当n=2时，最终结果2，a_int[0]=2，其他位0即可；

当n=3时，最终结果6，a_int[0]=6，其他位0即可；

.当n=4时，最终结果24，a_int[0]=24，此时需要十位和个位拆开，a_int[0]=4，a_int[1]=2，其他位0即可；

.当n=5时，最终结果120，此时a_int[0]=4，a_int[1]=2，其他位0，即120=24*5，4*5=20的个位0放在a_int[0]，而其十位需要进位，2*5=10为百位和十位，再加上刚才的进位，得百位十位为12，类似拆开a_int[1]=2，a_int[1]=1其他位0即可；

。。。

抽象分析，假设i!=1X2X3...Xi结果已经算出并保存在j个数据单元里即i!=a_int[j]a_int[j-1]...a_int[2]a_int[1]a_int[0],故输入i+1时，按照要求：

(i+1)!=1X2X3...XiX(i+1)，可写出此时乘法竖式：

对于个位a_int[0],a_int[0]*（i+1）所得数值取其个位保留在a_int[0].剩余的位作为进位；

对于十位位a_int[1],a_int[1]*（i+1）所得数值再加上个位的进位取其个位保留在a_int[1].剩余的位作为进位；

。。。

对于第k位a_int[k-1],a_int[k-1]*（i+1）所得数值取其个位保留在a_int[k-1].剩余的位作为进位；

.故算法如下：

void fun(){
     int i,j;
     memset(a_int,0,sizeof(a_int));
     a_int[0]=1;
     int n;
     int total,res;
     scanf("%d",&n);
     for(i=1;i<=n;i++)
     {res=0;total=0;
     for(j=0;j<M;j++)
     {
     total=a_int[j]*i+res;
     a_int[j]=total%10;
     res=total/10;
                     }
                      } 
for(j=M-1;j>=0;j--)if(a_int[j])break;//舍去无用的0
for(i=j;i>=0;i--)printf("%d",a_int[i]); 
printf("\n");
     }//阶乘的高精度运算