#include #define M_PI 3.14159265358979323846 // pi using namespace std; typedef long long ll; typedef unsigned long long ull; typedef pair P; typedef tuple t3; typedef tuple t4; #define rep(a,n) for(ll a = 0;a < n;a++) #define repi(a,n) for(ll a = n - 1;a >= 0;a--) template void chmax(T& reference, T value) { reference = max(reference, value); } template void chmaxmap(map& m, T key, T value) { if (m.count(key)) { chmax(m[key], value); } else { m[key] = value; } } template void chmin(T& reference, T value) { reference = min(reference, value); } #include using namespace atcoder; typedef modint1000000007 mint; constexpr ll mod = 1000000007; constexpr ll mpow(ll x, ll n) { ll ans = 1; x %= mod; while (n != 0) { if (n & 1) ans = ans * x % mod; x = x * x % mod; n = n >> 1; } return ans; } class Primes { private: vector Prime_Number; vector is_prime_; public: Primes(int N) { is_prime_.resize(N + 1, true); is_prime_[0] = is_prime_[1] = false; for (int i = 0; i < N + 1; i++) { if (is_prime_[i]) { Prime_Number.push_back(i); for (int j = 2 * i; j <= N; j += i) is_prime_[j] = false; } } } int operator[](int i) { return Prime_Number[i]; } int size() { return Prime_Number.size(); } int back() { return Prime_Number.back(); } bool isPrime(int q) { return is_prime_[q]; } }; class SDivisor { public: vector F; SDivisor(ll N) { for (ll i = 1; i * i <= N; i++) { if (N % i == 0) { F.push_back(i); if (i * i != N) F.push_back(N / i); } } sort(begin(F), end(F)); } int size() { return F.size(); } ll operator[](int k) { return F[k]; } const vector& factors() { return F; } }; class Divisor { private: vector F; vector> pfactorize; public: Divisor(ll N) { for (ll i = 1; i * i <= N; i++) { if (N % i == 0) { F.push_back(i); if (i * i != N) F.push_back(N / i); } } sort(begin(F), end(F)); Primes p((ll)sqrt(N) + 1); for (int i = 0; i < p.size(); i++) { pfactorize.emplace_back(p[i], 0); while (N % p[i] == 0) { N /= p[i]; pfactorize.back().second++; } if (pfactorize.back().second == 0) pfactorize.pop_back(); } if (N > 1) pfactorize.emplace_back(N, 1); } int size() { return F.size(); } const vector>& pfac() { return pfactorize; } ll operator[](int k) { return F[k]; } const vector& factors() { return F; } }; int main() { ll a, b; cin >> a >> b; if (a == b) { cout << "Yes" << endl; } else { if (a > b) swap(a, b); if (a == 1) { cout << "No" << endl; } else { ll c = 0; while (b % a == 0) { c++; b /= a; } if (b == 1 && c == a) { cout << "Yes" << endl; } else { cout << "No" << endl; } } } return 0; }