【模板】超级幂(欧拉降幂)

题意:求a ^ a ^ a ^ a…(b次)对mod取模

#include
#include
#include
using namespace std;
long long phi(long long x)
{
    long long res=x;
    for(long long i=2;i*i<=x;i++)
    {
        if(x%i==0)
        {
            res=res/i*(i-1);
            while(x%i==0) x/=i;
        }
    }
    if(x!=1) res=res/x*(x-1);
    return res;
}
long long pow(long long a,long long b,long long mod)
{
    long long res=1;
    a=a%mod;
    while(b)
    {
        if(b&1) res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
long long solve(long long a,long long b,long long mod)
{
    if(mod==1) return 0;
    if(b==0) return 1;
    long long t=phi(mod);
    long long res=solve(a,b-1,t);
    if(res>0&&res

你可能感兴趣的:(模板,数论)