精确求阶乘

《算法竞赛入门》：用较大的数组存放所求阶乘，在此定义的MAXN表示最大可求到数据长达3000位

//阶乘的精确值

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

#define MAXN 3000
int f[MAXN]; //3000位数据，从0开始为个、十、百

int main()
{
    int i, j, s, c, n;

    scanf("%d", &n);
    memset(f, 0, sizeof(f));
    f[0] = 1; //初始值

    for (i = 2; i <= n; i++) {
        c = 0;

        for (j = 0; j < MAXN; j++) {

            s = f[j]*i + c; //计算当i增长时与f的各位的乘积
            f[j] = s%10;    //取出个位
            c = s/10;       //取出除去个位的高位
        }
    }
    for (j = MAXN-1; j >= 0; j--) //截掉高位0
        if (f[j] != 0)
            break;

    for (j; j >= 0; j--)    //输出有效数位
        printf("%d", f[j]);
    //printf("\n");

    return 0;
}

测试：

输入：30

输出：265252859812191058636308480000000

^_^ THE END