【技巧】简单理解快速幂(求模)

2023-12-15 22:38

本文主要是介绍【技巧】简单理解快速幂(求模),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

【技巧】简单理解快速幂(求模)



      今天讲了一个特别有用的东西就是快速幂,为了弄懂这个百度了一上午。。。还是没咋明白。。。


      怕忘了就先开个博文记一下代码什么的。。。


      快速幂,就是更快速地计算一个数的次方的方法。传统方法求幂计算,数小了还好,数大了就容易超时。

      这个方法据说大部分比赛都不会超时,灰常地腻害呢~

      具体理论总是太高大上了,还是举栗子好吃,简单又粗暴!

      比如我们来算3的10次幂,把3乘10次脑袋就炸了,怎么算呢,这么算!


                              3*3*3*3*3*3*3*3*3*3…………………………(10个3相乘)

                            =(3*3)*(3*3)*(3*3)*(3*3)*(3*3)…………………②

                            =(3*3)^5

                            =((3*3)*(3*3))^2*(3*3)………………………………………③


      …啥?并没有感觉多好算?废话!人脑算起来当然难算,我们叫电脑来算啊~

      先来看①式,如果要电脑来算,10个3相乘,就要乘9次;对于②式,五个3×3(把3×3看成一个整体)相乘,就要乘4次;对于③,就是两个((3*3)*(3*3))相乘再多乘一个多出来的(3*3),只要乘3次就好了,对于电脑来讲,工作的循环次数越少就越省时间(这个时候我还没学时间复杂度呢,我就先这样理解了╮(╯_╰)╭)。

      这样就总结出来一个公式!


                                                    n^p   (p为偶数时)                                               n^p(p为奇数时)

                                                =(n^2)^(p/2)                                                        =((n^2)^(p/2))*n

                                                =((n^2)^2)^(p/2/2)                                =(((n^2)^2)^(p/2/2))*n

                                                .                                                                                                    .

                                                .                                                                                                    .

                                                .                                                                                                    .

                                                =(n^p)*(1)                                                                =这个没固定公式(因为p每次除2之后奇偶性不固定)

这个公式前提是不管p除多少个2商都是偶数


      而现实中情况更接近p为奇数的那种情况,,对于那种情况,变换也很简单,当p为奇数时,就把前面括号里的一堆东西(记为x)平方掉再乘以(p/2-0.5)(就是去尾法),注意还没完!还要再把去尾丢掉的一个x再乘上,就变成了  原式=(x^2)*(p/2-0.5)*x  当然程序里面如果p是整型变量就不用减去0.5了。


      综合上面的东西,可以得出快速计算a的p次幂的函数代码:


long long QuickPow(long long a,long long p)
{long long ans=1;while(p){if(p%2==1)//当p时奇数时,相当于往后面把那个少乘的x补乘上去{ans=ans*a;}p/=2;a*=a;}return ans;
}

          有时候会叫你求a的p次幂除以mod(mod只是一个数)的余数,这时候就要用到同余定理了,同余定理式子是这样的:

      (a*b)除以c的余数=(a除以c所得的余数)×(b除以c所得的余数),即(a*b)%c=(a%c)*(b%c),%是取余符号。


      这样就会得到一个引理:



      SO!

      我们的代码就可以这样写了:


__int64 quickpow(__int64 a,__int64 p,__int64 mod)
{__int64 ans=1;a=a%mod;while(p){if(p%2==1){ans=ans*a%mod;}p/=2;a=a*a%mod;}return ans%mod;
}
 

                  只是多往后面对mod取了个余罢了~

           而我们上面的代码,就是当mod=1时的情况。


                                              任务完成!



华丽分割-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=--=-=-=-=-=-=-=-=割分丽华


另附上学长给的代码:


//快速幂求模
#include<cstdio>
int quickpow(int n,int m,int mod)
{int ans=1,base=n;while(m){if(m&1){ans=(base*ans)%mod;}base=(base*base)%mod;m>>=1;printf("ans=%d base=%d m=%d\n",ans,base,m);}return ans;
}int main()
{int n,m,mod;while(~scanf("%d%d%d",&n,&m,&mod)){printf("%d\n",mod);printf("%d\n",quickpow(n,m,mod));}return 0;
}


       哎?有点不一样!?可以看到他给的函数中判断p(学长的函数中p是m)是不是奇数用了“m&1”,这个是按位与的意思,简单来说就是先把m转换为2进制,然后取最右边那一位,如果是1就说明m是奇数,是0说明m是偶数;还有p/2变成了m>>=1,这就是把二进制的m向右移了一位,把最右边的那一位给挤掉了,现在少了最右边的那一位,其实本质上还是把m给除了个2。。。

唉,果然只有那些写出来让人看不懂的代码才能达到装逼的效果。。。→_→

这篇关于【技巧】简单理解快速幂(求模)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/498156

相关文章

Ilya-AI分享的他在OpenAI学习到的15个提示工程技巧

Ilya(不是本人,claude AI)在社交媒体上分享了他在OpenAI学习到的15个Prompt撰写技巧。 以下是详细的内容: 提示精确化:在编写提示时,力求表达清晰准确。清楚地阐述任务需求和概念定义至关重要。例:不用"分析文本",而用"判断这段话的情感倾向:积极、消极还是中性"。 快速迭代:善于快速连续调整提示。熟练的提示工程师能够灵活地进行多轮优化。例:从"总结文章"到"用

认识、理解、分类——acm之搜索

普通搜索方法有两种:1、广度优先搜索;2、深度优先搜索; 更多搜索方法: 3、双向广度优先搜索; 4、启发式搜索(包括A*算法等); 搜索通常会用到的知识点:状态压缩(位压缩,利用hash思想压缩)。

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

hdu2289(简单二分)

虽说是简单二分,但是我还是wa死了  题意:已知圆台的体积,求高度 首先要知道圆台体积怎么求:设上下底的半径分别为r1,r2,高为h,V = PI*(r1*r1+r1*r2+r2*r2)*h/3 然后以h进行二分 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#includ

电脑桌面文件删除了怎么找回来?别急,快速恢复攻略在此

在日常使用电脑的过程中,我们经常会遇到这样的情况:一不小心,桌面上的某个重要文件被删除了。这时,大多数人可能会感到惊慌失措,不知所措。 其实,不必过于担心,因为有很多方法可以帮助我们找回被删除的桌面文件。下面,就让我们一起来了解一下这些恢复桌面文件的方法吧。 一、使用撤销操作 如果我们刚刚删除了桌面上的文件,并且还没有进行其他操作,那么可以尝试使用撤销操作来恢复文件。在键盘上同时按下“C

usaco 1.3 Prime Cryptarithm(简单哈希表暴搜剪枝)

思路: 1. 用一个 hash[ ] 数组存放输入的数字,令 hash[ tmp ]=1 。 2. 一个自定义函数 check( ) ,检查各位是否为输入的数字。 3. 暴搜。第一行数从 100到999,第二行数从 10到99。 4. 剪枝。 代码: /*ID: who jayLANG: C++TASK: crypt1*/#include<stdio.h>bool h

购买磨轮平衡机时应该注意什么问题和技巧

在购买磨轮平衡机时,您应该注意以下几个关键点: 平衡精度 平衡精度是衡量平衡机性能的核心指标,直接影响到不平衡量的检测与校准的准确性,从而决定磨轮的振动和噪声水平。高精度的平衡机能显著减少振动和噪声,提高磨削加工的精度。 转速范围 宽广的转速范围意味着平衡机能够处理更多种类的磨轮,适应不同的工作条件和规格要求。 振动监测能力 振动监测能力是评估平衡机性能的重要因素。通过传感器实时监

uva 10387 Billiard(简单几何)

题意是一个球从矩形的中点出发,告诉你小球与矩形两条边的碰撞次数与小球回到原点的时间,求小球出发时的角度和小球的速度。 简单的几何问题,小球每与竖边碰撞一次,向右扩展一个相同的矩形;每与横边碰撞一次,向上扩展一个相同的矩形。 可以发现,扩展矩形的路径和在当前矩形中的每一段路径相同,当小球回到出发点时,一条直线的路径刚好经过最后一个扩展矩形的中心点。 最后扩展的路径和横边竖边恰好组成一个直

【生成模型系列(初级)】嵌入(Embedding)方程——自然语言处理的数学灵魂【通俗理解】

【通俗理解】嵌入(Embedding)方程——自然语言处理的数学灵魂 关键词提炼 #嵌入方程 #自然语言处理 #词向量 #机器学习 #神经网络 #向量空间模型 #Siri #Google翻译 #AlexNet 第一节:嵌入方程的类比与核心概念【尽可能通俗】 嵌入方程可以被看作是自然语言处理中的“翻译机”,它将文本中的单词或短语转换成计算机能够理解的数学形式,即向量。 正如翻译机将一种语言

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu