数学知识 快速幂

快速幂

给定 n 组 a_i，b_i，p_i，对于每组数据，求出 $a_i^{bi}mod{p}_i$ 的值。

输入格式
第一行包含整数 n。

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

输出格式
对于每组数据，输出一个结果，表示 $a_i^{bi}mod{p}_i$ 的值。

每个结果占一行。

数据范围
1≤ n ≤100000,
1≤ a_i，b_i，p_i ≤2×10⁹
输入样例：

2
3 2 5
4 3 9

输出样例：

4
1

快速幂：快速求出 $a^{b}modp$ 的结果（时间复杂度 O(logb)）
基本思路：

1. 首先预处理出 $a^{2^0},a^{2^1},a^{2^2},\dots ,a^{2^{logb}}.$这b个数。
2. 想办法用 1 求出的数字表示 $a^b$ 。例如： $a^b=a^{2^i}\ast a^{2^j}\ast a^{2^k}...=a^{2^i+2^j+2^k...}$ 上述式子可以用二进制优化，即用二进制表示出来。

eg: $4^5=4^{(101{)}_2}=4^{2^0+2^2}=4^{2^0}\ast 4^{2^2}=4\ast 256$

#include<bits/stdc++.h>
using namespace std;
const int N = 100010;
typedef long long int LL;
LL n,a,b,p;
LL quick_pow(LL a,LL b,LL p)
{
    LL res=1;
    while(b)
    {
        int k=b&1;
        if(k) res=res*a%p;
        a=a*a%p; //预处理a
        b=b>>1;
    }
    return res;
}
int main()
{
    cin>>n;
    while(n--)
    {
        cin>>a>>b>>p;
        cout<<quick_pow(a,b,p)<<endl;
    }
    return 0;
}