NYOJ 大数阶乘

大数阶乘

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

难度： 3

描述

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

输入

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

输出

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

样例输入

50

样例输出

30414093201713378043612608166064768844377641568960512000000000000

思路：例如10！我们就先算出9！然后再将10个9！相加就可以转换成大数相加，9！又可以转换成9个8！相加，以此类推知道1！

 
#include<stdio.h>
#include<string.h>
int a[100000];
int main()
{
    int i,j,t;
    int n;
    int m;
    while(~scanf("%d",&n))
    {
        t=1;
        memset(a,0,sizeof(a));
        a[0]=1;
        for(i=1; i<=n; i++)
        {
            m=0;
            //printf("i=%d*********\n",i);
            for(j=0; j<t; j++)
            {
                //printf("####%d\n",a[j]);
                a[j]=a[j]*i+m;
                //printf("a[%d]=a[%d]*%d+%d=%d\n",j,j,i,m,a[j]);
                m=a[j]/10;
                //printf("m=a[%d]/10=%d\n",j,m);
                a[j]=a[j]%10;
                //printf("a[%d]=a[%d]%%10=%d\n",j,j,a[j]);
                if(m!=0&&t<=j+1)
                    t++;
                //printf("t=%d\n",t);
            }
        }
        for(i=t-1; i>=0; i--)
            printf("%d",a[i]);
        printf("\n");
    }
    return 0;
}        
将注释掉的显现出来就是整个数据的模拟过程。

例如 4！的模拟过程如下：

另一种简单的算方法：

 
#include<stdio.h>
#include<string.h>
const int maxn=20000;
int a[maxn];
int main()
{
	int n,i,j,s,c;
	scanf("%d",&n);
    memset(a,0,sizeof(a));
    a[0]=1;
    for(i=2;i<=n;i++)
	{c=0;
    for(j=0;j<=maxn;j++)
	{
    s=a[j]*i+c;
	a[j]=s%10;
	c=s/10;
	}
	}
	for(j=maxn;j>=0;j--) if(a[j]) break;
	for(i=j;i>=0;i--) printf("%d",a[i]);
	printf("\n");
	return 0;
}