标签：cdot res bmod 875 long int cdots 快速 AcWing

复杂度 $ O(log(k)) $ （k 是指数）

总体复杂度 $ log(2 \times 10^{9}) = 9 \times log(20) \approx 40 $

#include<iostream>

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;
        b >>= 1;
        a = (LL) a * a % p;
    }
    return res;
}
int main()
{
    int n;
    scanf("%d", &n);
    while (n --) {
        int a, b, p;
        scanf("%d %d %d", &a, &b, &p);
        printf("%lld\n", qmi(a, b, p));
    }
    return 0;
}

1. 快速幂的思路：
   ① $ b = c_1 \cdot2^{0} + c_2 \cdot 2^{1} + \cdots + c_k \cdot 2^{k} $
   $ a^{b} = a^{ c_1 \cdot2^{0} + c_2 \cdot 2^{1} + \cdots + c_k \cdot 2^{k} } = a^{ c_1 \cdot2^{0} } \cdot a^{ c_2 \cdot 2^{1} } \cdots a^{ c_k \cdot 2^{k} } $
   ② 由 ① ，只需找到 $ c_i = 1 $ 对应的乘积项，下标记为 $ i $ ~ $ j $ $ a^{b} \bmod p = ( \cdots (( a^{ c_i \cdot 2^{i} } \bmod p) \cdot a^{ c_{ i + 1 } \cdot 2^{i + 1} } \bmod p) \cdots) \cdot a^{ c_j \cdot 2^{j} } \bmod p $
2. 实现方式：
   每次让 $ b $ & 1 取出最低位，判断是否为 1，若是 1，就把这一位对应的乘积项代入上述公式计算，每次都让 $ b $ 右移一位，并把 $ a \times a \bmod p $ （强制类型转换为 long long，防止爆 int）

标签：cdot,res,bmod,875,long,int,cdots,快速,AcWing
来源： https://www.cnblogs.com/wKingYu/p/16248891.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。