2021-07-31

最新推荐文章于 2024-09-26 09:13:29 发布

兰德朔

最新推荐文章于 2024-09-26 09:13:29 发布

阅读量74

点赞数

分类专栏：算法文章标签：快速幂大整数运算模运算效率优化算法实现

本文链接：https://blog.csdn.net/qq_38229601/article/details/119270051

版权

算法专栏收录该内容

9 篇文章 0 订阅

订阅专栏

快速幂

核心思路：反复平方法

要求
$a^k\quad mod \quad m$
的结果，我们要先预计算出：
$a^{2^0}\quad mod \quad m\\ a^{2^1}\quad mod \quad m\\ a^{2^3}\quad mod \quad m\\ ...\\ a^{2^{log_2k}}\quad mod \quad m$
所以，原值
$a^k\quad mod \quad m= a^{x_1}\quad mod \quad m+a^{x_2}\quad mod \quad m+...+a^{x_t}\quad mod \quad m$
比如
$k = (110101)_2$
那么
$k = 2^0 + 2^2 + 2^4 + 2^5$
在预处理过程中，
$a^{2^{log_2k}}\quad mod \quad m = (a^{2^{log_2k - 1}}\quad mod \quad m)* (a^{2^{log_2k - 1}}\quad mod \quad m)\quad mod \quad m = (a^{2^{log_2k - 1}}\quad mod \quad m)^2\quad mod \quad m$
所以，我们预处理时，直接平方前一个数

给定 n 组 a_i,b_i,p_i，对于每组数据，求出 a^bi_i mod p_i 的值。

输入格式

第一行包含整数 nn。

接下来 n 行，每行包含三个整数 a_i,b_i,p_i。

输出格式

对于每组数据，输出一个结果，表示 a^bi_imodpi 的值。

每个结果占一行。

数据范围

1≤n≤100000
1≤a_i,b_i,p_i≤2×10⁹

输入样例：

2
3 2 5
4 3 9

输出样例：

4
1

#include<iostream>
#include<algorithm>
using namespace std;

typedef long long LL;

// a^k % p
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, k, p;
        scanf("%d%d%d", &a, &k, &p);
        printf("%d\n", qmi(a, k, p));
    }

    return 0;
}