C程序-快速幂

1、给定三个正整数:a,b,c(a<10^{9}, b<10^{6}, 1<m<10^{9}),求a^{b}%m%m。

时间复杂度为O(b)的代码:

typedef long long LL;

LL pow(LL a,LL b,LL m)
{
    LL ans = 1;
    for(int i=0;i<b;i++)
    {
        ans=ans*a%m;
    }
    return ans;
}

代码中使用long long 而不用int 的原因是防止两个int 变量相乘后溢出。

1、给定三个正整数:a,b,c(a<10^{9}, b<10^{18}, 1<m<10^{9}),求a^{b}%m%m。

O(b)的复杂度连支持b<10^{8}都已经很难了,更何况10^{18}

这里要使用快速幂的思想,它基于二分的思想,因此也称为二分幂。

  1. 如果b是奇数,那么有a^{b}%m=a*a^{b-1}
  2. 如果b是偶数,那么有a^{b}%m​​​​​​​=a^{b/2}*a^{b/2}

显然,b是奇数的情况下总可以在下一步转换为b是偶数的情况,而b是偶数的情况总可以在下一步转换为b/2的情况。这样,在log(b)级别次数转换后,就可以把b变为0,而任何正整数的0次方都是1。

快速幂的递归写法,时间复杂度为O(log b)

typedef long long LL
LL binaryPow(LL a,LL b,LL m)
{
    if(b == 0) return 1;//如果b为0,那么a^0=1
    
    if(b % 2 ==1)//b为奇数,转为b-1
        return a * binaryPow(a, b-1, m) % m;
    else//b为偶数,转为b/2
    {
        LL mul = binaryPow(a, b / 2, m);
        return mul * mul % m;
    }
}

上面的代码中,条件if(b%2==1)可以用if(b&1)代替,因为b&1进行位于操作,判断b的末位是否为1,这样写执行速度会快一点。

还要注意,当b%2==0时不要直接返回binaryPow(a, b/2, m)* binaryPow(a, b/2, m),这是因为调用了两次binaryPow()函数,导致时间复杂度变成O(2^{log (b)}​​​​​​​)=O(b);

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值