BZOJ 1052 覆盖问题 （二分 贪心）

1052: [HAOI2007]覆盖问题

Description
某人在山上种了N棵小树苗。冬天来了，温度急速下降，小树苗脆弱得不堪一击，于是树主人想用一些塑料薄膜把这些小树遮盖起来，经过一番长久的思考，他决定用3个L*L的正方形塑料薄膜将小树遮起来。我们不妨将山建立一个平面直角坐标系，设第i棵小树的坐标为（Xi,Yi），3个L*L的正方形的边要求平行与坐标轴，一个点如果在正方形的边界上，也算作被覆盖。当然，我们希望塑料薄膜面积越小越好，即求L最小值。

Input
第一行有一个正整数N，表示有多少棵树。接下来有N行，第i+1行有2个整数Xi,Yi，表示第i棵树的坐标，保证不会有2个树的坐标相同。

Output
一行，输出最小的L值。

Sample Input
4

0 1

0 -1

1 0

-1 0

Sample Output
1

HINT
100%的数据，N<=20000

思路：
对于这种题目，很容易想到二分，首先二分mid之后就只需要处理一个判定问题。
check的方法是找到能恰好围住所有树苗的矩形，塑料薄膜一定会覆盖矩形的某一个角，那么枚举四个角；对于剩余部分，我们用同样的办法放置第二个薄膜；第三个薄膜直接判断即可。
下面稍微证明一下，为什么塑料薄膜一定覆盖矩形的某一个角呢？
既然这个矩形恰好能围住所有点，那么矩形的四条边上都会点（如果点在矩形的顶点上，那么一定会占据一角）；贪心的想，放在角上可以覆盖两条边上的点，而不放在角上，只能覆盖一条边上的点，我们有矩形的四条边，但是只有三个膜，所以至少有一个会在角上。

#include <iostream>  
#include <cstdio>  
#include <cstring>  
#include <algorithm>  
#define Inf 2147483647  
#define LL long long  
using namespace std;

int n, mid, flag, vis[20005];

struct point{
    int x, y;
}a[20005];

void border(int &rt, int &lf, int &dn, int &up){
    for(int i=1; i<=n; i++)  
        if ( !vis[i] ){
            flag = 0;
            lf = min(lf, a[i].x);  
            rt = max(rt, a[i].x);  
            dn = min(dn, a[i].y);  
            up = max(up, a[i].y);  
        }  
}

void cover(int x, int y){  
    for(int i=1; i<=n; i++)  
        if(a[i].x<=x && a[i].x>=x-mid && a[i].y>=y && a[i].y<=y+mid)  
            vis[i]++;
}

void recover(int x, int y){  
    for(int i=1; i<=n; i++)  
        if(a[i].x<=x && a[i].x>=x-mid && a[i].y>=y && a[i].y<=y+mid)  
            vis[i]--;
}

void dfs(int now){  
    if( flag ) return;
    int rt, lf, dn, up; //栈中要存储的元素只能是局部变量 
    rt=-Inf, lf=Inf, dn=Inf, up=-Inf; flag = 1;
    border(rt, lf, dn, up);  
    if(flag == 1) return; 
    if(now == 3){
        if( max(rt-lf, up-dn) <= mid ) flag = 1;  
        return;
    }
    cover(rt, dn); dfs(now+1); recover(rt, dn);  
    cover(lf+mid, dn); dfs(now+1); recover(lf+mid, dn);  
    cover(rt, up-mid); dfs(now+1); recover(rt, up-mid);  
    cover(lf+mid, up-mid); dfs(now+1); recover(lf+mid, up-mid);  
}

int main(){  
    scanf("%d", &n);
    for(int i=1; i<=n; i++)  
        scanf("%d%d", &a[i].x, &a[i].y);  
    int ans = 0, L=0, R=1e9;  
    while(L <= R){
        mid = (L+R)>>1;  
        flag = 0;  
        dfs( 1 );  
        if( flag ) ans = mid, R = mid-1;  
        else L = mid+1;
        memset(vis, 0, sizeof(vis));  
    }
    printf("%d\n", ans);  
    return 0;
}