acwing基础课——快速幂

由数据范围反推算法复杂度以及算法内容 - AcWing

常用代码模板4——数学知识 - AcWing

基本思想：

     求一个数的n次时，我们的时间复杂度为O(n),当n特别大时，效率会很低可能超时，此时我们就需要运用到快速幂，将我们的时间复杂度降低到O(log n)。快速幂的核心思想是，我们预处理一些数据，然后将a^k用我们预处理的数据组合出来，即将a^k拆成若干个预处理的值，我们将k转化为二进制表示，就可以得到我们所需要的数据。

       那么 为什么快速幂的时间复杂度为O(log n)呢？因为我们首先要预处理出log k个数据，对于a的k次方，我们预处理a的2^0，a的2^1,....，a的2^log k，并且可以看出每一个数都是前一个数的平凡，我们直接累积就可以处理出来。

         然后，我们如何将a^k拆成我们预处理的值呢？我们可以将a^k拆成a的2^x1, .... ,a的2^xt的乘积，也等于a的2^x1 + 2^x2 + ... + 2^xt,这样我们就可以发现，我们将k转化成二进制，将对应1的预处理好的数据进行计算即可。这样我们只用花费O(log n)预处理数据，O(1)计算就可以得出答案。

         示例，4^5 mod 10, 我们预处理出 4 的 2^0, 4的2^1，4的2^2，然后将4^5拆解，5的二进制表示为101，即取出4 的 2^0和4的2^2相乘得到结果。

 875. 快速幂 - AcWing题库

给定 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×1e9

输入样例：

2
3 2 5
4 3 9

输出样例：

4
1

#include

using namespace std;

typedef long long LL;

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

int main()
{
    int n;
    cin >> n;
    while(n--)
    {
        int a, k, p;
        cin >> a >> k >> p;
        cout << qmi(a, k, p) << endl;
    }
    return 0;
}

876. 快速幂求逆元 - AcWing题库

         快速幂加上费马定理(假如p是质数，且gcd(a,p)=1（a和p互质），那么 a^(p-1)  ≡ 1（mod p），即 ( a^(p-1) )%p = 1)的运用，无解的情况为两数不互质，即a为p的倍数，否则我们一定可以通过费马定理构造出来一个解

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

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

输入格式

第一行包含整数 n。

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

输出格式

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

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

数据范围

1≤n≤1e5,
1≤ai,pi≤2∗1e9

输入样例：

3
4 3
8 5
6 3

输出样例：

1
2
impossible

#include
#include

using namespace std;

typedef long long LL;

int qmi(int a, int b,int p)
{
    LL res = 1;
    
    while(b)
    {
        if(b & 1) res = (LL) res * a % p;
        b >>= 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 == 0) cout << "impossible" << endl;
        else cout << qmi(a,p - 2,p) << endl;
    }
    
    return 0;
}