c++大数相乘

最新推荐文章于 2024-08-29 20:40:47 发布

真的不知道叫什么啊

最新推荐文章于 2024-08-29 20:40:47 发布

阅读量109

点赞数

分类专栏：算法文章标签： c++ 算法

本文链接：https://blog.csdn.net/m0_53237882/article/details/120146246

版权

算法专栏收录该内容

14 篇文章 0 订阅

订阅专栏

本文探讨了在处理大数乘法时，如何通过快速乘法技巧避免溢出问题，重点介绍了两种方法：二进制乘法和快速乘。通过实例展示了如何使用快速乘法函数ksc()来求解a*b对p的余数，并指出快速乘法在处理大规模数据时的优势。

摘要由CSDN通过智能技术生成

求 a 乘 b 对 p

取模的值。

输入格式

第一行输入整数a

，第二行输入整数b，第三行输入整数p

。

输出格式

输出一个整数，表示a*b mod p的值。

数据范围

1≤a,b,p≤1018

输入样例：

3
4
5

输出样例：

在写本题的时候，笔者最开始想到的是最原始的二进制乘法，这里初略介绍一下二进制乘法

二进制乘法

例：0010 * 1010，1010代表的是2 ^ 3 + 2 ^ 1, 我们将0010分别向左移动3位和移动1位，然后相加，就得出对应的结果

（参考文章：https://blog.csdn.net/wb175208/article/details/104278021?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522163093775216780357241029%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fall.%2522%257D&request_id=163093775216780357241029&biz_id=0&utm_medium=distribute.pc_search_result.none-task-blog-2~all~first_rank_v2~rank_v29-1-104278021.pc_search_result_cache&utm_term=%E4%BA%8C%E8%BF%9B%E5%88%B6%E4%B9%98%E6%B3%95&spm=1018.2226.3001.4187）

由此得出代码

#include<iostream>
#include<map>

using namespace std;

typedef long long LL;
LL a, b, mod;

LL muti(LL x, LL y)
{
    map<LL, int>bits;
    for(int i = 0; i < 64; i ++)
        bits.insert(pair<LL, int>(1 << i, i));//预处理2^n对应的n
    
    LL ans = 0;
    while(y > 0)
    {
        int bitshift = bits[y & (-y)];得出最后一位是二进制中从左到右的第几位
        ans = (ans % mod + (x << bitshift) % mod) % mod;//（a + b）% mod = (a % mod + b % mod) % mod
        y &= y - 1; //去掉最后一位
    }
    return ans;
}


int main()
{
    cin >> a >> b >> mod;
    cout << muti(a, b);
    return 0;
}

但是当测试数据太大的时候，会出现错误：

输入：

250182048980811753
413715569939057660
133223633696258584

输出：

75241317291714539

结果应该是

19308689043391716

这里的错误出现在x <<bitshift,x化为二进制数有47位，当bitshift足够大，用来存储左移后的x的匿名变量可能会将long long爆掉。在这里我们使用快速幂的亲戚——快速乘。

对于a*b，可以将b分解为2个二进制数x、y， a * b就化为 a * x + a * y ,再讲x转换为十进制x1，y1，就表示将x1个a和y1个a相加，这和快速幂相似，理解了快速幂的同学可能会更快理解快速乘。

对于的ac代码

#include<iostream>

using namespace std;

typedef long long LL;
LL a, b, mod;

LL ksc(LL x, LL y)
{
    LL sum = 0;
    while(y)
    {
        if(y & 1) sum = (sum + x) % mod;
        x = (x + x) % mod;
        y >>= 1;
    }
    return sum;
}

int main()
{
    cin >> a >> b >> mod;
    cout << ksc(a, b);
    return 0;
}