很久以前就听说这个大名鼎鼎的东西(暴力)了,现在才去写。其实很简单,对每一维排序,从左到右扫,维护一个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());
}
}
}