BZOJ4571: [Scoi2016]美味

BZOJ4571: [Scoi2016]美味

https://lydsy.com/JudgeOnline/problem.php?id=4571

分析:

• 比较好的一道题。
• 按位贪心，假设由已经确定的那些位组成的数为\(ans\)。
• 令\(c\)为\((ans+(1<<i))\;xor\;b\) 低于\(i\)的位都是\(0\)的结果。
• 那么合法的数的区间就是\([c-x,c+(1<<i)-1-x]\)。
• 每次主席树查区间是否存在一个数即可。

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
using namespace std;
#define N 200050
#define inf 100000
int n,a[N],m,root[N],cnt;
int ls[N*22],rs[N*22],siz[N*22];
void update(int l,int r,int x,int &p,int q) {
    p=++cnt; ls[p]=ls[q]; rs[p]=rs[q]; siz[p]=siz[q]+1;
    if(l==r) return ;
    int mid=(l+r)>>1;
    if(x<=mid) update(l,mid,x,ls[p],ls[q]);
    else update(mid+1,r,x,rs[p],rs[q]);
}
int query(int l,int r,int x,int y,int p,int q) {
    if(siz[p]==siz[q]) return 0;
    if(x<=l&&y>=r) return siz[p]-siz[q];
    int mid=(l+r)>>1,re=0;
    if(x<=mid) re+=query(l,mid,x,y,ls[p],ls[q]);
    if(y>mid) re+=query(mid+1,r,x,y,rs[p],rs[q]);
    return re;
}
int mask=(1<<18)-1;
int qian(int p,int d) {
    return d&(mask-((1<<p)-1));
}
int main() {
    scanf("%d%d",&n,&m);
    int i;
    for(i=1;i<=n;i++) scanf("%d",&a[i]),update(0,inf,a[i],root[i],root[i-1]);
    int ans;
    while(m--) {
        int b,x,l,r;
        scanf("%d%d%d%d",&b,&x,&l,&r);
        ans=0;
        for(i=17;i>=0;i--) {
            int t=ans+(1<<i);
            int c=qian(i,t^b);
            int mn=max(c-x,0),mx=min(c+(1<<i)-1-x,inf);
            if(mn<=mx&&query(0,inf,mn,mx,root[r],root[l-1])) {
                ans+=(1<<i);
            }
        }
        printf("%d\n",ans);
    }
}

转载于:https://www.cnblogs.com/suika/p/10093081.html