高精度模板

已于 2023-07-19 13:50:29 修改
阅读量198 收藏
点赞数
文章标签： 蓝桥杯 leetcode 算法 Powered by 金山文档
于 2023-03-18 15:11:29 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_58386305/article/details/129637167
版权

定义

高精度算法（High Accuracy Algorithm）是处理大数字的数学计算方法。在一般的科学计算中，会经常算到小数点后几百位或者更多，当然也可能是几千亿几百亿的大数字。一般这类数字我们统称为高精度数，高精度算法是用计算机对于超大数据的一种模拟加，减，乘，除，乘方，阶乘，开方等运算。对于非常庞大的数字无法在计算机中正常存储，于是，将这个数字拆开，拆成一位一位的，或者是四位四位的存储到一个数组中， 用一个数组去表示一个数字，这样这个数字就被称为是高精度数。高精度算法就是能处理高精度数各种运算的算法，但又因其特殊性，故从普通数的算法中分离，自成一家。

代码

本模板的数据输入

int main(){
    //读入
    string A,B;
    int b;
    vector<int> AI,BI;
    cin>>A>>B;
    for(int i = A.length()-1;i>=0;i--){
        AI.push_back(A[i]-48);
    } 
    for(int i = B.length()-1;i>=0;i--){
        BI.push_back(B[i]-48);
    } 
    // gjd_add 加法
    vector<int> C = gjd_add(AI,BI);
    for(int i = C.size()-1;i>=0;i--){
        cout<<C[i];
    }
    cout<<endl; 
    // gjd_sub 减法  
    C = gjd_sub(AI,BI);
    for(int i = C.size()-1;i>=0;i--){
        cout<<C[i];
    }
    cout<<endl; 
    // gjd_mul 乘法
    C = gjd_mul(AI,BI);
    for(int i = C.size()-1;i>=0;i--){
        cout<<C[i];
    }    
    cout<<endl;     
    // gjd_div 除法      
    vector<int> r; // 高精度除法返回的余数 
    C = gjd_div(AI,BI,r);
    for(int i = C.size()-1;i>=0;i--){
        cout<<C[i];
    } 
    cout<<endl; 
    for(int i = r.size()-1;i>=0;i--){
        cout<<r[i];
    } 
}

高精度数组的比大小

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> gjd_add(vector<int> &A, vector<int> &B){
    if (A.size() < B.size()) return gjd_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;
}

减法

vector<int> gjd_sub(vector<int> &A, vector<int> &B)
{     
    vector<int> C;
    if (!cmp(A,B)) 
    {
        C = gjd_sub(B,A);
        C[C.size() - 1] = -C[C.size() - 1];
        return 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;
    }
    while (C.size() > 1 && C.back() == 0) C.pop_back();
    return C;
}

乘法

vector<int> gjd_mul(vector<int> &A, vector<int> &B) {
    vector<int> C(A.size() + B.size() + 7, 0); // 初始化为 0，C的size可以大一点

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

    while (C.size() > 1 && C.back() == 0) C.pop_back(); // 必须要去前导 0，因为最高位很可能是 0
    return C;
}

除法

vector<int> gjd_div(vector<int> &A, vector<int> &B, vector<int> &r){
    vector<int> C;
    //判断A是否小于B 若小于则答案为0 求余数 
    if(!cmp(A,B)){
         C.push_back(0); 
         r.assign(A.begin(),A.end()); 
         return C; 
    }
    int j = B.size();
    r.assign(A.end()-j,A.end());
    while(j<=A.size()){
        int k=0;
        while(cmp(r,B)){
            vector<int> s = gjd_sub(r,B);
            r.clear();
            r.assign(s.begin(),s.end());
            k++;
        }
        C.push_back(k);
        if(j<A.size()) r.insert(r.begin(),A[A.size()-j-1]);
        if(r.size()>1&&r.back()==0) r.pop_back();
        j++;
    }
    reverse(C.begin(),C.end());
    while(C.size()>1&&C.back()==0) C.pop_back();
    return C;
}

最终代码（AC Code For 洛谷P1932 A+B A-B A*B A/B A%B Problem）：

#include<bits/stdc++.h>
using namespace std;
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> gjd_add(vector<int> &A, vector<int> &B)
{
    if (A.size() < B.size()) return gjd_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;
}
vector<int> gjd_mul(vector<int> &A, vector<int> &B) {
    vector<int> C(A.size() + B.size() + 7, 0); // 初始化为 0，C的size可以大一点

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

    while (C.size() > 1 && C.back() == 0) C.pop_back(); // 必须要去前导 0，因为最高位很可能是 0
    return C;
}
vector<int> gjd_sub(vector<int> &A, vector<int> &B)
{     
    vector<int> C;
    if (!cmp(A,B)) 
    {
        C = gjd_sub(B,A);
        C[C.size() - 1] = -C[C.size() - 1];
        return 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;
    }
    while (C.size() > 1 && C.back() == 0) C.pop_back();
    return C;
}

vector<int> gjd_div(vector<int> &A, vector<int> &B, vector<int> &r){
    vector<int> C;
    //判断A是否小于B 若小于则答案为0 求余数 
    if(!cmp(A,B)){
         C.push_back(0); 
         r.assign(A.begin(),A.end()); 
         return C; 
    }
    int j = B.size();
    r.assign(A.end()-j,A.end());
    while(j<=A.size()){
        int k=0;
        while(cmp(r,B)){
            vector<int> s = gjd_sub(r,B);
            r.clear();
            r.assign(s.begin(),s.end());
            k++;
        }
        C.push_back(k);
        if(j<A.size()) r.insert(r.begin(),A[A.size()-j-1]);
        if(r.size()>1&&r.back()==0) r.pop_back();
        j++;
    }
    reverse(C.begin(),C.end());
    while(C.size()>1&&C.back()==0) C.pop_back();
    return C;
}
int main(){
    //读入
    string A,B;
    int b;
    vector<int> AI,BI;
    cin>>A>>B;
    for(int i = A.length()-1;i>=0;i--){
        AI.push_back(A[i]-48);
    } 
    for(int i = B.length()-1;i>=0;i--){
        BI.push_back(B[i]-48);
    } 
    // gjd_add 加法
    vector<int> C = gjd_add(AI,BI);
    for(int i = C.size()-1;i>=0;i--){
        cout<<C[i];
    }
    cout<<endl; 
    // gjd_sub 减法  
    C = gjd_sub(AI,BI);
    for(int i = C.size()-1;i>=0;i--){
        cout<<C[i];
    }
    cout<<endl; 
    // gjd_mul 乘法
    C = gjd_mul(AI,BI);
    for(int i = C.size()-1;i>=0;i--){
        cout<<C[i];
    }    
    cout<<endl;     
    // gjd_div 除法      
    vector<int> r; // 高精度除法返回的余数 
    C = gjd_div(AI,BI,r);
    for(int i = C.size()-1;i>=0;i--){
        cout<<C[i];
    } 
    cout<<endl; 
    for(int i = r.size()-1;i>=0;i--){
        cout<<r[i];
    } 
}
爱吃葡萄的紫外线
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
C++高精度模板
03-05
支持（无限）高精度加、减、乘、除、模、大于、小于、大于等于、小于等于、位移、等于、不等于、++、--，cin>>，cout 总之，一切int类型的运算都支持
高精度（模板）
Coco_T的博客
10-05 551
今天写了一道需要高精度的模板题 结果高精乘写错一句，调了1个小时。。。STO 所以写了一遍常用的高精度，在这里再简单说一下高精度的思路：高精加模拟竖式计算 注意每次累加的时候都要加上余数d高精减有一个额外操作：判断两数大小，确定答案符号 第一标准：长度 第二标准每一位数的大小 注意：因为在存贮的时候是倒叙存储，所以在对比的时候也要从高位到低位倒叙对比 模拟竖式计算 借位的时候容易出错
EOJ 1200 聪明的张小强 大数
tjf的专栏
06-30 950
http://acm.cs.ecnu.edu.cn/problem.php?problemid=1200 题意：A+B 大数加法 #include #include #include #include using namespace std; //compare比较函数：相等返回0，大于返回1，小于返回-1 int compare(string str1,string str2) {
C++高精度模板（加减乘除模）
最新发布
02-04
C++竞赛高精度模板 优点：包含四则运算 抛弃传统数组，使用vector 缺点：乘法没有使用FFT算法，复杂度O(N^2) 在某些编辑器上可能无法使用资源中定义的常量，需要手动创建（CE的话可能是这个问题） 没有运算符...
高精度模板 c++实现
08-11
高精度c++模板，实测没有问题，可以放心使用，常数一般，代码可移植性和封装性较好
压位高精度模板，支持负数
04-25
自己写的压位高精度支持负数，保证无bug，单精高精运算时不会像其他人一样偷懒直接类型转换再运算，会很慢，可以+-*/%以及快速幂，压位大小可以自己改(不过输出流重载的地方也要改一下)
高精度(模板)
qq_44754132的博客
03-01 281
处理大数运算的时候都是用字符串存的，高级语言一般也都有大数类，分装成一个对象，实际上就是维护了一个字符串，运算就是两个字符串之间的操作，他的运算也有很多啊，毕竟数学运算就有很多，简单的有 大数加大数 大数减大数 大数乘小数 大数除小数 上述几种都是最常见也是最简单的，运算的代码只有几行，而且都是模拟字符串的方法实现。 实际上大小数可以统一起来，就像能用INT的地方用LONG LONG 必然是可...
为什么浮点类型变量不能直接比较？而是用一个范围？
doujianyoutiao的博客
10-24 1244
并不是不可以，而是不推荐，因为容易出错。从下面的例子中能看出（true打印出来是1，false打印出来是0），由于浮点数有精度的问题，数值结果很容易出现偏差。f2 - 81.2不等于0是因为，f2是float类型的，字面量81.2默认是double，所以由于81.2到float类型已经精度不准了，减去精准的81.2得到的结果肯定也不准了。f1/5=f3也是一样的道理，但由于f2和f3都是81.2转过去的，所以这两是相等的，f3跟81.2是不相等。==0也是一样的道理。如果要比较，结合数据类型和自己要求的精度
C++高精度模板实践：高效转换与操作
本文档主要介绍了如何在C++中实现高精度计算，特别是针对整数的高精度模板。作者提供了几个关键函数，如`convert_s`、`convert_b`、`convert_i`、`convert_b(LL)`、`comp`以及`addi`，它们用于处理大整数的表示、...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值