快速幂讲解

最新推荐文章于 2023-10-09 20:43:40 发布

unknown_coder

最新推荐文章于 2023-10-09 20:43:40 发布

阅读量1k

点赞数 3

分类专栏：快速幂文章标签：快速幂

本文链接：https://blog.csdn.net/ivy_uu/article/details/52219638

版权

快速幂专栏收录该内容

2 篇文章 0 订阅

订阅专栏

快速幂取模

用法：用于求解 a 的 b 次方，而b是一个非常大的数，用O（n）的复杂度会超时。那么就需要这个算法，注意它不但可以对数求次幂，而且可用于矩阵快速幂。

假如求 x ^ n 次方

我们可以把 n 表示为 2^k1 + 2k2 + 2^k3….，可以证明所有数都可以用前式来表示。（其实就是二进制表示数的原理）

那么 x^n = x^2^k1 * x^2^k2 * x^2^k3……

那么就可以利用二进制来加快计算速度了。

假如 x^23 , 23转化为二进制为 10111，即 x^23 = (x^16) * ( x^4) * ( x^2) * (x^1);
通过判断二进制x是否为1来判断是否加入运算中,(x>>1)每次向右移动一位，(x&1)判断是否为1，如果是1，则加入运算。如果不是，则x的幂次方继续乘方(1,2,4,8,16,32增长，作为待用乘法因子)。

那么是不是可以在O（logn）的复杂度求解。

代码：

#include<iostream>
#include<cstdio>
using namespace std;
int pw(int a,int n)
{
    int tmp=a,ans=1;
    while(n)//这里用a的1次方来试，就很容易写出
    {
        if(n&1)
        {
            ans*=tmp;
        }
        tmp*=tmp;
        n>>=1;
    }   
    return ans;
}
int main()
{
    int a,n; 
    cin>>a>>n;
    cout<<pw(a,n)<<endl;
}