#include<bits/stdc++.h>
#define y1 yy
using namespace std;
const int N=100010;
const int M=8000800;
const int inf=1000000000;
int t,x1[N],x2[N],y1[N],y2[N],S,T,res,n,m;
vector<int> mtr[2020][2020];
struct EDGE { int adj, w, next; } edge[M];
int top, gh[N], nrl[N];
void addedge(int x, int y, int w) {
    //printf("%d %d %d\n",x,y,w);
    edge[++top].adj = y;
    edge[top].w = w;
    edge[top].next = gh[x];
    gh[x] = top;
    edge[++top].adj = x;
    edge[top].w = 0;
    edge[top].next = gh[y];
    gh[y] = top;
}
int dist[N], q[N];
int bfs() {
    memset(dist, 0, sizeof(dist));
    q[1] = S; int head = 0, tail = 1; dist[S] = 1;
    while (head != tail) {
        int x = q[++head];
        for (int p=gh[x]; p; p=edge[p].next)
            if (edge[p].w && !dist[edge[p].adj]) {
                dist[edge[p].adj] = dist[x] + 1;
                q[++tail] = edge[p].adj;
            }
    }
    return dist[T];
}
int dinic(int x, int delta) {
    if (x==T) return delta;
    for (int& p=nrl[x]; p && delta; p=edge[p].next)
        if (edge[p].w && dist[x]+1 == dist[edge[p].adj]) {
            int dd = dinic(edge[p].adj, min(delta, edge[p].w));
            if (!dd) continue;
            edge[p].w -= dd;
            edge[p^1].w += dd;
            return dd;
        }
    return 0;
}
int work() {
    int ans = 0;
    while (bfs()) {
        memcpy(nrl, gh, sizeof(gh));
        int t; while (t = dinic(S, inf)) ans += t;
    }
    return ans;
}

int uni(int x,int y){
    int a,b,c,d;
    a=max(x1[x],x1[y]);
    b=max(y1[x],y1[y]);
    c=min(x2[x],x2[y]);
    d=min(y2[x],y2[y]);
    return max(0,c-a)*max(0,d-b);
}
int main(){
    //freopen("I.in","r",stdin);
    scanf("%d",&t);
    while (t--){
        top=1;
        scanf("%d%d",&m,&n);
        S=0;
        T=m+n+1;
        for (int i=1;i<=m;i++){
            scanf("%d%d%d%d",&x1[i],&y1[i],&x2[i],&y2[i]);
            addedge(S,i,1);
        }
        for (int i=1;i<=n;i++){
            scanf("%d%d%d%d",&x1[m+i],&y1[m+i],&x2[m+i],&y2[m+i]);
            mtr[x1[m+i]][y1[m+i]].push_back(m+i);
            addedge(i+m,T,1);
        }
        for (int i=1;i<=m;i++)
        for (int j=max(0,x1[i]-3);j<=x1[i]+3;j++)
        for (int k=max(0,y1[i]-3);k<=y1[i]+3;k++)
        for (int h=0;h<mtr[j][k].size();h++){
            if (2*uni(i,mtr[j][k][h])>=(x2[i]-x1[i])*(y2[i]-y1[i])) addedge(i,mtr[j][k][h],1);
        }
        res=work();
        printf("%d\n",res);
        memset(gh,0,sizeof(gh));
        memset(edge,0,sizeof(edge));
        for (int i=1;i<=n;i++) mtr[x1[m+i]][y1[m+i]].clear();
    }
}