快速幂——原理及实现

这篇文章讲一下快速幂的问题；

首先问一个简单的问题：
2³是几？
很简单啊，是不是？答案是8；
那么是怎么得来的呢？
222=8；
连续乘了3次2；

再看下边的问题：
2¹³=？
这个其实也很简单，不就是13个2相乘吗，连续计算13次就OK啦；

但是大家有没有想过，其实我们可以减少运算次数，更快的得到答案；
怎么做呢？

先看一下2¹⁶；
2¹⁶ = 2⁸ * 2⁸
对于上式，我们可以先运算8次得到2⁸，然后两个2⁸相乘得到2¹⁶；这样计算步骤由16个2连乘的16步降到了9步(8个2连乘，再平方)；
继续：2⁸ = 2⁴ * 2⁴;
计算几次？
5次！；
然后，又有2⁴ = 2² * 2²;
计算3次；
当然，2²就是计算1次啦；
综上：2¹⁶ 最少可以由4次计算得来： (((2²)²)²)² 注：16=2⁴；

再看2¹³怎样计算？？？
2¹³ = 2⁸ * 2⁴ * 2¹；
几次？3+3=6次；
怎么得到的呢？
在计算2⁸的过程中可以得到2⁴, 2¹,然后按个数相乘；
8, 4, 1是由13=2³ + 2² + 2¹ 得来；

大家可以发现，这样计算要比朴素的连乘计算要简单快速的多；不吗？你试试计算2¹⁰⁰⁰⁰,朴素的连乘计算要10000次；而上面的方法只要大约log₂10000次；

这就是快速幂算法；它把幂次运算由O(n)的复杂度简化到O(logn)；

下面推广到a^b;
首先令b = p₀2^k + p₁2^k-1 + … + p_k2⁰;
a^b = a^p₀2k * a^p₁2k-1 * … * a^p_k20;

下面就直接上代码啦：

typedef long long ll;
const ll mod=1e9+7;
ll _power(ll a, int b)//计算a^b;
{
	ll ans=1, res=a;
	while(b){
		if(b&1) ans=ans*res%mod;
		res=res*res%mod;
		b>>=1;
	}
	return ans%mod;
}