我发现高精度可以做出来很多愈发离谱的题,就像这个P1009 [NOIP1998 普及组] 阶乘之和。
题目描述
用高精度计算出 S = 1! + 2! + 3! + \cdots + n!S=1!+2!+3!+⋯+n!(n \le 50n≤50)。
其中 !
表示阶乘,定义为 n!=n\times (n-1)\times (n-2)\times \cdots \times 1n!=n×(n−1)×(n−2)×⋯×1。例如,5! = 5 \times 4 \times 3 \times 2 \times 1=1205!=5×4×3×2×1=120。
输入格式
一个正整数 nn。
输出格式
一个正整数 SS,表示计算结果。
输入输出样例
输入 #1复制
3
输出 #1复制
9
说明/提示
【数据范围】
对于 100 \%100% 的数据,1 \le n \le 501≤n≤50。
【其他说明】
注,《深入浅出基础篇》中使用本题作为例题,但是其数据范围只有 n \le 20n≤20,使用书中的代码无法通过本题。
这道题和要我们求阶乘和,程序也就自然在P1303 A*B Problem的基础上多了一大截。
你好程序:
#include<iostream>
#include<string>
#include<cmath>
using namespace std;
int cheng[5500];
int ans[5500];
int length = 1;
void gaojia();
void gaocheng(int x){
int yinzi = x;
for (int i = 0; i < length; i++)
cheng[i] = cheng[i] * yinzi;
for (int i = 0; i < length; i++){
cheng[i + 1] =cheng[i+1]+cheng[i] / 10;
cheng[i] = cheng[i] % 10;
}
if (cheng[length++] == 0)
length--;
gaojia();
}
void gaojia(){
for (int i = 0; i < length; i++)
ans[i] = ans[i] + cheng[i];
for (int i = 0; i < length; i++){
ans[i + 1] = ans[i+1]+ans[i] / 10;
ans[i] = ans[i] % 10;
}
if (ans[length++] == 0)
length--;
}
int main(){
int n;
cin >> n;
ans[0] = 0;
cheng[0] = 1;
for (int i = 1; i <= n; i++){
gaocheng(i);
}
for (int i = length - 1; i >= 0; i--)
cout << ans[i];
return 0;
}