bzoj4411 [Usaco2016 Feb]Load balancing（线段树+二分+贪心）

首先我们枚举横着一刀切在哪，然后考虑此时竖着一刀切在哪里最好。我们令横线上方的点记为0,下方的点记为1，竖线左边的点记为l,右边的点记为r，则随着竖线的右移，显然有$max(l_0,l_1)$逐渐变大，$max(r_0,r_1)$逐渐变小。
记$ans=max(max(l_0,l_1),max(r_0,r_1))$，则我们一定有ans先减小后增大，是个单峰函数，我们可以在线段树上二分一下找到这个峰值点，更新答案。

复杂度$O（nlogn）$

upd:这个单峰函数的最值可能在两端，需要取min。但是大概因为枚举了横线，就算此次求出的竖线不是当前横线最优的，那最优的也一定会在枚举某个横线时出现。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return  EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,num[N*10],ans=inf,X0,X1;
struct node{
    int x0,x1;
}tr[N*10*4];
struct P{
    int x,y;
    friend bool operator<(P a,P b){return a.y<b.y;}
}a[N];
inline void pushup(int p){
    tr[p].x0=tr[p<<1].x0+tr[p<<1|1].x0;
    tr[p].x1=tr[p<<1].x1+tr[p<<1|1].x1;
}
inline void build(int p,int l,int r){
    if(l==r){tr[p].x0=num[l];tr[p].x1=0;return;}
    int mid=l+r>>1;build(p<<1,l,mid);build(p<<1|1,mid+1,r);pushup(p);
}
inline void change(int p,int l,int r,int x){
    if(l==r){tr[p].x0--;tr[p].x1++;return;}
    int mid=l+r>>1;
    if(x<=mid) change(p<<1,l,mid,x);
    else change(p<<1|1,mid+1,r,x);pushup(p);
}
inline int ask(int p,int l,int r,int x0,int x1){//找峰值 
    if(l==r) return l;
    int mid=l+r>>1;int sum0=x0+tr[p<<1].x0,sum1=x1+tr[p<<1].x1;
    if(max(sum0,sum1)>=max(tr[1].x0-sum0,tr[1].x1-sum1))
        return ask(p<<1,l,mid,x0,x1);
    return ask(p<<1|1,mid+1,r,sum0,sum1);
}
inline void qsum(int p,int l,int r,int x,int y){
    if(x>y) return;
    if(x<=l&&r<=y){X0+=tr[p].x0;X1+=tr[p].x1;return;}
    int mid=l+r>>1;
    if(x<=mid) qsum(p<<1,l,mid,x,y);
    if(y>mid) qsum(p<<1|1,mid+1,r,x,y);
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();
    for(int i=1;i<=n;++i) a[i].x=read(),a[i].y=read(),num[a[i].x]++;
    sort(a+1,a+n+1);build(1,1,1e6);
    for(int i=1,j;i<=n;i=j){//枚举横向切在a[i].y上面 
        j=i;
        while(j<=n&&a[j].y==a[i].y) change(1,1,1e6,a[j].x),++j;
        int pos=ask(1,1,1e6,0,0);
        X0=0;X1=0;qsum(1,1,1e6,1,pos);
        ans=min(ans,max(X0,X1));
        X0=0;X1=0;qsum(1,1,1e6,1,pos-1);
        ans=min(ans,max(tr[1].x0-X0,tr[1].x1-X1));
    }printf("%d\n",ans);
    return 0;
}