阶乘之和——高精度算法

求阶乘之和，当测试数据过大时，超出了c++所能表示的最大范围，所以需要用高精度的算法求得结果。

其中用到的高精度加法和乘法，之前的博客有写过

https://blog.csdn.net/m0_66551858/article/details/128706501

这次用的高精度乘法稍微简洁一点，被乘数不需要采用数组形式，因为最多算到50的阶乘，也就是说被乘数最大为50。

每次计算阶乘不需要从1开始乘，因为在计算i的阶乘时，i-1的阶乘已经算好，再乘以i就可以，也就是说我们只需要两个数组，作为加法的两个加数，其中一个用来算阶乘，另一个存储加法最后的结果。

#include<iostream>
#include<cstdio>
using namespace std;
int n, a[101] = { 0 }, s[101] = { 0 };
void change(int x)//高精求积
{
	int g = 0;
	for (int i = 100; i >= 0; i--)
	{
		a[i] = a[i] * x + g;
		g = a[i] / 10;
		a[i] = a[i] % 10;
	}
}
void qh()//高精求和
{
	int g = 0;
	for (int i = 100; i >= 0; i--)
	{
		s[i] = s[i] + a[i] + g;
		g = s[i] / 10;
		s[i] = s[i] % 10;
	}
}
void sc()
{
	int w;
	for (int i = 0; i <= 100; i++)
	{
		if (s[i] != 0)
		{
			w = i;
			break;
		}
	}
	for (int i = w; i <= 100; i++)
		printf("%d", s[i]);
}
int main()
{
	cin>>n;
	s[100] = a[100] = 1;
	for (int i = 2; i <= n; i++)
	{
		change(i);
		qh();
	}
	sc();
	return 0;
}