P1034 [NOIP2002 提高组] 矩形覆盖

题目链接

题意：给定n个点的坐标和一个正整数k，这些点可以被k个矩形覆盖，问这k个矩形面积之和最小为多少（注意覆盖一个点或覆盖平行于坐标轴的直线上的点的矩形面积为0，并且矩形之间必须完全分开）

题解

由于k最多取到4，所以可以直接暴搜，这里dfs的思路是考虑每一个点，枚举所有矩形，分别计算将该点加入后的情况，获取最优解。

其实dfs的主体结构很简单，包括剪枝，即如果当前面积大于当前答案，直接返回；更新答案，即判断矩形不相交后，更新当前答案为当前面积；搜索，即枚举矩形，然后计算当前点在枚举到的矩形的情况，再继续搜索，注意搜完要回溯。

然后就是一些关于矩形的操作了。
首先是关于矩形的信息，我是记录了左上和右下的两点，也是比较常规的方法。
然后是判断矩形是否相交，这里没考虑两个矩形形成类似“十”的形状的情况，因为这种情况不会更优，所以只需要枚举矩形的四个顶点，只要有一个点在另一个矩形中，那就说明相交。
最后是最重要的加入点并计算增加面积的函数，那自然需要两个变量，即当前加入的点和当前矩形，这里将矩形的初始值设为-1，由于点的坐标都是非负数，所以当当前加入的矩形的初始值为-1时，说明该矩形并未加入任何点，那么直接修改矩形信息为当前加入的点并返回0即可，否则的话就需要判断一下当前加入的点和当前矩形的顶点关系，如果在矩形外围的话，那为了覆盖到这个点，我们需要将矩形信息修改到刚好覆盖该点（即该点在矩形的边上），这样能保证面积最小。然后再计算加入后和加入前矩形面积的差值并返回即可。

代码实现

#include <algorithm>
#include <iostream>
#include <stdio.h>
#include <iomanip>
#include <cstring>
#include <string>
#include <vector>
#include <queue>
#include <set>
#include <map>
#define ll long long
using namespace std;
int n, k;
int x[60], y[60];
int res = 1e9;
struct node{
    int lx, ly, rx, ry;
}mtr[5];
int add(int now, int last){         //加入点并返回增加面积
    if(mtr[last].lx == -1){         //如果当前矩形没有点，直接修改为当前加入点
        mtr[last].lx = mtr[last].rx = x[now];
        mtr[last].ly = mtr[last].ry = y[now];
        return 0;
    }
    int s1 = (mtr[last].rx - mtr[last].lx) * (mtr[last].ly - mtr[last].ry);  //否则扩展矩形到加入点
    if(mtr[last].lx > x[now]) mtr[last].lx = x[now];
    if(mtr[last].rx < x[now]) mtr[last].rx = x[now];
    if(mtr[last].ly < y[now]) mtr[last].ly = y[now];
    if(mtr[last].ry > y[now]) mtr[last].ry = y[now];
    int s2 = (mtr[last].rx - mtr[last].lx) * (mtr[last].ly - mtr[last].ry);
    return s2 - s1;                      //返回增加面积
}
bool rect(int id, int x, int y){         //判断点是否在矩形内
    if(mtr[id].ry <= y && y <= mtr[id].ly && mtr[id].lx <= x && x <= mtr[id].rx) return true;
    else return false;
}
bool check(){                           //判断是否有矩形相交
    for(int i = 1; i <= k; i++){
        for(int j = i + 1; j <= k; j++){
            if(rect(i, mtr[j].lx, mtr[j].ly)) return true;
            if(rect(i, mtr[j].lx, mtr[j].ry)) return true;
            if(rect(i, mtr[j].rx, mtr[j].ly)) return true;
            if(rect(i, mtr[j].rx, mtr[j].ry)) return true;
        }
    }
    return false;
}
void dfs(int now, int sum){         //搜索
    if(sum >= res) return;
    if(now == n + 1){
        if(!check()) res = min(res, sum);
        return;
    }
    for(int i = 1; i <= k; i++){
        node tmp = mtr[i];
        dfs(now + 1, sum + add(now, i));
        mtr[i] = tmp;
    }
}
int main(){
    cin >> n >> k;
    for(int i = 1; i <= n; i++) cin >> x[i] >> y[i];
    for(int i = 1; i <= k; i++) mtr[i].lx = mtr[i].ly = mtr[i].rx = mtr[i].ry = -1;  //初始化
    dfs(1, 0);
    cout << res;
    return 0;
}