1292765944的专栏

和POJ 3074一样 #include #define inf 1e8 #define M 16*16*16*16*16*4+100 #define N 16*16*16+100 int U[M],D[M],L[M],R[M],C[M],X[M];//x代表M所属的行，c代表M所属的列，U，D，L，R为一个元素的上下左右指针 int H[N],S[N],Q[N];//H是水平循环链表的头指针

建图还是坑了我一下午。。。。终于把数独搞定了(囧) 具体细节我加了注释 #include #define inf 1e8 #define M 9*9*9*9*9*4+100 #define N 9*9*9+100 int U[M],D[M],L[M],R[M],C[M],X[M];//x代表M所属的行，c代表M所属的列，U，D，L，R为一个元素的上下左右指针 int H[N],S[N],Q[

#include #define N 2011 #define M 102011 int U[M],D[M],L[M],R[M],C[M],X[M];//x代表M所属的行，c代表M所属的列，U，D，L，R为一个元素的上下左右指针 int H[N],S[N],Q[N];//H是水平循环链表的头指针，S代表每一列的元素个数，Q存储一个可能结果 int size,n,m; void remove(in

关键是建好图 把n*m的矩阵看成n*m个单位元，作为n*m列；每一个矩形一行。 问题即转化为从这些行中选择最少的一部分使每一列被覆盖且仅覆盖一次。 #include #define inf 10000000 #define N 1005 #define M 1024*505 int U[M],D[M],L[M],R[M],C[M];//C代表M所属的列，U，D，L，R为一个元素的上下左右指针

DLX重复覆盖对于当前矩阵的处理是，将当前列去掉，并将选作当前列的解的行能够覆盖到的列全部去掉。因为不需要每列仅由一个1去覆盖，因此不必要把能够覆盖某一列的所有行全部去掉。因此remove和resume函数的写法将会有所不同（具体可参考下面我写的代码） 第二个区别是，由于矩阵密度下降会变慢（因为去掉的少了），因此要加上一个强剪枝。这个剪枝利用的思想是A*搜索中的估价函数。即，对于当前的递归深度K