[NOIP1998 普及组] 阶乘之和

[NOIP1998 普及组] 阶乘之和

题目描述

用高精度计算出 $S=1!+2!+3!+\cdots +n!$（$n\le 50$）。

其中 ! 表示阶乘，例如：$5!=5\times 4\times 3\times 2\times 1$。

输入格式

一个正整数 $n$。

输出格式

一个正整数 $S$，表示计算结果。

样例 #1

样例输入 #1

3

样例输出 #1

9

提示

【数据范围】

对于 $100\%$ 的数据，$1\le n\le 50$。

【其他说明】

注，《深入浅出基础篇》中使用本题作为例题，但是其数据范围只有 $n\le 20$，使用书中的代码无法通过本题。

如果希望通过本题，请继续学习第八章高精度的知识。

思路

高精度加法、乘法。

代码

#include <bits/stdc++.h>

using namespace std;

int num[100] = {0};
int res[100] = {0};

void multi(int a)
{
    int add[100] = {0};
    for (int i = 1; i <= num[0]; i++)
    {
        add[i + 1] += (num[i] * a + add[i]) / 10;
        if (num[i + 1] || add[i + 1])
            num[0] = max(num[0], i + 1);
        num[i] = (num[i] * a + add[i]) % 10;
    }
}

void plus1()
{
    res[0] = max(res[0], num[0]);
    for (int i = 1; i <= num[0]; i++)
    {
        res[i + 1] += (num[i] + res[i]) / 10;
        if (res[i + 1] || num[i + 1])
            res[0] = max(res[0], i + 1);
        res[i] = (num[i] + res[i]) % 10;
    }
}

void fac(int n)
{
    memset(num, 0, sizeof(num));
    num[0] = 1;
    num[1] = 1;
    for (int i = 2; i <= n; i++)
    {
        multi(i);
    }
}

int main()
{
    int n;
    cin >> n;
    for (int i = 1; i <= n; i++)
    {
        fac(i);
        plus1();
    }
    for (int i = res[0]; i >= 1; i--)
        cout << res[i];
    return 0;
}