求10000以内n的阶乘

总时间限制: 5000ms 内存限制: 655360kB

描述

求10000以内n的阶乘。

输入只有一行输入，整数n（0<=n<=10000）。输出一行，即n!的值。样例输入

100

样例输出

93326215443944152681699238856266700490715968264381621468592963895217599993229915608941463976156518286253697920827223758251185210916864000000000000000000000000

来源http://noi.openjudge.cn/ch0106/14/

#include<iostream>
#include<cstring>
using namespace std;
int n,y,i;
int a[100000],ay[100000];

int main()
{
    cin>>n;
    a[0]=1;   //高精度数的位数 
    a[1]=1;   //高精度整数1（逆向存储，位数放在数组第0个元素） 
    for(y=1;y<=n;y++)
    {
        memset(ay,0,sizeof(ay));
        ay[0]=a[0];
        for (i=1;i<=a[0];i++)
        {
          ay[i]+=a[i]*y;
          ay[i+1]=ay[i]/10;  //向高位进位 
          ay[i]%=10;
        }
        while(ay[ay[0]+1]>0)  //向高位进位 
        {
            ay[ay[0]+2]=ay[ay[0]+1]/10;
            ay[ay[0]+1]%=10;
            ay[0]++;
        }
        for(i=1;i<=ay[0];i++) a[i]=ay[i];
        a[0]=ay[0];
    }
    for(i=a[0];i>=1;i--) cout<<a[i];
    cout<<endl;
    return 0;
}