首先声明,以下代码参考于AcWing.
高精度加法:
//c=A+B,A>=0,B>=0;
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;
}
高精度减法
//C=A-B,满足A>=B,A>=0,B>=0;
vector<int> sub(vector<int> &A,vector<int> &B)
{
vector<int>C;
for(int i=0,t=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();
return C;
}
高精度乘法
//C=A*b,A>=0,b>0
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;
}
return C;
}
高精度除法
//A/b==C...r,A>=0,b>0;
vector<int> div(vector<int>&A,int b,int &x)
{
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());
while(C.size()>1&&C.back()==0) C.pop_back();
return C:
}
3.1 高精度乘法
这个方法出自吴永辉老师。此代码简直让我拍手叫绝。
原理如下:
3 2 1 0 ——>数组a、b的下标
3 4 5 6 i ——>数组a[]
* 1 2 7 8 j ——>数组b[]
————————————
2 7 6 4 8
2 4 1 9 2
6 9 1 2
3 4 5 6
——————————————————
4 4 1 6 7 6 8 ——>数组c[]
6 5 4 3 2 1 0 i+j ——>数组c的下标
以上是俩个四位数相乘的竖式计算方法。可以看出,数的右面对齐,从低位向高位计算,计算结束后将一列结果相加即为答案。那么把俩个数从右向左依次标记为0、1、2…n,那么每一列的结果就是第一个数的下标为i的数与第二个数的下标为j的数相乘的结果,其存放在第i+j列。最终结果是每一列相加,就是i+j这一列所有数相加。所以可以用c[i+j] += a[i]*b[j]。
for(int i = 0;i < LA-1;i++)
for(int j = 0;j < LB-1;j++)
c[i+j] += a[i]*b[j];
for(int i = 0;i < LA+LB;i++)
if(c[i] >= 10){
c[i+1] += c[i]/10;
c[i] %= 10;
}