洛谷 P1034 矩形覆盖

最新推荐文章于 2024-09-15 10:34:21 发布

HJ921004

最新推荐文章于 2024-09-15 10:34:21 发布

阅读量273

点赞数

原文链接：http://www.cnblogs.com/cangT-Tlan/p/9911300.html

版权

题目描述

在平面上有 $\le 50n≤50)，每个点用一对整数坐标表示。例如：当 n=4n=4 时，44个点的坐标分另为：p_1p1(1,11,1)，p_2p2(2,22,2)，p_3p3(3,63,6)，P_4P4(0,70,7)，见图一。$

这些点可以用 $\le k \le 41≤k≤4)全部覆盖，矩形的边平行于坐标轴。当 k=2k=2 时，可用如图二的两个矩形 s_1,s_2s1,s2 覆盖，s_1,s_2s1,s2 面积和为44。问题是当nn个点坐标和kk给出后，怎样才能使得覆盖所有点的kk个矩形的面积之和为最小呢？约定：覆盖一个点的矩形面积为00；覆盖平行于坐标轴直线上点的矩形面积也为00。各个矩形必须完全分开（边线与顶点也都不能重合）。$

输入输出格式

输入格式：

$knkx_1 y_1x1y1x_2 y_2x2y2... ...$

$x_n y_nxnyn (0 \le x_i,y_i \le 5000≤xi,yi≤500)$

输出格式：

输出至屏幕。格式为：

输入输出样例

输入样例#1：复制

输出样例#1：复制

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int x[51],y[51];
int n,k,val,ans=0x7f7f7f7f;
struct nond{
    int l,r,u,d;
    bool flag;
}v[5];
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;    
}
bool jud(int i,int j){
    if(v[i].l<=v[j].l&&v[i].r>=v[j].l&&v[i].d>=v[j].d&&v[i].u<=v[j].d)    return true;
    if(v[i].l<=v[j].r&&v[i].r>=v[j].r&&v[i].d>=v[j].d&&v[i].u<=v[j].d)    return true;
    if(v[i].l<=v[j].l&&v[i].r>=v[j].l&&v[i].d>=v[j].u&&v[i].u<=v[j].u)    return true;
    if(v[i].l<=v[j].r&&v[i].r>=v[j].r&&v[i].d>=v[j].u&&v[i].u<=v[j].u)    return true;
    return false;
}
bool judge(){
    for(int i=1;i<=k;i++)
        if(v[i].flag)
        for(int j=1;j<i;j++)
            if(v[j].flag)
                if(jud(i,j))    return true;
    return false;
}
void dfs(int now){
    if(judge())    return ;
    val=0;
    for(int i=1;i<=k;i++)
        if(v[i].flag)
            val+=(v[i].r-v[i].l)*(v[i].d-v[i].u);
    if(val>ans)    return ;
    if(now==n+1){
        ans=val;
        return ;
    }
    for(int i=1;i<=k;i++)
        if(!v[i].flag){
            v[i].l=y[now];v[i].r=y[now];
            v[i].u=x[now];v[i].d=x[now];
            v[i].flag=1;
            dfs(now+1);
            v[i].flag=0;
        }
        else if(v[i].flag){
            int a=v[i].l,b=v[i].r,c=v[i].u,d=v[i].d;
            v[i].l=min(v[i].l,y[now]);
            v[i].r=max(v[i].r,y[now]);
            v[i].u=min(v[i].u,x[now]);
            v[i].d=max(v[i].d,x[now]);
            dfs(now+1);
            v[i].l=a;v[i].r=b;
            v[i].u=c;v[i].d=d;
        }
}
int main(){
    n=read();k=read();
    for(int i=1;i<=n;i++)
        scanf("%d%d",&x[i],&y[i]);
    dfs(1);
    cout<<ans;
}