ZOJ 3209 Treasure Map (DLX)

本文主要是介绍ZOJ 3209 Treasure Map (DLX)，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

题意：

思路：

DLX的精确覆盖
这道题让我明白了DLX的大体思路方向，就是把问题转化成一种类似二分图匹配的问题，想办法建立出行和列，使之产生可匹配的关系
对于这个题，每个矩阵都对应原大矩阵的一些面积，我们给原矩阵每块面积都标号，这样形成了对应关系，就可以用DLX精确覆盖了

错误及反思：

代码：

#include <iostream>
#include <stdio.h>
#include <string.h>
//精确覆盖问题的定义：给定一个由0-1组成的矩阵，是否能找到一个行的集合，使得集合中每一列都恰好包含一个1
const int MN = 1005;//最大行数
const int MM = 1005;//最大列数
const int MNN = 1e5 + 5 + MM; //最大点数struct DLX
{int n, m, si;//n行数m列数si目前有的节点数//十字链表组成部分int U[MNN], D[MNN], L[MNN], R[MNN], Row[MNN], Col[MNN];//第i个结点的U向上指针D下L左R右，所在位置Row行Col列int H[MN], S[MM]; //记录行的选择情况和列的覆盖情况int ansd, ans[MN];void init(int _n, int _m)  //初始化空表{n = _n;m = _m;for (int i = 0; i <= m; i++) //初始化第一横行（表头）{U[i] = D[i] = i;      //目前纵向的链是空的L[i] = i - 1;R[i] = i + 1;         //横向的连起来}R[m] = 0; L[0] = m;si = m;                 //目前用了前0~m个结点memset(S, 0, sizeof(S));memset(H, -1, sizeof(H));}void link(int r, int c)    //插入点(r,c){++S[Col[++si] = c];     //si++;Col[si]=c;S[c]++;Row[si] = r;//si该结点的行数为rD[si] = D[c];//向下指向c的下面的第一个结点U[D[c]] = si;//c的下面的第一个结点的上面为siU[si] = c;//si的上面为列指针D[c] = si;//列指针指向的第一个该列中的元素设为siif (H[r]<0)//如果第r行没有元素H[r] = L[si] = R[si] = si;else{R[si] = R[H[r]];//si的右边为行指针所指的右边第一个元素L[R[H[r]]] = si;//行指针所指的右边第一个元素的左侧为siL[si] = H[r];//si的左侧为行指针R[H[r]] = si;//行指针的右侧为si}}void remove(int c)        //列表中删掉c列{L[R[c]] = L[c];//表头操作  //c列头指针的右边的元素的左侧指向c列头指针左边的元素R[L[c]] = R[c];//c列头指针的左边的元素的右侧指向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[Col[j]];//该元素所在的列数目减一}}void resume(int c)        //恢复c列{for (int i = U[c]; i != c; i = U[i])//枚举该列元素for (int j = L[i]; j != i; j = L[j])//枚举该列元素所在的行++S[Col[U[D[j]] = D[U[j]] = j]];//D[U[j]]=j;U[D[j]]=j;S[Col[j]]++;L[R[c]] = R[L[c]] = c;//c列头指针左右相连}bool dance(int d) //选取了d行{if (ansd != -1 && ansd < d)return 0;if (R[0] == 0)//全部覆盖了{//全覆盖了之后的操作if(ansd==-1)ansd = d;else if (d < ansd) ansd = d;return 1;}int c = R[0];//表头结点指向的第一个列for (int i = R[0]; i != 0; i = R[i])//枚举列头指针if (S[i]<S[c])//找到列中元素个数最少的c = i;remove(c);//将该列删去for (int i = D[c]; i != c; i = D[i]){//枚举该列的元素ans[d] = Row[i];//记录该列元素的行for (int j = R[i]; j != i; j = R[j])remove(Col[j]);//将该列的某个元素的行上的元素所在的列都删去dance(d + 1);for (int j = L[i]; j != i; j = L[j])resume(Col[j]);}resume(c);return 0;}
}dlx;
inline int read() {char ch = getchar(); int x = 0, f = 1;while(ch < '0' || ch > '9') {if(ch == '-') f = -1;ch = getchar();} while('0' <= ch && ch <= '9') {x = x * 10 + ch - '0';ch = getchar();} return x * f;
}
int main()
{int T;T=read();while(T--){int n, m, k;n=read();m=read();k=read();dlx.init(k-1 ,n*m );for(int i=0;i<k;i++){int x1,x2,y1,y2;x1=read();y1=read();x2=read();y2=read();for(int j=x1;j<x2;j++)for(int l=y1;l<y2;l++)dlx.link(i,j*m+l+1);}dlx.ansd=-1;dlx.dance(0);printf("%d\n", dlx.ansd);}
}