大致题意: 有(n)家店,每个商品有一个标价。每天,都可能有某家商店进货,也可能有某人去购物。一个人在购物时,会于编号在区间([L_i,R_i])的商店里挑选一件进货(d_i)天以内的商品使得其标价与(x_i)的异或值最大,同时每家店有一个特殊商品在任何时候都可以选择。对于每名购物者求出这个最大异或值。
暴力
一个很暴力的想法,就是线段树套可持久化(Trie)树......
但这貌似是过不了的,因此需要优化。
线段树分治
这道题的正解是一个叫做线段树分治的神奇算法。
考虑对于一次询问,它询问的是一段时间区间内的答案。
则我们可以考虑开一棵线段树,然后把每个询问都扔到一个以时间为下标的线段树内的对应区间上。
接下来,考虑离线,在线段树上,用类似于整体二分的思想,维护每一个节点所对应的这段时间内的修改操作。
然后可以参照前面暴力的想法,用一棵可持久化(Trie)树。对于每一个节点,我们先清空可持久化(Trie)树,再把这个节点的修改对应到可持久化(Trie)树上。随后我们枚举询问,求出每个询问所询问的区间后,就可以求出每个询问的答案了。
注意一个询问可能扔到线段树可能会被分成多个,因此答案要取(max)。
代码
#include<bits/stdc++.h>
#define Tp template<typename Ty>
#define Ts template<typename Ty,typename... Ar>
#define Reg register
#define RI Reg int
#define Con const
#define CI Con int&
#define I inline
#define W while
#define N 100000
#define LN 20
#define max(x,y) ((x)>(y)?(x):(y))
#define Gmax(x,y) (x<(y)&&(x=(y)))
using namespace std;
int n,Et,Qt,ans[N+5];
struct Event {int i,v,t;I bool operator < (Con Event& o) Con {return i<o.i;}}e[N+5];
struct Query {int l,r,v,x,y,t;}q[N+5];
class FastIO
{
private:
#define FS 100000
#define tc() (A==B&&(B=(A=FI)+fread(FI,1,FS,stdin),A==B)?EOF:*A++)
#define pc(c) (C==E&&(clear(),0),*C++=c)
#define tn (x<<3)+(x<<1)
#define D isdigit(c=tc())
int T;char c,*A,*B,*C,*E,FI[FS],FO[FS],S[FS];
public:
I FastIO() {A=B=FI,C=FO,E=FO+FS;}
Tp I void read(Ty& x) {x=0;W(!D);W(x=tn+(c&15),D);}
Tp I void write(Ty x) {W(S[++T]=x%10+48,x/=10);W(T) pc(S[T--]);}
Tp I void writeln(Con Ty& x) {write(x),pc('
');}
I void clear() {fwrite(FO,1,C-FO,stdout),C=FO;}
}F;
class ChairmanTrie//可持久化Trie
{
private:
#define ST(x) O[rt1].S[x],O[rt2].S[x],v,d-1
int Nt,Vt,Rt[N+5];struct node {int V,S[2];I void Clear() {V=S[0]=S[1]=0;}}O[N*LN<<2];
I void Ins(CI rt1,int& rt2,CI v,CI d)//插入新数
{
O[rt2=++Nt]=O[rt1],++O[rt2].V;if(!~d) return;int t=(v>>d)&1;Ins(ST(t));
}
I int Qry(CI rt1,CI rt2,CI v,CI d)//求最大异或值
{
if(!rt2||!~d) return 0;int t=(v>>d)&1;
return O[O[rt1].S[t^1]].V^O[O[rt2].S[t^1]].V?Qry(ST(t^1))|(1<<d):Qry(ST(t));
}
public:
I void Insert(CI x) {++Vt,Ins(Rt[Vt-1],Rt[Vt],x,LN);}
I int Query(CI l,CI r,CI x) {return Qry(Rt[l-1],Rt[r],x,LN);}
I void Clear() {W(Nt) Rt[Nt--]=0;W(Vt) O[Vt--].Clear();}//清空
}C;
class SegmentTree
{
private:
#define PT CI l=1,CI r=Et,CI rt=1
#define LT l,mid,rt<<1
#define RT mid+1,r,rt<<1|1
#define pb push_back
int s[N+5];Event tl[N+5],tr[N+5];vector<int> P[N<<2];
public:
I void Ins(CI x,CI y,CI p,PT)//扔入询问
{
if(x<=l&&r<=y) return P[rt].pb(p);int mid=l+r>>1;//开vector存储一个节点上的询问
x<=mid&&(Ins(x,y,p,LT),0),y>mid&&(Ins(x,y,p,RT),0);
}
I void Solve(CI x,CI y,PT)//二分求答案
{
if(x>y) return;RI i,t,t1,t2,k=0;int mid=l+r>>1;vector<int>::iterator it;
for(C.Clear(),i=x;i<=y;++i) s[++k]=e[i].i,C.Insert(e[i].v);//将修改对应到可持久化Trie树上
for(it=P[rt].begin();it!=P[rt].end();++it)//枚举这个节点上的询问
t1=upper_bound(s+1,s+k+1,q[*it].l-1)-s,t2=upper_bound(s+1,s+k+1,q[*it].r)-s-1,//求出询问在可持久化Trie树上对应的区间
t=C.Query(t1,t2,q[*it].v),Gmax(ans[q[*it].t],t);//在可持久化Trie树上询问,更新答案
if(l==r) return;for(t1=t2=0,i=x;i<=y;++i) e[i].t<=mid?tl[++t1]=e[i]:tr[++t2]=e[i];
for(i=1;i<=t1;++i) e[x+i-1]=tl[i];for(i=1;i<=t2;++i) e[x+i+t1-1]=tr[i];
Solve(x,x+t1-1,LT),Solve(x+t1,y,RT);//递归操作
}
}S;
int main()
{
RI m,i,op,x;for(F.read(n),F.read(m),i=1;i<=n;++i) F.read(x),C.Insert(x);//对于特殊商品
W(m--) F.read(op),!op?(F.read(e[++Et].i),F.read(e[Et].v),e[Et].t=Et)://读入修改
(
F.read(q[++Qt].l),F.read(q[Qt].r),F.read(q[Qt].v),F.read(x),q[Qt].x=max(Et-x,0)+1,q[Qt].y=Et,//读入询问
ans[Qt]=C.Query(q[Qt].l,q[Qt].r,q[Qt].v),q[Qt].t=Qt//因为特殊商品随时都可以选,因此初始化答案
);
for(i=1;i<=Qt;++i) q[i].x<=q[i].y&&(S.Ins(q[i].x,q[i].y,i),0);sort(e+1,e+Et+1),S.Solve(1,Et);//把询问扔入线段树,离线二分求答案
for(i=1;i<=Qt;++i) F.writeln(ans[i]);return F.clear(),0;//输出答案
}