題目大意：

    給你一個$n*m$的矩形和 $p$ 個小矩形，求最少需要幾個小矩形可以精確覆蓋這個大矩形？？？

解題思路：

    明顯是最清晰的舞蹈鏈，那麼問題就在於建圖上，很多人接觸這一題應該都是剛學完模板不久，這裡要用 $p$ 個小矩形來填滿，求最少要幾個，那麼我們圖形中的“行”就變成了這 $p$ 個小矩形     行的問題解決了，列呢？？這裡因為大矩形是$n*m$，那麼我們將這整個圖分成$n*m$個小正方形，這個圖列就是$n*m$     換句話說，我們建的圖是行為 $p$，列為$n*m$的 $01$ 圖，為什麼是$n*m$呢，因為每個小矩形最多由這麼多個小正方形組成，相當於把每個矩形存放在了一行裡（二維變一維）

程式碼思路：

    注意，這裡與模板不同的是，模板只是找到了精確覆蓋就返回，本題要求最小覆蓋個數，不能帶有返回值

核心：初步掌握舞蹈鏈的思想並加以轉換

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

#define N 500010
#define M 1010
struct DancingLink{
    int n,s,ansd;//列數 節點總數 
    int S[M],A[M],H[M];//S[]該列節點總數  A[]答案  H[]行首指標 
    int L[N],R[N],U[N],D[N]; //L[],R[],U[],D[] 上下左右 
    int X[N],C[N];//X[] C[] 行列編號 
    void init(int n){//初始化 
        this->n=n;
        for(int i=0;i<=n;i++)
            U[i]=i,D[i]=i,L[i]=i-1,R[i]=i+1;
        R[n]=0,L[0]=n;s=n+1;
		memset(S,0,sizeof(S));
		memset(H,-1,sizeof(H));
    }
    void DelCol(int c){//刪除列 
        L[R[c]]=L[c];R[L[c]]=R[c];
        for(int i=D[c];i!=c;i=D[i])
            for(int j=R[i];j!=i;j=R[j])
                U[D[j]]=U[j],D[U[j]]=D[j],--S[C[j]];
    }
    void ResCol(int c){//恢復列 
        for(int i=U[c];i!=c;i=U[i])
            for(int j=L[i];j!=i;j=L[j])
                ++S[C[j]],U[D[j]]=j,D[U[j]]=j;
        L[R[c]]=c,R[L[c]]=c;
    }
    void AddNode(int r,int c){//新增節點 
        ++S[c],C[++s]=c,X[s]=r;
        D[s]=D[c],U[D[c]]=s,U[s]=c,D[c]=s;
        if(H[r]<0) H[r]=L[s]=R[s]=s;//行首節點
        else  R[s]=R[H[r]],L[R[H[r]]]=s,L[s]=H[r],R[H[r]]=s;
    }
    void dfs(int d){//深度，深搜遍歷 
    	if(ansd!=-1 && ansd<=d) return;
        if(!R[0]){
        	if(ansd==-1 || ansd>d)
            	ansd=d;
			return;
        }
        int c=R[0];
        for(int i=R[0];i;i=R[i]) if(S[i]<S[c]) c=i;
        DelCol(c);
        for(int i=D[c];i!=c;i=D[i]){
            A[d]=X[i];
            for(int j=R[i];j!=i;j=R[j]) DelCol(C[j]);
            dfs(d+1);
            for(int j=L[i];j!=i;j=L[j]) ResCol(C[j]);
        }
        ResCol(c);
    }

} dlx; 

int main(){
	int cas, n, m, p;
	scanf("%d", &cas);
	while(cas--) {
		scanf("%d%d%d", &n, &m, &p);
		dlx.init(n*m);
		for(int k=1; k<=p; k++) {
			int x1, x2, y1, y2;
			scanf("%d%d%d%d", &x1, &y1, &x2, &y2);
			for(int i=x1+1; i<=x2; i++)
				for(int j=y1+1; j<=y2; j++)
					dlx.AddNode(k, m*(i-1)+j);
		}
		dlx.ansd=-1;
		dlx.dfs(0);
		printf("%d\n",dlx.ansd);
	}
}