题目描述
火星上的一条商业街里按照商店的编号1，2 ，…，n ，依次排列着n个商店。商店里出售的琳琅满目的商品中，每种商品都用一个非负整数val来标价。每个商店每天都有可能进一些新商品，其标价可能与已有商品相同。
火星人在这条商业街购物时，通常会逛这条商业街某一段路上的所有商店，譬如说商店编号在区间[L,R]中的商店，从中挑选1件自己最喜欢的商品。每个火星人对商品的喜好标准各不相同。通常每个火星人都有一个自己的喜好密码x。对每种标价为val的商品，喜好密码为x的火星人对这种商品的喜好程度与val异或x的值成正比。也就是说，val xor x的值越大，他就越喜欢该商品。每个火星人的购物卡在所有商店中只能购买最近d天内（含当天）进货的商品。另外，每个商店都有一种特殊商品不受进货日期限制，每位火星人在任何时刻都可以选择该特殊商品。每个商店中每种商品都能保证供应，不存在商品缺货的问题。
对于给定的按时间顺序排列的事件，计算每个购物的火星人的在本次购物活动中最喜欢的商品，即输出val xor x的最大值。这里所说的按时间顺序排列的事件是指以下2种事件：
事件0，用三个整数0,s,v，表示编号为s的商店在当日新进一种标价为v 的商品。
事件1，用5个整数1,L,R,x,d，表示一位火星人当日在编号为L到R的商店购买d天内的商品，该火星人的喜好密码为x。
题解
考虑如果没有时间限制，我们可以选取的区间为连续一段，那样就是简单的按位贪心，可以用可持久化trie维护。
现在有了时间限制，我们考虑分治，每个询问所覆盖的都是连续一段区间，而且修改之间是独立的，我们可以把每个询问拆成log个挂在按照时间建立的线段树上，然后按时间在线段树上分治。
我们可以每次都重构trie，复杂度是对的。
代码
#include<iostream>
#include<cstdio>
#include<vector>
#include<algorithm>
#define N 210002
using namespace std;
vector<int>vec[N<<2];
int len=18,tot,inv[20],T[N],ans[N],n,m,tim,top;
inline int rd(){
int x=0;char c=getchar();bool f=0;
while(!isdigit(c)){if(c==‘-‘)f=1;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
return f?-x:x;
}
struct TRIE{
int ch[N*22][2],size[N*22];
inline void ins(int &now,int pre,int x,int deep){
now=++tot;ch[now][0]=ch[pre][0];ch[now][1]=ch[pre][1];size[now]=size[pre]+1;
if(deep<0)return;
if(x&(1<<deep))ins(ch[now][1],ch[pre][1],x,deep-1);
else ins(ch[now][0],ch[pre][0],x,deep-1);
}
inline int query(int now,int pre,int x,int deep){
if(deep<0)return 0;
int o=(x&(1<<deep))!=0,num=size[ch[now][!o]]-size[ch[pre][!o]];
if(num)return inv[deep]+query(ch[now][!o],ch[pre][!o],x,deep-1);
else return query(ch[now][o],ch[pre][o],x,deep-1);
}
}tr;
struct node{
int pos,val,tim;
bool operator <(const node &b)const{return pos<b.pos;}
}co[N],q1[N],q2[N];
struct nod{int l,r,x,st,en;}q[N];
void calc(int cnt,int l,int r){
int top=0;tot=0;
for(int i=l;i<=r;++i){
++top;
tr.ins(T[top],T[top-1],co[i].val,len);
}
for(int i=0;i<vec[cnt].size();++i){
int id=vec[cnt][i];node x;
x.pos=q[id].l-1;
int L=upper_bound(co+l,co+r+1,x)-co-l;//!!!!
x.pos=q[id].r;
int R=upper_bound(co+l,co+r+1,x)-co-l;
ans[id]=max(ans[id],tr.query(T[R],T[L],q[id].x,len));
}
}
void upd(int cnt,int l,int r,int L,int R,int x){
if(L>R)return;
if(l>=L&&r<=R){vec[cnt].push_back(x);return;}
int mid=(l+r)>>1;
if(mid>=L)upd(cnt<<1,l,mid,L,R,x);
if(mid<R)upd(cnt<<1|1,mid+1,r,L,R,x);
}
void solve(int cnt,int l,int r,int L,int R){
if(l>r||L>R)return;
int mid=(l+r)>>1;
calc(cnt,L,R);
int o=0,p=0;
for(int i=L;i<=R;++i){
if(co[i].tim<=mid)q1[++o]=co[i];
else q2[++p]=co[i];
}
for(int i=1;i<=o;++i)co[L+i-1]=q1[i];
for(int i=1;i<=p;++i)co[L+o+i-1]=q2[i];
if(l!=r)solve(cnt<<1,l,mid,L,L+o-1);
solve(cnt<<1|1,mid+1,r,L+o,R);
}
int main(){
n=rd();m=rd();int l,r,x,d,v,opt;
for(int i=1;i<=n;++i)x=rd(),tr.ins(T[i],T[i-1],x,len);
inv[0]=1;
for(int i=1;i<=len;++i)inv[i]=inv[i-1]<<1;
for(int i=1;i<=m;++i){
opt=rd();
if(opt){
l=rd();r=rd();x=rd();d=rd();
q[++top]=nod{l,r,x,max(1,tim-d+1),tim};
ans[top]=tr.query(T[r],T[l-1],x,len);
}
else{
tim++;x=rd();v=rd();co[tim]=node{x,v,tim};
}
}
sort(co+1,co+tim+1);
for(int i=1;i<=top;++i)upd(1,1,tim,q[i].st,q[i].en,i);
solve(1,1,tim,1,tim);
for(int i=1;i<=top;++i)printf("%d\n",ans[i]);
return 0;
}
原文地址：https://www.cnblogs.com/ZH-comld/p/10146883.html