快速幂取余

应用:快速计算出a的n次方对mod取余

对于计算a的n次方,暴力的方法时间复杂度是O(n),对于n>1e8的情况下已经不能使用了

我们知道a的n次方可以写成:a^n=a^{b1}*a^{b2}*\cdot \cdot \cdot *a^{bk}(其中b1+b2+\cdot \cdot \cdot +bk=n

所以只要想办法找出b数组即可;我们考虑到一个数使用二进制表示的特点,对于一个数的二进制表示为1的位次分别是第i1、i2、i3、……ik,则n = 2^{i1-1}+2^{i2-1}+\cdot \cdot \cdot +2^{ik-1}  ,所以我们只需要从低位到高位依次找到二进制位是1,每找到一个为1的位就将结果更新。具体代码实现还有一些技巧,只需短短几行代码就能很快的解决这个问题

例如:n=13,其二进制表示为13=001101=2^3+2^2+2^0,b1=2^0、b2=2^2、b3=2^3

ans=a^{b1}*a^{b2}*a^{b3}=a^{0}*a^{4}*a^{8}

代码实现:

#include <bits/stdc++.h>
using namespace std;
#define maxn 21
typedef long long int ll;
const ll mod=1e9+7;
ll quick_mi(ll a,ll b)//a的b次方对mod取余
{
    ll tp=1;
    while(b)
    {
        if(b&1)//这一位是1,也即表示此时b是奇数
            tp=(tp*a)%mod;
        a=(a*a)%mod;
        b>>=1;//向右移一位,也即b/=2
    }
    return tp;
}
int main()
{
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值