[FJOI2015]火星商店问题 --线段树分治+可持久化trie

题意：https://www.luogu.com.cn/problem/P4585

询问有两种限制

1：商店区间限制

2：时间区间限制

思路：

像这样两个限制我们考虑线段树二分（应该是这么考虑吧。。。

线段树二分我的理解其实就是整体二分只不过在线段树上

线段树二分有时候也利用dfs的顺序解题，比如LOJ121（我也就做了这几题=-=

我们把时间放到线段树上就不需要考虑时间了

试问：不需要考虑时间，给你一个询问还不好做吗，二分啊，那n个询问就是整体二分

就没了。。

#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
#include <cstdio>//sprintf islower isupper
#include <cstdlib>//malloc  exit strcat itoa system("cls")
#include <iostream>//pair
#include <fstream>//freopen("C:\Users\13606\Desktop\Input.txt","r",stdin);
#include <bitset>
//#include <map>
//#include<unordered_map>
#include <vector>
#include <stack>
#include <set>
#include <string.h>//strstr substr strcat
#include <string>
#include <time.h>// srand(((unsigned)time(NULL))); Seed n=rand()%10 - 0~9;
#include <cmath>
#include <deque>
#include <queue>//priority_queue<int, vector<int>, greater<int> > q;//less
#include <vector>//emplace_back
//#include <math.h>
#include <cassert>
#include <iomanip>
//#include <windows.h>//reverse(a,a+len);// ~ ! ~ ! floor
#include <algorithm>//sort + unique : sz=unique(b+1,b+n+1)-(b+1);+nth_element(first, nth, last, compare)
using namespace std;//next_permutation(a+1,a+1+n);//prev_permutation
//******************
clock_t __START,__END;
double __TOTALTIME;
void _MS(){__START=clock();}
void _ME(){__END=clock();__TOTALTIME=(double)(__END-__START)/CLOCKS_PER_SEC;cout<<"Time: "<<__TOTALTIME<<" s"<<endl;}
//***********************
#define rint register int
#define fo(a,b,c) for(rint a=b;a<=c;++a)
#define fr(a,b,c) for(rint a=b;a>=c;--a)
#define mem(a,b) memset(a,b,sizeof(a))
#define pr printf
#define sc scanf
#define ls rt<<1
#define rs rt<<1|1
typedef pair<int,int> PII;
typedef vector<int> VI;
typedef unsigned long long ull;
typedef long long ll;
typedef double db;
const db E=2.718281828;
const db PI=acos(-1.0);
const ll INF=(1LL<<60);
const int inf=(1<<30);
const db ESP=1e-9;
const int mod=(int)1e9+7;
const int N=(int)1e5+10;

int ans[N];

class TRIE
{
public:
    int tot;
    int root[N*20];
    int tr[N*20][2];
    int sz[N*20];
    void Init()
    {
        tot=0;
        sz[0]=root[0]=0;
    }
    void insert(int &x,int u,int val)
    {
        int now;now=x=++tot;
        for(int i=17;i>=0;i--)
        {
            int d=(val&(1<<i))!=0;
            tr[now][d^1]=tr[u][d^1];tr[now][d]=++tot;
            now=tr[now][d];u=tr[u][d];
            sz[now]=sz[u]+1;
        }
    }
    int query(int l,int r,int val)
    {
        int res=0;
        for(int i=17;i>=0;i--)
        {
            int d=(val&(1<<i))!=0;
            if(sz[tr[r][d^1]]-sz[tr[l][d^1]]>0)
                l=tr[l][d^1],r=tr[r][d^1],res+=(1<<i);
            else l=tr[l][d],r=tr[r][d];
        }
        return res;
    }
}TR;

vector<vector<int> >Seg(N*4);
void update(int L,int R,int val,int l,int r,int rt)//查询区间和
{
    if(L>R)return;
    if(L<=l&&r<=R)
    {
        Seg[rt].push_back(val);
        return;
    }
    int mid=(l+r)>>1;
    if(L<=mid)update(L,R,val,l,mid,rt<<1);
    if(R> mid)update(L,R,val,mid+1,r,rt<<1|1);
}
struct node
{
    int day,shop,v;
    friend bool operator<(node a,node b)
    {
        return a.shop<b.shop;
    }
}goods[N];int cg;
struct node0{int l,r,x,L,R;}q[N];int cq;

int temp[N];
void cal(int rt,int gol,int gor)
{
    int cnt=0;
    TR.Init();
    for(int i=gol;i<=gor;++i)temp[++cnt]=goods[i].shop,TR.insert(TR.root[cnt],TR.root[cnt-1],goods[i].v);
    for(int i=0,sz=Seg[rt].size();i<sz;++i)
    {
        int now=Seg[rt][i];
        int l=lower_bound(temp+1,temp+1+cnt,q[now].l)-temp;
        int r=upper_bound(temp+1,temp+1+cnt,q[now].r)-temp-1;
        ans[now]=max(ans[now],TR.query(TR.root[l-1],TR.root[r],q[now].x));
    }
}
node t1[N],t2[N];
void solve(int l,int r,int gol,int gor,int rt)
{
    if(gol>gor)return;
    int cnt1=0,cnt2=0;
    int mid=(l+r)>>1;
    cal(rt,gol,gor);
    if(l==r)return;
    //左右放  类似整体二分
    for(int i=gol;i<=gor;i++)
    {
        if(goods[i].day<=mid)t1[++cnt1]=goods[i];
        else t2[++cnt2]=goods[i];
    }
    //必须按顺序摆放 符合单调性  因为你要二分呀艹
    for(int i=1;i<=cnt1;i++)goods[i+gol-1]=t1[i];
    for(int i=1;i<=cnt2;i++)goods[i+cnt1+gol-1]=t2[i];
    solve(l,mid,gol,cnt1+gol-1,ls);
    solve(mid+1,r,gor-cnt2+1,gor,rs);
}

int main()
{
    int n,m;
    sc("%d%d",&n,&m);
    TR.Init();
    for(int i=1,temp;i<=n;++i)sc("%d",&temp),TR.insert(TR.root[i],TR.root[i-1],temp);
    int op;
    for(int i=1;i<=m;++i)
    {
        int l,r,shop,v,d,x;
        sc("%d",&op);
        if(op==0)sc("%d%d",&shop,&v),goods[++cg]={cg,shop,v};
        else
        {
            sc("%d%d%d%d",&l,&r,&x,&d);
            q[++cq]={l,r,x,max(1,cg-d+1),cg};
            ans[cq]=TR.query(TR.root[l-1],TR.root[r],x);
        }
    }
    for(int i=1;i<=cq;++i)update(q[i].L,q[i].R,i,1,cg,1);
    sort(goods+1,goods+1+cg);
    solve(1,cg,1,cg,1);
    for(int i=1;i<=cq;++i)pr("%d
",ans[i]);
    return 0;
}

/**************************************************************************************/