ZOJ-3209___Treasure Map —— DLX精确覆盖

题目链接：点我啊╭(╯^╰)╮

题目大意：

    给你一个$n\ast m$的矩形和 $p$ 个小矩形，求最少需要几个小矩形可以精确覆盖这个大矩形？？？

解题思路：

    明显是最清晰的舞蹈链，那么问题就在于建图上，很多人接触这一题应该都是刚学完模板不久，这里要用 $p$ 个小矩形来填满，求最少要几个，那么我们图形中的“行”就变成了这 $p$ 个小矩形
    行的问题解决了，列呢？？这里因为大矩形是$n\ast m$，那么我们将这整个图分成$n\ast m$个小正方形，这个图列就是$n\ast m$
    换句话说，我们建的图是行为 $p$，列为$n\ast m$的 $01$ 图，为什么是$n\ast m$呢，因为每个小矩形最多由这么多个小正方形组成，相当于把每个矩形存放在了一行里（二维变一维）

代码思路：

    注意，这里与模板不同的是，模板只是找到了精确覆盖就返回，本题要求最小覆盖个数，不能带有返回值

核心：初步掌握舞蹈链的思想并加以转换

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

struct DancingLink{
	const static int N = 500010;
	const static int M = 1010;
    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);
	}
}