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

时间:2023-03-09 04:05:05
bzoj4571/luogu3293 美味 (主席树+贪心)

首先想到建出可持久化trie树然后在上面贪心,但是它加了一个数所以不能这么做

但依然可以贪心,仿照上面那个的过程,如果设y是在第i位上^b是1的数(前面的位数已经贪好了),我只要在[l,r]范围内能有[y-x,y+(1<<i)-x-1)]的数,那这位异或出来就是可以是1的

 #include<bits/stdc++.h>
#define pa pair<int,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn=2e5+,maxp=maxn*; inline ll rd(){
ll x=;char c=getchar();int neg=;
while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
while(c>=''&&c<='') x=x*+c-'',c=getchar();
return x*neg;
} int root[maxn],N,M,L=;
int ch[maxp][],v[maxp],pct; void insert(int pre,int &p,int l,int r,int x){
p=++pct;
v[p]=v[pre]+;
if(l<r){
int m=l+r>>;
if(x<=m) ch[p][]=ch[pre][],insert(ch[pre][],ch[p][],l,m,x);
else ch[p][]=ch[pre][],insert(ch[pre][],ch[p][],m+,r,x);
}
} int query(int pre,int p,int l,int r,int x,int y){
if(x<=l&&r<=y) return v[p]-v[pre];
int m=l+r>>,re=;
if(x<=m) re=query(ch[pre][],ch[p][],l,m,x,y);
if(y>=m+) re+=query(ch[pre][],ch[p][],m+,r,x,y);
return re;
} inline int solve(int b,int x,int pre,int p){
int n=;
for(int i=;i>=;i--){
int y=n|((((b>>i)&)^)<<i);
// printf("%d %d\n",n,y);
if(query(pre,p,,L,max(,y-x),min(L,y+(<<i)-x-))) n=y;
else if((b>>i)&) n|=(<<i);
}
return n^b;
} int main(){
//freopen("","r",stdin);
int i,j,k;
N=rd();M=rd();
for(i=;i<=N;i++){
int a=rd();
insert(root[i-],root[i],,L,a);
}
for(i=;i<=M;i++){
int b=rd(),x=rd(),l=rd(),r=rd();
printf("%d\n",solve(b,x,root[l-],root[r]));
}
return ;
}