大整数阶乘的计算

#include <stdio.h>

#define mo 10000
int ans[1000001];
//这里高精度整数存储格式是：ans[0]存储位数，低位在前，高位在后，每一个ans[i]存储四位整数

void mul(int *,int);
void print_arr(int *);
void fact(int n);

int main()
{
    int n,i;
    scanf("%d",&n);
    printf("%d!=",n);
    fact(n);
    print_arr(ans);/**/

    /*ans[0]=1;
    ans[1]=1000;
    n=15;
    mul(ans,n);
    print_arr(ans);
    for(i=0;i<=ans[0];i++) printf("%d ",ans[i]);*/

    return 0;
}
void mul(int * a,int c)
{
    int i,x=0;
    for (i=1;i<=a[0];i++)
    {
        a[i]=a[i]*c+x;
        x=a[i]/mo;  //进位
        a[i]=a[i]%mo;
    }
    while(x>0)      //向更高位进位
    {
        a[++a[0]]=x%mo;
        x=x/mo;
    }
}

void print_arr(int * a)
{
    int i;

    //单独输出高位。避免在循环体内用%04d的格式输出从而产生无意义的前缀0。
    printf("%d",a[a[0]]);

    for (i=a[0]-1;i>0;i--)
        printf("%04d",a[i]);
    putchar('\n');
}

void fact(int n)
{
    if (n==1)
    {
        ans[0]=1;
        ans[1]=1;
        return ;
    }
    fact(n-1);   //先求（n-1）的阶乘
    mul(ans,n);  //然后求n的阶乘
}