ACwing算法备战蓝桥杯——Day23——快速幂与逆元

模板：

思路：

对于暴力法，要进行k次res*a%p这一操作。但是k可以通过二进制法凑出来，对于int类型的k最多进行31次操作就可以凑出。

并且k是a的阶数，每次k的变化可以转换成a进行平方；(通过二进制凑数的方法很常见)

主要内容	参数	作用
函数pmi（）	int a , int k , int p	在将求（a^k mod p）的时间复杂度从O(k)降到O(log k);

求 m^k mod p，时间复杂度 O(logk)。

int qmi(int m, int k, int p)
{
    int res = 1 % p, t = m;
    while (k)
    {
        if (k&1) res = res * t % p;
        t = t * t % p;//每次都会乘，所以也叫反复平方法
        k >>= 1;
    }
    return res;
}

要注意的是，数据范围是否会爆int.

模板题1：

给定 n 组 ai,bi,pi，对于每组数据，求出 abiimodpi 的值。

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

接下来 n 行，每行包含三个整数 ai,bi,pi。

输出格式
对于每组数据，输出一个结果，表示 abiimodpi 的值。

每个结果占一行。

数据范围
1≤n≤100000,
1≤ai,bi,pi≤2×109
输入样例：
2
3 2 5
4 3 9
输出样例：
4
1

代码：
#include <iostream>
#include <cmath>
#include <vector>
using namespace std;

typedef long long LL;

const int N = 1e5+10;

//反复平方
LL qmi(int a,int b,int p){
    
    LL res=1;
    
    while(b){
        
        if(b&1) res=res*a%p;
        
        a=a*(LL)a%p;//每次都乘,相当于预处理
        
        b>>=1;
        
    }
    
    return res;
}

int main(){
    int n;
    cin>>n;
    
    while(n--){
        int a,b,p;
        scanf("%d%d%d",&a,&b,&p);
        
        printf("%lld\n",qmi(a,b,p));
    }
    
    return 0;
}

模板题2：

对任意正整数  N ，计算 XNmod233333 的值。

输入格式
共一行，两个整数 X 和 N。

输出格式
共一行，一个整数，表示 XNmod233333 的值。

数据范围
1≤X,N≤109
输入样例：
2 5
输出样例：
32

代码：
#include <iostream>

using namespace std;

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

int main(){
    int x,n;
    
    cin>>x>>n;
    
    cout<<qmi(x,n,233333)<<endl;
    
    return 0;
}

快速幂求逆元：

逆元的作用：

可以将除法问题转换成乘法问题。

求a/b余数的问题转换成求a*(b的逆元的问题），mod上一个m

思路：

1.逆元的性质：b的逆元*b 模上 m =1

2.费马定理：对于两个互质的数b,m，b的m-1次方模上m等于1

根据1,2两条信息得,b的逆元就是b的m-2次方（m是质数，m-2大于等于0）

所以求用快速幂求b的m-2次方mod m即可

例题：

给定 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 为质数时，b的m−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 <iostream>
#include <algorithm>

using namespace std;

typedef long long LL;


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

//非法情况：因为规定了p是质数，所以如果a与p之间存在倍数关系，只能是p|a，因为p是质数，没有别的约数了。在这里，不存在逆元就只有a是p的倍数这种情况

int main()
{
    int n;
    scanf("%d", &n);
    while (n -- )
    {
        int a, p;
        scanf("%d%d", &a, &p);
        
        if (a % p == 0) puts("impossible");//当p|a时，两者不存在乘法逆元
        
        else printf("%lld\n", qmi(a, p - 2, p));//  
    }

    return 0;
}