题目描述
用高精度计算出 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;
}