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

时间限制1.00s         内存限制125.00MB         难易度：普及−

【题目描述】

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

其中 ! 表示阶乘，定义为  n!=n×(n−1)×(n−2)×⋯×1。例如， 5!=5×4×3×2×1=120。

【输入格式】

一个正整数 n。

【输出格式】

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

【输入输出样例】

输入 #1

3

输出 #1

9

【说明/提示】

【数据范围】

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

【其他说明】

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

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

NOIP1998 普及组 第二题

【算法分析】

这道题需要用高精度来算，如果n=50，S就有65位，远超long long 的10^18 的数据大小。本题我们可以将数化为数组，每个数组存储一位，一位一位乘，一位一位加，一位一位输出，加和乘可以模拟竖式计算，注意处理进位。

【参考代码】

#include<bits/stdc++.h>
using namespace std;
int a[1000],b[1000];
int main()
{
 	int n;
 	cin>>n;
 	a[1]=1;
 	for(int i=1;i<=n;i++)
	{
		int cnt=0;
		for(int j=1;j<=1000;j++) //乘(注意进位)。
		{
			a[j]=a[j]*i+cnt;
			cnt=a[j]/10;
			a[j]=a[j]%10;
		}
		cnt=0;
		for(int j=1;j<=1000;j++) //加(注意进位)。
		{
			b[j]=b[j]+a[j]+cnt;
			cnt=b[j]/10;
			b[j]=b[j]%10;
		}
 	 }
	 int k;
	 for(int i=1000;i>=0;i--) //因为最高位不可能为0,所以从后往前找第一个不为0的数即最高位。
	 	if(b[i]!=0)
	 	{
	 		k=i;
	 		break;
 		}
	 for(int i=k;i>=1;i--)
	 	cout<<b[i];
    return 0;
}