快速幂是非常实用的一种工具,也是一项很基础的技能
对于幂运算,我们可以这么理解:
设有b,p,k三个数值,求出 b^p mod k 的值
首先,看到这个问题,脑子中便想到了暴力,一重循环,一边进行乘法运算,一边模运算,爽歪歪
但是当数值更大的时候,会怎么样呢? 用暴力为正解吊唁
线性如果也会超时呢?
得想出更快的解法了 分治大法好
对于b^p,我们何尝不能想象为 b^(p/2) * b^(p/2) 呢?
由此我们又可以分为 b^(p/4) * b^(p/4) * b^(p/4) * b^(p/4)
总结以上,便可以得到:b^p = b^(p/n) * b^(p/n) * … *b^(p/n) (运算个数为log2为底n的对数次)
这何尝不是一种分解,分解,再分解的运算呢?
好吧,这就是分治。 分解的过程,我们便可以用递归来解决
使用递归的话,首先重要的便是边界了
对此,我们可以抉择出两种最特殊,最简单的边界
if (p == 0) return 1; //所有数的0次方都是1
if (p == 1) return b; //所有数的一次方无非是原数
当然,我们现在已经解决了最简的子问题了,便是要向上“升华”一下了——求部分的幂值
由于是递归的方法,便是可以直接的调用来层层分解了
但是最重要的,我们要搞清楚,我们分解的是什么?——是p,是b的指数
对于不同性质的p我们是要有不同的解决策略的
我们可以将p分成奇数和偶数来解决
对于比较不同得到奇数,我们进行不同的运算
姑且设t为b的一个局部的,进行拆分过的(运算完成所有偶数部分的)幂值,每一个奇数,对于2的取余的值,便都是1,那么,我们可以想象目前运算的,是最后一步,那么当前的值,无非是 ttb
之后,便是偶数的情况了
偶数时是最简易的情况了,可以直接为 t*t
把以上的汉语用程序语言表达出来
便是
if (p % 2 == 0) return t * t; //偶数情况
else
return t * t * b; //奇数情况
但是,我们还有一步不能省略,便是重要的取模了,
用指数级的大整数进行运算的话,算到中间还没有取模的时候就已经爆了
在这里分享一个取模的公式:
b * b mod k = (b mod k) * (b mod k) mod k
完整的代码如下
#include <bits/stdc++.h>
using namespace std;
long long b,p,k;
long long counting(long long b,long long p,long long k)
{
if (p == 0) return 1 % k; //边界
if (p == 1) return b % k;
long long t = counting(b,p/2,k); //算局部幂值
if (p % 2 == 0) return (t % k) * (t % k) % k; //综合“升华”
else
return (t % k) * (t % k) % k * (b % k) % k;
}
int main()
{
cin>>b>>p>>k;
cout<<counting(b,p,k);
}