阶乘运算——ACM

大数阶乘

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

难度：3

 

描述

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

 

输入

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

输出

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

样例输入

50

样例输出

30414093201713378043612608166064768844377641568960512000000000000

看到大数运算首先想到的是运用数组在表示大数，以前是用过int型数组，但是考虑到这将会占用大量内存。所以决定尝试 char型数组，char和int的转化很简单，(int)(n - 48)和(char)(c + 48)。
然后就是解决大数的乘法问题，乘法中的难点在进位。想起来很简单，可实际上程序实现起来非常复杂，差不多用了我一下午时间，真是惭愧。下面就是程序：

#include <stdio.h>
#include <string.h>
#include <stdlib.h>

#define NUM 20000

int main(char argc, char**argv)
{
    char a[NUM];
    char t[NUM];
    char nstr[];
    int i, n, oldj = ,len = ;
    scanf("%d", &n);

    memset(a,,NUM*sizeof(char));
    a[] = '';
    memset(t,,NUM*sizeof(char));
    t[] = '';

    for(i = ; i <= n; i++)
    {
        int count = , lastAdd = , j, c = ;
        int ii = i;
        // 拷贝字符串数组
        strcpy(a, t);
        // 整数n的转化
        while(ii)
        {
            ii/=;
            count++;
        }
        sprintf(nstr, "%d", i);

        for(c = ; c < count - ; c++)
        {
            if(a[oldj + c] == '\0')
                a[oldj + c] = '';
        }

        for(j = ; j < NUM; j++)
        {
            int remain;
            int nn = , k;

            if(a[j] == '\0')
                break;

            for(k = ; k <= j && k < count; k++)
            {
                nn += (int)(nstr[count--k] - ) * (int)(a[j - k] - );
            }
            nn += lastAdd;
            lastAdd = nn / ;
            remain = nn % ;
            if(lastAdd !=  && a[j + ] == '\0')
                a[j + ] = '';
            t[j] = (char)(remain + );
        }
        oldj = j;
    }

    for(i = oldj-; i >= ; i--)
    {
        printf("%c", t[i]);
    }
    printf("\n");

    system("pause");
}

提交代码后发现还有更优的算法，而且在时间和空间上比我的程序好好几倍。瞬间感觉自己弱爆了。

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

基本思路基本一致，只不过它用的是int型数组，但是为什么我的程序运行速度更慢，也许是字符和数字之间的转化太多造成的。弄巧成拙了。。