【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 = ; x = ;
        char c = getchar();
        for    (; !isdigit(c); c = getchar()) { if (c == '-') f = -f; }
        for (; isdigit(c); c = getchar()) x = (x << ) + (x << ) + c - '';
        x *= f;
}
template <typename T> inline void write(T x) {
        if (x < ) { putchar('-'); x = -x; }
        if (x > ) write(x/);
        putchar(x%+'');
}
template <typename T> inline void writeln(T x) {
        write(x);
        puts("");
}

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

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

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

int main() {

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

        return ;
}