快速幂

C++

#include<bits/stdc++.h>
using namespace std;
long long a,b,mod;
long long qmod(long long n,long long m){
	long long ans=1;
	while(m){
		if(m&1) ans=ans*n%mod;
		n=n*n%mod;
		m>>=1;
	}
	return ans;
}
int main(){
    cin>>a>>b>>mod;
    cout<<a<<"^"<<b<<" mod "<<mod<<"="<<qmod(a,b);
	return 0;
}

py

a,b,p=input().split()
a,b,p=int(a),int(b),int(p)
m=b
num=1
n=a
while(m):
    if(m%2==1):
        num=num*n%p
    n=n*n%p
    m=m//2
print(f"{a}^{b} mod {p}={num}")