bzoj1935 [Shoi2007]园丁的烦恼

bzoj1935 [Shoi2007]园丁的烦恼

有N个点坐标为(xi,yi)，M次询问，询问(a,b)-(c,d)的矩形内有多少点。

0≤n≤500000，1≤m≤500000，0≤xi，yi≤10000000

看完数据范围傻眼系列。

做法：

离线处理

因为这个范围肯定不能把x y都离散，

所以只把点和询问的y坐标放在一起都离散化，然后一起按x坐标排序，

再开一个树状数组，求出每个需要的点的二维前缀和。

不理解的话拿画图板画一下就好了。

细节见代码

#include<bits/stdc++.h>
#define MAXN 500005
#define lowbit(a) (a&-a)
using namespace std;
inline int read(){
    int x=,t=;char c;c=getchar();
    while(c<''||c>''){if(c=='-')t=-;c=getchar();}
    while(c>=''&&c<=''){x=x*+c-'',c=getchar();}
    return x*t;
}
int N,M,Y[*MAXN],cnt,tx[MAXN],ty[MAXN],sum[MAXN][],
xa[MAXN],ya[MAXN],xb[MAXN],yb[MAXN],c[MAXN*],siz;
struct Work{
    int x,y,id,f;
    bool operator <(const Work z)const{
        return x<z.x||(x==z.x&&f<z.f);
    }
}a[*MAXN];
inline void add(int pos,int val){
    for(;pos<=siz;pos+=lowbit(pos))
        c[pos]+=val;
}
inline int query(int pos){
    int s=;
    for(;pos>;pos-=lowbit(pos))
        s+=c[pos];
    return s;
}
int main()
{
    N=read(),M=read();
    for(int i=;i<=N;i++){
        tx[i]=read()+,ty[i]=read();
        Y[++cnt]=ty[i];
    }
    for(int i=;i<=M;i++){
        xa[i]=read()+,ya[i]=read(),xb[i]=read()+,yb[i]=read();
        Y[++cnt]=ya[i],Y[++cnt]=yb[i];
    }
    sort(Y+,Y++cnt);
    siz=unique(Y+,Y++cnt)-Y-;    cnt=;
    for(int i=;i<=N;i++){
        ty[i]=lower_bound(Y+,Y++siz,ty[i])-Y;       //把点和询问放在一起离散化，然后排序。
        a[++cnt]=(Work){tx[i],ty[i],,};
    }
    for(int i=;i<=M;i++){
        ya[i]=lower_bound(Y+,Y++siz,ya[i])-Y;
        yb[i]=lower_bound(Y+,Y++siz,yb[i])-Y;
        a[++cnt]=(Work){xb[i],  yb[i],  i,};         //把需要计算矩阵和需要的四个点加进操作数组中，二位前缀和不用解释吧。
        a[++cnt]=(Work){xa[i]-,ya[i]-,i,};
        a[++cnt]=(Work){xa[i]-,yb[i],i,};
        a[++cnt]=(Work){xb[i],ya[i]-,i,};
    }
    sort(a+,a++cnt);
    for(int i=;i<=cnt;i++){
        if(a[i].id)　　　　　　　　　　　　　　　　　　　  　//遇到点把他加进树状数组
            sum[a[i].id][a[i].f]=query(a[i].y);
        else
            add(a[i].y,);
    }
    for(int i=;i<=M;i++)
        printf("%d\n",sum[i][]+sum[i][]-sum[i][]-sum[i][]);
    return ;
}

附个离散写法

sort -> unique -> lower_bound