首页 > 代码库 > 【BZOJ】1052: [HAOI2007]覆盖问题(贪心)

【BZOJ】1052: [HAOI2007]覆盖问题(贪心)

http://www.lydsy.com/JudgeOnline/problem.php?id=1052

首先膜拜题解orz,表示只能想到二分。。。

贪心就是每一次找到一个最小的能包围所有点的矩阵,然后枚举四个角放正方形。

不会证QAQ

然后被sb错吊打了半小时。。。

又是自以为的开全局变量没想到局部改变。。。

#include <cstdio>#include <cstring>#include <cmath>#include <string>#include <iostream>#include <algorithm>#include <queue>#include <set>#include <map>using namespace std;typedef long long ll;#define pii pair<int, int>#define mkpii make_pair<int, int>#define pdi pair<double, int>#define mkpdi make_pair<double, int>#define pli pair<ll, int>#define mkpli make_pair<ll, int>#define rep(i, n) for(int i=0; i<(n); ++i)#define for1(i,a,n) for(int i=(a);i<=(n);++i)#define for2(i,a,n) for(int i=(a);i<(n);++i)#define for3(i,a,n) for(int i=(a);i>=(n);--i)#define for4(i,a,n) for(int i=(a);i>(n);--i)#define CC(i,a) memset(i,a,sizeof(i))#define read(a) a=getint()#define print(a) printf("%d", a)#define dbg(x) cout << (#x) << " = " << (x) << endl#define error(x) (!(x)?puts("error"):0)#define printarr2(a, b, c) for1(_, 1, b) { for1(__, 1, c) cout << a[_][__]; cout << endl; }#define printarr1(a, b) for1(_, 1, b) cout << a[_] << ‘\t‘; cout << endlinline const int getint() { int r=0, k=1; char c=getchar(); for(; c<‘0‘||c>‘9‘; c=getchar()) if(c==‘-‘) k=-1; for(; c>=‘0‘&&c<=‘9‘; c=getchar()) r=r*10+c-‘0‘; return k*r; }inline const int max(const int &a, const int &b) { return a>b?a:b; }inline const int min(const int &a, const int &b) { return a<b?a:b; }const int N=40005, oo=~0u>>2;int n, vis[N];struct dat { int x, y; }a[N];void ifind(int &x1, int &y1, int &x2, int &y2) {	x1=y1=oo; x2=y2=-oo;	for1(i, 1, n) if(!vis[i]) x1=min(x1, a[i].x), x2=max(x2, a[i].x), y1=min(y1, a[i].y), y2=max(y2, a[i].y);}void fix(const int &x1, const int &y1, const int &x2, const int &y2, const int &dep) { for1(i, 1, n) if(!vis[i] && a[i].x>=x1 && a[i].x<=x2 && a[i].y>=y1 && a[i].y<=y2) vis[i]=dep; }bool find(int dep, const int &L) {	int x[2], y[2];	ifind(x[0], y[0], x[1], y[1]);	if(max(x[1]-x[0], y[1]-y[0])<=L) return true;	if(dep==3) return false;	rep(i, 2) rep(j, 2) {		if(i==0) {			if(j==0) fix(x[i], y[j], x[i]+L, y[j]+L, dep);			else fix(x[i], y[j]-L, x[i]+L, y[j], dep);		}		else {			if(j==0) fix(x[i]-L, y[j], x[i], y[j]+L, dep);			else fix(x[i]-L, y[j]-L, x[i], y[j], dep);		}		//printf("dep:%d, (%d,%d), (%d,%d)\n", dep, x[0], y[0], x[1], y[1]);		//for1(i, 1, n) printf("%d ", vis[i]); puts("");		if(find(dep+1, L)) return true;		for1(i, 1, n) if(vis[i]==dep) vis[i]=0;	}	return false;}bool check(int L) {	memset(vis, 0, sizeof(int)*(n+2)); //dbg(L);	return find(1, L);}int main() {	read(n);	for1(i, 1, n) read(a[i].x), read(a[i].y);	int x[2], y[2];	ifind(x[0], y[0], x[1], y[1]);	int l=1, r=max(x[1]-x[0], y[1]-y[0]);	while(l<=r) {		int mid=(l+r)>>1;		if(check(mid)) r=mid-1;		else l=mid+1;	}	print(r+1);	return 0;}

  

 


 

 

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

 

Source

【BZOJ】1052: [HAOI2007]覆盖问题(贪心)