【算法基础——第五讲】高精度

ฅ(๑˙o˙๑)ฅ 大家好, 欢迎大家光临我的博客：面向阿尼亚学习
算法学习笔记系列持续更新中~

---

---

一、前言

由于在C++中没有处理大整数的类，我们需要用字符串string或者c++STL容器，来处理大整数的加减乘除

---

二、高精度的介绍

高精度算法，属于处理大数字的数学计算方法。在一般的科学计算中，会经常算到小数点后几百位或者更多，当然也可能是几千亿几百亿的大数字。一般这类数字我们统称为高精度数，高精度算法是用计算机对于超大数据的一种模拟加，减，乘，除，乘方，阶乘，开方等运算。

三、高精度的运用

1. 高精度加法模板

大整数的存储：由于进位的缘故，我们需要将整数逆序的读入vector容器当中
大整数的计算：模拟人工计算，从末位开始加减，用取余的方式进行进位。如果最后还有余数，则最后一位需进1
如下：

// C = A + B, A >= 0, B >= 0
// big number
#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, 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);
    for (int i = C.size() - 1; i >= 0; i--) printf("%d", C[i]);
    return 0;
}

2.高精度减法模板

主要的难点在于借位与进位
如下：

// C = A - B, 满足A >= B, A >= 0, B >= 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;
}
// A-B，在A>=B的前提下
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;
    }
    // 避免出现类似001的情况
    while (C.size() > 1 && C.back() == 0) C.pop_back(); //去掉前导0
    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--) printf("%d", C[i]);
    } else {
        auto C = sub(B, A);
        printf("-");
        for (int i = C.size() - 1; i >= 0; i--) printf("%d", C[i]);
    }
    return 0;
}
}

3.高精度乘低精度

注意如t还未变为0时，需要重复执行
如下:

// A / b = C ... r, A >= 0, b > 0
#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;
   }
    return C;
}
int main() {
    string a;
    int b;
    vector<int> A;
    cin >> a >> b;

    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--) printf("%d", C[i]);

    return 0;
}

4.高精度除以低精度

余数需要同时传入r的引用
C最后要进行一次倒置，来保证与其它四则运算输出兼容
如下:

// A / b = C ... r, A >= 0, b > 0
#include <iostream>
#include <vector>
#include <algorithm>

using namespace std;

// A / b，余数r，商为C
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());
    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 = div(A, b, r);

    for (int i = C.size()-1; i >= 0; i--) printf("%d", C[i]);
    cout << endl << r << endl;

    return 0;
}

---

最后

莫言真理无穷尽，寸进自有寸进欢