快速幂实现

#include <bits/stdc++.h>
#define int long long//a^b%p
using namespace std;
int qmi(int a,int b,int p){
	int res=1;
 	while(b){// 13=2^8+2^4+2^1  二进制：1101
		 if(b&1)res=res*a%p;//如果是1，则乘a
		 a=a*a%p;//a^2,a^4,a^8....
		 b>>=1;//110  11   1
	 }
 return res;
}
signed main(){
	 int a,b,p;cin>>a>>b>>p;
	 cout<<qmi(a,b,p);
}

同模定理

(a + b) %c =( a%c + b%c ) % c;//加法
(a - b) %c =( a%c - b%c ) % c;//减法
(a - b) %c =(a - b) %p + p) % p;//负数在模p后,变成正数
(a * b) %c =( a%c * b%c ) % c;//乘法，除法不适用