【BZOJ1052】覆盖问题（贪心）

【BZOJ1052】覆盖问题（贪心）

题面

BZOJ
洛谷

题解

这题好神仙啊。
很明显可以看出来要二分一个边长。
那么如何\(check\)呢？
我们把所有点用一个最小矩形覆盖，
那么必定每个边界上都至少存在一个点，
但是我们有\(4\)个边界，但是只有\(3\)个矩形，
意味着至少有一个矩形卡住了两个边界，
那么我们递归处理，每次枚举卡在了哪个角上就好了。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define ll long long
#define MAX 20020
inline int read()
{
    int x=0;bool t=false;char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=true,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return t?-x:x;
}
int n;
struct Node{int x,y;}p[MAX];
int vis[MAX];
bool operator<(Node a,Node b)
{
    if(a.x!=b.x)return a.x<b.x;
    return a.y<b.y;
}
void cmin(int &x,int y){x=min(x,y);}
void cmax(int &x,int y){x=max(x,y);}
void Cover(int x1,int x2,int y1,int y2,int id)
{
    for(int i=1;i<=n;++i)
        if(!vis[i]&&x1<=p[i].x&&p[i].x<=x2&&y1<=p[i].y&&p[i].y<=y2)
            vis[i]=id;
}
void UnCover(int id){for(int i=1;i<=n;++i)if(vis[i]==id)vis[i]=0;}
bool dfs(int tot,int L)
{
    int x[2],y[2];x[0]=y[0]=2e9;x[1]=y[1]=-2e9;
    for(int i=1;i<=n;++i)
        if(!vis[i])cmin(x[0],p[i].x),cmax(x[1],p[i].x),cmin(y[0],p[i].y),cmax(y[1],p[i].y);
    if(max(x[1]-x[0],y[1]-y[0])<=L)return true;
    if(tot==3)return false;
    for(int i=0;i<2;++i)
        for(int j=0;j<2;++j)
        {
            if(i==0)
            {
                if(j==0)Cover(x[0],x[0]+L,y[0],y[0]+L,tot);
                else Cover(x[0],x[0]+L,y[1]-L,y[1],tot);
            }
            else
            {
                if(j==0)Cover(x[1]-L,x[1],y[0],y[0]+L,tot);
                else Cover(x[1]-L,x[1],y[1]-L,y[1],tot);
            }
            if(dfs(tot+1,L))return true;
            UnCover(tot);
        }
    return false;
}
bool check(int L)
{
    memset(vis,0,sizeof(vis));
    return dfs(1,L);
}
int main()
{
    n=read();
    for(int i=1;i<=n;++i)p[i].x=read(),p[i].y=read();
    int l=0,r=2e9,ret=2e9;
    while(l<=r)
    {
        int mid=((ll)l+r)/2;
        if(check(mid))r=mid-1,ret=mid;
        else l=mid+1;
    }
    printf("%d\n",ret);
}

转载于:https://www.cnblogs.com/cjyyb/p/9444110.html