高精度算法

高精度

高精度适用情况

大数的长度最长为 10^6位

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

高精度储存

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

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-HGw3PIPX-1641802022731)(…/…/…/图库/1590506552684.png)]

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>
    vector<int> 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; // vector C 储存和
    int t = 0;
    for (int i = 0; i < A.size(); i ++ )
    {
        t += A[i];//进位累加
        if (i < B.size()) t += B[i]; // 数B没了
        C.push_back(t % 10);
        t /= 10;
    }

    if (t) C.push_back(t);//最高位如果进位
    return C;
}

例题

https://www.acwing.com/problem/content/submission/code_detail/1607027/

#include <bits/stdc++.h>
using namespace std;
// 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; // vector C 储存和
    int t = 0;
    for (int i = 0; i < A.size(); i ++ )
    {
        t += A[i];//进位累加
        if (i < B.size()) t += B[i]; // 数B没了
        C.push_back(t % 10);
        t /= 10;
    }

    if (t) C.push_back(t);//最高位如果进位
    return C;
}
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>
    vector<int> C = add(A, B);
    // 倒序输出
    for (int i = C.size() - 1; i >= 0; i--) printf("%d", C[i]);
    return 0;
}

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

A-B中两种情况:

image-20210519194622060

A和B对位相减的话,两种情况:

image-20210519194706568

  • 从低位到高位计算 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、答案有可b能是负数
    ​ 2、最后存储在 vector 中的字符串可能会存在前导 0
#include <iostream>
#include <vector>
using namespace std;

const int N = 1000010;
//判断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;
}
//去前导零
void trimZero(vector<int> &A)
{
    while (A.back() == 0 && A.size() > 1) A.pop_back();
}

vector<int> sub(vector<int> &A, vector<int> &B)
{
    vector<int> C;//C为A-B的差
    int t = 0;
    for (int i = 0; i < A.size(); i++)
    {
        t = A[i] - t;//减去借位
        if (i < B.size()) t -= B[i]; //前两行实际是 t = A[i] - t - B[i]
        C.push_back((t + 10) % 10);//加上借位,如果t为负,那需要借位的
        if (t < 0) t = 1; //向前借位
        else t = 0;
    }
    trimZero(C);//去除前导0

    return C;
}

int main()
{
    string a, b;
    cin >> a >> b;
    vector<int> A, B, C;
    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)) C = sub(A, B); //若B比A大,将A-B变为B-A,加上负号
    else {
        C = sub(B, A);
        printf("-");
    }
    for (int i = C.size() - 1; i >= 0; i--) cout << C[i];//倒序输出

    return 0;
}

高精度乘法

以123 x 12 为例子

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

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-dgIdxush-1641802022737)(…/…/…/图库/1590450263226.png)]

高精度乘以高精度



#include <iostream>
#include <vector>

using namespace std;

vector<int> mul(vector<int> &A, vector<int> &B) {
    // 初始化为 0,且999*99最多 5 位
    vector<int> C(A.size() + B.size(), 0); 

    for (int i = 0; i < A.size(); i++)
        for (int j = 0; j < B.size(); j++)
            C[i + j] += A[i] * B[j];

    int t = 0;
    for (int i = 0; i < C.size(); i++) { 
        // i = C.size() - 1时 t 一定小于 10
        t += C[i];
        C[i] = t % 10;
        t /= 10;
    }
	// 必须要去前导 0,因为最高位很可能是 0
    while (C.size() > 1 && C.back() == 0) C.pop_back(); 
    return C;
}

int main() {
    string a, b;
    cin >> a >> b; // a = "1222323", b = "2323423423"

    vector<int> 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');

    vector<int> C = mul(A, B);

    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) {
    vector <int> C;

    int t = 0;
    for (int i = 0; i < A.size(); i ++) {
        t += A[i] * b;       // 例如:t + A[i] * b = 7218
        C.push_back(t % 10); // 只取个位 8
        t /= 10;             // 721 看作 进位
    }

    while (t) {            // 处理最后剩余的 t
        C.push_back(t % 10);
        t /= 10;
    }

    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');

   vector <int> C = mul(A, b);

    for (int i = C.size() - 1; i >= 0; i --) {
        cout << C[i];
    }

    return 0;
}

高精度除法

需要带一个余数 r
下一次除的被除数就是 r*10 + A[i]
r = (r*10 + A[i])% b
C[i] =(r*10 + A[i])/ b
从高位开始除

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-iXXvXAe5-1641802022738)(…/…/…/图库/1590451533408.png)]

// A / b = C ... r, A >= 0, b > 0
vector<int> div(vector<int>& A, int& b, int& r){ // r引用
    vector<int> C;
    r = 0;//r为余数
    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;
}
 
// 高精度除法
#include <iostream>
#include<vector>
#include<algorithm>
using namespace std;
// A被除数 b为除数 r为余数
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];// 余数*10+当前这一位
        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; //余数
    vector<int> C = div(A, b, r);
    for (int i = C.size() - 1; i >= 0; i--) {
        cout << C[i];
    }
    cout << endl << r << endl;
    return 0;
}

  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值