算法--高精度计算

算法基础系列

---

前言

这里的高精度计算(大整数计算)只涉及以下内容(设A,B为大整数,a为较大数)
A+B，A-B，Axa，A/a

思想

数据存储

利用数组，将大整数一位位存入数组中，低位到高位依次存入数组中

Q：为什么数组的从低位开始存储？
A：方便进位，在数据的末尾加位置容易

计算方法——模拟竖式计算

回忆小学的竖式计算，用代码模拟其过程

高精度加法

用一个字母t 表示是否进位

高精度减法

前提：保证A >= B (若 B > A，按A - B 计算，最后加负号)
分两步，第一步判断Ai - Bi 情况
第二步，若大于等于0,则正常减去Ai - Bi -t
若 小于0，则向上一位借位 即：Ai - Bi + 10 -t

高精度乘法

和前面思路类似，按位计算，用一个字母t表示进位
t += A[i] * b
t % 10 表示当前位上数字
t /= 10 表示进位

高精度除法

除法是从高位开始计算，结果要翻转一下
余数r
用r * 10 + A[i] 表示每一位的被除数
然后用 r / b 表示每一位商
r % b 表示每一位余数

代码模板

用Vector存数据

加法

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

减法

vector<int> sub(vector<int> &A, vector<int> &B)
{
    vector<int> C;
    // A > B
    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;
}

乘法

vector<int> mul(vector<int> &A, int b)
{
    vector<int> C;
    int t = 0; // 进位
    if(b == 0) // 特判一下 b 是否为0
    {
        C.push_back(0);
        return C;
    }
    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;
}

除法

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

题

高精度加法

传送门

代码：

#include <bits/stdc++.h>

using namespace std;

const int N = 1e6 + 10;

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

高精度减法

传送门
代码：

#include <bits/stdc++.h>

using namespace std;

bool cmp(vector<int> &A, vector<int> &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;
    // A > B
    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;
}

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);
        cout << '-';
        for (int i = C.size() - 1; i >= 0; i--)
            printf("%d", C[i]);
    }
    return 0;
}

高精度乘法

传送门

代码：

#include <bits/stdc++.h>

using namespace std;

const int N = 1e5 + 10;

vector<int> mul(vector<int> &A, int b)
{
    vector<int> C;
    int t = 0; // 进位
    if(b == 0) // 特判一下 b 是否为0
    {
        C.push_back(0);
        return C;
    }
    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;
}

高精度除法

传送门

#include <bits/stdc++.h>

using namespace std;

const int N = 1e5 + 10;

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;

    vector<int> A;

    cin >> a >> b;

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