Xors on Segments Codeforces

http://codeforces.com/problemset/problem/620/F

此题是莫队，但是不能用一般的莫队做，因为是最优化问题，没有办法在删除元素的时候维护答案。

---

这题的方法（好像还有个名字叫"回滚莫队“，说明：https://blog.csdn.net/maverickfw/article/details/72988286）：

官方题解：http://codeforces.com/blog/entry/22936

设完成一次"加入贡献"的操作的复杂度是O(p)

首先定块大小sz，定每个坐标应该属于的块。

对于左右端点在同一块内的询问，暴力处理即可。每一个询问需要复杂度O(p*sz)

对于其他询问，按左端点所属的块分类。

分类完后，对于同一类的询问，按右端点从小到大排序。对于每一类的询问分开处理。

对于每一类，先把"当前维护区间"的左端点定为 属于该块的坐标 中最靠右的坐标（这样的话，该类询问的左端点都小于等于这个点且离这个点距离不超过sz），右端点定为左端点-1（为了让区间为空）。还要清空当前答案。

对于该类中每一个询问，先照常移动右端点到与询问的右端点相同，然后记下此时的答案，然后把左端点移动到与询问的左端点相同，得到该询问的答案并记录；然后把左端点的移动还原，把答案还原即可（这样就避开了维护答案）

对于每一类，右端点最多有O(n)次移动，而有n/sz类，这一部分复杂度是O(p*n*n/sz)；对于每一个询问，额外还要O(p*sz)的复杂度进行左端点的调整

因此，总复杂度是O(p*(n*n/sz+m*sz))，当sz=sqrt(n*n/m)时最小

---

解决了莫队的问题，还要解决异或最大值维护的问题。这个是用01字典树做的，就不写了。。。看官方题解吧

（大概要支持在01字典树中查询：集合中，对于所有"第二权值"小于等于/大于等于特定值的元素，查询与另一给定值异或的最大值，挺奇怪的）

错误记录：

1.trie没有开垃圾回收，内存不够

2.计算块大小时，没有注意块大小可能变成0

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<set>
#include<queue>
#include<vector>
using namespace std;
int f[2000100];
int n,m,sz,sz1;
//int ttt;
int lft[40];
namespace Trie
{
const int l2n=20;
int sz[4000100],ch[4000100][2],dat[4000100];
multiset<int> ss[4001000];
int x,kkk;
queue<int> qqq;
int getnode()
{
    //yttt=max(ttt,int(2000100-qqq.size()));
    int t=qqq.front();qqq.pop();
    return t;
}
void delnode(int x)
{
    qqq.push(x);
}
void deltree(int x)
{
    if(!x)    return;
    deltree(ch[x][0]);
    deltree(ch[x][1]);
    ss[x].clear();ch[x][0]=ch[x][1]=dat[x]=sz[x]=0;
    delnode(x);
}
void init()
{
    int i;
    for(i=1;i<4000100;i++)    qqq.push(i);
    lft[0]=1;
    for(i=1;i<=l2n;i++)    lft[i]=lft[i-1]<<1;
}
int getdat1(int x)    {return x?dat[x]:0x3f3f3f3f;}
void _ins1(int p,int &num)
{
    if(!num)    num=getnode(),dat[num]=0x3f3f3f3f;
    sz[num]++;
    if(p>=0)    _ins1(p-1,ch[num][!!(x&lft[p])]),dat[num]=min(getdat1(ch[num][0]),getdat1(ch[num][1]));
    else    ss[num].insert(kkk),dat[num]=*ss[num].begin();
}
void ins1(int d,int kk,int &num)    {x=d;kkk=kk;_ins1(l2n-1,num);}
//ins1维护插入kk的最小值
void _ins2(int p,int &num)
{
    if(!num)    num=getnode();
    sz[num]++;
    if(p>=0)    _ins2(p-1,ch[num][!!(x&lft[p])]),dat[num]=max(dat[ch[num][0]],dat[ch[num][1]]);
    else    ss[num].insert(kkk),dat[num]=*ss[num].rbegin();
}
void ins2(int d,int kk,int &num)    {x=d;kkk=kk;_ins2(l2n-1,num);}
//ins2维护插入kk的最大值


void _era1(int p,int &num)
{
    sz[num]--;
    if(p>=0)    _era1(p-1,ch[num][!!(x&lft[p])]),dat[num]=min(getdat1(ch[num][0]),getdat1(ch[num][1]));
    else    ss[num].erase(ss[num].find(kkk)),dat[num]=ss[num].empty()?0x3f3f3f3f:*ss[num].begin();
}
void era1(int d,int kk,int &num)    {x=d;kkk=kk;_era1(l2n-1,num);}
void _era2(int p,int &num)
{
    sz[num]--;
    if(p>=0)    _era2(p-1,ch[num][!!(x&lft[p])]),dat[num]=max(dat[ch[num][0]],dat[ch[num][1]]);
    else    ss[num].erase(ss[num].find(kkk)),dat[num]=ss[num].empty()?0:*ss[num].rbegin();
}
void era2(int d,int kk,int &num)    {x=d;kkk=kk;_era2(l2n-1,num);}


int que1(int x,int kk,int num)//que1查询维护值小于等于kk的数中与x的xor最大值
{
    int ans=0,i;bool t;
    for(i=l2n-1;i>=0;--i)
    {
        t=x&lft[i];
        if(sz[ch[num][!t]]&&dat[ch[num][!t]]<=kk)    ans|=lft[i],num=ch[num][!t];
        else    num=ch[num][t];
    }
    return ans;
}
int que2(int x,int kk,int num)//que2查询维护值大于等于kk的数中与x的xor最大值
{
    int ans=0,i;bool t;
    for(i=l2n-1;i>=0;--i)
    {
        t=x&lft[i];
        if(sz[ch[num][!t]]&&dat[ch[num][!t]]>=kk)    ans|=lft[i],num=ch[num][!t];
        else    num=ch[num][t];
    }
    return ans;
}
}
using Trie::ins1;using Trie::ins2;
using Trie::era1;using Trie::era2;
using Trie::que1;using Trie::que2;
using Trie::deltree;
struct Q
{
    int l,r,num;
    Q(){}
    Q(int a,int b,int c):l(a),r(b),num(c){}
};
int a[50100],ans[5010];
int rt1,rt2;
int anss;
void add1(int p)//加入p的贡献
{
    ins1(f[p-1],p,rt1);ins2(f[p],p,rt2);
    anss=max(anss,que1(f[p],p,rt1));
    //printf("a%d
",f[p]^f[p-1]);
    anss=max(anss,que2(f[p-1],p,rt2));
}
void del1(int p)
{
    era1(f[p-1],p,rt1);era2(f[p],p,rt2);
}
vector<Q> q[1010];
bool cmp(const Q &a,const Q &b)    {return a.r<b.r;}
int main()
{
    int i,j,j2,l,r,tans;Trie::init();
    for(i=1;i<=2000000;i++)    f[i]=f[i-1]^i;
    scanf("%d%d",&n,&m);sz=max(1,int(sqrt(double(n)/m*n)));sz1=(n-1)/sz;
    for(i=1;i<=n;i++)    scanf("%d",&a[i]);
    for(i=1;i<=m;i++)
    {
        scanf("%d%d",&l,&r);
        if((l-1)/sz==(r-1)/sz)
        {
            deltree(rt1);deltree(rt2);
            rt1=rt2=anss=0;
            for(j=l;j<=r;j++)    add1(a[j]);
            ans[i]=anss;
            //printf("a%d %d
",i,ans[i]);
        }
        else
            q[(l-1)/sz].push_back(Q(l,r,i));
    }
    for(i=0;i<=sz1;i++)
    {
        sort(q[i].begin(),q[i].end(),cmp);
        l=(i+1)*sz;r=l-1;rt1=rt2=anss=0;
        for(j=0;j<q[i].size();j++)
        {
            while(r<q[i][j].r)    add1(a[++r]);
            tans=anss;
            for(j2=l-1;j2>=q[i][j].l;j2--)    add1(a[j2]);
            ans[q[i][j].num]=anss;
            for(j2=l-1;j2>=q[i][j].l;j2--)    del1(a[j2]);
            anss=tans;
        }
        deltree(rt1);deltree(rt2);
    }
    for(i=1;i<=m;i++)    printf("%d
",ans[i]);
    //printf("a%d
",ttt);
    return 0;
}