bzoj 1935 Tree 园丁的烦恼

题目大意：

一些点，每次查询一个矩形内有多少个点

思路：

因为空间太大

所以不能用什么二维树状数组

需要把这些点和所有查询的矩阵的左下和右上离线下来

先离散化

然后每个子矩阵像二维前缀和那样查询

按照x升序加入点，对于矩阵的点查询

 #include<iostream>
 #include<cstdio>
 #include<cmath>
 #include<cstdlib>
 #include<cstring>
 #include<algorithm>
 #include<vector>
 #include<queue>
 #define inf 2139062143
 #define ll long long
 #define MAXN 500100
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(!isdigit(ch)) {if(ch=='-') f=-;ch=getchar();}
     while(isdigit(ch)) {x=x*+ch-'';ch=getchar();}
     return x*f;
 }
 int n,c[MAXN*],vis[MAXN],m,cnt,k,ans[MAXN];
 struct data
 {
     int x,y,pos,yy;
 }g[MAXN*],a[MAXN*];
 bool cmp1(data a,data b)
 {
     if(a.x<b.x) return ;
     if(a.x>b.x) return ;
     if(a.x==b.x)
     {
         if(a.pos<=n&&b.pos>n) return ;
         if(a.pos>n&&b.pos<=n) return ;
         return a.y<b.y;
     }
 }
 bool cmp2(data a,data b) {return a.y<b.y||(a.y==b.y&&a.x<b.x);}
 int lowbit(int x) {return x&(-x);}
 void add(int x) {for(int i=x;i<=k;i+=lowbit(i)) c[i]++;}
 int query(int x) {int res=;for(int i=x;i;i-=lowbit(i)) res+=c[i];return res;}
 int main()
 {
     n=read(),m=read();
     for(int i=;i<=n;i++)
         a[i].pos=g[i].pos=i,g[i].x=read(),g[i].y=read();
     for(int i=n+;i<=n+*m;i++)
         a[i].pos=g[i].pos=i,g[i].x=read(),g[i].y=read();
     sort(g+,g+n+*m+,cmp2);
     g[].x=g[].y=-,cnt=;
     //离散化---------------------------------------------------------
     for(int i=;i<=n+*m;i++)
     {
         if(g[i].y!=g[i-].y) a[g[i].pos].y=++cnt;
         else a[g[i].pos].y=cnt;
     }
     k=cnt+;sort(g+,g+n+*m+,cmp1);cnt=;
     for(int i=;i<=n+*m;i++)
     {
         if(g[i].x!=g[i-].x) a[g[i].pos].x=++cnt;
         else a[g[i].pos].x=cnt;
     }
     //对于每个子矩阵的点需要记录另一个点的纵坐标---------------------
     for(int i=;i<=*m;i+=)
         {--a[i+n].x,a[i+n].yy=a[i++n].y,a[i++n].yy=a[i+n].y;}
     sort(a+,a+n+m*+,cmp1);
     for(int i=;i<=n+*m;i++)
     {
         if(a[i].pos>n)//如果是查询
         {
             int h=(a[i].pos-n+)>>;vis[h]++;
             if(vis[h]==) ans[h]=query(a[i].y-)-query(a[i].yy);//左下角的点
             if(vis[h]==) ans[h]+=query(a[i].y)-query(a[i].yy-);//右上角的点
         }
         else add(a[i].y);//不是查询的点
     }
     for(int i=;i<=m;i++) printf("%d\n",ans[i]);
 }