关于快速幂

对于一个数求幂,我们有下面的公式

当n为偶数,a^n=a^{n/2}a^{n/2}

当n为奇数,a^n=aa^{n/2}a^{n/2}

这里的除法是c++里面的整数除法。

快速幂的代码:

float power(const float x, const int k){
    float result = 1, base = x;
    int p = (k > 0) ? k : -k;
    while(p){
        if(p & 1) result *= base;
        base *= base;
        p = (p >> 1);
    }
    return (k < 0) ? 1 / result :
        (k > 0) ? result : 1;
}

一开始看这个代码,想不明白为什么指数幂倒着除,而底数正着往上乘。想了半天,最后搜维基百科发现快速幂的另一个写法:

当n为偶数a^n=(a^2)^{n/2}

当n为奇数a^n=a(a^2)^{n/2}

也就是说,当前的幂底下,result要乘上一个当前的底数,然后更新当前的底数为a^2.

最后,由于所有的数最后除法都会到1或者-1,所以最后的值会乘一个当前的底数base,最后跳出循环输出。

这里,result就是记录奇数情况下,多出来的那个底数,base就是不断地更新平方后的底数。

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值