BZOJ4411: [Usaco2016 Feb]Load balancing

题目大意：给定n个点，请你确定一条横线和一条竖线，使得平面划分成的四个区域内点最大值最小

首先枚举横着的那条线，然后考虑在这条线上方随着竖线的平移，上面两个区域的max函数是一个单峰函数，下面同理，所以将两个函数取max之后还是一个单峰函数，那么我们就可以用二分或者三分来求这个峰，这个答案就是当横线在这里的最小值

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 100010
using namespace std;
struct ppp{int x,y;}a[N];
bool cmp(ppp x,ppp y){return x.x<y.x;}
int c[2][N*10],nn;
void change(int k,int x,int v)
{
	for(;x<=nn;x+=x&-x)
	c[k][x]+=v;
}
int check(int k,int x)
{
	int ret=0;
	for(;x;x-=x&-x)
	ret+=c[k][x];
	return ret;
}
int main()
{
	int n;
	scanf("%d",&n);
	int i,j,x,y;
	for(i=1;i<=n;i++)
	{
		scanf("%d%d",&a[i].x,&a[i].y);
		nn=max(nn,a[i].y);
	}
	for(i=1;i<=n;i++)
	change(1,a[i].y,1);
	sort(a+1,a+n+1,cmp);
	int l,r,mid;
	int ans=707185547;
	for(i=1;i<=n;i=j)
	{
		j=i;
		while(a[j].x==a[i].x&&j<=n)
		{
			change(1,a[j].y,-1);
			change(0,a[j].y,1);
			j++;
		}
		l=1;r=nn;
		int t1,t2;
		int tt1=check(0,nn),tt2=check(1,nn);
		while(l<r)
		{
			mid=(l+r)>>1;
			t1=check(0,mid),t2=check(1,mid);
			if(max(t1,t2)>=max(tt1-t1,tt2-t2))
			r=mid;
			else l=mid+1;
		}
		t1=check(0,l),t2=check(1,l);
		ans=min(ans,max(max(t1,t2),max(tt1-t1,tt2-t2)));
	}
	printf("%d",ans);
}