zju 3209 dancing links 求取最小行数

题目可以将每一个格子都看做是一列，每一个矩形作为1行，将所有格子进行标号，在当前矩形中的格子对应行的标号为列，将这个点加入到十字链表中

最后用dlx求解精确覆盖即可，dance()过程中记得剪枝

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <climits>

using namespace std;
#define N 1005
#define M 505
#define MAXNODE 500100
const int INF = INT_MAX;

struct DLX{
    int n,m,size;
    int L[MAXNODE] , R[MAXNODE] , U[MAXNODE] , D[MAXNODE];
    int col[MAXNODE] , row[MAXNODE] , first[M];
    int cnt_col[N];
    int ans;

    void init(int _n , int _m)
    {
        n=_n , m=_m , size = _m , ans=INF;
        for(int i=0 ; i<=m ; i++){
            U[i] = D[i] = i;
            L[i] = i-1 , R[i] = i+1;
        }
        L[0] = m , R[m] = 0;
        for(int i=1 ; i<=n ; i++) first[i] = -1;
        for(int i=1 ; i<=m ; i++) cnt_col[i] = 0;
    }

    void link(int r , int c)
    {
        ++size;
        D[size] = D[c] , U[D[c]] = size;
        U[size] = c , D[c] = size;
        cnt_col[c]++;

        if(first[r]<0){
            first[r] = size;
            L[size]= R[size] = size;
        }
        else{
            R[size] = R[first[r]] , L[R[first[r]]] = size;
            L[size] = first[r] , R[first[r]] = size;
        }
        row[size]=r , col[size]=c;
    }

    void Remove(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]){
                D[U[j]] = D[j] , U[D[j]] = U[j];
                --cnt_col[col[j]];
            }
        }
    }

    void Resume(int c)
    {
        for(int i=U[c] ; i!=c ; i=U[i]){
            for(int j=L[i] ; j!=i ; j=L[j]){
                U[D[j]] = D[U[j]] = j;
                ++cnt_col[col[j]];
            }
        }
        R[L[c]] = L[R[c]] = c;
    }

    void Dance(int d)
    {
        if(d>=ans) return;
        if(!R[0]){
            ans = min(ans , d);
            return;
        }
        int st = R[0];
        for(int i=st ; i!=0 ; i=R[i])
            if(cnt_col[st]>cnt_col[i]) st = i;

        Remove(st);
        for(int i=D[st] ; i!=st ; i=D[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(st);
    }
}dlx;

int main()
{
   // freopen("a.in" , "r" , stdin);
    int T;
    scanf("%d" , &T);
    while(T--)
    {
        int n,m,p;
        scanf("%d%d%d" , &n , &m , &p);
        dlx.init(p , n*m);
        for(int i=1 ; i<=p ; i++){
            int x1,y1,x2,y2;
            scanf("%d%d%d%d" , &x1 , &y1 , &x2 , &y2);
            for(int j=x1+1 ; j<=x2 ; j++){
                for(int k=y1+1 ; k<=y2 ; k++){
                    dlx.link(i , (k-1)*n+j);
                }
            }
        }
        dlx.Dance(0);
        if(dlx.ans == INF) puts("-1");
        else printf("%d
" , dlx.ans);
    }
    return 0;
}