快速幂 快速乘原理讲解(模板)
目录
1 问题描述
2 原因分析
3 解决方法
4 快速幂讲解
5 快速乘讲解
6 完整代码
7 References
1 问题描述
我们发现,在 int 型下使用 pow 函数求 ,结果为 124 。
2 原因分析
pow 函数的返回值为 double型 ,因浮点数长度问题,存在截断误差。
3 解决方法
将变量定义为 double 型
有没有更快求幂的方法?
4 快速幂讲解
假设我们要求 ,按照朴素算法就是把 a 连乘 b 次,这样一来时间复杂度是 O(b) ,即是 O(n) 级别。但快速幂能做到 O(logn) 的复杂度。
快速幂:
对于二进制的位运算,我们需要用到 "&" 与 ">>" 运算符,详见 位运算符的应用
先上实现快速幂运算的具体代码:
long long ksm(long long a, long long b) {long long ans = 1, base = a;while(b != 0) {if(b & 1 != 0) {ans *= base;}base *= base;b >>= 1;}return ans;
}
其中 “b & 1” 指取 b 的二进制数的最末位,如 11 的二进制数为 1011 ,第一次循环,取的是最右边的 “1” ,以此类推。
而 “b >>= 1” 等效于 b = b >> 1,即右移 1 位,删去最低位。
以 a ^ 11 为例:
b 的二进制数为 1011,二进制从右向左算,但乘出来的顺序是 ,是从左向右的。我们不断的让 目的是累乘,以便随时对 ans 做出贡献。
要理解 这一步:因为 == ,下一步再乘,就是 () * () == ,然后同理 () * () == ,由此可以做到 → → → → → .......指数正好是 。再看上面的例子, = ,这三项就可以完美解决了,快速幂就是这样。
如还有不明白的地方,建议手动模拟代码的运行过程。
5 快速乘讲解
我们知道,在计算机中做加法运算会比乘法快得多(参考模电中的加法器),做乘法运算往往会溢出,即使用 long long 类型也拯救不了。因此需要寻找一种能高效完成乘法运算且不会溢出的算法,这就是快速乘算法。
快速乘与快速幂原理相似,也是将运算转换为二进制处理:
以 a * 11 为例: 11 的二进制数为 1011,则有
就是把快速幂中的 * 号改为+号
long long ksc(long long a, long long b) {long long ans = 0;while(b != 0) {if(b & 1 != 0) {ans += a;}a += a;b >>= 1;}return ans;
}
此版本的复杂度和快速幂一样,也是 O(logn) 。如果需要特别卡常数,可以去了解 O(1) 版本的快速乘。
6 完整代码
为了防止溢出,一般快速幂和快速乘的算法会在 mod 下运用,下面给出取模运算代码。
#include <bits/stdc++.h>
using namespace std;typedef long long ll;
const ll mod = 1e7;//a ^ b
ll ksm(ll a, ll b, ll mod) {ll ans = 1, base = a;while(b != 0) {if(b & 1 != 0) {ans = (ans * base) % mod;}base = (base * base) % mod;b >>= 1;}return ans;
}//a * b
ll ksc(ll a, ll b, ll mod) {ll ans = 0;while(b != 0) {if(b & 1 != 0) {ans = (ans + a) % mod;}a = (a + a) % mod;b >>= 1;}return ans;
}int main() {cout << "5 ^ 3 = " << ksm(5, 3, mod) << endl;cout << "345352 * 11 = " << ksc(345352, 11, mod) << endl;return 0;
}
运算结果:
7 References
- https://www.cnblogs.com/CXCXCXC/p/4641812.html
以上,有问题欢迎指正!
快速幂 快速乘原理讲解(模板)相关推荐
- 快速幂算法的原理及实现
幂运算,即次方运算,例如计算 的值即是幂运算,在实现的时候我们往往是这样写的: int __pow(int a,int b){int ans = 1;while(b--){ans *= a;}retu ...
- 【牛客每日一题】4.16 逆序对 ( 数学 , 排列组合 ,快速幂 , 快速乘 )
[每日一题]逆序对 链接:https://ac.nowcoder.com/acm/problem/14731 来源:牛客网 题目描述 求所有长度为n的01串中满足如下条件的二元组个数: 设第i位和第j ...
- 快速幂----快速求解底数的n次幂
目录 一.快速幂 1.问题的引入 2.快速幂的介绍 3.核心思想 4.代码实现 二.Pow(x, n) 1.题目描述 2.问题分析 3.代码实现 三.猴子碰撞的方法数 1.题目描述 2.问题分析 3. ...
- 【HDU - 5187】zhx's contest (快速幂+ 快速乘,模板)
题干: 2018百度之星复赛晋级名单出炉(增加20%晋级名额)~ zhx's contest Time Limit: 2000/1000 MS (Java/Others) Memory Limi ...
- 快速幂||取余运算【模板】(洛谷P1226题题解,Java语言描述)
题目要求 P1226题目链接 分析 标准的快速幂取模算法板子,之前这个算法我在这篇文章中讲过了:<快速幂算法详解&&快速幂取模算法详解>. 这里选择使用比较简单的API实现 ...
- HDU4549(矩阵快速幂+快速幂)
f(n)=a^f(n-1) + b^f(n-2):计算矩阵部分用矩阵快速幂:计算a的幂次和b的幂次用快速幂. #include<iostream> #include<algorith ...
- 快速幂(快速幂取余)
引入公式 (a*b) %c = ((a % c)*(b % c)) %c 普通求幂的解法 public static int pow(int x,int n) {int result = 1;for ...
- 数学--数论--快速幂--最大公约数--位运算模板
ACM常用模板合集 //位运算求解最大公约数 long long gcd(long long a,long long b) { if(a<b) return gcd(b,a);if(b==0) ...
- 【BZOJ-4522】密钥破解 数论 + 模拟 ( Pollard_Rho分解 + Exgcd求逆元 + 快速幂 + 快速乘)...
4522: [Cqoi2016]密钥破解 Time Limit: 10 Sec Memory Limit: 512 MB Submit: 290 Solved: 148 [Submit][Stat ...
最新文章
- 比尔盖茨宣布离开微软董事会:昔日全球首富致力于改变世界
- 浅析C#中的套接字编程
- 洛谷 2820 局域网
- buaaoo_fourth_assignment
- 机器学习(四)——SVM(2)
- Qt学习笔记-QSS装饰控件
- CV Code|计算机视觉开源周报20200503期
- Ubuntu 20.04 LTS 开发周期的重要任务:移除 Python 2
- 查找重复文件_重复文件查找和磁盘整理工具:Tidy Up
- PE文件格式学习之PE头移位
- 【Vue2.0】—Vue脚手架配置代理(二十二)
- 实体店想多赚钱就要学会互联网思维
- 5G 是否能让国产手机回到群雄割据时代?
- 《游戏改变企业》一一第 2 章 实时的经济信息反馈(为什么游戏很重要) 游戏改变企业...
- 计算机三级数据库要安装什么软件,计算机三级数据库技术题库app
- win10与win7系统之间文件共享
- Python开发网站步骤
- JAVA基础-多线程中锁机制
- 数据传输加密技术分析
- leetcode每日一题--雀巢原理;抽屉算法;Quorum机制;分布式应用