高精度--阶乘之和（南昌理工学院ACM寒假集训队）

身为蒟蒻的我学习了很久，有点小经验就分享给大家啦！
这道例题呢，来自洛谷P1009 [NOIP1998 普及组] 阶乘之和
大家可以先看看例题哦，对于初学高精度的还是很有帮助的！

1. 我们呢首先来分析一下题目，题目意思呢，超级简单哦，嘻嘻，给你一个int型的数n，要求计算n的阶乘之和（n<=50)
2. 进入分析之前呢，我想分享一下我第一次的想法：
   用一个for循环把每次得到的阶乘加起来，然后用一个递归函数把每次的i的阶乘给算出来，哈哈哈哈。简单吧，上代码

#include<iostream>
using namespace std;
int i,sum=0,n;
int jie(int x)
{
   if(x==1)return 1;
   return x*jie(x-1);
}
int main()
{
   cin>>n;
for(i=1;i<=n;i++)sum+=jie(n);
cout<<sum<<endl;
   return 0;
}

嗯不错，我第一下试了一下1，第二次试了一下3，感觉没问题了，结果提交了！
可提交上去，我哭了，0分！也就是一次测试都没有过，我就试了一下大的数，试了30；
结果发现等于一个负数！

这是为什么呢？因为30的阶乘超过了int型的数据范围，所以输出的必然是错误的！

题目给的是1~50的范围，那么想想最大50的时候，这个数据其实已经超出了long long 的范围了。。。

用整型还是解决不了问题，那么就得用高精度了

什么是高精度呢？高就是大的意思，精度就是这个数据的范围，可以超过10的1000000000次方的位数
注意我说的是位数！什么概念，这个数有100000000000这么多位！
那么我们可以用string来保存，不过在这呢，我推荐大家用vector容器装这个数，为什么呢？
因为不管是读入这个数字无论多大都能用a.push_back(int b)来一直向数组a的后面加数，还是这个数组的长度a.size()都很方便！
好，这样的话，我们就开始今天的题解了！

1. 思路呢，和之前差不多，先定义一个int n；用来读入需要计算的阶乘之和的数字；
2. 定义vector leijia; vector sum;
3. 第一个（外循环）for就是计算每一个 i 从1到n的阶乘的和；
4. 第二个（内循环）for就是计算 i 的阶乘；
5. 最后呢，就是输出这个结果数组了；

int main()
{
	cin >> n;

	vector<int> leijia;//存累加的结果

	for (int j = 1; j <= n; j++)
	{
		vector<int> sum;//存阶乘的结果
		sum.push_back(1);
		for (int i = 1; i < j; i++)//注意是i<j，我写的时候把j放下面了
		{
			sum = mul(sum, i + 1);//当前j的阶乘，从i=1开始和后面的数i+1相乘
		}
		leijia = jia(leijia, sum);
	}

	for (int i = leijia.size() - 1; i >= 0; i--)cout << leijia[i];//在这为什么要逆序输出呢，原因就是计算的时候低位在数组最开始的位置！
	cout << endl;
	return  0;
}

先看看mul函数计算阶乘，详细在注释上面

vector<int> mul(vector<int> &sum, int b)//开头是返回int型vector数组，后面呢在sum前面加一个&，为引用
{
	int t=0;//处理进位
	vector<int> jie;//jie是计算后的结果
	for (int i = 0; i < sum.size()||t; i++)//只要t还有进位就继续
	{
		if(i<sum.size())t+=sum[i]*b;//每一个进位为当前数sum[i]*b加上上一次进位
		jie.push_back(t % 10);//将得到的数进入数组里
		t /= 10;//进位除掉个位数
	}
	return jie;
}

再来看看两位数的相加

vector<int> jia(vector<int> &a, vector<int> &b)
{
	int t = 0;//处理进位
	vector<int>jia;
	for (int i = 0; i < a.size() || i < b.size(); i++)
	{
		if(i<a.size())t += a[i];//只要数组a没计算完，就加上
		if (i < b.size())t += b[i];//只要数组b没计算完，就加上
		jia.push_back(t % 10);
		t /= 10;//进位数等于当前和的/10
	}
	if (t)a.push_back(1);//计算到最后还有进位的话，那就进位吧
	return jia;
}

完整代码

#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;
int n;
vector<int> jia(vector<int> &a, vector<int> &b)
{
	int t = 0;//处理进位
	vector<int>jia;
	for (int i = 0; i < a.size() || i < b.size(); i++)
	{
		if(i<a.size())t += a[i];
		if (i < b.size())t += b[i];
		jia.push_back(t % 10);
		t /= 10;
	}
	if (t)a.push_back(1);
	return jia;
}
vector<int> mul(vector<int>& sum, int b)
{
	int t = 0;
	vector<int> jie;
	for (int i = 0; i < sum.size() || t; i++)
	{
		if (i < sum.size())t += sum[i] * b;
		jie.push_back(t % 10);
		t /= 10;
	}
	return jie;
}
int main()
{
	cin >> n;

	vector<int> leijia;//存累加的结果

	for (int j = 1; j <= n; j++)
	{
		vector<int> sum;//存阶乘的结果
		sum.push_back(1);
		for (int i = 1; i < j; i++)
		{
			sum = mul(sum, i + 1);//当前j的阶乘

		}
		leijia = jia(leijia, sum);
	}

	for (int i = leijia.size() - 1; i >= 0; i--)cout << leijia[i];
	cout << endl;
	return  0;
}

哎，写完不容易啊，毕竟我是个蒟蒻；
希望对大家有帮助哦！