bzoj4137 [FJOI2015]火星商店问题

　　比较容易想到的做法是线段树套字典树，修改操作时在字典树上经过的节点维护一个最近被访问过的时间，这样询问操作只经过满足时间条件的节点，时间复杂度O(NlogN^2)但是因为线段树每个节点都要套个字典树，这样的话空间是不够的，不过由于可以离线处理，我们可以先把每个修改和询问操作所访问的线段树节点保存下来，然后一个个节点去做，这样的话空间复杂度就可以保证了。

代码

#include<cstdio>
#include<set>
#include<vector>
#define N 1000010
#define M 10000010
using namespace std;
int l[N],r[N],n,m,i,a,typ,b,L,R,ti,ans[N],tot,p[N],cnt,flag;
int s[M][2],t[M];
vector<pair<int,int> > vec[N];
vector<int> id[N];
void build(int x,int a,int b)
{
    int m;
    l[x]=a;r[x]=b;
    if (b-a>1)
    {
        m=(a+b)>>1;
        build(2*x,a,m);
        build(2*x+1,m,b);
    }
}
void insert(int x,int a,int b,pair<int,int> c,int typ)
{
    int m;
    if (typ==0)
    {
        vec[x].push_back(c);
        id[x].push_back(typ);
    }
    if ((a<=l[x])&&(r[x]<=b))
    {
        if (typ)
        {
            vec[x].push_back(c);
            id[x].push_back(typ);
        }
        return; 
    }
    m=(l[x]+r[x])>>1;
    if (a<m) insert(2*x,a,b,c,typ);
    if (m<b) insert(2*x+1,a,b,c,typ);
} 
void cl(int x)
{
    int i;
    for (i=1;i<=17;i++)
    {
        p[17-i+1]=x%2;
        x=x/2;
    }
}
void change(int x,int y)
{
    int now,i;
    cl(x);
    now=1;
    for (i=1;i<=17;i++)
    {
        if (s[now][p[i]]==0) 
            s[now][p[i]]=++tot;
        now=s[now][p[i]];
        t[now]=max(t[now],y);
    }
}
int query(int x,int y)
{
    int now,i,ans=0;
    cl(x);
    now=1;
    for (i=1;i<=17;i++)
    if ((s[now][1-p[i]])&&(t[s[now][1-p[i]]]>=y))
    {
        now=s[now][1-p[i]];
        ans=ans+(1<<(17-i));
    }
    else
    now=s[now][p[i]];
    return ans;
}
void gao(int x)
{
    int i,a,b;
    for (i=1;i<=tot;i++)
    s[i][0]=s[i][1]=t[i]=0;
    tot=1;
    for (i=0;i<vec[x].size();i++)
    {
        a=vec[x][i].first;
        b=vec[x][i].second;
        if (id[x][i]==0)
        change(a,b);
        else
        ans[id[x][i]]=max(ans[id[x][i]],query(a,b));
    }
    if (r[x]-l[x]==1)
    return;
    gao(2*x);
    gao(2*x+1);
} 
int main()
{
    scanf("%d%d",&n,&m);
    build(1,0,n);
    for (i=1;i<=n;i++)
    {
        scanf("%d",&a); 
        insert(1,i-1,i,make_pair(a,m+1),0);
    }
    ti=1;
    for (i=1;i<=m;i++)
    {
        scanf("%d",&typ);
        if (typ==1)
        {
            scanf("%d%d%d%d",&L,&R,&a,&b);
            insert(1,L-1,R,make_pair(a,ti-b+1),++cnt);
        }
        else
        {
            ti++;
            scanf("%d%d",&a,&b);
            insert(1,a-1,a,make_pair(b,ti),0);
        }
    }
    gao(1);
    for (i=1;i<=cnt;i++)
    printf("%d
",ans[i]);
}