super_log

最新推荐文章于 2019-09-03 20:10:14 发布

木白星枝

最新推荐文章于 2019-09-03 20:10:14 发布

阅读量141

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/qq_42759120/article/details/100522201

版权

数论专栏收录该内容

2 篇文章 0 订阅

订阅专栏

题目来源：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;
}