高精度运算合集，加减乘除，快速幂，详细代码，OJ链接

最新推荐文章于 2024-02-04 20:47:10 发布
最新推荐文章于 2024-02-04 20:47:10 发布
阅读量1k 收藏 21
点赞数 16
分类专栏： 数据结构与算法 文章标签： c++ 算法 数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/EQUINOX1/article/details/135788594
版权

文章目录

零、前言

高精度运算是某些题目涉及大数值运算且范围超出语言内置类型允许范围时采取的处理手段，原理就是小学列竖式计算的代码化，比较易懂，本文介绍高精度加法、减法、乘法、除法、以及高精度快速幂，附模板即OJ链接练手。

一、加法

高精度加法步骤

  • 倒序存储加数，低位放数组前面，高位放数组后面(类似小端存储)
  • 从低位到高位，累加，进位，求余
  • 答案数组从末尾到开头即为高位——低位

P1601 A+B

P1601 A+B Problem（高精） - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
#define IOTIE ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);

#define ll long long
const int N = 505;
int A[N]{0}, B[N]{0}, C[N]{0};
int la, lb, lc;

void add()
{
	//累加、进位、取余
    for (int i = 0; i < lc; i++)
        C[i] += A[i] + B[i], C[i + 1] += C[i] / 10, C[i] %= 10;

    if (C[lc]) //末位进位，位长+1
        lc++;
}

int main()
{
    IOTIE
    // freopen("in.txt", "r", stdin);
    string a, b;
    cin >> a >> b, la = a.size(), lb = b.size(), lc = max(la, lb);
    for (int i = 0; i < la; i++)
        A[la - i - 1] = (a[i] ^ 48);
    for (int i = 0; i < lb; i++)
        B[lb - i - 1] = (b[i] ^ 48);
    add();
    for (int i = lc - 1; i >= 0; i--)
        cout << C[i];
    return 0;
}

二、减法

高精度减法步骤

  • 倒序存储加数，低位放数组前面，高位放数组后面(类似小端存储)
  • 被减数A，减数B，如果A < B，交换A，B，输出负号
  • 从低位到高位，逐位求差，借位，存差
  • 答案数组从末尾到开头即为高位——低位

P2142 高精度减法

P2142 高精度减法 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
#define IOTIE ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);

#define ll long long
const int N = 10100;
int A[N]{0}, B[N]{0}, C[N]{0};
int la, lb, lc;

bool cmp()
{
    if (la != lb)
        return la > lb;
    for (int i = la - 1; i >= 0; i--)
        if (A[i] != B[i])
            return A[i] > B[i];
    return true;
}

void sub()
{
    for (int i = 0; i < lc; i++)
    {
        if (A[i] < B[i])//借位
            A[i] += 10, A[i + 1]--;
        C[i] = A[i] - B[i];
    }
    while (lc && !C[lc]) //前导0
        lc--;
}

int main()
{
    IOTIE
    // freopen("in.txt", "r", stdin);
    string a, b;
    cin >> a >> b, la = a.size(), lb = b.size(), lc = max(la, lb);
    for (int i = 0; i < la; i++)
        A[la - i - 1] = (a[i] ^ 48);
    for (int i = 0; i < lb; i++)
        B[lb - i - 1] = (b[i] ^ 48);
    if (!cmp())
        swap(A, B), cout << '-';
    sub();
    for (int i = lc; i >= 0; i--)
        cout << C[i];
    return 0;
}

三、乘法

高精度乘法步骤

  • 倒序存储加数，低位放数组前面，高位放数组后面(类似小端存储)
  • 从低位到高位，累加乘积，进位，求余
  • 答案数组从末尾到开头即为高位——低位

P1303 A*B

P1303 A*B Problem - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
#define IOTIE ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);

#define ll long long
const int N = 4010;
int A[N]{0}, B[N]{0}, C[N]{0};
int la, lb, lc;

void mul()
{
    for (int i = 0; i < la; i++)
        for (int j = 0; j < lb; j++)
            C[i + j] += A[i] * B[j], C[i + j + 1] += C[i + j] / 10, C[i + j] %= 10;
    while (lc && !C[lc])//前导0
        lc--;
}

int main()
{
    IOTIE
    // freopen("in.txt", "r", stdin);
    string a, b;
    cin >> a >> b, la = a.size(), lb = b.size(), lc = la + lb;
    for (int i = 0; i < la; i++)
        A[la - i - 1] = (a[i] ^ 48);
    for (int i = 0; i < lb; i++)
        B[lb - i - 1] = (b[i] ^ 48);

    mul();
    for (int i = lc; i >= 0; i--)
        cout << C[i];
    return 0;
}

四、除法

高精度除法步骤

  • 倒序存储加数，低位放数组前面，高位放数组后面(类似小端存储)
  • 从高位到低位，当前被除数，存商，求余数
  • 答案数组从末尾到开头即为高位——低位

P1480 A/B

P1480 A/B Problem - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
#define IOTIE ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);

#define ll long long
const int N = 10100;
int A[N]{0}, b, C[N]{0};
int la, lc;

void div()
{
    ll num = 0;
    for (int i = la - 1; i >= 0; i--)
        num = (num << 3) + (num << 1) + A[i], C[i] = num / b, num %= b;
    while (lc && !C[lc])
        lc--;
}

int main()
{
    IOTIE
    // freopen("in.txt", "r", stdin);
    string a;
    cin >> a >> b, la = a.size(), lc = la;
    for (int i = 0; i < la; i++)
        A[la - i - 1] = (a[i] ^ 48);

    div();
    for (int i = lc; i >= 0; i--)
        cout << C[i];
    return 0;
}

五、高精度快速幂

关于二分快速幂：二分快速幂和快读快写模板【C++】-CSDN博客

其实就是把普通二分快速幂中结果和底数相乘亦即底数自乘的部分换成高精度乘法即可，这里不给出代码只需要替换普通二分快速幂中的乘法为高精度乘法即可，我们直接看一道例题来练手。

麦森数

[P1045 NOIP2003 普及组] 麦森数 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

其实就是一个高精度快速幂裸题，由于题目只要求我们最后输出后500位，所以我们高精度乘法的时候只需要进行后500位的计算即可

那么位数如何求呢？

显然存在k ，10^k <= 2^p - 1 < 10^(k + 1)

显然2^p - 1 的位数就是k + 1，而2^p - 1和 2^p的十进制位数相同，那么有

10^k = [2 ^ p]，k = p * log10 (2)

我们的位数就是 p * log10 (2) + 1，log10可以用库函数

代码如下：

#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <cstdio>
#include <vector>
#include <algorithm>
using namespace std;
#define IOTIE ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
#define sc scanf
#define int long long
const int N = 500;
typedef vector<int> VI;
int p;
VI ans(N + 1), a(N + 1);

VI mul(VI &x, VI &y)
{
    VI ret(N + 1);
    for (int i = 0; i < N; i++)
        for (int j = 0; j + i < N; j++)
            ret[i + j] += x[i] * y[j], ret[i + j + 1] += ret[i + j] / 10, ret[i + j] %= 10;

    return ret;
}
void quickpower()
{
    ans[0] = 1, a[0] = 2;
    while (p)
    {
        if (p & 1)
            ans = mul(ans, a);
        a = mul(a, a), p >>= 1;
    }
    ans[0]--;
}
signed main()
{
    IOTIE
    // freopen("in.txt", "r", stdin);
    cin >> p;
    cout << floor(p * log10(2) + 1) << '\n';
    quickpower();
    for (int i = 499; i >= 0; i--)
    {
        cout << ans[i];
        if ((500 - i) % 50 == 0)
            cout << '\n';
    }
    return 0;
}
EQUINOX1
关注
  • 16
    点赞
  • 踩
  • 21
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
快速幂——高精度求幂
ivan_9的博客
01-31 4089
文章目录前言一、快速幂（Fast Exponentiation）的定义：二、快速幂原理：三、常规求幂：四、简单快速求幂：四、递归快速求幂：五、位运算快速求幂：六、高精度快速求幂：六、python高精度快速求幂：七、实战演习：总结 前言 本文讲述快速幂的原理，以及用法 一、快速幂（Fast Exponentiation）的定义： 定义：快速求，取base为底数的exp次幂，即求：baseexp； 时间复杂度： O(log₂N) 二、快速幂原理： 思想：每一步都把指数分成两半，而相应的底数做平方运算。不仅能把非
高精度 大数加法 乘法 除法 幂 模板
MyHeaven7的博客
08-17 547
转载请注明出处：http://blog.csdn.net/u012860063?viewmode=contents 一：加法 1、普通两个大数相加 代码如下： [cpp] view plain copy  print? #include    #include    #include    void fan(ch
POJ 1001——高精度快速幂
qq_43513855的博客
02-24 918
POJ 1001——高精度快速幂 本题主要用到两个东西，一个是高精度，一个是快速幂。 虽然是浮点数的乘法，但是实际上可以转化为整数的乘法，然后再用高精度乘法，高精度算法本身不难，但是在没重载运算符的情况下使用数组还是挺麻烦的，不少地方的初始化都需要特别注意。 关于数组的初始化： 数组名在函数内部时，它的sizeof数组名获得的是数组本身的长度，而数组作为函数形参的时候，sizeof获得的是这个指针本身的长度，这一点特别坑，这会导致memset函数出现问题。 关于快速幂： 快速幂是一种快速求出数字的幂的方法，
高精度求幂取模（附C++、python3代码）
沉默的大多数
02-08 1501
求x的n次方对mod（1e9+7）取模，当n巨大（）时，连快速幂取模也不行。此时用下面的方法： C++： #include &lt;bits/stdc++.h&gt; using namespace std; const int mod=1e9+7; int main(){ char n[100001]; int x; cin&gt;&gt;x&gt;&gt;n;/...
高精度、快速幂与快速排序
最新发布
m0_66115819的博客
02-04 622
高精度、快速幂、快速排序
大数高精度运算库
10-20
MapM库由Arthur O'Dwyer开发，提供包括加减乘除、取模、指数、对数、平方根等在内的多种大数运算功能。它的设计目标是能够在各种操作系统环境下运行，包括老旧的DOS系统以及Windows系列（如XP、2000、NT和9x）。 在...
C++实现的大整数加减乘除
02-28
C++实现的大整数加减乘法，已经测试了，并且在oj平台上测试通过
高精度幂次方计算 C++实现
12-27
在编程领域，高精度计算是处理大整数乘法、除法、加法、减法以及更复杂的数学运算的一种技术。对于题目“高精度幂次方计算 C++实现”，我们聚焦于实现一个C++程序，它能计算任意两个大整数的乘方，即求解形如`a^b`的...
C++实现oj算法代码-贪心算法.zip
10-28
贪心算法是一种在每一步选择中都采取当前状态下最好或...总之，这个压缩包提供了C++实现的oj算法代码，主要关注贪心算法。通过研究这些代码，可以加深对贪心算法的理解，提升C++编程技能，并为解决实际问题打下基础。
OI、OJ、ACM、PAT、CSP 题目常用代码模板
02-23
在编程竞赛领域，如OI（信息学奥赛）、OJ（在线判题系统）、ACM（国际大学生程序设计竞赛）、PAT（全国计算机等级考试）、CSP（中国计算机学会认证）等，参赛者需要掌握一系列高效的代码模板来快速解决各类问题。...
C++ 高精度算法（并附实例）
WIND_NAN的博客
07-26 5078
C++ 高精度算法（并附实例）
C++基础算法①——高精度加减法计算
博樽
03-14 6813
高精度导论、高精度+低精度、高精度+高精度、高精度减法。
高精度（加减乘除幂次）
weixin_72223753的博客
03-02 2123
len_a也直接当作ans的长度用了（主要是后面len_a就用不到了懒得多定义），这个有个不好的地方就是你要是后续还想对a进行什么别的操作就很不方便了，当然对于这个短小的代码就无伤大雅了，甚至还十分短小精悍方便很很。（如上述竖式计算中3+8=11保留一进一，在进行更高位2+1的计算中要加上进的这一位1）（表现在竖式计算中中间区段的存在，i+j有好几种组合，如1=1+0/1+1）对于超大数据的一种模拟加，减，乘，除，乘方，阶乘，开方等运算。
【算法】高精度幂
Dpi168的博客
02-04 1572
在高精度王国中，有个格外简单的算法，高精度幂，幂是什么？幂就是求次方，我敢说，高精度幂低精度，只要不考虑超时，比高精度加法还简单，因为就是不停地加自己，位数相同，高精度加法还得考虑位数不相同的，那么高精度幂该如何实现呢？请看下面： #include<iostream> #include<cstdio> #include<cmath> #include<cstring> using namespace std; string a="2",b; int c; in
【高精度】高精度减法（详细代码）
qq_40802813的博客
12-02 1323
【题目描述】 求两个大的正整数相减的差。 【输入】 共2行，第1行是被减数a，第2行是减数b(a > b)。每个大整数不超过200位，不会有多余的前导零。 【输出】 一行，即所求的差。 【输入样例】 9999999999999999999999999999999999999 9999999999999 【输出样例】 9999999999999999999999990000000000000
C++实现-高精度减法
store_的博客
04-14 385
if(a[i]<b[i]){//2-4---若现在这一位不够减，则向后一位借一。//若被减数小于减数，交换两个的值。//记录差中不是0的数字的下标。//不能直接使用L1，因为有可能交换了s1和s2的值。//长度相等 223 234 ASCII。//长度--12 123。//高精度减法的运算。
高精度减法 (C++)
Freddy1208的博客
05-28 374
输入两个大整数a，b，输出他们的差。(1<=a,b<=10^100)两个大非负整数a,b。(1<=a,b<=10^100)输出a减b的差。这道题依然是高精度，如果你会做高精度加法，那这题会变得很简单，不会也没关系，听作者慢慢分析。首先，仍然是令人头疼输入问题，用字符串型输入就行了，然后是怎么做减法，举个例子：111-19 1 1 1
高精度减法（C语言实现）
treble_rock的博客
11-05 246
高精度的本质是将数字以字符串的形式读入，然后将每一位分别存放入`int`数组中，通过模拟每一位的运算过程，来实现最终的运算效果。
求2的N次幂（C++）解决高精度运算
qq_63320529的博客
06-22 3674
计算高精度的2的N次方数字。
高精快速幂python
09-26
高精度快速幂是一种在计算大整数的幂时，使用快速幂算法的方法。Python中可以使用内置函数`pow()`来进行高精度快速幂的计算。以下是一个示例代码： ```python n, m = map(int, input().split()) sum = pow(m, n, 998244353) print(sum) ``` 在这个代码中，我们首先通过`map()`函数将输入的两个整数赋值给变量`n`和`m`。然后，使用`pow()`函数进行高精度快速幂的计算，其中第一个参数是底数m，第二个参数是指数n，第三个参数是模数998244353。最后，将计算结果打印出来。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

EQUINOX1

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值