UVALive 7141 BombX

离散化，线段树。$2014$年$ACM/ICPC$亚洲区域赛上海站$D$题。

可以处理出炸任意相邻的$h$行能消灭的点的数量，以及炸任意相邻的$w$列能消灭的点的数量，分别用$py[i]$和$px[i]$记。

然后可以枚举炸哪个相邻的$h$行，这相邻的$h$行中有些位置可能有点在，所以有一些位置的$px$值是不可取的，要将这些$px$删去之后找一个最大值$m$，利用$m+py[i]$更新答案。一个点会导致一段连续的$px$不可取，所以可以用线段树区间更新将这一段区间的值减去一个比$n$大的值，使得这一段区间的值都是负的。

#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
#include<iostream>
using namespace std;

const int maxn=3e5+;
int n,w,h;
struct Point
{
    int x,y;
    void read() { scanf("%d%d",&x,&y); }
}p[maxn];
vector<int>x,y,in[maxn],out[maxn];
int px[maxn],py[maxn],W,H,ans;

long long s[maxn*],f[maxn*];

void build(int l,int r,int rt)
{
    f[rt]=;
    if(l==r) { s[rt]=px[l]; return ; }

    int m=(l+r)/;
    build(l,m,*rt);
    build(m+,r,*rt+);

    s[rt]=max(s[*rt],s[*rt+]);
}

void pushDown(int rt)
{
    if(f[rt]==) return ;

    s[*rt]+=f[rt], f[*rt]+=f[rt];
    s[*rt+]+=f[rt], f[*rt+]+=f[rt];
    f[rt]=;
}

void update(int L,int R,int val,int l,int r,int rt)
{
    if(L<=l&&r<=R) { s[rt]+=val, f[rt]+=val; return ; }

    int m=(l+r)/;
    pushDown(rt);

    if(L<=m) update(L,R,val,l,m,*rt);
    if(R>m) update(L,R,val,m+,r,*rt+);

    s[rt]=max(s[*rt],s[*rt+]);
}

int get(int g,bool flag)
{
    if(flag==) return lower_bound(x.begin(),x.end(),g)-x.begin()+;
    return lower_bound(y.begin(),y.end(),g)-y.begin()+;
}

int main()
{
    int T,cas=; scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&n,&w,&h);
        for(int i=;i<=n;i++) p[i].read();  x.clear(); y.clear();
        for(int i=;i<=n;i++)
        {
            y.push_back(p[i].y); y.push_back(p[i].y-h+); y.push_back(p[i].y+);
            x.push_back(p[i].x); x.push_back(p[i].x-w+); x.push_back(p[i].x+);
        }
        sort(x.begin(),x.end()); x.erase(unique(x.begin(),x.end()),x.end());
        sort(y.begin(),y.end()); y.erase(unique(y.begin(),y.end()),y.end());

        memset(px,,sizeof px); memset(py,,sizeof py);

        W=x.size(); H=y.size();
        for(int i=;i<=H;i++) { in[i].clear(); out[i].clear(); }

        for(int i=;i<=n;i++)
        {
            int xx=get(p[i].x-w+,), yy=get(p[i].y-h+,);
            int xxx=get(p[i].x+,), yyy=get(p[i].y+,);
            px[xx]++; px[xxx]--; py[yy]++; py[yyy]--;
            in[yy].push_back(i); out[yyy-].push_back(i);
        }

        ans=;
        for(int i=;i<=W;i++) px[i]=px[i]+px[i-], ans=max(ans,px[i]);
        for(int i=;i<=H;i++) py[i]=py[i]+py[i-], ans=max(ans,py[i]);

        build(,W,);

        for(int i=;i<=H;i++)
        {
            for(int j=;j<in[i].size();j++)
            {
                int ll=get(p[in[i][j]].x-w+,), rr=get(p[in[i][j]].x,);
                update(ll,rr,-n,,W,);
            }

            if(s[]>) ans=max(ans,py[i]+(int)s[]);

            for(int j=;j<out[i].size();j++)
            {
                int ll=get(p[out[i][j]].x-w+,), rr=get(p[out[i][j]].x,);
                update(ll,rr,n,,W,);
            }
        }
        printf("Case #%d: %d\n",cas++,ans);
    }
    return ;
}