#include <bits/stdc++.h>
using namespace std;
long long cal(long long a, long long n, long long p) {
  if (a * n < p) return max(a, p - a * n);
  long long auxn = a * n / p;
  if (a * n % p < p - p % a - a) --auxn;
  return cal(min(a - p % a, p % a), auxn, a);
}
void solve() {
  int t;
  scanf("%d", &t);
  while (t--) {
    long long a, n, p, h, res;
    scanf("%I64d%I64d%I64d%I64d", &a, &n, &p, &h);
    a %= p;
    if (a * n < p)
      res = a;
    else
      res = cal(a, n, p);
    if (res > h)
      puts("NO");
    else
      puts("YES");
  }
}
int main() {
  solve();
  return 0;
}
