【题解·洛谷】P1009 [NOIP1998 普及组] 阶乘之和

#温馨提示：此文章使用的是MD编辑器，在小程序上查看可能会出现乱码，敬请谅解
#本题解有部分文字来源于P1009 [NOIP1998 普及组] 阶乘之和 题解的内容，如果侵犯了您的权益，请联系我删除。

题目介绍

题目描述

本题难度：普及-
算法标签：数学 高精度

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

其中 ! 表示阶乘，定义为 $n!=n\times (n-1)\times (n-2)\times \cdots \times 1$。例如，$5!=5\times 4\times 3\times 2\times 1=120$。

输入格式

一个正整数 $n$。

输出格式

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

样例 #1

样例输入 #1

$\text{3}$

样例输出 #1

$\text{9}$

提示

【数据范围】

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

【其他说明】

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

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

NOIP1998 普及组 第二题

题解代码及思路

本题是一道考察高精度知识的题目。
由于数据范围过大，

当我们输入 $50$ 时，
会输出 $31035053229546199656252032972759319953190362094566672920420940313$

所以我们需要使用高精度来完成。
其他解析见代码注释。

AC代码如下

#include<bits/stdc++.h>
using namespace std;

int a[100000], n, i, y, xy[100000], s[100000]; // 定义高精度数组 a 和 s，用于存储阶乘值和阶乘和

void add() { // 高精度加法函数：将当前阶乘值 a 加到阶乘和 s 上
    int i;
    memset(xy, 0, sizeof(xy)); // 清零临时数组 xy，并初始化第一个元素为 max(s[0], a[0])
    xy[0] = max(s[0], a[0]);
    for (i = 1; i <= xy[0]; i++) { // 进行每一位的相加，并处理进位
        xy[i] += s[i] + a[i];
        xy[i + 1] = xy[i] / 10;
        xy[i] %= 10;
    }
    while (xy[xy[0] + 1] > 0) { // 处理可能存在的更高位进位
        xy[xy[0] + 2] = xy[xy[0] + 1] / 10;
        xy[xy[0] + 1] %= 10;
        xy[0]++;
    }
    s[0] = xy[0]; // 更新阶乘和数组 s 的内容
    for (i = 1; i <= xy[0]; i++) s[i] = xy[i];
}

int main() {
    cin >> n; // 输入要求计算阶乘之和的数 n
    a[0] = 1; // 初始化 a 数组，a[0] 表示阶乘值的长度，a[1] 表示 1! 的值，即 1
    a[1] = 1;
    s[0] = 1; // 初始化阶乘和数组 s，s[0] 表示阶乘和的长度，s[1] 表示初始时只有 1!
    s[1] = 0;
    for (y = 1; y <= n; y++) { // 计算从 1 到 n 的阶乘值并累加到阶乘和数组 s 中
        memset(xy, 0, sizeof(xy)); // 计算 y! 的每一位，并存入临时数组 xy
        xy[0] = a[0];
        for (i = 1; i <= a[0]; i++) {
            xy[i] += a[i] * y;
            xy[i + 1] = xy[i] / 10;
            xy[i] %= 10;
        }
        while (xy[xy[0] + 1] > 0) { // 处理高位进位，更新 a 数组的内容
            xy[xy[0] + 2] = xy[xy[0] + 1] / 10;
            xy[xy[0] + 1] %= 10;
            xy[0]++;
        }
        for (i = 1; i <= xy[0]; i++) a[i] = xy[i];
        a[0] = xy[0];
        add(); // 将新计算出的阶乘值累加到阶乘和数组 s 上
    }
    for (i = s[0]; i >= 1; i--) cout << s[i]; // 输出阶乘之和的结果
    cout << endl;
	return 0;
}