#include #include #include using namespace std; class Xorshift { private: unsigned x, y, z, w; public: Xorshift(unsigned seed=88675123, int loop=50){ x = 123456789; y = 362436069; z = 521288629; w = seed; while(--loop >= 0){ (*this)(); } } unsigned operator()(){ unsigned t=(x^(x<<11)); x=y; y=z; z=w; return w=(w^(w>>19))^(t^(t>>8)); } unsigned operator()(unsigned size){ return (*this)() % size; } unsigned long long get64(){ unsigned long long a = (*this)(); unsigned long long b = (*this)(); return (a << 32) | b; } unsigned long long get64(unsigned long long size){ return get64() % size; } template void shufflte(vector& v){ unsigned n = v.size(); for(unsigned i=0; i 0){ if(b & 1){ ret += a; if(ret > mod) ret -= mod; } a <<= 1; if(a > mod) a -= mod; b >>= 1; } return ret; } unsigned long long modpow(unsigned long long a, unsigned long long b, unsigned long long mod) { unsigned long long ret = 1; unsigned long long tmp = a; while(b > 0){ if(b & 1) ret = modmul(ret, tmp, mod); tmp = modmul(tmp, tmp, mod); b >>= 1; } return ret; } bool millerRabinPrimalityTest(long long n, int loop=30) { if(n < 2) return false; long long d = n - 1; int s = 0; while(d % 2 == 0){ ++ s; d /= 2; } while(--loop >= 0){ long long a = xorshift.get64(n-1) + 1; long long x = modpow(a, d, n); if(x != 1){ bool isComposite = true; for(int r=0; r 6e18 / base){ return 2e18; } } return ans; } bool testcase(long N){ if(N == 2){ return false; }else if(N % 2 == 0){ return true; } for(long i = 2; i < N; i *= 2){ for(int j = 1; j < 100; j++){ long left = 1, right = 2e18; // [left, right] while(left < right){ long middle = (left + right) / 2; long power_num = power(middle, j); if(power_num > N - i){ right = middle - 1; }else if(power_num < N - i){ left = middle + 1; }else if(millerRabinPrimalityTest(middle)){ return true; }else{ break; } } } } return false; } int main(){ int Q; cin >> Q; for(int i = 0; i < Q; i++){ long N; cin >> N; if(testcase(N)){ cout << "Yes" << endl; }else{ cout << "No" << endl; } } }