概念:对两个超级大的数进行计算时,int、long long等数据类型的数据范围显然不够使用时,我们就需要引入高精度算法。
本质:高精度算法本质上时用字符串模拟数字进行计算。其步骤与人工演算过程一致。
最常用的高精度算法有四种:
- 大整数相加(整数位数约为10^6)
- 大整数相减(整数位数约为10^6)
- 大整数乘小整数(大整数的位数<=10^6,小整数的数值<=10^9)
- 大整数除以小整数
一、vector容器
在学习高精度算法之前,我们需要了解C++中vector容器(可以理解为动态数组)的基本操作:
1.头文件:
#include<vector>
2.初始化:
vector<int> a;
vector<int> a(n); //定义一个长度为n的vector数组
vector<int> a(n,m);//定义一个长度为n,每个元素为m的vector数组
vector<int> a[10]; //定义10个vector数组
3.vector常用内置函数:
a.size() //返回vector容器中元素的个数
a.empty() //返回a是否为空,若为空,返回ture
//以上两种函数是所有STL容器都支持的函数,时间复杂度为O(1)
a.clear() //清空
a.front() //返回a的第一个元素
a.back() //返回a最后一个元素
a[i] //返回a的第i个元素,当且仅当a存在
a.push_back() //在末尾插入一个数
a.pop_back() //删除最后一个数
a==b //b也为vector类型。a和b的比较操作还有!= >= <= 等
4.赋值:
- 常用的赋值操作:
vector<int> a;
for(int i=0;i<10;i++)a.push_back(i);
- 常见的错误赋值操作
vector<int>a;
for(int i=0;i<10;i++) a[i]=i;
//下标只能用来获取已经存在的元素
5.访问:
for(int i=0;i<10;i++)a.push_back(i);
//第一种遍历:下标访问
for(int i=0;i<a.size();i++)cout<<a[i]<<' ';
//第二种遍历:迭代器访问
for(vector<int>::iterator i=a.begin();i!=a.end();i++)cout<<*i<<' ';
//第三种遍历:auto关键字
for(auto i=a.begin();i!=a.end();i++)cout<<*i<<' ';
或者
for(auto i:a) cout<<i<<' ';
6.常用的函数
#include<algorithm>
//对a中的从a.begin()(包括它)到a.end()(不包括它)的元素进行从小到大排列
sort(a.begin(),a.end());
//对a中的从a.begin()(包括它)到a.end()(不包括它)的元素倒置,但不排列,如a中元素为1,3,2,4,倒置后为4,2,3,1
reverse(a.begin(),a.end());
//把a中的从a.begin()(包括它)到a.end()(不包括它)的元素复制到b中,从b.begin()+1的位置(包括它)开始复制,覆盖掉原有元素
copy(a.begin(),a.end(),b.begin()+1);
//在a中的从a.begin()(包括它)到a.end()(不包括它)的元素中查找10,若存在返回其在向量中的位置
find(a.begin(),a.end(),10);
本条出自C++_vector操作_会敲代码的地质汪的博客-CSDN博客
二、高精度加法
模板题:
注意点:
1、大整数的存储
(1)我们将大整数的每一位存到数组里面,那么下标为0的位置应该存储个位还是最高位?
在加法计算中,我们从低位到高位开始计算。当计算到最高位的时候。可能出现进位的情况。例如500+500=1000,于是我们需要开辟一个位置,存储新的最高位。假设下标为0的位置存储最高位,那么在开辟新位置存储新的最高位时,我们需要把数组中所有的数都往后移一位,显然这很麻烦;而如果下标为0的位置存储个位,那么数组的末尾存储最高位,对数组的末尾进行元素的增添十分容易。
所以,大整数在存储时,下标为0的位置存储的是个位数字。
(2)如何将大整数存储到数组中?
由于大整数远远超出int,甚至long long的数据范围,所以我们选择用字符串来进行大整数的输入。由(1)可以知道,大整数下标为0的位置存储个位数字,所以在给数组赋值时,必须从字符串的末尾开始。
例如定义字符串 s ,输入整数“123”,那么s[0]=‘1’,s[1]=‘2’,s[2]=‘3’;定义数组A,用字符串s给A赋值,那么A[0]=s[2]=‘3’-‘0’=3,A[1]=s[1]=‘2’-‘0’=2,A[2]=s[0]=‘1’-‘0’=1。
2、大整数的运算
在纸上模拟过程:
![](https://img-blog.csdnimg.cn/96c70fbb68c24da78f7606bb72578c85.png)
![](https://img-blog.csdnimg.cn/d2e66a8a886f423e901cb9be0afffced.png)
代码实现:
#include<iostream>
#include<vector>
using namespace std;
vector<int> add(vector<int> &A,vector<int> &B){
vector<int> C;
int t=0; //用t存储当前计算的值
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);
t/=10;
}
if(t)C.push_back(t); //如果t不为0,则进一位
return C;
}
int main(){
//输入大整数
string a,b;
cin>>a>>b;
//将大整数按位存储到数组中
vector<int> A,B,C;
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');
C=add(A,B);
for(int i=C.size()-1;i>=0;i--)cout<<C[i];
cout<<endl;
return 0;
}
三、高精度减法
模板题:
注意点:
1.当我们在计算减法的过程中,可能会出现两种情况:较大的数减去较小的数;较小的数减去较大的数。当我们计算较小的数减去较大的数时,例如11-22,往往会先计算22-11=11,然后再将结果前面加上一个符号‘-’,最终得到的 ‘-11’就是所求的答案。那么高精度减法也是一样。
因此高精度减法相对于其他高精度算法,要多一个步骤:判断A和B的大小。在sub函数中,一律用较大的整数减去较小的整数。
2.在运算过程中(也就是在sub函数中),可能会出现最高位为0的情况。例如5553-5523=0030,所以我们需要去除答案中处于最高位的0,直到最高位数为非零数值为止。
3.高精度减法的运算模拟:
![](https://img-blog.csdnimg.cn/7a7ee3430d574c8889c89869856ca2b0.png)
![](https://img-blog.csdnimg.cn/596fdee9b23b49bb84dd0c7482dc2561.png)
![](https://img-blog.csdnimg.cn/959556b3535b470aaa35279818b4d352.png)
代码实现:
#include<iostream>
#include<vector>
using namespace std;
//比较A和B的大小,如果A>=B,那么返回true
bool cmp(vector<int> &A,vector<int> &B){
//先比较A和B的位数。如果A的最高位高于B的最高位,那么A>B,返回true
if(A.size()!=B.size())return A.size()>B.size();
//从A、B的最高位开始比较
for(int i=A.size()-1;i>=0;i--)
if(A[i]!=B[i])return A[i]>B[i];
//如果A=B,那么返回true
return true;
}
//计算A-B的值C
vector<int> sub(vector<int> &A,vector<int> &B){
vector<int> C;
int t=0; //存储当前位数的计算结果
for(int i=0;i<A.size();i++){
t+=A[i];
if(i<B.size())t-=B[i];
//这两行计算t=a[i]-b[i]+t
C.push_back((t+10)%10); //+10是考虑到t为负数时的情况
if(t<0)t=-1; //如果t<0,那么需要借位,即高位数-1
else t=0; //如果t>=0,那么无需借位
}
//去除答案中处于最高位的0
while(C.size()>1&&C.back()==0)C.pop_back();
return C;
}
int main(){
//输入大整数
string a,b;
cin>>a>>b;
vector<int> A,B,C;
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');
//比较A和B的大小
if(cmp(A,B))C=sub(A,B);
else{
C=sub(B,A); //如果A<B,则将B减去A
cout<<'-'; //并在结果前添加负号
}
//输出结果
for(int i=C.size()-1;i>=0;i--)cout<<C[i];
return 0;
}
四、高精度乘法
模板题:
注意点:
1.需要判断小整数为0的情况
2.高精度乘法的运算模拟:
![](https://img-blog.csdnimg.cn/5c611601d4c64ff99daef48915edacc4.png)
![](https://img-blog.csdnimg.cn/bef82230551d4bb897a562934574ad14.png)
![](https://img-blog.csdnimg.cn/7c79d3b4aa5f4e18bf1ae94451ed9d35.png)
代码实现:
#include<iostream>
#include<vector>
using namespace std;
vector<int> mul(vector<int> &A,int b){
vector<int> C;
int t=0;
for(int i=0;i<A.size()||t;i++){
if(i<A.size())t+=A[i]*b;
C.push_back(t%10);
t/=10;
}
//当b为0时,不断弹出最高位,直到数组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,C;
for(int i=a.size()-1;i>=0;i--)A.push_back(a[i]-'0');
C=mul(A,b);
for(int i=C.size()-1;i>=0;i--)cout<<C[i];
return 0;
}
五、高精度除法
模板题:
注意点:
1.我们计算除法时,习惯上从高位开始计算,这与其他高精度算法不同
2.在运算过程中(也就是在div函数中),可能会出现最高位为0的情况。例如526/13=040,所以我们需要去除答案中处于最高位的0,直到最高位数为非零数值为止。
3.高精度除法的运算模拟:
![](https://img-blog.csdnimg.cn/44b6100d0c084f46b141018c4fc1e2d9.png)
![](https://img-blog.csdnimg.cn/e24af5b2ec4649cea70018552cbd5406.png)
代码实现:
#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
vector<int> div(vector<int> &A,int b,int &r){ //r为余数
vector<int> C;
for(int i=A.size()-1;i>=0;i--){ //从位数高到低的顺序遍历
r=r*10+A[i];
C.push_back(r/b); //C数组下标越大,位数越低
r%=b;
}
reverse(C.begin(),C.end()); //翻转C,使得C数组下标越大,位数越高
while(C.size()>1&&C.back()==0)C.pop_back(); //去除最高位的0
return C;
}
int main(){
string a;
int b,r=0;
cin>>a>>b;
vector<int> A,C;
for(int i=a.size()-1;i>=0;i--)A.push_back(a[i]-'0'); //A数组下标越大,位数越高
C=div(A,b,r);
for(int i=C.size()-1;i>=0;i--)cout<<C[i];
cout<<endl;
cout<<r;
return 0;
}