P1009 [NOIP1998 普及组] 阶乘之和

题目描述

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

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

输入格式

一个正整数 nn。

输出格式

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

输入输出样例

输入 #1复制

3

输出 #1复制

9

说明/提示

【数据范围】

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

【其他说明】

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

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

#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include <iostream>
#pragma warning(disable:4996) 
using namespace std;
typedef long long ll;
//总结 高精乘+高精加
int a[100][10000];//用二维数组存放每个阶层
int sum[10000];//高精加
int getlen(int now)//得出now阶层的位数
{
    int i = 0;
    int len = 10000;//不知道他有多长，反正长点没错
    for (i = len; i >= 1; i--)
    {
        if (a[now][len] == 0)
        {
            len--;
        }
        else
        {
            break;
        }
    }
    return len;
}
void mul(int now, int b)//函数表示求  now层*b  = now+1层
{
    int i = 0;
    int len = 0;
    len = getlen(now);//now阶层的位数
    int q = 0;//进位
    for (i = 1; i <= len+2; i++)//可能进位b为50多加两位，考虑到进位
    {
        a[now + 1][i] = a[now][i] * b + q;
        q = a[now + 1][i] / 10;
        a[now + 1][i] %= 10;
    }
}
void sum1(int n)//n层  高精加
{
    int i = 0;
    int j = 0;
    for (i = 1; i <= n; i++)//1-n的和
    {
        for (j = 1; j <= getlen(i); j++)//1-len位相加
        {
            if (sum[j] + a[i][j] < 10)
            {
                sum[j] += a[i][j];
            }
            else
            {
                sum[j] += a[i][j];
                sum[j] -= 10;
                sum[j + 1]++;
            }
        }
    }
}
int main()
{
    int n = 0;
    cin >> n;
    a[1][1] = 1;
    int i = 0;
    for (i = 2; i <= 50; i++)//i表示填入第i层的数
    {
        mul(i - 1, i);//第i层为  i-1层的数 *i
    }
    sum1(n);//求和
    int f = 0;
    for (i = 10000; i >= 1; i--)
    {
        if (sum[i] != 0)
        {
            f = 1;
        }
        if (f == 1)
        {
            printf("%d", sum[i]);
        }
    }
    return 0;
}