[bzoj4411] [Usaco2016 Feb]Load balancing

　　先离散化一下（也可以不用

　　枚举横坐标，用线段树维护两边纵坐标上的节点数。

　　每次在线段树上二分。。。（感觉似乎树状数组也行？

 #include<cstdio>
 #include<iostream>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 const int maxn=;
 const int mxnode=;
 struct poi{
     int x,y;
 }a[maxn];
 struct yy{
     int y,id;
 }y[maxn];int cnty;
 int sz[][mxnode],lnum[],rnum[];
 int i,j,k,n,m,ans,size;

 int ra;char rx;
 inline int read(){
     rx=getchar(),ra=;
     while(rx<''||rx>'')rx=getchar();
     while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
 }
 bool cmpy(yy a,yy b){return a.y<b.y;}
 bool cmppoi(poi a,poi b){return a.x<b.x;}
 inline int min(int a,int b){return a<b?a:b;}
 inline int max(int a,int b){return a>b?a:b;}
 inline int query(){
     int p=,lc,rc;bool flag=;
     lnum[]=lnum[]=rnum[]=rnum[]=;
     while(){
         lc=p<<,rc=lc|;
         if(p>size)lc=rc=p,flag=;
         if(max(lnum[]+sz[][lc],lnum[]+sz[][lc]) < max(rnum[]+sz[][rc],rnum[]+sz[][rc]))
             lnum[]+=sz[][lc],lnum[]+=sz[][lc],p=rc;
         else rnum[]+=sz[][rc],rnum[]+=sz[][rc],p=lc;
         //printf("  %d %d  %d %d  %d\n",lnum[0],lnum[1],rnum[0],rnum[1],p);
         if(flag)
             return max(max(lnum[],lnum[]),max(rnum[],rnum[]));
     }

 }
 int main(){
     n=read();
     for(i=;i<=n;i++)a[i].x=read(),a[i].y=y[i].y=read(),y[i].id=i;
     sort(y+,y++n,cmpy);
     for(i=;i<=n;i++)
         cnty+=y[i].y!=y[i-].y,a[y[i].id].y=cnty;
     for(size=;size<cnty;size<<=);size--;//printf("  %d\n",size);
     sort(a+,a++n,cmppoi);

     for(i=;i<=n;i++)
         ++sz[][a[i].y+size];
     for(i=size;i;i--)sz[][i]=sz[][i<<]+sz[][i<<|];

     int now=;ans=query();
     for(i=;i<=n;i++)if(a[i].x!=a[i+].x)
     {
     //  printf("try:%d\n",a[i].x);
         for(;now<=i;now++)
             for(j=a[now].y+size;j;j>>=)++sz[][j],--sz[][j];
         ans=min(ans,query());
     }
     printf("%d\n",ans);
     return ;
 }