java位运算求幂,程序员必学:快速幂算法

前阵子,有小伙伴在我B站的算法教程底下留言

824807511b1ebb4366a25abf56bec1fb.png

小伙伴们有任何疑问或者希望我解说任何内容,都可以在我的小我私家B站或民众号(xmg_mj)留言哦,我会尽我最大能力、只管抽时间去写文章\录视频来回应人人。

关于快速幂

实在快速幂相关的问题,是加入算法竞赛(NOI、ACM等)的小伙伴必须要掌握的一小块基础内容。固然,就算你不设计加入算法竞赛,小我私家以为只要你是一名程序员,就必须要掌握快速幂算法。

在《计算机程序设计艺术》一书中就有提到快速幂算法,此书的英文名是The Art of Computer Programming,简称TAOCP。

5d8309c3878e14df9332926ad1c05ce7.png

TAOCP出自Donald Ervin Knuth先辈之手。Knuth先辈是在计算机领域成就颇丰的着名科学家,是著名的KMP算法的发明人之一,在1974年获得“计算机领域的诺贝尔奖”:图灵奖(昔时他才36岁)。现在TAOCP已经出书了第1、2、3、4A卷,凭据设计,另有第4B、5、6、7卷未出书。第一卷首发于1968年,Knuth先辈今年是82岁高寿,听说他设计在105岁之前完成这部巨著。

78ce8be0f1a2f8d0b428a9341ac2efa5.png

关于TAOCP,微软创始人Bill Gates曾说过

If you think you’re a really good programmer… read (Knuth’s) Art of Computer Programming… You should definitely send me a resume if you can read the whole thing.

也许意思是:若是你以为自己是一位异常优异的程序员,那就应该阅读Knuth的TAOCP;若是你能读懂全部内容,可以直接给我发送一份简历。听说Knuth先辈的言辞加倍犀利:看不懂就别当程序员了!不外TAOCP对于新手来说,阅读难度简直对照大,书中的所有示例都使用了Knuth先辈自创的MIX汇编语言。

阅读本文之前的提醒

今天就抽闲写一篇文章来解说一下经典的快速幂算法哈。不外要想彻底看懂本文,有几个前提条件

熟悉算法中的2个基础观点:时间复杂度、空间复杂度

若是你压根没听过这2个观点,说明你的算法基础完全为0,真的没有在开顽笑!

可以向民众号发送复杂度获取相关教程

熟悉二进制和十进制的转换

若是连这个都不熟悉的话,那你的编程基础就真的需要好好补补啦

可以向民众号发送进制获取相关教程

熟悉常见的位运算操作

n & 1的效果是n最低二进制位的值,也可以用于判断n的奇偶性

求正整数n / 2,可以用位运算取代:n >> 1

若是不明白上述操作的原理,可以向民众号发送位运算获取相关教程

什么是幂(Power)?

众所周知,x的n次幂,是指x的n次方,也就是n个x相乘,好比2的4次幂就是2 * 2 * 2 * 2。

为了简化形貌,后面x的n次幂,我就简化为x ^ n(本文中的 ^ 并不是按位异或的意思)

那若何通过编程求幂?假设只思量x、n是整数且n大于即是0的情形,最容易想到的方式如下所示(这里接纳的编程语言是Java,但没有涉及Java特殊的语法。以是就算你没用过Java,也可以看懂)

int power(int x, int n) {

int result = 1;

while (n-- > 0) {

result *= x;

}

return result;

}

很显然,这种方式的时间复杂度是O(n)、空间复杂度是O(1)

什么是快速幂?

所谓快速幂,就是用效率更高(时间复杂度更低)的方式求幂,可以将时间复杂度优化至O(logn)。这里先容2种求解方式:递归、非递归

递归

e4d9f98f403df62830795b82716bfcf8.png

凭据上图中的等式,不难写出以下代码

int fastPower(int x, int n) {

if (n == 0) return 1;

int result = fastPower(x, n >> 1);

result *= result;

return (n & 1) == 0 ? result : result * x;

}

这个方式的时间、空间复杂度都是O(logn)。

那若何剖析出这个方式的复杂度呢?

若是你的算法功底对照微弱,可以代入特定值作一个也许的剖析,好比当n为16时,方式的递归挪用历程如下图所示

80377a9a53b915bb129d43c15d565b6b.png

不难看出,每次挪用时,n的规模都减半,以是时间和空间复杂度都是O(logn)

若是你的算法功底还行,那就可以用更专业的方式去剖析它的复杂度(没有一定的算法基础,可能会看不懂)

这实在是典型的应用分治计谋的算法

假设T(n)是数据规模为n时的时间复杂度,不难得出递推式:T(n) = T(n / 2) + O(1)

最后凭据递推式 + 主定理(Master Theorem)可以直接得出结论:T(n) = O(logn)

非递归

我们以求3 ^ 21为例子,来剖析一下非递归的代码应该怎么写。

首先21的二进制形式是10101

93c44fc713826fcd00633d09294261bd.png

cbc415ccdac64dda3ba20576eba91a5d.png

不难得出以下结论

3 ^ n(n为2、4、8、16)都可以由3 ^ 1累乘出来

每一个3 ^ n都有对应的二进制位

3 ^ 1对应二进制位的值是1,实在是二进制10101的最后1位

3 ^ 2对应二进制位的值是0,实在是二进制1010的最后1位

3 ^ 4对应二进制位的值是1,实在是二进制101的最后1位

3 ^ 8对应二进制位的值是0,实在是二进制10的最后1位

3 ^ 16对应二进制位的值是1,实在是二进制1的最后1位

若是3 ^ n对应二进制位的值是0,就不用乘进最终效果

好比3 ^ (8 * 0)、3 ^ (2 * 0)

由于它们最终的值都是3 ^ 0,也就是1

若是3 ^ n对应二进制位的值是1,就需要乘进最终效果

好比3 ^ (16 * 1)、3 ^ (4 * 1)、3 ^ (1 * 1)

以是,综合以上种种结论,可以总结出以下解题步骤

行使3 ^ 1,不停累乘出3 ^ n(n为2、4、8、16)

每当累乘出一个3 ^ n,就查看其对应二进制位的值是1照样0,来决议是否要将它乘进最终效果

int fastPower(int x, int n) {

int result = 1;

while (n != 0) {

if ((n & 1) == 1) {

result *= x;

}

x *= x;

n >>= 1;

}

return result;

}

代入3和21,fastPower(3, 21)的执行流程如下

第1轮while循环

第4行代码

n的二进制是10101(十进制是21)

x = 3 ^ 1, 其对应二进制位的值是1(n的最后一个二进制位)

以是需要执行第5行代码:将x乘进最终效果

result = 3 ^ 1

第7行代码

x = (3 ^ 1) * (3 ^ 1) = 3 ^ 2

第8行代码

n右移1位,其二进制变成了1010(对应的十进制是啥?不重要!!!)

第2轮while循环

第4行代码

n的二进制是1010

x = 3 ^ 2, 其对应二进制位的值是0(n的最后一个二进制位)

以是不需要执行第5行代码:不需要将x乘进最终效果

result = 3 ^ 1

第7行代码

x = (3 ^ 2) * (3 ^ 2) = 3 ^ 4

第8行代码

n右移1位,其二进制变成了101(对应的十进制是啥?不重要!!!)

第3轮while循环

第4行代码

n的二进制是101

x = 3 ^ 4, 其对应二进制位的值是1(n的最后一个二进制位)

以是需要执行第5行代码:将x乘进最终效果

result = (3 ^ 1) * (3 ^ 4)

第7行代码

x = (3 ^ 4) * (3 ^ 4) = (3 ^ 8)

第8行代码

n右移1位,其二进制变成了10(对应的十进制是啥?不重要!!!)

第4轮while循环

第4行代码

n的二进制是10

x = 3 ^ 8, 其对应二进制位的值是0(n的最后一个二进制位)

以是不需要执行第5行代码:不需要将x乘进最终效果

result = (3 ^ 1) * (3 ^ 4)

第7行代码

x = (3 ^ 8) * (3 ^ 8) = 3 ^ 16

第8行代码

n右移1位,其二进制变成了1(对应的十进制是啥?不重要!!!)

第5轮while循环

第4行代码

n的二进制是1

x = 3 ^ 16, 其对应二进制位的值是1(n的最后一个二进制位)

以是需要执行第5行代码:将x乘进最终效果

result = (3 ^ 1) * (3 ^ 4) * (3 ^ 16)

第7行代码

x = (3 ^ 16) * (3 ^ 16) = 3 ^ 32

第8行代码

n右移1位,其二进制变成了0

最后

由于n = 0,以是退出while循环

最终result = (3 ^ 1) * (3 ^ 4) * (3 ^ 16)

复杂度剖析

每执行一次while的循环体,n >>= 1, 会导致n的值减半

以是时间复杂度:O(logn)、空间复杂度:O(1)

Leetcode

Leetcode上的第50号题50. Pow(x, n),恰好就可以用今天解说的快速幂算法。以下是我的代码实现

// 递归

public double myPow(double x, int n) {

if (n == 0) return 1;

if (n == -1) return 1 / x;

double half = myPow(x, n >> 1);

half *= half;

return ((n & 1) == 1) ? half * x : half;

}

// 非递归

public double myPow(double x, int n) {

long y = (n < 0) ? -((long) n) : n;

double result = 1.0;

while (y > 0) {

if ((y & 1) == 1) {

result *= x;

}

x *= x;

y >>= 1;

}

return (n < 0) ? (1 / result) : result;

}

需要提醒的是

这里我用的编程语言是Java,人人可以凭据自己熟悉的编程语言,对一些语法细节作出响应的调整

Leetcode上的n可能是个负数,以是上面的代码针对负数的情形作了一些处置

更多快速幂相关的问题

时间有限,这篇文章就先说到这了哈。给小伙伴们留2个快速幂相关的问题,有空的话,可以去研究一下

使用矩阵快速幂求斐波那契数列

请设计一个算法求x的y次幂模z的效果:(x ^ y) % z

假设x、y都可能是很大的整数(y大于即是0,z不即是0)

若是你稀奇希望我写点什么方面的内容,也可以留言建议,谢谢。迎接关注

原文链接:https://www.cnblogs.com/mjios/p/12690097.html

本站声明:网站内容来源于网络,若有侵权,请联系我们,我们将及时处置。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/572738.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

java过滤器

过滤器 1、Filter工作原理&#xff08;执行流程&#xff09; 当客户端发出Web资源的请求时&#xff0c;Web服务器根据应用程序配置文件设置的过滤规则进行检查&#xff0c;若客户请求满足过滤规则&#xff0c;则对客户请求&#xff0f;响应进行拦截&#xff0c;对请求头…

假设检验

假设检验分参数假设和非参数假设。 假设 先假设原假设H0&#xff0c;对应的反面叫做备择假设H1。SAS一般沿用的规则是NEYMAN和PEARSON提出的&#xff1a;在控制犯第一类错误的原则下&#xff0c;是犯第二类错误的概率尽量小&#xff08;即&#xff0c;原假设受到保护&#xff0…

jeesite在eclipse中部署

1&#xff1a;下载下来最新版本的jeesite&#xff0c;首先要在本地安装好maven运行环境 2&#xff1a;运行 bin/eclipse.bat 生成工程文件并下载jar依赖包 如果需要修改默认项目名&#xff0c;请打开pom.xml修改第7行artifactId&#xff0c;然后再执行eclipse.bat文件 3&…

php抖音关注列表,网页PHP抖音批量取消关注JS代码

在很多时候我们抖音关注的人太多&#xff0c;想要批量取消但是一个一个点击太麻烦了&#xff0c;如何解放双手批量取消关注呢&#xff1f;今天分享一段JS代码&#xff0c;可在线批量取消关注&#xff0c;需电脑才能操作。首先打开抖音创作服务平台登录&#xff1a;https://crea…

python 灰度改二值,Python Image 对验证码进行灰度和二值法处理 去掉边框

def binarizing(img, threashold):img img.convert("L") # 转灰度pixdata img.load()w, h img.sizefor y in range(h):for x in range(w):if pixdata[x, y] < threashold:pixdata[x, y] 0else:pixdata[x, y] 255return imgdef removeFrame(img, width)::para…

前端js模糊搜索(模糊查询)

1.html结构&#xff1a; <label for"searchShop" class"clear pos-a" style"top:17px;"> <input type"text" id"searchShop" placeholder"场所搜索"> <input type"but…

php代码审计思路,代码审计思路讨论

原标题&#xff1a;代码审计思路讨论前言首先感谢手电筒的建议&#xff0c;建议我从Bluecms来开始学习代码审计。感谢桑桑格&#xff0c;风哥的问题解答。文章主要是写下我自己在代码审计的时候所思所想。当然并不是太过严谨&#xff0c;也有可能存在错误。这里主要是写出了审计…

Dubbo源码分析:ThreadPool

定义了通过URL对象作为参数获取Executor对象的get&#xff25;xecutor方法。所有实现ThreadPool接口的类都是基于ThreadPoolExecuotr对象来实现的。 类图 转载于:https://www.cnblogs.com/wspgbw/p/6409933.html

MySQL存储过程及触发器

一、存储过程 存储过程的基本格式如下&#xff1a; -- 声明结束符 -- 创建存储过程 DELIMITER $ -- 声明存储过程的结束符 CREATE PROCEDURE pro_test() --存储过程名称(参数列表) BEGIN-- 可以写多个sql语句; -- sql语句流程控制SELECT * FROM employee; END $ -- 结…

destoon 短信发送函数及短信接口修改

// $DT在common.inc.php中定义&#xff0c; $CACHE cache_read(module.php); $DT $CACHE[dt]; 从缓存里读取网站配置信息。 //$db 数据库类&#xff0c;$DT 全局配置&#xff0c;$DT[sms] 开启手机短信&#xff0c;$DT_TIME系统时间&#xff0c;$DT_IP 当前ip地址&#xff0…

安卓界面组件-文本框

前面三篇文章本质上是围绕着View类进行的。View是安卓UI的基础类&#xff0c;我们的安卓开发“千里之行”是从View开始的。 安卓界面UI有大量的组件&#xff0c;组件的继承和间接继承于View。有一类组件很基本&#xff0c;这就是TextView文本框&#xff0c;作用是显示文本。在T…

redis linux 运行依赖,在Ubuntu系统下安装Redis Desktop Manager的方法

以下为你介绍在Ubuntu16.04系统下安装Redis可视化工具Redis Desktop Manager。不管是哪个Redis Desktop Manager版本&#xff0c;都可以采用下面的安装方式进行。Redis Desktop Manager安装方法1.首先在终端中启动redis服务执行命令&#xff1a;sudo service redis start 与 ps…

bzoj4631

4631: 踩气球 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 260 Solved: 133[Submit][Status][Discuss]Description 六一儿童节到了&#xff0c; SHUXK 被迫陪着M个熊孩子玩一个无聊的游戏&#xff1a;有N个盒子从左到右排成一排&#xff0c;第i个盒子里装着Ai个气球。SH…

hpcp5225设置linux网络,惠普HP color laserjet CP5225打印机驱动

软件标签&#xff1a;惠普CP5225驱动是为同名打印机推出的驱动程序&#xff0c;维持着打印机的正常工作&#xff0c;小编提供了驱动下载&#xff0c;若是你的惠普CP5225打印机在使用过程中出现不能打印或者无法被电脑识别等问题&#xff0c;快来本站下载该驱动来解决问题吧!安装…

Dubbo底层采用Socket进行通信详解

由于Dubbo底层采用Socket进行通信&#xff0c;自己对通信理理论也不是很清楚&#xff0c;所以顺便把通信的知识也学习一下。 n 通信理论 计算机与外界的信息交换称为通信。基本的通信方法有并行通信和串行通信两种。 1.一组信息&#xff08;通常是字节&#xff09;的各位数据被…

linux 网络协议栈变化,ZZ Linux网络协议栈学习

最近学习linux内核网络协议栈&#xff0c;把数据包接收流程大致理了一下&#xff0c;前面也看了瀚海书香兄的总结&#xff0c;感觉总结的比我精炼&#xff0c;抓住了主干&#xff0c;是一目了然的那种我的这篇本来是自己看得&#xff0c;因此把我自己学习中一些遇到的问题写了出…

linux安装各种文件格式,Embeded linux中的各类文件系统

Cramfs&#xff1a;(优点&#xff1a;只读&#xff0c;目的&#xff1a;不可更改文件、防入侵)1.下载cramfs-1.1源码工具http://files.cnblogs.com/files/pokerface/cramfs-1.1.tar.gz2.解压tar xvzf cramfs-1.1.tar.gz -C /home/liuzhengwu/tools3.编译make4.得到COPYING cra…

初始Angularjs2

一、为什么使用Angularjs2 Angularjs2全面使用了模块化、组件化的思想&#xff0c;它有如下的特性&#xff1a; 1&#xff09;模块化 在ng2的应用中&#xff0c;所有的系统功能都是模块化的&#xff0c;开发者只需要按需导入模块化的好处在于可以当应用加载时&#xff0c;是按需…

linux中输入ls出现蓝色的点,linux上ls的蓝色太深,怎么处理

linux下ls的蓝色太深&#xff0c;怎么办&#xff1f;在linux下使用ls命令时&#xff0c;会将目录以蓝色显示。在某些显示器上&#xff0c;黑底蓝字&#xff0c;看起来相当费力&#xff0c;以至于我有时候不得不使用dir命令&#xff0c;全部显示成白色。先看一个图&#xff1a;看…

Discuz常见小问题-如何取消登陆发帖验证码

1 正常情况下&#xff0c;用户点击登录之后&#xff0c;需要填写验证码 2 进入后台&#xff0c;点击防灌水&#xff0c;验证设置&#xff0c;然后下面的各个选项可以设置是否启用验证码。 转载于:https://www.cnblogs.com/acetaohai123/p/6504754.html