[模板] 大数阶乘模板

最新推荐文章于 2021-01-31 20:43:28 发布

原创

最新推荐文章于 2021-01-31 20:43:28 发布 · 1.5k 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#ACM #数论 #模板

本文分享了一种用于计算大数阶乘的模板，适用于ACM竞赛和数论问题。作者在忙碌的项目工作中抽出时间编写，以保持编程技能的活跃。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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

#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