NYOJ.28.大数阶乘

28.大数阶乘

时间限制：3000 ms | 内存限制：65535 KB
难度：3
描述
我们都知道如何计算一个数的阶乘，可是，如果这个数很大呢，我们该如何去计算它并输出它？
输入
输入一个整数m(0<m<=5000)
输出
输出m的阶乘，并在输出结束之后输入一个换行符
样例输入
50样例输出
30414093201713378043612608166064768844377641568960512000000000000

#include <stdio.h>
#define E 100000
//大数与小数乘法，digit代表位数
void mul(int m[],int n,int &digit){
    int i,r=0;//r记录进位，个位无进位，故初始化为0
    for(i=0;i<digit;i++){
        r+=m[i]*n;
        m[i]=r%E;
        r/=E;
    }
    while(r){//判断相乘后大数长度是否有所增加并进行相关操作
        m[digit++]=r%E;
        r/=E;
    }
}
int main(){
    int f[3300],d,i,n;
    while(~scanf("%d",&n)){
        for(i=2,d=f[0]=1;i<=n;i++)
            mul(f,i,d);
        printf("%d",f[--d]);//首先输出大数最高位，因为最高位不需要补0
        while(d--)
            printf("%05d",f[d]);//一个数组元素存储5位数字，故输出的时候要顾及到0的位数
        printf("\n");
    }
    return 0;
}

使用了int型数组进行大数的存储，一个数组元素存储5位数字，故输出的时候要顾及到0的位数