NYOJ 28 大数阶乘

大数阶乘

时间限制：3000 ms  |  内存限制：65535 KB

难度：3

 

描述

我们都知道如何计算一个数的阶乘，可是，如果这个数很大呢，我们该如何去计算它并输出它？

 

输入

输入一个整数m(0<m<=5000)

输出

输出m的阶乘，并在输出结束之后输入一个换行符

样例输入

50

样例输出

30414093201713378043612608166064768844377641568960512000000000000

模拟题：

#include<stdio.h>
#include<string.h>
#define N 1000010
int a[N];
int imitation(int n)
{
    int i;
    for(i=0;a[i]>=0;i++)
    {
        a[i]*=n;
    }
    for(i=0;a[i]>=0;i++)
    {
        if(a[i]>=10)
        {
            if(a[i+1]<0)
                a[i+1]=0;
            a[i+1]+=a[i]/10;
            a[i]%=10;
        }

    }
    return i+1;

}
int main()
{
    int i,m,flag;
    memset(a,-1,N*sizeof(int));
    scanf("%d",&m);
    a[0]=1;
    for(i=2;i<=m;i++)
    {
        flag=imitation(i);
    }
    for(i=flag;i>=0;i--)
    {
        if(a[i]>=0)
            printf("%d",a[i]);
    }
    printf("
");
    return 0;
}

 
#include<iostream>
#include<iomanip>
using namespace std;
//存储20000以内的阶乘
int a[15470];
int main()
{
    //freopen("1.txt","r",stdin);
    //freopen("2.txt","w",stdout);
    int n;
    cin>>n;
    a[1]=1;
    a[0]=1;
    int up;
    for(int i=2;i<=n;++i)
    {
        up=0;
        for(int j=1;j<=a[0];++j)   //各个位相乘
        {
            a[j] *=i;
            a[j] +=up;
            up=a[j]/100000;
            a[j] %=100000;
        }
        if(up!=0)
        {
            a[0]++;
            a[a[0]]=up;
        }
    }
    if(a[0]==1) cout<<a[1];
        else
        {
            cout<<a[a[0]];
            for (int i=a[0]-1;i>0;i--)
            {
                cout<<setfill('0')<<setw(5)<<a[i];
            }
        }
}