《DFS》矩形覆盖

 

Problem F. 矩形覆盖

 

时间限制 1000 ms
内存限制 128 MB

题目描述

  在平面上有 n 个点(n < = 50),每个点用一对整数坐标表示。例如:当 n=4 时,4个点的坐标分另为:p1(1,1),p2(2,2),p3(3,6),P4(0,7)。见图一。

A_zjzj
  这些点可以用 k 个矩形(1< =k< =4)全部覆盖,矩形的边平行于坐标轴。当 k=2 时,可用如图二的两个矩形 sl,s2 覆盖,s1,s2 面积和为 4。问题是当 n 个点坐标和 k 给出后,怎样才能使得覆盖所有点的 k 个矩形的面积之和为最小呢。约定:覆盖一个点的矩形面积为 0;覆盖平行于坐标轴直线上点的矩形面积也为0。各个矩形必须完全分开(边线与顶点也都不能重合)。

输入数据

格式为
n k
xl y1
x2 y2
... ...
xn yn (0< =xi,yi< =500)

输出数据

一个整数,即满足条件的最小的矩形面积之和。

样例输入

4 2
1 1
2 2
3 6
0 7

样例输出

4

思路解析:

n<=50,k<=4, 额,暴力可以吗?妹妹。

对n个点进行dfs搜索,

难点有2:

1.如何判断此长方形和其他长方形是否相交。

2. 如何更新点。

对于第一个问题:

建立isOk函数: a>b  是x轴的坐标, d>c是y轴的坐标。

bool isOk(int index)
{
    for(int i=0;i<rtNum;i++)
    {
        if(index==i)
            continue;
        //矩形四个顶点(b,c) (b,d) (a,c) (a,d) 点只要有一个在其他矩形中就不可以。
        if(rt[index].b>=rt[i].b&&rt[index].b<=rt[i].a&&rt[index].c>=rt[i].c&&rt[index].c<=rt[i].d) return false;
        if(rt[index].b>=rt[i].b&&rt[index].b<=rt[i].a&&rt[index].d>=rt[i].c&&rt[index].d<=rt[i].d) return false;
        if(rt[index].a>=rt[i].b&&rt[index].a<=rt[i].a&&rt[index].c>=rt[i].c&&rt[index].c<=rt[i].d) return false;
        if(rt[index].a>=rt[i].b&&rt[index].a<=rt[i].a&&rt[index].d>=rt[i].c&&rt[index].d<=rt[i].d) return false;
    }
    return true;
}

对于第二个问题:

每个点可以放在任何一个已建立的长方形中,同时更新长方形

如果长方形数<给定的长方形数k,还可以新建立长方形,并更新

当然最后都要回溯。

void dfs(int index)
{
    int sum=0;
    for(int i=0;i<rtNum;i++)
        sum+=area(rt[i]);
    if(sum>=ans)//剪枝
        return;

    if(index==n)
    {
        ans=min(sum,ans);
        return;
    }

    for(int i=0;i<rtNum;i++)
    {
        // a,b,c,d  a>b d>c  a,b x轴方向,  c d y轴方向
        int tempA=rt[i].a;
        int tempB=rt[i].b;
        int tempC=rt[i].c;
        int tempD=rt[i].d;

        rt[i].a=max(p[index].x,rt[i].a);
        rt[i].b=min(p[index].x,rt[i].b);
        rt[i].c=min(p[index].y,rt[i].c);
        rt[i].d=max(p[index].y,rt[i].d);
        if(isOk(i))
            dfs(index+1);

        rt[i].a=tempA;//回溯
        rt[i].b=tempB;
        rt[i].c=tempC;
        rt[i].d=tempD;
    }
    if(rtNum<k)//增加矩阵放在后面是为了反正出现 如一行点,然后一个矩阵就可以解决了,而且面积为0,  先用少的矩阵尝试解决,再用多的矩阵解决,从而可以更好的覆盖所有情况
    {
        //初始化
        rt[rtNum].a=p[index].x;
        rt[rtNum].b=p[index].x;
        rt[rtNum].c=p[index].y;
        rt[rtNum].d=p[index].y;

        rtNum++;

        dfs(index+1);

        rtNum--;//回溯

    }
}

提示一点:

题目的意思是指k个长方体不一定用完。

AC代码:

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<string>
#include<stack>
#include<queue>
#include<set>
#include<map>
#include<cmath>
#define INF 0x7FFFFFFF //int 最大值
#define nMax 50
#define kMax 4
using namespace std;

int n,k;
int ans=INF;//答案
int rtNum=0;//答案矩阵数量

struct Point{
    int x,y;
};
Point p[nMax+5];


struct Rectangle{
    int a,b,c,d;// (a,c)右下顶点,(b,d)左上顶点。
};
Rectangle rt[kMax+5];

int area(Rectangle rect)
{
    return (rect.a-rect.b)*(rect.d-rect.c);
}

bool isOk(int index)
{
    for(int i=0;i<rtNum;i++)
    {
        if(index==i)
            continue;
        //矩形四个顶点(b,c) (b,d) (a,c) (a,d) 点只要有一个在其他矩形中就不可以。
        if(rt[index].b>=rt[i].b&&rt[index].b<=rt[i].a&&rt[index].c>=rt[i].c&&rt[index].c<=rt[i].d) return false;
        if(rt[index].b>=rt[i].b&&rt[index].b<=rt[i].a&&rt[index].d>=rt[i].c&&rt[index].d<=rt[i].d) return false;
        if(rt[index].a>=rt[i].b&&rt[index].a<=rt[i].a&&rt[index].c>=rt[i].c&&rt[index].c<=rt[i].d) return false;
        if(rt[index].a>=rt[i].b&&rt[index].a<=rt[i].a&&rt[index].d>=rt[i].c&&rt[index].d<=rt[i].d) return false;
    }
    return true;
}


void dfs(int index)
{
    int sum=0;
    for(int i=0;i<rtNum;i++)
        sum+=area(rt[i]);
    if(sum>=ans)//剪枝
        return;

    if(index==n)
    {
        ans=min(sum,ans);
        return;
    }

    for(int i=0;i<rtNum;i++)
    {
        // a,b,c,d  a>b d>c  a,b x轴方向,  c d y轴方向
        int tempA=rt[i].a;
        int tempB=rt[i].b;
        int tempC=rt[i].c;
        int tempD=rt[i].d;

        rt[i].a=max(p[index].x,rt[i].a);
        rt[i].b=min(p[index].x,rt[i].b);
        rt[i].c=min(p[index].y,rt[i].c);
        rt[i].d=max(p[index].y,rt[i].d);
        if(isOk(i))
            dfs(index+1);

        rt[i].a=tempA;//回溯
        rt[i].b=tempB;
        rt[i].c=tempC;
        rt[i].d=tempD;
    }
    if(rtNum<k)//增加矩阵放在后面是为了反正出现 如一行点,然后一个矩阵就可以解决了,而且面积为0,  先用少的矩阵尝试解决,再用多的矩阵解决,从而可以更好的覆盖所有情况
    {
        //初始化
        rt[rtNum].a=p[index].x;
        rt[rtNum].b=p[index].x;
        rt[rtNum].c=p[index].y;
        rt[rtNum].d=p[index].y;

        rtNum++;

        dfs(index+1);

        rtNum--;//回溯

    }
}

int main()
{
    ios::sync_with_stdio(false);
    cin>>n>>k;
    for(int i=0;i<n;i++)
        cin>>p[i].x>>p[i].y;

    ans=INF;
    rtNum=0;
    dfs(0);

    cout<<ans<<endl;

    return 0;
}

再提一点:

几何知识真的得加强了,其次很难想象,这种题竟然被出到贪心算法题里了,不懂0.0

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值