#pragma GCC optimize("O3")
//#pragma GCC target("avx2")
//#pragma GCC optimize("unroll-loops")

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef vector<int> vi;
typedef vector<ll> vll;
typedef vector<ld> vld;
typedef vector<bool> vb;
typedef pair<ll,ll> Pll;
typedef pair<int,int> Pin;

ll INF = 1e16;
int inf = 1e9;

#define ALL(x) (x).begin(), (x).end()
#define FOR(i, m, n) for (ll i = (m); i < (n); ++i)
#define REVFOR(i, m, n) for (ll i = ((n) - 1); i >= (m); --i)
#define REP(i, n) FOR(i, 0, n)
#define REVREP(i, n) REVFOR(i, 0, n)
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define eb emplace_back
#define bcnt __builtin_popcountll

// inputs
#define VLLIN(x, n) vll x(n); REP(i, n) cin >> x[i];
#define VIIN(x, n) vi x(n); REP(i, n) cin >> x[i];
#define VLDIN(x, n) vld x(n); REP(i, n) cin >> x[i];
#define SIN(s, n) string s; cin >> s;

#ifdef LOCAL
#include <prettyprint.hpp>
#define debug(...)  cerr << "[" << #__VA_ARGS__ << "]: ", d_err(__VA_ARGS__);
#else
#define debug(...) 83;
#endif

void d_err() {
    cerr << endl;
}

template <typename H, typename... T>
void d_err(H h, T... t) {
    cerr << h << " ";
    d_err(t...);
}

template <typename T>
void print(T x) {
    cout << x << "\n";
}

template <typename T>
void print(vector<T>& x) {
    int N = x.size();
    REP(i, N) {
        if (i > 0) cout << " ";
        cout << x[i];
    }
    cout << "\n";
}

// return primes between [2, n]
void SOE(ll n, vll& p) {
    vll s(n + 5);
    REP(i, s.size()) s[i] = i;

    p.eb(2);
    for (ll j = 4; j <= n; j += 2) {
        s[j] = 2;
    }

    for (ll i = 3; i <= n; i += 2){
        if (s[i] != i) continue;
        p.pb(i);
        for (ll j = i+i; j <= n; j += i) {
            s[j] = i;
        }
    }
}

int main(){
    cin.tie(0);
    ios_base::sync_with_stdio(false);
    cout << fixed << setprecision(20);

    ll x, y, a, b;
    cin >> x >> a >> y >> b;

    vll prime;
    SOE(1000000, prime);

    map<ll, int> ca, cb;
    for (auto &p: prime) {
        while(x % p == 0) {
            ca[p]++;
            x /= p;
        }

        while(y % p == 0) {
            cb[p]++;
            y /= p;
        }
    }

    if (x > 1) ca[x]++;
    if (y > 1) cb[y]++;

    for (auto &e: cb) {
        if (ca[e.fi] * a < cb[e.fi] * b) {
            print("No");
            return 0;
        }
    }

    print("Yes");
}