《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）。见图一。

　　这些点可以用 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