码迷,mamicode.com
首页 > 其他好文 > 详细

快速幂和快速乘

时间:2020-06-07 14:36:24      阅读:54      评论:0      收藏:0      [点我收藏+]

标签:container   判断   快速   计算   toc   div   算法分析   star   快速乘   

快速幂和快速乘

1. 算法分析

1.1 快速幂

计算a ^ k % p
把k拆成二进制表示形式,比如k等于5的时候,k = (101)2 = c1 * 1 + c2 * 0 + c3 * 1
预处理a ^ c1, a ^ c2, ..., a ^ ct
这样计算a^k%p时,答案即为res = (a^c1) * (a^c3)

1.2 快速乘

a * (2 ^ 0) = (2 ^ 0) * a;
a * (2 ^ 1) = (2 ^ 1) * a;
a * (2 ^ 2) = (2 ^ 2) * a;
a * (2 ^ 3) = (2 ^ 3) * a;
...
a * (2 ^ k) = (2 ^ k) * a , res1 = (2 ^ k) * a;
a * (2 ^ (k + 1)) = (2 ^ k + 2 ^ k) * a, res2 = (res1 + (2 ^ k)) * a;

可以看出2的(k+1)次是在2的k次的基础上再加上2的k次
比如求 a * 13 % p = a*(1 * 1 + 2 * 0 + 4 * 1 + 8 * 1) % p
即第1为乘1,第2为乘0,第3为乘1,第4为乘1,要判断每一位是否为1,如果为1,才需要

2. 板子

  1. 快速幂
LL qmi(LL a, LL k, LL p)  
{
    LL res = 1 % p;  // res记录答案, 模上p是为了防止k为0,p为1的特殊情况
    while(k)  // 只要还有剩下位数
    {
        if (k & 1) res = (LL)res * a % p;  // 判断最后一位是否为1,如果为1就乘上a,模上p, 乘法时可能爆int,所以变成long long
        k >>= 1;  // 右移一位
        a = (LL) a * a % p;  // 当前a等于上一次的a平方,取模,平方时可能爆int,所以变成long long
    }
    return res;
}
  1. 快速乘
// 快速乘(只能计算到a = 2^18, b = 2 ^ 18)
LL qmul(LL a, LL k, LL b)
{
    LL res = 0;
    while (k)
    {
        if (k & 1) res = (res + a) % b;
        a = (a + a) % b;
        k >>= 1;
    }
    return res;
}

快速幂和快速乘

标签:container   判断   快速   计算   toc   div   算法分析   star   快速乘   

原文地址:https://www.cnblogs.com/spciay/p/13060329.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!