c++中的高精度算法

一.引入：

在计算中，我们经常会遇见许多高位数字计算，此时，因为int 与 long long 都有一个数字范围的限制，如果计算的数字过大，可能会爆掉，因此，我们需要一个工具来计算这些高位数字的运算，而这个工具就是高精度。

准备工作：

首先，我们要有一台电脑，然后，我们需要下载一个c语言编译器，最后，需要学习一下vector构建动态数组。

二.高精度加法：

这个高精度的加法特别的简单，就类似一个小学生算加法的竖式运算，给个草图回忆一下：

我们可以看到，这个加法是依次从最低位开始运算的，逢十进一， 但是这样还不够，我们需要一个东西来将每个数的每个位的数字依次存入，所以需要两个数组，并且不是平常的int或double类型的数组，我们可以用字符型的数组。

为什么？因为字符数组每个空间只能存放一个字符，并且字符和数字可以转换，只需一个-‘0’就可以做到。

下面是代码的实现：

#include<bits/stdc++.h>
using namespace std;
vector<int>add(vector<int>&A,vector<int>&B)//引用传递是因为让函数运行时快一点。
{
	vector<int>c;//设置出口数组。
	int t=0;//设置进位
	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];//防止某个数组没了还在加
		c.push_back(t%10);//取多于10的余数填入低位。
		t/=10;//存进位
	}
	if(t)c.push_back(t);
	return c;
}
int main()
{
	string a,b;
	vector<int>A,B;
	cin>>a>>b;//输入要计算的两个数
	for(int i=a.size()-1;i>=0;i--)A.push_back(a[i]-'0');//数字从个位到最高位存入，方便计算。
	for(int i=b.size()-1;i>=0;i--)B.push_back(b[i]-'0');
	auto c=add(A,B);//auto 就是vector<int>c,c++11的auto可以自动识别要命名的变量类型
	for(int i=c.size()-1;i>=00;i--)cout<<c[i];
	return 0;
}

 注意：输出时要从高位到低位输出，因为c数组原来是从低到高填入数字的，输出时要逆着输出。

下面的代码都是和这个类似的。 

三.高精度减法：

减法和加法类似，不过是进位改成借位，首先要让数组的数字减掉借位。半径有借有还。

下面是草图模拟：

 可以看到，个位的借位为零，此时0-0还是零，但是0减不了2了，所以需要向十位借一，此时个位加10，减2为8，但是十位的借位为1，0-1=-1，-1减1减不了，向高位借一，此时十位为9，9-1为8。到最后一位，此时借位为1，2-1=1，1-0=1，所以最高位为1，此时借位已还，所以借位成0了。

代码实现：

#include<bits/stdc++.h>
using namespace std;
vector<int>sub(vector<int>&A,vector<int>&B)
{
	int t;
	vector<int>c;//设置出口
	for(int i=0,t=0;i<A.size();i++)
	{
		t=A[i]-t;
		if(i<B.size())t-=B[i];//算减法有2种情况，减得了和减不了，前者t大于0，后者t<0;
		c.push_back((t+10)%10);//将两种情况合并
		if(t<0)t=1;//小于零说明借位了，所以t为1；之后下一位要减t；
		else t=0;//大于等于零说明没有借位，所以t为0；
	}
	while(c.size()>1&&c.back()==0)c.pop_back();//注意有前导零，如222-222=000，这一步是去前导0
	return c;
}
bool cmp(vector<int>&A,vector<int>&B)
{
	if(A.size() != B.size()) return A.size() > B.size();//判断A和B的大小
	    for(int i = A.size(); i >= 0; i--)
	        if(A[i] != B[i])
	            return A[i] > B[i];
	    return true;
}
int main()
{
	string a,b;
	vector<int>A,B;
	cin>>a>>b;
	for(int i=a.size()-1;i>=0;i--)A.push_back(a[i]-'0');
	for(int i=b.size()-1;i>=0;i--)B.push_back(b[i]-'0');
	if(cmp(A,B))//A>B的情况
	{
		auto c=sub(A,B);
		for(int i=c.size()-1;i>=0;i--)cout<<c[i];
	}
	else//A<B的情况，此时为-（B-A）
	{
		cout<<'-';
		auto c=sub(B,A);
		for(int i=c.size()-1;i>=0;i--)cout<<c[i];
	}
}

四：高精度的乘法：

乘法较上面就有点改变了，上面都是两个字符数组，乘法是一个字符数组，一个整型的数。我认为这样会比较好计算。。。。

下面是草图模拟：

 首先，个位的数先乘22，进2余2，则2直接填入，之后十位乘22，然后加进位的2，得24，进2余4，所以4直接填入，最后的高位乘22后加进位的2，得24，直接填入。

可以看出，有点像加法，但是又不完全像，下面我们来代码实现一下：

#include<bits/stdc++.h>
using namespace std;
vector<int>mul(vector<int>&A,int b)
{
	int t=0;
	vector<int>c;//老规矩设出口
	for(int i=0;i<A.size()||t;i++)//将计算完之后t的去向同时合并了
	{
		if(i<A.size())t+=A[i]*b;//判断A数组是否运算完了；
		c.push_back(t%10);
		t/=10;
	}
	/*将t单独出个步骤：
		if(t)c.push_back(t);
	*/
	while(c.size()>1&&c.back()==0)c.pop_back();//还是和减法一样去前导0
	return c;
}
int main()
{
	string a;
	int b;
	cin>>a>>b;
	vector<int>A;
	for(int i=a.size()-1;i>=0;i--)A.push_back(a[i]-'0');
	auto c=mul(A,b);
	for(int i=c.size()-1;i>=0;i--)cout<<c[i];
	return 0;
}

五：高精度的除法：

到最难的一个东西来了,这玩意和之前的东西都不同，在运算时上面的都是从低位到高位，而这个东西却是从高位到低位，并且还有余数要考虑

咱们先搞个草图模拟模拟：

 

首先是从高位开始，0乘10为0，加上第一位为1，1/4=0，余数为1，然后到下一位，余数乘10加上下一位的数，余数为12，12/4=3，余数为0，然后乘10加下一位8，余数为8，8/4=2，余数又为0，0乘10为零，加上最后一位数0还是零，0/4=0，所以得出 最终结果为0320，余数为0，所以代码实现时要记得去掉前导零。

下面是代码实现：

#include<bits/stdc++.h>
using namespace std;
vector<int>biv(vector<int>&A,int b,int&r)//与之前的不同，我们需要计算余数，所以直接传余数地址
{
	vector<int>c;
	r=0;//开始时余数为0
	for(int i=A.size()-1;i>=0;i--)//注意从高位开始除
	{
		r=r*10+A[i];//开始算每个步骤下的余数
		c.push_back(r/b);//计算除法
		r%=b;//算该步骤余数
	}
	reverse(c.begin(),c.end());//这个函数是将c数组的数字进行反转，方便去前导0；
	while(c.size()>1&&c.back()==0)c.pop_back();
	return c;
}
int main()
{
	string a;
	int b;
	cin>>a>>b;
	vector<int>A;
	for(int i=a.size()-1;i>=0;i--)A.push_back(a[i]-'0');
	int r;
	auto c=biv(A,b,r);
	for(int i=c.size()-1;i>=0;i--)cout<<c[i];
	cout<<endl<<r;
	return 0;
}

总结：

虽然这些都有点复杂和蠢，但是对于刚开始接触算法的萌新来说（比如我）还是很好理解的，对于以上的代码，无一不用到了STL里的函数，所以，建议去学习一下STL里的函数，这玩意儿对于写代码还是很好用的。

最后，希望本篇文章对您有帮助，虽然很简陋。如果有更好的思路，可以的话希望能提出，毕竟是学术性的交流，虽然作者很弱。。。