#include <bits/stdc++.h>
using namespace std;
long long max_dis(long long a, long long p, long long n) {
  long long m = p / a;
  long long n1 = a * n / p;
  if (n1 && a * n % p < (m - 1) * a) n1--;
  if (!n1) return a * n / p ? a : max(a, p - a * n);
  return max_dis(min(p % a, p - p % a), a, n1);
}
bool gao(long long a, long long n, long long p, long long h) {
  return min(max_dis(a % p, p, n), a % p) <= h;
}
int main() {
  int T;
  scanf("%d", &T);
  while (T--) {
    int a, n, p, h;
    scanf("%d%d%d%d", &a, &n, &p, &h);
    puts(gao(a, n, p, h) ? "YES" : "NO");
  }
}
