[BZOJ3884]上帝与集合的正确用法 欧拉函数

求ANS=2^2^2^2^2…mod p的值。
设ANS=2^n，假设2^k||p，那么2^(n-k)=ooo (mod p/2^k) ，ANS=ooo*2^k<满足欧拉定理的条件，于是2^((n-k)%φ(p/2^k))=ooo (mod p/2^k) 转化为求n%φ(p/2^k)，因为n和ANS都是相同形式，递归即可。边界是p==1是返回0。
代码：

#include<iostream>
#include<cstdio>
#define ll long long
using namespace std;
int num,pri[1000010],phi[10000010];
bool flag[10000010];
void getpri(int n)
{
    flag[1]=1;phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!flag[i]) {pri[++num]=i;phi[i]=i-1;}
        for(int j=1;j<=num&&i*pri[j]<=n;j++) 
        {
            flag[i*pri[j]]=1;
            if(i%pri[j]==0) {phi[i*pri[j]]=phi[i]*pri[j];break;}
            else phi[i*pri[j]]=phi[i]*(pri[j]-1);
        }
    }
}
ll ksm(ll a,ll b,int mod){ll re=1; while(b){if(b&1)re=re*a%mod;a=a*a%mod;b>>=1;} return re;}
ll solve(int x)
{
    ll tmp=0;
    while(!(x&1)){x>>=1;tmp++;}
    if(x==1) return 0;
    return ksm(2,((solve(phi[x])-tmp)%phi[x]+phi[x])%phi[x],x)*(1ll<<tmp);
}

int main()
{
    getpri(10000000);
    int ca;
    scanf("%d",&ca);
    while(ca--)
    {
        int mod;
        scanf("%d",&mod);
        printf("%lld\n",solve(mod));
    }
}