875. 快速幂
文章目录
- 题目描述
- 输入格式:
- 输出格式:
- 数据范围
- 输入样例
- 输出样例
- 方法:快速幂
- 解题思路
- 代码
- 复杂度分析:
题目描述
给定 n n n 组 a i , b i , p i a_i,b_i,p_i ai,bi,pi,对于每组数据,求出 a i b i m o d p i a_i^{b_i}mod\ p_i aibimod pi 的值。
输入格式:
第一行包含整数 n。
接下来 n 行,每行包含三个整数 a i , b i , p i a_i,b_i,p_i ai,bi,pi。
输出格式:
对于每组数据,输出一个结果,表示 a i b i m o d p i a_i^{b_i}mod\ p_i aibimod pi 的值。
每个结果占一行。
数据范围
- 1 ≤ N ≤ 1 0 5 1≤N≤10^5 1≤N≤105
- 1 ≤ a i , b i , p i ≤ 2 × 1 0 9 1≤a_i,b_i,p_i≤2\times10^9 1≤ai,bi,pi≤2×109
输入样例
2
3 2 5
4 3 9
输出样例
4
1
方法:快速幂
解题思路
我们用 k 替换 b。
最朴素的做法是,循环 k 次,每次乘 a 再模上 p.
由于这里的 k 比较大,只能另觅他法了。
快速幂的思想就是把 a 的 k 次幂看作是多个 a 的 2 的多少次幂相乘。
每一次循环,先判断 k 的最后一位是否为 1,如果为 1,则需要把
r
e
s
∗
a
m
o
d
p
res * a \ mod \ p
res∗a mod p;
然后将 k 右移 1 位,a 要平方再 mod p
Tips
- 因为这里的数据范围比较大,往往数论的题目都要定义为 long long,不然会爆 int.
代码
#include <iostream>
using namespace std;
typedef long long LL;
int n;
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() {
cin >> n;
while(n--) {
int a, k, p;
scanf("%d%d%d", &a, &k, &p);
printf("%d\n", qmi(a, k, p));
}
return 0;
}
复杂度分析:
- 时间复杂度: O ( l o g 2 k ) O(log_2k) O(log2k)
- 空间复杂度: O ( 1 ) O(1) O(1)