[转载] 数值的整数次方

题目:实现函数 double Power(double base、int exponent)。求base的exponent次方,不得使用库函数,
同时不需要考虑大数问题。

最朴素的求幂方法

也就是平常使用pow函数,最简单的实现就是一直累乘,可以得到这样的代码:

复制代码

 1 int Pow(int a,int b){
 2 
 3     int ans = 1;
 4 
 5     for(int i = 0;i < b;i++){
 6 
 7         ans *= a;
 8 
 9     }
10 
11     return ans;
12 
13 }

复制代码
可以看到,算法的时间复杂度是O(n)。为了降低时间复杂度,我们可以使用快速幂算法,将时间复杂度降低到O(logn),n是幂。

快速幂:

首先,快速幂的目的就是做到快速求幂,假设我们要求a^b,

 假设我们要求a^b,那么其实b是可以拆成二进制的,该二进制数第i位的权为2^(i-1),例如当b==11时

                             a11=a(2^0+2^1+2^3)
  11的二进制是1011,11 = 2³×1 + 2²×0 + 2¹×1 + 2º×1,因此,我们将a¹¹转化为算 a2^0*a2^1*a2^3,也就是a1*a2*a8
,看出来快的多了吧原来算11次,现在算三次
那么怎么算呢
可以考虑成根据二进制的权值来求解的。那么在关于位运算的部分,我们可以逐位获取b的位,碰到0,就累乘,
碰到1,就将累乘的值并且将乘到答案。由此可以得到代码:
代码如下:




```cpp
int poww(int a, int b) {
    int ans = 1, base = a;
    while (b != 0) {
        if (b & 1 != 0)
{
            ans *= base;
}
            base *= base;
            b >>= 1;
    }
    return ans;
}

按 Ctrl+C 复制代码
以b==11为例,b=>1011,二进制从右向左算,但乘出来的顺序是 a(20)*a(21)a(23),是从左向右的。我们不断的让base=base目的即是累乘,以便随时对ans做出贡献。

其中要理解base*=base这一步:因为 basebase==base2,下一步再乘,就是base2base2==base4,然后同理 base4base4=base8,由此可以做到base–>base2–>base4–>base8–>base16–>base32…指数正是 2^i ,再看上面的例子,a¹¹= a1a2*a8,这三项就可以完美解决了,快速幂就是这样。

🔺快速幂取模

根据之前的博客同余定理,我们知道

(ab)%m = ((a%m)(b%m))%m;

其实快速幂取模也是用到这个

那么根据上面的定理可以推导出另一个定理:

(a^b) mod c = (a * a * a…)%c = ((a%c)(a%c)(a%c)*…)%c = (a%c)^b %c;

这就是快速幂取模

代码如下:

复制代码

 1 int pow_mod(int a ,int b)
 2 {
 3     int ans = 1 ;
 4     int base = a % c;
 5         while(b>0)
 6     {
 7           if(b&1!=0)
 8             ans = (ans *base)%c;
 9     }
10      base = (base*base)%c;
11     b >>= 1;
12     return ans;
13 }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值