洛谷 P1009 [NOIP1998 普及组] 阶乘之和

[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

3

样例输出 #1

9

提示

【数据范围】

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

【其他说明】

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

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

NOIP1998 普及组 第二题

因为结果可能会特别的大，所以使用高精度的做法，思路是把数拆分开放进数组里，对每一位a[j]进行乘法或者加法，如果结果大于9就进位到a[j+1]，这里我用变量cr来保存进位数，遍历数组并进行运算，我们就得到了想要的那个阶乘数字，然后使用相同的方法进行高精度加法。最后倒序输出数组b[]，因为我们是从低到高存储结果的

通过某种玄学 数学手段（chatgpt）判断出结果是个65位的数字，所以最大数组设置为65

#include <bits/stdc++.h>
using namespace std;
int main()
{
	int a[65]={0},b[65]={0},i,j,n,cr=0;
	cin>>n;
	a[0]=1,b[0]=1;

	for(i = 2;i<=n;i++)
	{
		for(j = 0;j<65;j++)
		{
			a[j]=a[j]*i;
			a[j]=a[j]+cr;
			cr = 0;
			if(a[j]>9)
			{
				cr = a[j]/10;
				a[j]=a[j]%10;
			}			
		}
		
		for(j = 0;j<65;j++)
		{
			b[j]=b[j]+a[j];
			b[j]=b[j]+cr;
			cr = 0;
			if(b[j]>9)
			{
				cr = b[j]/10;
				b[j]=b[j]%10;
			}
		}
	}
	for(j=65;b[j]==0;j--);
	for(;j>=0;j--)
	{
		cout<<b[j];
	}
	
	return 0;
}