快速幂
前提:
- 任何一个数k都可以由二进制表示,即
k = a l o g k 2 l o g k + . . . + a 1 2 1 + a 0 2 0 k = a_{logk}2^{logk} + ...+a_12^1 + a_02^0 k=alogk2logk+...+a121+a020
其中 a l o g k a l o g ( k ) − 1 . . . . a 1 a 0 a_{logk}a_{log(k)-1}....a_1a_0 alogkalog(k)−1....a1a0就对应了k的二进制!! - 一个数的幂可以由二分法拆成:
x k = x a 0 2 0 + a 1 2 1 + . . . + a l o g k 2 l o g k = ( x 2 0 ) a 0 ∗ ( x 2 1 ) a 1 ∗ . . . ( x 2 l o g k ) a l o g k x^k = x^{a_02^0 +a_1 2^1 +...+a_{logk}2^{logk}} = (x^{2^0})^{a^0}*(x^{2^1})^{a^1}*...(x^{2^{logk}})^{a^{logk}} xk=xa020+a121+...+alogk2logk=(x20)a0∗(x21)a1∗...(x2logk)alogk
所以不妨将 x 2 i x^{2^i} x2i全部预处理出来,至于要不要乘以该数,只需判断k对应的二进制位是否为1:
典型例题:
剑指offer 27. 数值的整数次方:
实现函数double Power(double base, int exponent),求base的 exponent次方。
不得使用库函数,同时不需要考虑大数问题。
注意:
不会出现底数和指数同为0的情况
当底数为0时,指数一定为正
样例1
输入:10 ,2
输出:100
样例2
输入:10 ,-2
输出:0.01
有个地方绕了本人一下,如果没有影响到读者,就当废话, x 2 0 x^{2^0} x20对应 x x x而非1
double Power(double base, int n) {
typedef long long LL;
double res = 1;
bool is_minus = n < 1; //判断正负
for(LL k = abs(LL(n));k ; k>>=1){
if(k&1)res*= base; //若当前k的位数为1,则乘以当前的base,第一次base对应
base *= base; //更新下一个x^2^i
}
if(is_minus) res = 1.0/res;
return res;
}