【4.4 快速幂详解及快速幂求逆元】

$更好的阅读体验$

4.4.1 快速幂详解

概念

• 快速求出$a^k\mathrm{m}\mathrm{o}\mathrm{d}p$的结果

思想

• 预处理出$a^{2^0},a^{2^1},a^{2^2}\dots a^{2^{lo{g}_2^k}}$的结果
• 则使得$k=2^{p_1}+2^{p_2}+\cdots +2^{p_i}$
• 即：$a^k=a^{2^{p_1}}\times a^{2^{p_2}}\times \cdots \times a^{2^{p_i}}$
• 对于$a^{2^0}\times a^{2^0}=a^{2^1},a^{2^1}\times a^{2^1}=a^{2^2}$，即 $a^{2^{p_i}}=a^{2^{p_{i-1}}}\times a^{2^{p_{i-1}}}$
• 综上所述，在操作时记录$a^{p_i}$的值，和累乘的结果
• 将$k$化为二进制表示，按位>>操作，若当前位是$1$，则对当前累乘的结果$\times a^{p_i}\mathrm{m}\mathrm{o}\mathrm{d}p$
• 每次对$k$进行>>操作后，更新$a^{p_{i+1}}=a^{p_i}\times a^{p_i}\mathrm{m}\mathrm{o}\mathrm{d}p$
• 当二进制下的$k$无法再>>时，累乘结果即为答案

模板

typedef long long LL;

LL qmi(LL a,LL k,LL p){  //计算 a^k % p 的结果
    
    LL res=1;  //记录累乘结果
    
    while(k){
        
        if(k&1) res=res*a%p;  //k&1得到当前位，若为1则累乘a^pi
        a=a*a%p;  //更新a^pi
        k>>=1;  //右移1位
    }
    
    return res;
    
}

---

4.4.2 快速幂求逆元

---

概念

• 同余：设$m$为正整数，$a$和$b$是整数，若$m|a-b$，则称$a$模$m$同余于$b$，或$a$与$b$模$m$同余，记作$a\equiv b(modm)$
• 若$ab\equiv 1(modm)$，则称$b$为$a$的模$m$逆，记作$a^{-1}(modm)$或$a^{-1}$

注意

• $a$的模$m$逆存在 $\iff$ $a$与$m$互质

• 当$m$为质数时，用费马小定理求

• 当$m$不为质数时，用扩展欧几里得算法求

思想

• 利用快速幂实现$m$为质数时用费马小定理求逆元

• 费马小定理：设$p$为素数，且$a$与$p$互质，则$a^{p^{-1}}\equiv 1(modp)$

• $$\begin{array}{rl}{a}^{p^{-1}}\equiv 1(modp)\to & a\times a^{p^{-2}}\equiv 1(modp)\\ & a\times b\equiv 1(modp)\\ & 即：b=a^{p^{-2}}\end{array}$$

---

模板例题 876. 快速幂求逆元

原题链接

描述

给定 n 组 ai,pi，其中 pi 是质数，求 ai 模 pi 的乘法逆元，若逆元不存在则输出 impossible。

注意：请返回在 0∼p−1 之间的逆元。

乘法逆元的定义
若整数 b，m 互质，并且对于任意的整数 a，如果满足 b|a，则存在一个整数 x，使得 a/b≡a×x(modm)，则称 x 为 b 的模 m 乘法逆元，记为 b−1(modm)。
b 存在乘法逆元的充要条件是 b 与模数 m 互质。当模数 m 为质数时，bm−2 即为 b 的乘法逆元。

输入格式
第一行包含整数 n。

接下来 n 行，每行包含一个数组 ai,pi，数据保证 pi 是质数。

输出格式
输出共 n 行，每组数据输出一个结果，每个结果占一行。

若 ai 模 pi 的乘法逆元存在，则输出一个整数，表示逆元，否则输出 impossible。

数据范围
1≤n≤105,
1≤ai,pi≤2∗109
输入样例：

3
4 3
8 5
6 3

输出样例：

1
2
impossible

代码

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

typedef long long LL;

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

int main(){
    
    int n;
    
    cin>>n;
    
    while(n--){
        
        int a,p;
        
        cin>>a>>p;
        
        if(a%p==0) cout<<"impossible"<<endl;  //a与p不互质则说明无逆元
        else cout<<qmi(a,p-2,p)<<endl;
        
    }
    
    return 0;
        
}

---