快速幂、快速幂取模的分析与代码实现

2023-05-16

写在前面

在网上搜了相关内容,感觉写的都不是特别详细,也没有人讲,只能自己理解了。下面会写一下这3个算法的分析与实现。当然都是基于自己的理解。因为博主搜了很多博客都是没有详细的解释,数学渣一脸懵逼啊。所以关于原理的解释如果有错误请一定要评论我改正!

朴素的求幂算法

也就是平常使用pow函数,最简单的实现就是一直累乘,可以得到这样的代码:

int Pow(int a,int b){
    int ans = 1;
    for(int i = 0;i < b;i++){
        ans *= a;
    }
    return ans;
}

可以看到,算法的时间复杂度是O(n)。为了降低时间复杂度,我们可以使用快速幂算法,将时间复杂度降低到O(logn),n是幂。

快速幂

关于快速幂,博主的理解是使用位运算。下面是数学证明:

关于a^b,举一个实际的例子——2^10。

那么对于6而言,如果我们将10变成二进制,那么就是:1010,如果变成加权的情况可以得到表达式:

0*2^0 + 1*2^1 + 0*2^2 + 1*2^3

代入原来的2^10可以得到表达式:

2^(0*2^0 + 1*2^1 + 0*2^2 + 1*2^3)

然后再拆开并化简这个表达式,可以得到:

2^(2^1) * 2^(2^3)

也就是说,我们在求解2^10的时候,可以考虑成根据二进制的权值来求解的。那么在关于位运算的部分,我们可以逐位获取b的位,碰到0,就累乘,碰到1,就将累乘的值乘到答案。由此可以得到代码:

int Pow(int a,int b){
    int ans = 1;
    int base = a;
    while(b){
        if(b & 1) ans *= base;
        base *= base;
        b >>= 1;
    }
    return ans;
}

 

关于位运算的部分,这里只简单提一下,更多详细内容请百度。if语句中的b & 1,也就是将b与1按位与,那么1的二进制是1,也就是说b除了最后一位,其他位都和0相与,那么得到的值一定都是0。而最后一位与上1,那么原来位实际上是不变的。这样我们就获得了b的最后一位的值。while循环中b >>= 1,是移位运算,将b向右移动1位。因为每次我们都是要获取最后一位,看是0还是1,那么每次移动1位,下一次在碰见if语句,就得到原来最后一位的前一位了。

核心代码是下面这一部分:

if(b & 1) ans *= base;
base *= base;
b >>= 1;

博主理解这里理解半天Orz,这里其实是对上面一句话的模拟:“2^(2^1) * 2^(2^3)”。实际操作一下吧:

可以看到,出现0的时候,我们就累乘,得到的是幂的累积值。出现1的时候,我们将累积的幂值进行累乘。或者说碰见一个1,就出现一次底数2。

如果仍然没有太理解,debug一下会有很多帮助。

快速幂取模

这一部分就跟数论关系很大了。取模也是数论问题中经常出现的。那么对于幂来取模,如果我们直接用模运算实际上是速度很慢的(因为试除法)。所以我们不妨在求快速幂的时候添加一些内容,从而得到结果。这个算法需要了解一下数论的一个定理:

(a*b) mod c = ((a mod c)*(b mod c)) mod c

那么根据上面的定理可以推导出另一个定理:

(a^b) mod c = (a mod c)^b mod c

具体的证明这里不再赘述,主要是看第二个定理,恰好符合我们的小标题——快速幂取模。我们可以在求快速幂的时候,通过对底数取模的方式,不断缩小底数的规模。那么我们在上面快速幂的基础上,添加取模,就可以完成整个操作。

int pow_mod(int a,int b,int c){
    int ans = 1;
    int base = a%c;
    while(b){
        if(b & 1) ans = (ans*base)%c;
        base = (base*base)%c;
        b >>= 1;
    }
    return ans;
}

如果能理解上面的快速幂算法,那么这个也会比较好理解了。定理里面,底数是要有一次取模运算的。这里我们在给base赋值的时候就运行了一次。那么对于后面的一次取模,我们实际上利用了分配率,即:

(a*b) mod c = ((a mod c)*(b mod c)) mod c

我们求幂的本质仍然是求积。所以每次我们对base或者ans进行运算的时候,都必须使用一次分配率,所以都要mod c。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

快速幂、快速幂取模的分析与代码实现 的相关文章

  • 【Linux】Ubuntu 代理配置

    apt get 设置代理 proxy 方法 方法一 xff1a 这是一种临时的手段 xff0c 如果你仅仅是暂时需要通过http代理使用apt get xff0c 你可以使用这种方法 在使用 apt get 之前 xff0c 在终端中输入以
  • 百度之星之E:C++ 与Java

    E C 43 43 与Java 时间限制 2000ms 内存限制 65536kB 描述 在百度之星的贴吧里面 xff0c Java的爱好者和C 43 43 的爱好者总是能为这两种语言哪个更好争论上几个小时 Java的爱好者会说他们的程序更加
  • 并查集详解

    并查集是我暑假从高手那里学到的一招 xff0c 觉得真是太精妙的设计了 以前我无法解决的一类问题竟然可以用如此简单高效的方法搞定 不分享出来真是对不起party了 xff08 party xff1a 我靠 xff0c 关我嘛事啊 xff1f
  • ubuntu18.04 开启ssh远程服务

    1 查看ssh服务是否已经开启 说明 xff1a 1 ssh agent 指的是ubuntu的ssh服务的客户端 xff0c 用于该ubuntu远程连接其它Linux主机 如果没有ssh agent的话 xff0c 该ubuntu主机也无法
  • Python必备知识之“if __name__ == ‘__main__‘:”

    在学习Python的过程中经常会看到 if name 61 61 39 main 39 这行代码 xff0c 那么这行代码的作用究竟是什么呢 xff1f if name 61 61 39 main 39 这行代码的主要作用是调试某个模块的正
  • Windows Server 网络连接由公用网络改为专用网络

    主题 xff1a Windows Server 网络连接由公用网络改为专用网络 关键字 xff1a 问题描述 xff1a Windows Server 2012 r2 启动后网络连接被识别为公用网络 xff0c 导致远程桌面等服务无法使用
  • 关于书籍(WPF及其它)

    原文 xff1a On Books WPF and Otherwise 有人让我去看coding horror comparison xff0c 这篇文章来至于Charles Petzold和Adam Nathan的书籍 xff0c 是关于
  • pip,pip安装源

    介绍 Python在使用pip安装第三方包 第三方功能库的时候 xff0c pip3 pip install xxx走的是国外源 xff0c 有点慢 我们可以采用国内源加快下载的速度 常用pip源 xff1b 豆瓣 xff1a https
  • 安装Anaconda时安装路径错误,提示Directory" xxx is not empty ,please choose a different location."问题的解决方案

    错误如下图所示 重新选择路径 xff0c 选择平时安装的盘 xff0c 然后手动输入Anaconda xff0c 即可正常安装 xff08 在这一步之前一定要删除卸载 先前安装产生的文件夹 xff09 进QQ群 xff08 77980901
  • vue项目引入PWA(vue-cli4)

    1 概念 PWA 全称为 Progressive Web App xff0c 中文译为渐进式 Web APP 其目的是通过各种 Web 技术实现与原生 App 相近的用户体验 也就是说 xff0c 只要你使用浏览器 xff0c 就可以实现免
  • Linux远程管理协议(RFB、RDP、Telnet和SSH)

    提到远程管理 xff0c 通常指的是远程管理服务器 xff0c 而非个人计算机 个人计算机可以随时拿来用 xff0c 服务器通常放置在机房中 xff0c 用户无法直接接触到服务器硬件 xff0c 只能采用远程管理的方式 远程管理 xff0c
  • Python第三方库(模块)下载和安装(使用pip命令)

    Python第三方库是由社区开发者编写的代码包 xff0c 用于增强Python的功能和提供各种特定的功能 通常 xff0c 这些库被打包为模块 xff0c 可以通过使用Python包管理工具pip来下载和安装 以下是使用pip下载和安装P
  • 计蒜客T1098 大整数加法

    求两个不超过 200 位的非负整数的和 输入格式 有两行 xff0c 每行是一个不超过 200 位的非负整数 xff0c 可能有多余的前导 0 输出格式 一行 xff0c 即相加后的结果 结果里不能有多余的前导 0 xff0c 即如果结果是
  • Linux系统学习(三)Linux系统管理

    用户和组管理 1 配置文件 passwd文件 位置 xff1a etc passwd xff1b 对任何用户可读 作用 xff1a 用于保存各用户的账户信息 shadow文件 位置 xff1a etc shadow xff1b 只对root
  • HTTP Host 头攻击 -- 学习笔记

    目录 1 HTTP Host头攻击 2 HTTP Host头的作用 3 什么是HTTP Host头攻击 4 如何发掘HTTP Host头攻击 修改Host值 添加重复的Host头 使用绝对路径的URL 添加缩进或换行 注入覆盖Host头的字
  • Linux 网络流量监控工具

    Linux 网络流量监控 Linux 网络流量监控是捕获和分析企业的 Linux 网络流量的过程 为什么要监控 Linux 网络流量 深入了解网络流量对于测量和管理带宽使用情况非常重要 分析 Linux 网络流量有助于识别带宽瓶颈 最高用量
  • 【解决“Authentication is required to create a color profile/managed device“】

    解决Ubantu系统 34 Authentication is required to create a color profile managed device 34 问题 xff1a 在Windows下使用远程桌面连接到工作站的Uban
  • 漫谈微信开放平台一(小程序服务器url设置)

    点击查看文档 这里的是需要用开放平台设置特约商户的域名 两种域名 xff0c 1 服务器域名 2 业务域名 两种域名设置方案相似 xff0c 我以服务器域名设置为例 需要注意 xff1a 1设置的域名需要在开放平台进行设置 xff08 注意
  • linux / ubuntu / 添加和查看环境变量的方法

    一 添加 1 export 指令 export PATH 61 PATH home xiaoming Doc 将 home xiaoming Doc 放到了名为 PATH 的环境变量的后面 或者 export PATH 61 home xi
  • 504 Gateway Time-out原因及解决方法

    1 今天在webpack通过proxy开发的时候 xff0c 接口时正常的 xff0c 但是上到测试机 xff0c 通过nginx转发的时候 xff0c 就会出现504 Gatway time out 思路 1 xff0c 初步判断时ngi

随机推荐