【BZOJ 3884】 上帝与集合的正确用法

【题目链接】

           点击打开链接

【算法】

          通过欧拉拓展定理，列出递推公式

 【代码】

             

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

ll T,N;
map<ll,ll> M;

template <typename T> inline void read(T &x) {
        ll f = 1; x = 0;
        char c = getchar();
        for    (; !isdigit(c); c = getchar()) { if (c == '-') f = -f; }
        for (; isdigit(c); c = getchar()) x = (x << 3) + (x << 1) + c - '0';
        x *= f;
}
template <typename T> inline void write(T x) {
        if (x < 0) { putchar('-'); x = -x; }    
        if (x > 9) write(x/10);
        putchar(x%10+'0'); 
}
template <typename T> inline void writeln(T x) {
        write(x);
        puts("");    
}

ll power(ll a,ll n,ll P) {
        ll res;
        if (n == 0) return 1;
        if (n == 1) return a % P;
        res = power(a,n>>1,P);
        res = res * res % P;
        if (n & 1) res = res * a % P;
        return res;     
}

ll phi(ll x) {
        ll i,ret=x;
        for (i = 2; i <= sqrt(x); i++) {
                if (x % i == 0) {
                        while (x % i == 0) x /= i;
                        ret = ret / i * (i - 1);
                }
        }    
        if (x > 1) ret = ret / x * (x - 1);
        return ret; 
} 

inline ll calc(ll n) {
        if (M.count(n)) return M[n];
        return M[n] = power(2,calc(phi(n))+phi(n),n);     
}

int main() {
    
        M[1] = 0;
        read(T);
        while (T--) {
                read(N);
                writeln(calc(N));    
        }
        
        return 0;
}