关于这个问题,它的核心就是 (a*b)%n =(a%n*b%n)%n
那么下面给出自己写的两个算法:
ll modexp(ll a, ll x, ll n)
{
ll ret=1;
ll temp=a;
while(b)
{
if(x&0x1) ret=ret*temp%n;
temp=temp*temp%n;
x>>=1;
}
return ret;
}
//递归计算,注意一些边界条件
void expmod(int a,int b,int n,int &ans)
{
if(a==0){ans=0;return;}
else if(b==0){ans=1;return;}
else if(b==1)
{
ans=a%n;
return ;
}
else
{
if(b%2==0)
{
expmod(a,b/2,n,ans);
ans=ans*ans%n;
}
else
{
expmod(a,b/2,n,ans);
int temp;
temp=ans*(a%n)%n;
ans=ans*temp%n;
}
}
}