bzoj4571 [Scoi2016]美味(主席树+按位贪心)

考虑如果没有+x,那么直接可持久化Trie+按位贪心即可。
有+x怎么办呢?我们考虑在Trie树上跑的过程,实际上每次选择时就是看一个区间内有没有数。那么我现在就是要把这个区间左移x位即可。用主席树模拟Trie树的这一过程即可。
复杂度 O(nlogn+mlog2n)

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m=0,rt[N<<1],a[N<<1],owo=0;
struct node{
    int lc,rc,x;
}tr[N*2*20];
inline void ins(int &p,int l,int r,int x){
    tr[++owo]=tr[p];p=owo;tr[p].x++;
    if(l==r) return;int mid=l+r>>1;
    if(x<=mid) ins(tr[p].lc,l,mid,x);
    else ins(tr[p].rc,mid+1,r,x);
}
inline int ask(int p1,int p2,int l,int r,int x,int y){
    if(x>y) return 0;if(x<=l&&r<=y) return tr[p2].x-tr[p1].x;
    int mid=l+r>>1,res=0;
    if(x<=mid) res+=ask(tr[p1].lc,tr[p2].lc,l,mid,x,y);
    if(y>mid) res+=ask(tr[p1].rc,tr[p2].rc,mid+1,r,x,y);return res;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();int ofo=read();
    for(int i=1;i<=n;++i) a[i]=read(),m=max(m,a[i]);
    for(int i=1;i<=n;++i) rt[i]=rt[i-1],ins(rt[i],0,m,a[i]);
    while(ofo--){
        int b=read(),x=read(),l=read(),r=read(),ans=0;
        for(int i=17;i>=0;--i){
            int y=ans|(((b>>i&1)^1)<<i);
            if(ask(rt[l-1],rt[r],0,m,max(y-x,0),min(y+(1<<i)-1-x,m))) ans=y;
            else ans|=(b>>i&1)<<i;
        }printf("%d\n",ans^b);
    }return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值