当前位置: 首页 > ds >正文

快速幂算法详解:从暴力到优雅的数学优化

文章目录

  • 一、朴素幂运算的问题
  • 二、快速幂的数学原理
  • 三、快速幂的递归实现
  • 四、快速幂的迭代实现
  • 五、模运算下的快速幂
  • 六、快速幂的应用场景
  • 七、总结

快速幂是一种高效计算幂运算的算法,能够将时间复杂度从朴素的 O (n) 降低到 O (log n)。本文将深入探讨快速幂的原理、实现和应用场景。

一、朴素幂运算的问题

计算 a^n 最直接的方法是循环 n 次:

long long power(long long a, long long n) {long long result = 1;for (int i = 0; i < n; i++) {result *= a;}return result;
}

这种方法的时间复杂度为 O (n),当 n 非常大时(如 10^9),计算效率极低,甚至可能超时。

二、快速幂的数学原理

快速幂的核心思想是利用指数的二进制分解。例如,计算 3^13:

  1. 将指数 13 转换为二进制:13 = 1101 (2) = 8 + 4 + 1
  2. 3^13 = 3^(8+4+1) = 3^8 × 3^4 × 3^1

这样,我们只需要计算 3^1, 3^2, 3^4, 3^8 这几个值,然后将指数二进制表示中对应位为 1 的项相乘即可。

三、快速幂的递归实现

递归实现快速幂更加直观:

long long quickPower(long long a, long long n) {if (n == 0) return 1;if (n % 2 == 1) return a * quickPower(a, n - 1);else {long long temp = quickPower(a, n / 2);return temp * temp;}
}

递归的思路是:

  • 如果 n 为 0,返回 1
  • 如果 n 为奇数,分解为 a × a^(n-1)
  • 如果 n 为偶数,分解为 (a^(n/2))^2

四、快速幂的迭代实现

迭代实现更加高效,避免了递归带来的函数调用开销:

long long quickPower(long long a, long long n) {long long result = 1;while (n > 0) {if (n & 1) result *= a;  // 如果当前位为1,累乘到结果a *= a;  // 底数平方n >>= 1;  // 指数右移一位}return result;
}

迭代的核心逻辑是:

  1. 初始化结果为 1
  2. 循环处理指数的每一位
  3. 如果当前位为 1,将当前底数乘入结果
  4. 底数平方,指数右移

五、模运算下的快速幂

在实际应用中,幂运算的结果往往非常大,需要对结果取模:

long long quickPower(long long a, long long n, long long mod) {long long result = 1;a %= mod;  // 防止初始值过大while (n > 0) {if (n & 1) result = (result * a) % mod;a = (a * a) % mod;n >>= 1;}return result;
}

六、快速幂的应用场景

  1. 密码学:RSA 算法中大量使用模幂运算
  2. 数论问题:如计算大指数的余数
  3. 动态规划:状态转移方程中可能涉及幂运算
  4. 矩阵快速幂:计算递推数列的高效方法

七、总结

快速幂算法通过利用指数的二进制分解,将幂运算的时间复杂度从 O (n) 优化到 O (log n),是一种非常高效的算法。迭代实现避免了递归调用的开销,是实际应用中的首选。在处理大数问题时,模运算下的快速幂尤为重要。

http://www.xdnf.cn/news/14248.html

相关文章:

  • Python脚本开发入门:从基础到进阶技巧
  • SpringBoot ​@ControllerAdvice 处理异常
  • 鸿蒙app 开发中 如何 看 app 页面的ui结构
  • JS 数组转Object和Map
  • PHP基础-运算符
  • 【62 Pandas+Pyecharts | 智联招聘大数据岗位数据分析可视化】
  • 如何VMware虚拟机扩容磁盘,有很详细图文
  • Blazor Web Assembly - 使用Power Automate Desktop来跟踪一下Blazor页面的内存使用情况
  • 动态规划:求最长回文子串
  • OpenMMlab导出MaskFormer/Mask2Former实例分割模型并用onnxruntime和tensorrt推理
  • DB2连接池监控与挂起连接释放指南
  • Win32OpenSSL工具下载地址
  • Electron截取响应体
  • @Validation 的自定义校验实现, Spring Boot 和 java
  • 实现网页中嵌入B站视频播放器:解决high_quality=1 失效的问题
  • struct stat结构体
  • NY230NY233美光固态闪存NY237NY246
  • 【Transformer拆解】-2. 位置编码(Positional Encoding)
  • 一个密码实现库crypto-work
  • Pandas数据工程深度解析
  • 四数之和-力扣
  • XSS (Reflected)-反射型XSS
  • 晶振常见封装工艺及其特点
  • 深入讲解 Ollama 的源码
  • 【Java多线程从青铜到王者】定时器的原理和实现(十一)
  • Spring依赖注入源码学习:基于XML配置的DI源码解析
  • PGCP:用于比较基因组学的植物基因组综合数据库-文献精读144
  • 信息学奥赛一本通 1543:【例 3】与众不同
  • ubuntu之坑(十四)——安装FFmpeg进行本地视频推流(在海思平台上运行)
  • UVM同步的方法