高精度渐加减乘除运算
在一般的科学计算中,会经常算到小数点后几百位或者更多,当然也可能是几千亿几百亿的大数字。一般这类数字我们统称为高精度数,高精度算法是用计算机对于超大数据的一种模拟加、减、乘、除、乘方、阶乘、开方等运算。 对于一个很大的数字N >= 10^ 100,很显然这样的数字无法在计算机中正常存储。于是, 我们想到了办法,将这个数字拆开,拆成一位一位的或者是四位四位的存储到一个数组中,用一个数组去表示一个数字。这样这个数字就被称谓是高精度数。 对于高精度数,也要像平常数一样做加减乘除以及乘方的运算,于是就有了高精度算法。
由于计算机输入计算结果的精度通常受到计算机的限制,如:在双精度方式下,计算机最多只能输出16位有效数字,如果超过16位,则只能按浮点形式输出,另外,一般计算机实数表示的范围为1038,如果超过这个范围,计算机就无法表示了。但是我们可以通过一些简单的办法来解决这个问题。这就是我们要说的高精度计算机。
高精度加法
题目:给定两个正整数(不含前导 0),计算它们的和。
输入格式
共两行,每行包含一个整数。
输出格式
共一行,包含所求的和
数据范围
1≤整数长度≤100000
为了方便进位,将输入的数字倒序存储到数组中。代码如下:
#include<iostream>
#include<string>
#include<vector>
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);
t/=10;
}
if(t) C.push_back(t);//处理进位
return C;
}*/
//方法二:省去了一部判断的步骤
vector<int> add(vector<int> &A,vector<int> &B)
{
if(A.size()<B.size())//提前判断大小
return add(B,A);
vector<int> C;
int t=0;
for(int i=0;i<A.size();i++)
{
t+=A[i];
if(i<B.size()) t+=B[i];
C.push_back(t%10);
t/=10;
}
if(t) C.push_back(t);
return C;
}
int main()
{
vector<int> A,B;
string 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);
for(int i=C.size()-1;i>=0;i--)
{
cout<<C[i];
}
return 0;
}
高精度减法
给定两个正整数(不含前导 0),计算它们的差,计算结果可能为负数。
输入格式
共两行,每行包含一个整数。
输出格式
共一行,包含所求的差。
数据范围
1≤整数长度≤105
不同于加法,高精度减法需要考虑减数和被减数的大小关系,如果被减数大于减数,正常计算,如果被减数小于减数,换位置进行计算,再在结果输出前边加上负号,代码如下:
#include<iostream>
#include<string>
#include<vector>
using namespace std;
//判断是否有A>=B
bool cmp(vector<int> &A,vector<int> &B)
{
//判断是否A>B
if(A.size()!=B.size()) return A.size()>B.size();
for(int i=A.size()-1;i>=0;i--)
{
if(A[i]!=B[i])
return A[i]>B[i];
}
//如果都没有A>B那么就是A=B
return true;
}
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]-t;
if(i<B.size()) t-=B[i];
C.push_back((t+10)%10);
if(t<0) t=1;
else t=0;
}
//去除前导零
while(C.size()>1&&C.back()==0) C.pop_back();//C.back()表示返回末尾元素的引用
return C;
}
int main()
{
vector<int> A,B;
string 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');
vector<int> C;
if(cmp(A,B))
{
C=sub(A,B);
}
else
{
C=sub(B,A);
cout<<"-";
}
for(int i=C.size()-1;i>=0;i--)
{
cout<<C[i];
}
cout<<endl;
return 0;
}
高精度乘法
给定两个非负整数(不含前导 0) A 和 B,请你计算 A×B 的值。
输入格式
共两行,第一行包含整数 A,第二行包含整数 B。
输出格式
共一行,包含 A×B 的值。
数据范围
1≤A的长度≤100000,
0≤B≤10000
首先看数据范围,A是高精度,B是低精度(高进度乘高精度是会死人的),所以B用整型表示即可不需要转换。
int的取值范围为:-2^31 ---- 2^31-1 ,即:-2147483648 - 2147483647
乘法的计算是分别拿大整数的每一位去和B相乘,保留个位加入到C中,然后相乘结果/10,进入到下次循环,代码如下:
#include<iostream>
#include<string>
#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;
}
while(C.size()>1 && C.back()==0)
{
C.pop_back();
}
return C;
}
int main()
{
vector<int> A;
string a;
int b;
cin>>a>>b;
for(int i=a.size()-1;i>=0;i--) A.push_back(a[i]-'0');
vector<int> C=mul(A,b);
for(int i=C.size()-1;i>=0;i--)
{
cout<<C[i];
}
cout<<endl;
return 0;
}
高精度除法:
给定两个非负整数(不含前导 0) A,B,请你计算 A/B 的商和余数。
输入格式
共两行,第一行包含整数 A,第二行包含整数 B。
输出格式
共两行,第一行输出所求的商,第二行输出所求余数。
数据范围
1≤A的长度≤100000,
1≤B≤10000,
B 一定不为 0
同高精度乘法,是高精度除以低精度,与前三种运算不同的是,需要传入额外的参数r表示余数
另外,由于除法是从高位开始算的所以在反转之后,i要从A.size()-1开始计算,比如123除以11,首先计算1除以11,余数为1,除得结果为0,将0加入到C中,更新余数r,进入下一次循环…代码如下
#include<iostream>
#include<string>
#include<vector>
#include<algorithm>
using namespace std;
vector<int> div(vector<int> &A,int b,int &r)//这里使用引用为了提高效率
{
vector<int> C;
r=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());//reverse()对指定区域进行翻转
while(C.size()>1&&C.back()==0) C.pop_back();//处理前导零
return C;
}
int main()
{
vector<int> A;
string a;
int b;
int r;//表示所求余数
cin>>a>>b;
for(int i=a.size()-1;i>=0;i--) A.push_back(a[i]-'0');
auto C=div(A,b,r);
for(int i=C.size()-1;i>=0;i--) cout<<C[i];
cout<<endl;
cout<<r<<endl;
return 0;
}