高精度的四则运算

高精度

• 高精度四则运算通过对数位分别储存运算，并维护一个用于处理进位的变量，从而在不造成溢出的前提下实现大数四则运算。
• Java和Python中虽然不需要考虑大整数四则运算的溢出问题，但这种算法的思想很值得借鉴。

高精度储存

• 从高位对应着 数据中表示也是高位。
• 当然输出的时候需要倒叙输出

int main()
{
    // 使用字符串读入
    string a, b;
    vector<int> A, B;
    cin >> a >> b; //a = "123456"
    // 使用vector逆序读入,变成整数需要减去偏移量0
    for (int i = a.size() - 1; i >= 0; i--) A.push_back(a[i] - '0');//A = [6,5,4,3,2,1]
    for (int i = b.size() - 1; i >= 0; i--) B.push_back(b[i] - '0');
    // 相当于vector<int>
    auto C = add(A, B);
    // 倒序输出
    for (int i = C.size() - 1; i >= 0; i--) printf("%d", C[i]);
    return 0;
}

高精度加法

很简单。就是模拟竖式相加。注意：

• t为进位
• 从低到高位计算 c[i] = (A[i] + B[i] + t) % 10, 进位 t = (A[i] + B[i] + t) / 10; 计算完之后如果 t大于0，记得进位加到 c数组最后。

// 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;
}

高精度减法，A - B(注意需要保证 A >= B)

• 从低位到高位计算 c[i] = (A[i] - B[i] - t) >= 0 ? (A[i] - B[i] - t) % 10 : ((A[i] - B[i] - t) + 10) % 10, 借位 t = (A[i] - B[i] - t) >= 0 ? 0 : 1
• 注意 : 1、答案有可能是负数
  ​ 2、最后存储在 vector 中的字符串可能会存在前导 0

// C = A - B, 满足 A >= B, A >= 0, B >= 0
bool cmp(vector<int>& A, vector<int>& B){ // 判断A是否>=B
    if(A.size() != B.size()) // 先比较A和B的位数
        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;
    int t = 0;
    for(int i = 0; i < A.size(); i ++){
        t = A[i] - t;
        if(i < B.size()) // 如果B里面还有数
            t -= B[i];
        C.push_back((t + 10) % 10); // 若t < 0, 就应该取(t + 10), 否则，取t, 合并起来就是((t + 10) % 10)
        if(t < 0) // 如果有借位
            t = 1;
        else
            t = 0;
    }
    while(C.size() > 1 && C.back() == 0) // 处理前导0
        C.pop_back();
    return C;
}

高精度乘法

以123 x 12 为例子

• 123的每一位都乘以 12
• c为每一位要输出的数字
• t为进位
• 可能有前导0,需要去除

// 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;
    }

    while (C.size() > 1 && C.back() == 0) C.pop_back();// 去除前导0

    return C;
}

高精度除法

• 需要带一个余数 r
• 下一次除的被除数就是 r*10 + Ai
• 从高位开始除
  

// A / b = C ... r, A >= 0, b > 0
vector<int> div(vector<int>& A, int& b, int& r){ // 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) // 去除前导0
        C.pop_back();
    return C;
}