#include using namespace std; #define rep(i, n) for (int i = 0; i < (int)n; i++) using ll = long long int; int mod_pow(int x, int n, int mod) { int ret = 1; while (n > 0) { if ((n & 0x01) != 0) ret = (long long)ret * x % mod; x = (long long)x * x % mod; n >>= 1; } return ret; } bool miller_rabin(ll n) { // is prime if (n == 2) return true; if (n == 1 || n % 2 == 0) return false; ll s = 0; while (((n-1) & -(n-1)) > (1 << s)) s++; int d = (n-1) / ((n-1) & -(n-1)); // n-1 = 2^s*d vector tests = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37}; // n - 1 = a^s * d for (int a : tests) { if (a == n) continue; ll x = mod_pow(a, d, n); // x = a^d(mod n) if (x == 1) continue; int r = 0; while (x != n-1) { // r = 0 -> s x = mod_pow(x, 2, n); // x = (a^d)^2(mod n) = a^(2^(rd))(mod n) r++; if (x == 1 || r == s) // x^2 = 1(mod n) -> not prime number return false; } } return true; } int main() { ll N; cin >> N; vector tests = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37}; for (int n : tests) { ll p = pow(2, n) - 1; if (miller_rabin(p)) { ll q = pow(2, n-1) * p; if (N == q) { cout << "Yes" << endl; return 0; } } } cout << "No" << endl; return 0; }