【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
;
return M
= power(2,calc(phi(n))+phi(n),n);
}

int main() {

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

return 0;
}