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