#define _USE_MATH_DEFINES #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; unsigned xorshift() { static unsigned x=123456789, y=362436069, z=521288629, w=88675123; unsigned t=(x^(x<<11)); x=y; y=z; z=w; return w=(w^(w>>19))^(t^(t>>8)); } unsigned long long xorshift64() { unsigned long long a = xorshift(); a <<= 32; a |= xorshift(); return a; } long long modmul(long long a, long long b, long long mod) { long long ret = 0; for(int i=62; i>=0; --i){ ret *= 2; if(a & (1LL << i)) ret += b; ret %= mod; } return ret; } long long modpow(long long a, long long b, long long mod) { long long ret = 1; 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=50) { long long d = n - 1; int s = 0; while(d % 2 == 0){ ++ s; d /= 2; } while(--loop >= 0){ long long a = xorshift64() % (n - 1) + 1; long long x = modpow(a, d, n); if(x != 1){ bool isComposite = true; for(int r=0; r> q; while(--q >= 0){ long long n; cin >> n; if(solve(n)) cout << "Yes" << endl; else cout << "No" << endl; } return 0; }