快速冪算法

求a^b%m的值。
当b是奇数时，有a ^ b = a * a ^(b-1)
当b是偶数时，那么有a ^ b= a ^ (b/2) * a ^ (b/2)
根据这两个条件写递归式。

#include<iostream>
using namespace std;
#define ll long long
ll a,b,m;
ll func(ll a,ll b,ll m){
 if(b==0)return 1;
 if(b&1)return a*func(a,b-1,m)%m;
 return func(a,b/2,m)*func(a,b/2,m)%m;
}
int main(){
cin>>a>>b>>m;
cout<<func(a,b,m)<<endl;
return 0;
}

快速幂的迭代写法
对于 a ^ b来说，若果把 b 写成2 进制，那么b 就可以写成若干二次幂之和，如13 的二进制 1101，于是3 号位 、2号位、0号位就都是1，那么就可以得到13 = 2^3 + 2^2 + 2^1 = 8 + 4 + 1。所以a ^13 = a^8 * a^4 * a^1。

通过同样的推导，我们可以把任意的a ^ b 表示成 a ^ (2 ^ k)……、a ^ 8、a^ 4、a^ 2、a^ 1中若干的乘积。若果二进制的i号位为1.那么想中的a^ (2^ i)就被选中。于是可以得到计算a^ b的大致思路：令i 从0到k枚举b的二进制的每一位，如果为1 那就累计a^(2 ^i)。注意
a^ (2^ k)……、a^ 8、a^ 4、a^ 2、a^ 1前一项总是等于后一项的平方（正是因为如此，每对b的二进制判断一位数，底数a就要平方）。具体步骤。
（1）初始令ans = 1,用来存放累积的结果。
（2）判断b的二进制末尾是否为1 ，（及判断 b&1 是否为 1），也可以理解为判断b 是否为奇数。如果是的话，令ans乘上a的值。
（3）令a平方，并使b右移一位，（也可以理解为，b/2）
（4）只要b 大于0，就返回（2）。

#include
using namespace std;
#define ll long long
ll a,b,m;
int main(){
cin>>a>>b>>m;
int ans=1;
while(b){
if(b&1)ans=ans*a%m;
a=a*a%m;
b=b>>1;
}
cout<<ans<<endl;
return 0;
}