洛谷 P4585 [FJOI2015]火星商店问题

（勿看，仅作笔记）

bzoj权限题。。。

https://www.luogu.org/problemnew/show/P4585

对于特殊商品，直接可持久化trie处理一下即可

剩下的，想了一段时间cdq，但是没想出来。。。应该是不行的

事实上，如果询问的不是最大值，而是一些满足[l,r]的答案等于[1,r]的答案-[1,l-1]的答案的东西，那么的确可以每个询问拆成两个直接cdq。。。

但是这题就不能。。不过可以线段树分治，这是基于[l,r]的答案可以被分成多个线段树上区间（这些区间的并等于[l,r]）的答案的最大值（要对时间分治，[l,r]指对于某个询问合法的时间区间）

此题应该也可以线段树套可持久化trie做，但是最大的问题是垃圾回收。。。估计要写基于引用计数的垃圾回收？或者根本写不了？但是用分治的话只要每次solve之后把这次solve中用到的点全部回收掉就行了

跟题解学了个技巧：把询问和修改分开放进两个数组，方便处理

貌似那个vector也是可以去掉的，常数可以变小

使用vector是因为一个区间询问可能同时被分进两个子区间

那么只要处理出应该被分进左子区间的，然后solve左子区间，然后处理出应该被分进右子区间的，然后solve右子区间就行了；不一定要同时把应该分进左右子区间的处理出来

曾经错误：空间只算了后面solve部分所用字典树的空间，无视了特殊商品所占用的空间

 #include<cstdio>
 #include<algorithm>
 #include<vector>
 #include<cstring>
 #include<queue>
 #include<map>
 #define pb push_back
 using namespace std;
 namespace T
 {
 const int l2n=;
 int lft[];
 int sz[],ch[][];
 int x;
 int st[],top;
 void init()
 {
     int i;lft[]=;
     for(i=;i<=l2n;i++)  lft[i]=lft[i-]<<;
     for(i=;i<;i++)    st[++top]=i;
 }
 int getnode()
 {
     int t=st[top--];sz[t]=ch[t][]=ch[t][]=;
     return t;
 }
 void delnode(int x)    {st[++top]=x;}
 inline void cp(int &num)
 {
     int t=num;num=getnode();sz[num]=sz[t];ch[num][]=ch[t][];ch[num][]=ch[t][];
 }
 void _ins(int p,int &num)
 {
     cp(num);sz[num]++;
     if(p>=)    _ins(p-,ch[num][!!(x&lft[p])]);
 }
 void ins(int d,int &num)    {x=d;_ins(l2n-,num);}
 int que(int x,int r1,int r2)
 {
     int ans=,i;bool t;
     for(i=l2n-;i>=;i--)
     {
         t=x&lft[i];
         if(sz[ch[r2][!t]]-sz[ch[r1][!t]])   ans|=lft[i],r2=ch[r2][!t],r1=ch[r1][!t];
         else    r2=ch[r2][t],r1=ch[r1][t];
     }
     return ans;
 }
 }
 int n,m;
 int rt[];
 int ans[];
 struct Q
 {
     bool type;
     int L,R,x,l,r,num;
 };
 bool c1(const Q &a,const Q &b)    {return a.L<b.L;}
 vector<Q> qq0,qq1;
 void solve(const vector<Q> &q,const vector<Q> &c,int l,int r)
 {
     int i;
     int qz=q.size(),cz=c.size();
     map<int,int> rtt;rtt[]=;int rt1,rt2=,tmem=T::top;
     for(i=;i<cz;i++)
     {
         T::ins(c[i].x,rt2);rtt[c[i].L]=rt2;
     }
     for(i=;i<qz;i++)
     {
         if(q[i].l<=l&&r<=q[i].r)
         {
             rt1=(--rtt.upper_bound(q[i].L-))->second;
             rt2=(--rtt.upper_bound(q[i].R))->second;
             ans[q[i].num]=max(ans[q[i].num],T::que(q[i].x,rt1,rt2));
         }
     }
     T::top=tmem;
     if(l==r)    return;
     int mid=l+(r-l)/;vector<Q> q1,q2,c1,c2;
     for(i=;i<qz;i++)
     {
         if(!(q[i].l<=l&&r<=q[i].r))
         {
             if(q[i].l<=mid)    q1.pb(q[i]);
             if(mid<q[i].r)    q2.pb(q[i]);
         }
     }
     for(i=;i<cz;i++)
     {
         if(c[i].l<=mid)    c1.pb(c[i]);
         else    c2.pb(c[i]);
     }
     solve(q1,c1,l,mid);
     solve(q2,c2,mid+,r);
 }
 bool type[];
 int main()
 {
     int i,t,a,b,c,d,idx,dd=;
     T::init();
     scanf("%d%d",&n,&m);
     for(i=;i<=n;i++)
     {
         scanf("%d",&t);
         rt[i]=rt[i-];T::ins(t,rt[i]);
     }
     for(i=;i<=m;i++)
     {
         scanf("%d",&idx);
         if(idx==)
         {
             scanf("%d%d",&a,&b);
             qq1.pb((Q){,a,,b,++dd,,i});
         }
         else
         {
             type[i]=;
             scanf("%d%d%d%d",&a,&b,&c,&d);
             ans[i]=T::que(c,rt[a-],rt[b]);
             if(d!=)    qq0.pb((Q){,a,b,c,max(,dd-d+),dd,i});
         }
     }
     sort(qq1.begin(),qq1.end(),c1);
     solve(qq0,qq1,,dd);
     for(i=;i<=m;i++)
         if(type[i])
             printf("%d\n",ans[i]);
     return ;
 }