BZOJ 1052 二分+贪心

想到了一部分，没想全，还是看了题解。。

也不知道为什么对。。。

首先，求出所有点的4个边界值形成的一个矩形，第一个正方形的一个边界一定与这个矩形的4个角中的一个重合，枚举4次即可，
然后再找到剩下的点中的边界，重复一遍上面的操作，最后判断一下一个正方形是否可以覆盖剩余的所有矩形

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>

#define N 22222
#define INF 2147483647

using namespace std;
/*
首先，求出所有点的4个边界值形成的一个矩形，第一个正方形的一个边界一定与这个矩形的4个角中的一个重合，枚举4次即可，
然后再找到剩下的点中的边界，重复一遍上面的操作，最后判断一下一个正方形是否可以覆盖剩余的所有矩形
*/
struct P
{
    int x,y;
}p[N],p1[N],p2[N];

int n;

inline void read()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d%d",&p[i].x,&p[i].y);
}

inline bool check(int len)
{
    if(n==0) return true;
    int sd[4][4];
    sd[1][1]=INF; sd[1][2]=-INF;
    sd[2][1]=INF; sd[2][2]=-INF;
    for(int i=1;i<=n;i++)
    {
        sd[1][1]=min(sd[1][1],p[i].x); sd[1][2]=max(sd[1][2],p[i].x);
        sd[2][1]=min(sd[2][1],p[i].y); sd[2][2]=max(sd[2][2],p[i].y);
    }
    for(int i=1;i<=2;i++)
        for(int j=1;j<=2;j++)
        {
            int m=0;
            for(int k=1;k<=n;k++)
                if(abs(p[k].x-sd[1][i])>len||abs(p[k].y-sd[2][j])>len) p1[++m]=p[k];
            if(m==0) return true;
            
            int sp[4][4];
            sp[1][1]=INF; sp[1][2]=-INF;
            sp[2][1]=INF; sp[2][2]=-INF;
            for(int k=1;k<=m;k++)
            {
                sp[1][1]=min(sp[1][1],p1[k].x); sp[1][2]=max(sp[1][2],p1[k].x);
                sp[2][1]=min(sp[2][1],p1[k].y); sp[2][2]=max(sp[2][2],p1[k].y);
            }
            for(int ii=1;ii<=2;ii++)
                for(int jj=1;jj<=2;jj++)
                {
                    int s=0;
                    for(int kk=1;kk<=m;kk++)
                        if(abs(p1[kk].x-sp[1][ii])>len||abs(p1[kk].y-sp[2][jj])>len) p2[++s]=p1[kk];
                    if(s==0) return true;
                    
                    int sq[4][4];
                    sq[1][1]=INF; sq[1][2]=-INF;
                    sq[2][1]=INF; sq[2][2]=-INF;
                    for(int kk=1;kk<=s;kk++)
                    {
                        sq[1][1]=min(sq[1][1],p2[kk].x); sq[1][2]=max(sq[1][2],p2[kk].x);
                        sq[2][1]=min(sq[2][1],p2[kk].y); sq[2][2]=max(sq[2][2],p2[kk].y);
                    }
                    if(sq[2][2]-sq[2][1]<=len&&sq[1][2]-sq[1][1]<=len) return true;
                }
        }
    return false;
}

inline void go()
{
    int l=0,r=2000000000,mid,ans;
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(check(mid)) ans=mid,r=mid-1;
        else l=mid+1;
    }
    printf("%d\n",ans);
}

int main()
{
    read(),go();
    return 0;
}