bzoj4571/luogu3293 美味 (主席树+贪心)

 1 #include<bits/stdc++.h>
 2 #define pa pair<int,int>
 3 #define CLR(a,x) memset(a,x,sizeof(a))
 4 using namespace std;
 5 typedef long long ll;
 6 const int maxn=2e5+10,maxp=maxn*20;
 7 
 8 inline ll rd(){
 9     ll x=0;char c=getchar();int neg=1;
10     while(c<‘0‘||c>‘9‘){if(c==‘-‘) neg=-1;c=getchar();}
11     while(c>=‘0‘&&c<=‘9‘) x=x*10+c-‘0‘,c=getchar();
12     return x*neg;
13 }
14 
15 int root[maxn],N,M,L=100000;
16 int ch[maxp][2],v[maxp],pct;
17 
18 void insert(int pre,int &p,int l,int r,int x){
19     p=++pct;
20     v[p]=v[pre]+1;
21     if(l<r){
22         int m=l+r>>1;
23         if(x<=m) ch[p][1]=ch[pre][1],insert(ch[pre][0],ch[p][0],l,m,x);
24         else ch[p][0]=ch[pre][0],insert(ch[pre][1],ch[p][1],m+1,r,x);
25     }
26 }
27 
28 int query(int pre,int p,int l,int r,int x,int y){
29     if(x<=l&&r<=y) return v[p]-v[pre];
30     int m=l+r>>1,re=0;
31     if(x<=m) re=query(ch[pre][0],ch[p][0],l,m,x,y);
32     if(y>=m+1) re+=query(ch[pre][1],ch[p][1],m+1,r,x,y);
33     return re;
34 }
35 
36 inline int solve(int b,int x,int pre,int p){
37     int n=0;
38     for(int i=17;i>=0;i--){
39         int y=n|((((b>>i)&1)^1)<<i);
40         // printf("%d %d\n",n,y);
41         if(query(pre,p,0,L,max(0,y-x),min(L,y+(1<<i)-x-1))) n=y;
42         else if((b>>i)&1) n|=(1<<i);
43     }
44     return n^b;
45 }
46 
47 int main(){
48     //freopen("","r",stdin);
49     int i,j,k;
50     N=rd();M=rd();
51     for(i=1;i<=N;i++){
52         int a=rd();
53         insert(root[i-1],root[i],0,L,a);
54     }
55     for(i=1;i<=M;i++){
56         int b=rd(),x=rd(),l=rd(),r=rd();
57         printf("%d\n",solve(b,x,root[l-1],root[r]));
58     }
59     return 0;
60 }