#快速幂 ——快速幂求逆元 ~20.8.28

最新推荐文章于 2024-03-15 12:51:25 发布

Y.YL

最新推荐文章于 2024-03-15 12:51:25 发布

阅读量268

点赞数 1

分类专栏：学习心得

。

本文链接：https://blog.csdn.net/weixin_45916959/article/details/108285595

版权

学习心得专栏收录该内容

31 篇文章 0 订阅

订阅专栏

文章目录

快速幂

快速幂

一、快速幂

AcWing 875. 快速幂
在这里插入图片描述

思路：反复平方法

将k拆分成若干个2的n次方相加的形式，并用预处理好的a^{2的0次方~2的logk次方}进行组合相乘得到a^k。
第一步：预处理出 a^2的0次方 ~2的logk次方^ mod P的结果，一共 logk 个。
预处理时，a^2的k次方 = a^{2的(k-1)次方}* a^{2的(k-1)次方} = a^{2的(k-1)次方}的平方
第二步：将a^k组合出来

int qmi(int a, int k, int p){
    int res = 1;
    while(k){
        if(k & 1) res = (ll)res * a % p;//如果k的二进制位有1那么就乘上当前的a
        k >>= 1;//k向右移一位
        a = (ll)a * a % p;//将a变成下一步的a
    }
    return res;
}

完整代码

#include<bits/stdc++.h>
using namespace std;
int n;
typedef long long ll;
int qmi(int a, int k, int p){
    int res = 1;
    while(k){
        if(k & 1) res = (ll)res * a % p;
        k >>= 1;
        a = (ll)a * a % p;
    }
    return res;
}
int main(){
    cin >> n;
    while(n --){
        int a, k, p;
        cin >> a >> k >> p;
        cout << qmi(a, k, p) << endl;
    }
    return 0;
}

二、快速幂求逆元

AcWing 876. 快速幂求逆元
在这里插入图片描述

思路

……
这里输入的a是题目定义里面的b,p是定义里的m。
所以求a的逆元的话，只需满足 a，p互质，
然后根据公式算出 a^p-2 mod p 的值就好。
.
……
总的来说，直接求 (a/b)%m可能会出错，于是想办法用一个x来代替b^-1的作用，使得(a *x)mod p同余于(a/b)%m

……辅助说明

1、转自 @无糖卡布奇诺的文章《逆元（关于除法取模）》

逆元辅助说明1
2、来自 @Dunyb 的文章《逆元》
辅助说明2

完整代码

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;

int qmi(int a, int k, int p){
    int res = 1;
    while(k){
        if(k & 1) res = (ll)res * a %  p;
        k >>= 1;
        a = (ll) a * a % p;
    }
    return res;
}

int main(){
    int n;
    cin >> n;
    while(n --){
        int a, p;
        cin >> a >> p;
        
        if(a % p) {//p已经是质数了，那么如果a还不能整除p，就说明它们互质
           int res = qmi(a, p - 2, p);
           cout << res << endl;
        }
        else cout << "impossible" << endl;
    }
    return 0;
}