bzoj千题计划180：bzoj4411: [Usaco2016 Feb]Load balancing

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

用树状数组维护扫描线

一个树状数组维护扫描线之上的y<=i点，另一个维护扫描线之下y<=i的点

将点按x排好序，开始全部插入扫描线之下的树状数组

枚举x这一条线，线上的在第一个树状数组里加上，第二个树状数组里减去

最大值是一个单峰函数

可以用二分或三分

二分的话，哪边大往哪边移

#include<cstdio>
#include<iostream>
#include<algorithm>

using namespace std;

#define N 100001

#define lowbit(x) x&-x

int m;

int c[][N*];

struct node
{
    int x,y;
}e[N];

void read(int &x)
{
    x=; char c=getchar();
    while(!isdigit(c)) c=getchar();
    while(isdigit(c)) { x=x*+c-''; c=getchar(); }
}

bool cmp(node p,node q)
{
    return p.x<q.x;
}

void add(int t,int x,int w)
{
    while(x<=m)
    {
        c[t][x]+=w;
        x+=lowbit(x);
    }
}

int query(int t,int x)
{
    int sum=;
    while(x)
    {
        sum+=c[t][x];
        x-=lowbit(x);
    }
    return sum;
}

int main()
{
    int n;
    read(n);
    for(int i=;i<=n;++i)
    {
        read(e[i].x);
        read(e[i].y);
        m=max(m,e[i].y);
    }
    for(int i=;i<=n;++i) add(,e[i].y,);
    sort(e+,e+n+,cmp);
    int j;
    int ans=n;
    int l,r,mid,tmp;
    int sum0,sum1,tot0,tot1;
    for(int i=;i<=n;i=j)
    {
        j=i;
        while(j<=n && e[j].x==e[i].x)
        {
            add(,e[j].y,-);
            add(,e[j].y,);
            j++;
        }
        l=; r=m;
        tot0=query(,m);
        tot1=query(,m);
        tmp=;
        while(l<=r)
        {
            mid=l+r>>;
            sum0=query(,mid);
            sum1=query(,mid);
            if(max(sum0,sum1)>=max(tot0-sum0,tot1-sum1))
            {
                tmp=mid;
                r=mid-;
            }
            else l=mid+;
        }
        sum0=query(,tmp);
        sum1=query(,tmp);
        ans=min(ans,max(max(sum0,sum1),max(tot0-sum0,tot1-sum1)));
    }
    cout<<ans;
}