快速幂

31 篇文章 0 订阅
5 篇文章 0 订阅

快速幂是非常实用的一种工具,也是一项很基础的技能

对于幂运算,我们可以这么理解:
设有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);
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值