题目来源:The Preliminary Contest for ICPC Asia Nanjing 2019
这个题目,,,
求((a^a)^a)^a.........b个a次幂 对 m取模
用到的知识点是扩展欧拉定理及降幂
递归他的幂次数,每次他的幂次数都要取模,模数也不相同,快速幂求的时候要判断情况
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e5+10;
long long a,b,m;
long long Eular(long long n) //计算n的欧拉值
{
long long i;
long long ans=n;
for(i=2; i*i<=n; ++i)
{
if(n%i==0) //i是质数
{
ans=ans-ans/i;
while(n%i==0) //除掉所有的i
{
n/=i;
}
}
}
if(n>1)
{
ans=ans-ans/n;
}
return ans;
}
long long Suan(long long a,long long x,long long mod)
{
long long ans=1,sum=0;
while(x)
{
if(x&1)
{
ans=ans*a;
if(ans>=mod)
{
ans%=mod;
sum=1;
}
}
x>>=1;
if(x==0)
break;
a=a*a;
if(a>=mod)
{
a%=mod;
sum=1;
}
}
if(sum)
return ans%mod+mod;
return ans%mod;
}
long long solve(int l,int r,long long mod)
{
if(l==r)
{
return a%mod;
}
if(mod==1)
{
if(a<mod)
return a%mod;
else
{
return a%mod+mod;
}
}
long long phi=Eular(mod);
return Suan(a,solve(l+1,r,phi),mod);
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
scanf("%lld%lld%lld",&a,&b,&m);
if(b==0)
{
printf("%d\n",1%m);
}
else
{
printf("%lld\n",solve(1,b,m)%m);
}
}
return 0;
}