数学知识：快速幂（快速幂求逆元）

ACWing 875. 快速幂

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

输入格式
第一行包含整数n。
接下来n行，每行包含三个整数ai,bi,pi。

输出格式
对于每组数据，输出一个结果，表示ai^bi mod pi的值。
每个结果占一行。

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

Code:

#include <iostream>
#include <algorithm>
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
        k >>= 1;//k的末位删掉
        a = (LL) a * a % p;
    }
    return res;
}

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

AcWing 876. 快速幂求逆元

给定n组ai,pi，其中pi是质数,求ai模pi的乘法逆元，若逆元不存在则输出impossible。
注意：请返回在0∼p−1之间的逆元。

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

输入格式
第一行包含整数n。
接下来n行，每行包含一个数组ai,pi，数据保证pi是质数。

输出格式
输出共n行，每组数据输出一个结果，每个结果占一行。
若ai模pi的乘法逆元存在，则输出一个整数，表示逆元，否则输出impossible。

数据范围
1≤n≤105,
1≤ai,pi≤2∗10⁹
输入样例：
3
4 3
8 5
6 3
输出样例：
1
2
impossible

快速幂求逆元的思路：
逆元的性质：若b的逆元是b^-1, 则有 b * b^-1 ≡ 1(mod m);
费马小定理：如果m是一个质数，且整数a不是m的倍数，则有 a^m-1≡ 1(mod m).
由费马小定理可得 a * a^m-2 ≡ 1(mod m),所以 a 的逆元为a^m-2.（条件是a与m互质）

Code:

#include <iostream>
#include <algorithm>
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;
    scanf("%d", &n);
    while(n --)
    {
        int a, p;
        scanf("%d%d", &a, &p);
        
        int res = qmi(a, p - 2, p);
        if(a % p)  printf("%d\n", res);
        else puts("impossible");
    }
    return 0;
}