hdu1042

题目链接：hdu1042    tzc1213

代码一：转自 http://blog.pfan.cn/yzfy/25193.html ，有所改动。

#include<stdio.h>
#define N 10005    
int main()
{
    int t,a,b;
    int n;
    while(scanf("%d",&n)!=EOF)
    {
         t=2,a=1,b=0;
         if(n==0||n==1){printf("1\n");continue;}
         long s[N]={1,1};
        for(;a<=*s||(++t<=n?(b=0,a=1):0);(*s==a++&&b)?(*s)++:0)
            s[a]=(b+=s[a]*t)%10000,b/=10000;
        for(printf("%d",s[*s]);--*s>0;)printf("%04d",s[*s]);
        printf("\n");
    }   
    return 0;
}

代码二：转自http://blog.pfan.cn/yzfy/24430.html ，有所改动。

//以10进制输出结果，本程序计算1000!用时不到50ms
#include <stdio.h>
#include <conio.h>
#include <time.h>  //用于计时
#define MaxNum    10005  //阶乘最大的数(最高位数)
//如果要算更大的n!请自行修改此值

void LargeNumberTimes(long *num,long &nMax,long nTimes)
{
    long z1,z2,z3=0;
    for(z1=0;z1<=nMax;z1++)
    {
        if((z2=num[z1]*nTimes+z3)>=10000)
        {
            z3=z2/10000,z2%=10000;
            if(z1==nMax)nMax++;
        }
        else z3=0;
        num[z1]=z2;
    }
}

int main()
{    //雨中飞燕之作;
    long n;
    while(scanf("%d",&n)!=EOF)
    {    
        //输入n，求n!
        long n1,nt,nMax=0,n0=0,n5=0,nc2=-10,nMod5,num[MaxNum]={1};
        long t=clock();//开始计时
        for(nMod5=n1=2;n1<=n;n1++,nMod5++)
        {
            if(nMod5==10)n0++,nc2++,nt=n1/10,nMod5=0;else
            if(nMod5==5)
            {
                n0++,nt=n1/5;
                while(nt%5==0 && nt)n0++,nt/=5;
            }
            else nt=n1;
            while(n0>=nc2 && (nt&1)==0 && nt)nc2++,(nt>>=1);
            if(nt>1)LargeNumberTimes(num,nMax,nt);
        }
        if((nc2+=10)>n0)LargeNumberTimes(num,nMax,(1<<(nc2-n0)));
        t=clock()-t;//结束计时
     
        printf("%d",num[nMax]);//输出计算结果
        for(n1=nMax-1;n1>=0;n1--)printf("%04d",num[n1]);
        for(n1=0;n1<n0;n1++)printf("0");printf("\n");
    }
    return 0;
}