扩展欧拉定理
注意φ(m)的求法 和细节处理
#include<bits/stdc++.h>
using namespace std;
#define int long long
int a,m;
void MOD(int& x) {
x=(x%m+m)%m;
}
string s;
int qpow(int a,int b) {
int tmp=1;
while(b) {
if(b&1)tmp*=a,MOD(tmp);
b>>=1;
a=a*a,MOD(a);
}
return tmp;
}
signed main() {
cin>>a>>m;
MOD(a);
int mm=m;//calculate φ(m)
int phi=1;
for(int i=2; i*i<=mm; i++) {
if(mm%i==0) {
phi*=(i-1);
mm/=i;
while(mm%i==0) {
phi*=i;
mm/=i;
}
}
}
if(mm>1)phi*=(mm-1);
int bm=0;
char ch=getchar();
while(ch<'0'||ch>'9')ch=getchar();
int flag=0;
while(ch>='0'&&ch<='9') {
bm=(int)10*bm+ch-'0';
if(bm>phi)flag=1,bm=(bm%phi+phi)%phi;
ch=getchar();
}
if(bm>=phi)flag=1,bm=(bm%phi+phi)%phi;
if(flag)bm+=phi;
int ans=qpow(a,bm);
MOD(ans);
cout<<ans<<"\n";
return 1+1==3;
}