bitset做多维偏序

很久以前就听说这个大名鼎鼎的东西（暴力）了，现在才去写。其实很简单，对每一维排序，从左到右扫，维护一个bitset表示这位之前的某个点是否出现（即出现则该位为1否则0），查询某个数时则在每一维的排序完数组二分找到最右边<=这个数这一维权值的位置，取出该位的bitset，将每一维的bitset与起来即可得到答案，每次查询n/w，但这样空间要维数*n^2/w，因此可以分块增加查找时的一些时间换空间。

链接：五维偏序

http://hihocoder.com/problemset/problem/1236

#include<bits/stdc++.h>
#define Bst bitset<50010>
using namespace std;
const int N=50010,bsz=250;

template<class T>
void rd(T &x)
{
    char c=getchar();x=0;bool f=0;
    while(!isdigit(c))f|=(c=='-'),c=getchar();
    while(isdigit(c))x=x*10+c-48,c=getchar();
    if(f)x=-x;
}

struct gg{
    int a[5],id;
    void in()
    {rd(a[0]),rd(a[1]),rd(a[2]),rd(a[3]),rd(a[4]);}
}nd[N];
int n,m,bg[260],blk,cmp_id,mp[5][N],wh[5][N],lm[260],rm[260],bl[N];
Bst pre[5][260],bst;

bool cmp(gg x,gg y)
{return x.a[cmp_id]<y.a[cmp_id];}

Bst ask(int id,int x)
{
    static Bst res;res.reset();
    int ps=upper_bound(mp[id]+1,mp[id]+n+1,x)-mp[id],tmp;
    while(mp[id][ps]>x||ps>n)--ps;
    if(!ps)return res;
    tmp=bl[ps]-1;
    if(tmp!=0)res|=pre[id][tmp];
    for(int i=rm[tmp]+1;i<=ps;i++)res.set(wh[id][i]);
    return res;
}

int main()
{
    int T,ans,Q;
    rd(T);
    while(T--)
    {
        rd(n),rd(m),ans=0;
        for(int i=1;i<=n;i++)
            nd[i].in(),nd[i].id=i;
        for(int id=0;id<5;id++)
        {
            cmp_id=id;
            sort(nd+1,nd+n+1,cmp);
            for(int i=1;i<=n;i++)
                mp[id][i]=nd[i].a[id],wh[id][i]=nd[i].id;
            bst.reset();
            for(int i=1,j=1;i<=n;i+=bsz,j++)
            {
                for(int k=0;k<bsz&&i+k<=n;k++)
                    bst.set(wh[id][i+k]),bl[i+k]=j;
                pre[id][j]=bst,lm[j]=i,rm[j]=min(i+bsz-1,n);
            }
        }
        gg tmp;int ps;
        rd(Q);
        while(Q--)
        {
            tmp.in();
            bst.set();
            for(int i=0;i<5;i++)
                bst&=ask(i,tmp.a[i]^ans);
            printf("%d\n",ans=bst.count());
        }
    }
}