快速幂，龟速乘，光速乘模板

最新推荐文章于 2023-09-06 19:57:00 发布

_小可爱

最新推荐文章于 2023-09-06 19:57:00 发布

阅读量489

点赞数 3

分类专栏：算法刷题文章标签：算法

本文链接：https://blog.csdn.net/m0_51156601/article/details/124181855

版权

算法刷题专栏收录该内容

6 篇文章 1 订阅

订阅专栏

本文解析了三种高效求解大整数幂运算的算法：快速幂的光速乘法解法，利用长整型和浮点数计算的巧妙方法；龟速乘法，适用于逐位相乘；以及常规的快速幂实现。通过对比，讲解了它们在处理10^18范围内的模运算时的效率提升。

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

快速幂

求 $a^b mod p (1 \le a, b, p \le 10^{18})$

ll quick_power(ll a, ll b, ll p) { // ll 是long long类型
    ll ans = 1 % p;
    for(; b; b >>= 1) {
        if(b & 1) ans = ans * a % p;
        a = a * a % p;
    }
    return ans;
}

龟速乘（解法一）

求 $\ mod p (1 \le a, b, p \le 10^{18})$

ll quick_mul(ll a, ll b, ll p) { // ll 是long long类型
    ll ans = 0;
    for(; b; b >>= 1) {
        if(b & 1) ans = (ans + a) % p;
        a = a * 2 % p;
    }
    return ans;
}

光速乘（解法二）

求 $\ mod p (1 \le a, b, p \le 10^{18})$

ull mul(ull a, ull b, ull p) { // ull 是 unsigned long long 类型，ll 是long long类型
    a %= p, b %= p;
    ull c = (long double)a * b / p;
    ull x = a * b, y = c * p;
    ll ans = (long long)(x % p) - (long long)(y % p);
    if(ans < 0) ans += p;
    return ans;
}