P1009 阶乘之和 【高精度】

题目描述

用高精度计算出S=1!+2!+3!+…+n!(n≤50)S=1!+2!+3!+…+n! (n≤50)S=1!+2!+3!+…+n!(n≤50)

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

一个正整数NNN。
输出格式

一个正整数SSS，表示计算结果。
输入输出样例
输入 #1

3

输出 #1

9

说明/提示

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

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

#include <iostream>
#include <cstring>
using namespace std;

const int maxn = 1e5 + 10;
struct bign{
	int len;
	int d[maxn];
	bign() {
		len = 0;
		memset(d, 0, sizeof(d));
	}
};

bign mul(bign a, int b) {
	bign res;
	int carry = 0;
	for(int i = 0; i < a.len; i++) {
		int tem = a.d[i] * b + carry;
		res.d[res.len++] = tem % 10;
		carry = tem / 10;
	}
	while(carry)	res.d[res.len++] = carry % 10, carry /= 10;
	while(res.len > 1 && res.d[res.len - 1] == 0)	res.len--;
	return res;
}

bign add(bign a, bign b) {
	int carry = 0;
	bign ans;
	for(int i = 0; i < max(a.len, b.len); i++) {
		int tem = a.d[i] + b.d[i] + carry;
		ans.d[ans.len++] = tem % 10;
		carry = tem / 10;
	}
	while(carry)	ans.d[ans.len++] = carry % 10, carry /= 10;
	return ans;
}

void output(bign x) {
	for(int i = x.len - 1; i >= 0; i--)	cout << x.d[i];
	cout << endl;
}

int main() {
	int n;
	cin >> n;
	bign num, ans;
	num.len = 1;
	num.d[0] = 1;
	for(int i = 1; i <= n; i++) {
		num = mul(num, i);
		ans = add(ans, num);
	}
	output(ans);
	return 0;
}