zoj 3209 Treasure Map(精确覆盖)

#include <iostream>
#include <string>
#include <stdio.h>
#include <cmath>
#include <algorithm>
#include <cstring>
#define INF 1000000000
#define MAXN 500005
using namespace std;

int n,m,p,tot,col[MAXN],u[MAXN],d[MAXN],l[MAXN],r[MAXN],sum[1005],ans;

int add(int c)
{
    ++tot;
    u[tot]=d[tot]=l[tot]=r[tot]=tot;
    col[tot]=c;
    sum[c]++;
    
    return tot;
}

void linklr(int head,int x)
{
    l[x]=l[head];
    r[l[head]]=x;
    l[head]=x;
    r[x]=head;
}
void linkud(int head,int x)
{
    u[x]=u[head];
    d[u[head]]=x;
    u[head]=x;
    d[x]=head;
}

void del(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];
         sum[col[j]]--;    
     }
}

void res(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]]=j;
         d[u[j]]=j;
         sum[col[j]]++;    
     }
     l[r[c]]=c;
     r[l[c]]=c;
}

void dfs(int num)
{
    int c,i,j;
    
    if (!r[0])//如果列为空的话
    {
         ans=min(ans,num);//有很多种解决方案 选中其中步数最小的
         return ;
    }
    c=r[0];//得到第一列
    for (i=r[0];i!=0;i=r[i]) //枚举还存在的每一列 选中该列相面相同数量最少的进行操作
        if (sum[i]<sum[c]) c=i;
    del(c);//先删除该列所影响的
    for (i=d[c];i!=c;i=d[i])//枚举该列下面所有存在的行
    {
        for (j=r[i];j!=i;j=r[j])//枚举相同行 中的所有列 删除
        {
             del(col[j]);
        }
        dfs(num+1) ;//传递删除了几行 当前加1
        for (j=l[i];j!=i;j=l[j]) res(col[j]);//然后恢复该行的
    }
    res(c);    //恢复该列
}

int main()
{
    int t0,x1,y1,x2,y2,cur,h,fst;
    
    scanf("%d",&t0);
    while (t0--)
    {
        scanf("%d%d%d",&n,&m,&p);//n*m大小的矩形  p个矩形块
        ans=INF;
        tot=-1;
        memset(sum,0,sizeof(sum));
        for (int t=0;t<n*m+1;t++)//枚举n*m个单位 主要是所有列名建立连接r,l,u,d(前n*m都是列)
        {
             add(t);
             linklr(0,t);//r[0] 是所有列的头部指向
        }
        
        for (int t=0;t<p;t++)//为n个块做连接
        {
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);//得到左上和右下角的点
            fst=0;
            for (int i=x1;i<x2;i++)//遍历这块面积中的每一个点 建立连接
             for (int j=y1;j<y2;j++)
             {
                 cur=add(i*m+j+1); //在该列下面重复添加一次  即该数量sum记录了列下面相同的个数和自己
                 if (fst++) linklr(h,cur); else h=cur;//如5→7→9
                 linkud(i*m+j+1,cur);//列 作为头部连接该点
             }
        }
        dfs(0);//删除了0行
        if (ans!=INF) printf("%d\n",ans);else printf("-1\n");
    }    
    return 0;
}