P3293 [SCOI2016]美味(主席树+最大异或合）

传送门

题意:给你n个数,每次给你一个l到r区间，还有一个bi和xi，每次输出l到r中使得bi XOR (aj+xi)值最大，输出最大值

题解：首先是最大异或合，最大异或合的求法是基于贪心的思想，来源于：10000XOR00000>10000XOR11111，所以对于我们要找的最大值，我们将所有的数拆成二进制，从高位开始选取，如果有使当前位异或为1的值就直接选取，由于只是贪心了当前位所以对于序列的值来说就会出现一个范围：

        if(bi[a]==0){
            ansl=er[a]+ans-xi;
            ansr=er[a+1]-1+ans-xi;
        }
        else{
            ansl=ans-xi;
            ansr=er[a]-1+ans-xi;
        }

将序列的值换成二进制插入主席树后对应区间里面是否存在这个值，在主席树种就表现为当前记录的这一位的个数是否大于0

AC代码：

#include<stdio.h>
#include<vector>
#include<string.h>
#include<algorithm>
#include<iostream>
#define ll long long
using namespace std;
const int maxn=2e5+5;
//const ll inf=9223372036854775800;
const int inf=1e5+5;
int ls[maxn*20],rs[maxn*20],num[maxn*20],ti[maxn],ai[maxn],tot;
int bi[20],xi,ans,p;
int er[20];
void add(int &now,int lst,int l,int r,int z){
    now=++tot;
    ls[now]=ls[lst];
    rs[now]=rs[lst];
    num[now]=num[lst]+1;
    if(l==r)
        return ;
    int mid=(l+r)>>1;
    if(z<=mid)
        add(ls[now],ls[lst],l,mid,z);
    else
        add(rs[now],rs[lst],mid+1,r,z);
}
int query(int now,int lst,int l,int r,int l1,int r1){
    if(num[now]-num[lst]==0)
        return 0;
    if(l1<=l&&r1>=r)
        return num[now]-num[lst];
    int mid=(l+r)>>1,num1=0;
    if(mid>=l1)
        num1+=query(ls[now],ls[lst],l,mid,l1,r1);
    if(num1>=1)
        return num1;
    if(mid+1<=r1)
        num1+=query(rs[now],rs[lst],mid+1,r,l1,r1);
    return num1;
}
void get( ){
    for(int j=17;j>=0;j--)
        bi[j+1]=((p>>j)&1);
}
void work(int now,int lst){
    //printf("%d %d\n",now,lst);
    get();
    int ansl,ansr,ansz=0;
    for(int a=18;a>=1;a--){
        if(bi[a]==0){
            ansl=er[a]+ans-xi;
            ansr=er[a+1]-1+ans-xi;
        }
        else{
            ansl=ans-xi;
            ansr=er[a]-1+ans-xi;
        }
        ansl=max(ansl,0);
        ansr=min(ansr,inf);
          int num1=query(now,lst,0,inf,ansl,ansr);
            if(num1&&bi[a]==0)
                ans+=er[a];
            else if(num1==0&&bi[a]==1)
                ans+=er[a];
    }
}
int main( )
{
    er[0]=0;
    er[1]=1;
    for(int a=2;a<=19;a++)
        er[a]=er[a-1]*2;
    int n,m;
    scanf("%d%d",&n,&m);
    for(int a=1;a<=n;a++)
        scanf("%d",&ai[a]);
    for(int a=1;a<=n;a++)
        add(ti[a],ti[a-1],0,inf,ai[a]);
    int li,ri;
    for(int a=1;a<=m;a++){
        ans=0;
        scanf("%d%d%d%d",&p,&xi,&li,&ri);
        work(ti[ri],ti[li-1]);
        printf("%d\n",ans^p);
    }
}