[FJOI2015]火星商店问题(分治+可持久化)

题目描述

火星上的一条商业街里按照商店的编号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;
}