10000的阶乘

#include<iostream>
#define N 10000//设阶乘数为10000
#define SIZE 40000//预留40000位保存结果
using namespace std;

void mul(int *n1,int n2)//两数相乘
{
    for(int i=0;i!=SIZE;++i)
        n1[i]*=n2;
    for(int i=0;i!=SIZE;++i){
        n1[i+1]+=(n1[i]/10);
        n1[i]%=10;
    }
}
int main()
{
    int *num=new int[SIZE];
    memset(num,0,SIZE);
    num[0]=1;

    //N的阶乘
    for(int i=1;i!=N+1;++i){
        mul(num,i);
    }

    //输出时排除开头的0
    for(int i=SIZE-1;i!=-1;--i){
        if(0!=num[i]){
            for(int j=i;j!=-1;--j){
                cout<<num[j];
            }
            break;
        }
    }
    delete[] num;
    return 0;
}

 

  这个实现版本需要预留空间，这就要我先预算一下会有多少位了。运算时间35秒左右，测试机器为P3 733。