[noi1773]function

以统计x坐标的数量为例：x为下标建一棵线段树，然后对每一个区间按照y坐标建一棵可持久化线段树（每一个x只保留最大的一个y），询问时，二分找到这个区间内最大的y以前的点并统计，复杂度为$o(nlog^{2}n)$

还有一种做法是bitset+分块，预处理出：1.第i个块到第j个块的bitdet；2.每一个块的前缀的bitset；3.每一个块的后缀的bitset，根据这个就可以快速求出对应区间的bitset，然后通过左移和右移去除区间外的点，求和即可，时间复杂度$o(n^{2}/32+n\sqrt{n})$

#include<bits/stdc++.h>
using namespace std;
#define N 50005
#define y1 y11
#define T(p) p,1,1,m[p]
#define L (k<<1)
#define R (L+1)
#define mid (l+r>>1)
struct ji{
    int x,y;
    bool operator < (const ji &k)const{
        return y<k.y;
    }
}a[N];
struct tr{
    int k,ls,rs;
}f[N*1000];
vector<int>rt[2][N<<2]; 
vector<ji>v[2][N<<2];
int V,m[2],n,q,p,x1,y1,x2,y2,ans1,ans2,lst[N],X[N],Y[N];
bool cmp(ji x,ji y){
    return x.x<y.x;
}
void update(int &k,int l,int r,int x,int y){
    f[++V]=f[k];
    k=V;
    f[k].k+=y;
    if (l==r)return;
    if (x<=mid)update(f[k].ls,l,mid,x,y);
    else update(f[k].rs,mid+1,r,x,y);
}
int query(int k,int l,int r,int x,int y){
    if ((!k)||(l>y)||(x>r))return 0;
    if ((x<=l)&&(r<=y))return f[k].k;
    return query(f[k].ls,l,mid,x,y)+query(f[k].rs,mid+1,r,x,y);
}
void push(int p,int k,int l,int r,int x,int y){
    v[p][k].push_back(ji{x,y});
    if (l==r)return;
    if (x<=mid)push(p,L,l,mid,x,y);
    else push(p,R,mid+1,r,x,y);
}
void build(int p,int k,int l,int r){
    if (l<r){
        build(p,L,l,mid);
        build(p,R,mid+1,r);
    }
    sort(v[p][k].begin(),v[p][k].end());
    int las=0,root=0;
    for(int i=0;i<v[p][k].size();i++){
        update(root,1,m[p^1],v[p][k][i].y,1);
        if (lst[v[p][k][i].x])update(root,1,m[p^1],lst[v[p][k][i].x],-1);
        lst[v[p][k][i].x]=v[p][k][i].y;
        if ((i+1==v[p][k].size())||(v[p][k][i+1].y!=v[p][k][las].y)){
            while (las++<=i)rt[p][k].push_back(root);
            las--;
        }
    }
    for(int i=0;i<v[p][k].size();i++)lst[v[p][k][i].x]=0;
}
int query(int p,int k,int l,int r,int x1,int x2,int y1,int y2){
    if ((x1>r)||(l>x2))return 0;
    if ((x1<=l)&&(r<=x2)){
        int x=lower_bound(v[p][k].begin(),v[p][k].end(),ji{0,y2+1})-v[p][k].begin()-1;
        if (x==-1)return 0;
        return query(rt[p][k][x],1,m[p^1],y1,y2);
    }
    return query(p,L,l,mid,x1,x2,y1,y2)+query(p,R,mid+1,r,x1,x2,y1,y2);
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%d%d",&a[i].x,&a[i].y);
    sort(a+1,a+n+1,cmp);
    for(int i=1;i<=n;i++)X[i]=a[i].x;
    sort(a+1,a+n+1);
    for(int i=1;i<=n;i++)Y[i]=a[i].y;
    m[0]=unique(X+1,X+n+1)-X-1;
    m[1]=unique(Y+1,Y+n+1)-Y-1;
    for(int i=1;i<=n;i++){
        a[i].x=lower_bound(X+1,X+m[0]+1,a[i].x)-X;
        a[i].y=lower_bound(Y+1,Y+m[1]+1,a[i].y)-Y;
    }
    for(int i=1;i<=n;i++)
    scanf("%d%d",&q,&p);
    for(int i=1;i<=q;i++){
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        x1=lower_bound(X+1,X+m[0]+1,x1^(p*ans1))-X;
        y1=lower_bound(Y+1,Y+m[1]+1,y1^(p*ans2))-Y;
        x2=upper_bound(X+1,X+m[0]+1,x2^(p*ans1))-X-1;
        y2=upper_bound(Y+1,Y+m[1]+1,y2^(p*ans2))-Y-1;
    }
}