#ifndef ONLINE_JUDGE #define _GLIBCXX_DEBUG #endif #include using namespace std; #define pass (void)0 #define INF (1<<30)-1 #define INFLL (1LL<<60)-1 #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define repr(i, n) for (int i = (int)(n) - 1; i >= 0; i--) #define rep2(i, a, b) for (int i = (int)(a); i < (int)(b); i++) #define repr2(i, a, b) for (int i = (int)(b) - 1; i >= (int)(a); i--) #define all(x) (x).begin(), (x).end() #define rall(x) (x).rbegin(), (x).rend() #define sz(x) ((int)(x).size()) #define YesNo(cond) cout << ((cond) ? "Yes\n" : "No\n") #define YESNO(cond) cout << ((cond) ? "YES\n" : "NO\n") using ll = long long; using pii = pair; using pll = pair; using vi = vector; using vl = vector; using vvi = vector; using vvl = vector; template void print(const T& value) { cout << value << "\n"; } template void print(const vector& vec) { for (auto& v : vec) cout << v << " "; cout << "\n"; } template void input(vector& vec) { for (auto& v : vec) cin >> v; }; template bool chmin(T& a, const T& b) { if (a > b) { a = b; return true; } return false; } template bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } return false; } #include using namespace atcoder; using mint = modint998244353; using vm = vector; ll MOD; long long pow2(long long x, long long n) { long long ret = 1; while (n > 0) { if (n & 1) ret = ret * x % MOD; // n の最下位bitが 1 ならば x^(2^i) をかける x = x * x % MOD; n >>= 1; // n を1bit 左にずらす } return ret; } ll inverse(ll n, ll d) { return n % MOD * pow2(d, -1) % MOD; } ll geometric_progression(ll a, ll r, ll n) { return a*((pow2(r, n)-1)%MOD)%MOD*inverse(1, r-1)%MOD; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(10); ll N; cin >> N >> MOD; if (N == 0) { YesNo(false); return 0; } if (MOD == 2) { YesNo(N%2 == 1); return 0; } if ((N+1)%MOD == 0 || geometric_progression(1, MOD-1, N+1) == 0) { YesNo(true); return 0; } random_device rd; mt19937 rng(rd()); clock_t start = clock(); while ((double)(clock()-start) <= 1.95) { uniform_int_distribution dist(2, MOD); ll n = dist(rng); if (geometric_progression(1, n, N+1) == 0) { YesNo(true); return 0; } } YesNo(false); }