【BZOJ 3884】上帝与集合的正确用法【欧拉降幂】

题目：

做法:
欧拉降幂：$a^x\equiv a^{x \ mod \ \phi(p) + \phi(p)} (mod \ p)\ (x \geq p)$
证明贴个地址：地址虽然似乎是爬虫搞出来的但是再也找不到了原出处了.
有了欧拉降幂,可以设$f(p)=2^{2^{2^{2^{...}}}}\ mod\ p$
带入欧拉降幂:

$$f(p)=2^{2^{2^{...}}mod\phi(p)+\phi(p)}mod\ p$$
$$因为2^{2^{2^{2^{...}}}}\geq p显然成立$$
注意到现在的指数和 $f$有相同的结构，即:
$$f(p)=2^{f(\phi(p))+\phi(p)}\ mod\ p$$
由此我们可以递归计算 $f$,并将$f(1)=0$（显然任意数 $mod\ 1=0$）作为边界即可.
时间复杂度似乎并不显然,但是注意到递归次数就是使得 $\phi(\phi(\phi(...\phi(p))))=1$的迭代次数.
而有以下结论: $$若p是偶数，则\phi(p)\leq p\frac{2-1}{2}=\frac{p}{2}$$
$$若p是奇数,则设p的质因子为q_i,所以\phi(p)=p\frac{q_1-1}{q_1}\frac{q_2-1}{q_2}...\frac{q_n-1}{q_n}.注意到p是奇数,则q_i为奇数，则q_i-1为偶数,故\phi(p)为偶数,转化为上一个情况$$
所以,每两次递归调用时 $p$至少减半，即调用次数为$O(\log p)$.加上欧拉函数的计算，总时间复杂度为 $O(\sqrt p \log p)$.

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>

using namespace std;

namespace MYS
{
    const int MAXN = 10000005;
    typedef long long LL;

    int phi_table[MAXN];

    int phi(int x)
    {
        if(phi_table[x]) return phi_table[x];
        int &ret = phi_table[x], upper = sqrt(x + 0.5);
        ret = x;
        for(int i=2;i<=upper;i++) {
            if(x % i == 0) {
                ret /= i;
                ret *= i-1;
                while(x % i == 0) x /= i;
            }
        }
        if(x > 1) { ret /= x; ret *= x-1; }
        return ret;
    }

    LL qpow(LL a, LL n, int p)
    {
        LL ret = 1, tmp = a % p;
        while(n) {
            if(n & 1) ret = ret * tmp % p;
            tmp = tmp * tmp % p;
            n >>= 1;
        }
        return ret;
    }   

    LL f(int p)
    {
        if(p == 1) return 0;
        int phip = phi(p);
        return qpow(2, f(phip) + phip, p);
    }

    void work()
    {
        int p;
        scanf("%d", &p);
        printf("%d\n", (int)f(p));
    }

    void Main()
    {
        int T;
        scanf("%d", &T);
        while(T--) work();
    }
}


int main()
{
    freopen("in.txt", "r", stdin);
    MYS::Main();
    return 0;
}