[模板] 大数阶乘模板

最近项目工程量大，好像没什么空码代码，抽空写个简单大数阶乘，证明我还没有彻底荒废。。。

#include<iostream>
#include<cmath>
using namespace std;

int bit[3000];

void update(int hbit)  // 进位
{
    for(int i = 0; i <= hbit; i++)
    {
        if(bit[i] > 9)
        {
            bit[i+1] += bit[i] / 10;
            bit[i] %= 10;
        }
    }
}

int fact(int n) // 返回n!结果的位数
{
    double sum = 0;
    for(int i = 1; i <= n; i++)
        sum += log10(i);
    int digit = (int)sum+1;

    for(int i = 0; i <= digit; i++)
        bit[i] = 0;
    bit[0] = 1;
    int hbit = 1;
    for(int i = 2; i <= n; i++)
    {
        int j = digit;
        while(bit[j] == 0)
            hbit = --j;

        for(int j = 0; j <= hbit; j++)
            bit[j] *= i;
        update(hbit);
    }
    int j = digit;
    while(bit[j] == 0)
        hbit = --j;
    return hbit;
}

int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        int hbit = fact(n);
        while(hbit >= 0)
            printf("%d",bit[hbit--]);
        printf("n");
    }
    return 0;
}