#include <cstdio>

bool cx[110000];
bool cy[110000];
int sx[110000],sy[110000];

int main() {
    int T;
    scanf("%d",&T);
    while(T--) {
        int n,m,K,Q;
        scanf("%d%d%d%d",&n,&m,&K,&Q);
        for(int i=0;i<=n;++i) cx[i] = false;
        for(int i=0;i<=m;++i) cy[i] = false;
        for(int i=0;i<K;++i) {
            int x,y;
            scanf("%d%d",&x,&y);
            cx[x] = true;
            cy[y] = true;
        }
        for(int i=1;i<=n;++i) {
            sx[i] = sx[i-1] + (cx[i] ? 1 : 0);
        }
        for(int i=1;i<=m;++i) {
            sy[i] = sy[i-1] + (cy[i] ? 1 : 0);
        }
        for(int i=0;i<Q;++i) {
            int x1,x2,y1,y2;
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            int xx = sx[x2] - sx[x1-1];
            int yy = sy[y2] - sy[y1-1];
            if( xx == x2-x1+1 || yy == y2-y1+1 ) {
                puts("Yes");
            }
            else {
                puts("No");
            }
        }
    }
    return 0;
}