基础算法之高精度

高精度

高精度算法用于长度为10^6数的加减乘除运算
为什么要用高精度算法呢？
int类型是四个字节，32位.
所以int的最大值位2^32-1=2147483647，大约是20亿，数据范围是 0~10 ^ 9
long类型是八个字节，64位
所以long的取值范围是-2^63 ~ 2 ^63-1,最大可以达到10 ^18.
例如阶乘运算，最后的结果可能会非常大，数的长度达到几十位，即使是最long long 也会爆，
因此对于大整数运算，我们可以用高精度算法。
高精度算法的本质还是模拟我们小学加减乘除运算的过程

高精度加法

1.大整数的存储
对于大整数，普通的long long也无法存储，所以我们可以用字符串的形式来存储大整数。
2.模拟加法运算的过程

#include <iostream>
#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(1);
    return C;
}
int main(){
    string a;
    string b;
    vector<int>A;
    vector<int>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;
}

高精度减法

和高精度加法一样，模拟减法过程

#include <iostream>
#include <vector>
using namespace std;
//判断A和B谁大谁小
bool cmp(vector<int>&A,vector<int>&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];
        }
    }
    return true;
}
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;
        //判断一下B[i]是否存在，因为B的位数可能比A的位数要小
        if(i<B.size()) {
            t-=B[i];
        }
        //这里的t有两种情况，要么大于0，要么小于0
        //(t+10)%10 相当于将两种情况合并
        c.push_back((t+10)%10);
        //这里根据t的值，判断要不要进位
        if(t<0) t=1;
        else  t=0;
    }
    //去掉前导0 
    while(c.size()>1 && c.back()==0) {
        c.pop_back();
    }
    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');
    }
    if(cmp(A,B)){
        auto C = sub(A,B);
        for(int i=C.size()-1;i>=0;i--){
            cout<<C[i];
        }
    }else {
        cout<<"-";
        auto C = sub(B,A);
        for(int i=C.size()-1;i>=0;i--){
            cout<<C[i];
        }
    }
    return 0;
}

高精度乘法

#include <iostream>
#include <vector>
using namespace std;
vector<int> mul(vector<int>& A,int b){
    //t表示进位
    int t=0;
    vector<int>c;
    //当i还没循环完或者进位还没有处理完时，一直执行循环
    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的时候，123*0=000，就要去除前导零
    while(c.size()>1 && c.back()==0) c.pop_back();
    return c;
}
int main(){
    string a;
    int b;
    cin>>a;
    cin>>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;
};

高精度除法

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
// A / d. . . r  
vector<int> div(vector<int>&A,int d,int &r){
    vector<int>c;//c为商
    for(int i=A.size()-1;i>=0;i--){
        r=r*10+A[i];
        c.push_back(r/d);
        r%=d;
    }
    //因为是高位在前，低位在后，所以要将它反转过来以便去除前导0
    reverse(c.begin(),c.end());
    while(c.size()>1 && c.back()==0){
        c.pop_back();
    }
    return c;
}
int main(){
    string a;
    int d;
    cin>>a;
    cin>>d;
    vector<int>A;
    for(int i=a.size()-1;i>=0;i--){
        A.push_back(a[i]-'0');
    }
    // r是余数
    int r=0;
    auto C = div(A,d,r);
    for(int i=C.size()-1;i>=0;i--){
        cout<<C[i];
    }
    cout<<endl<<r<<endl;
    return 0;
}