bzoj4888 [Tjoi2017]异或和(按位+树状数组)

我们首先做出前缀和s[i],那么一个连续和就可以表示成s[i]-s[j],0<=j< i
我们考虑s[i]的第k位-s[j]之后为1的有多少个,如果奇数个就对答案有贡献。

如果s[i]的第k位为1,s[j]需要满足第k位为0且前k-1位小于s[i]的前k-1位,或者第k位为1,且前k-1位大于s[i]的前k-1位。如果s[i]的第k位为0同理。
我们可以用两个权值BIT来维护。
复杂度 O(nlognlogw)

#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,s[N],bin[21],c[2][N*10],ans=0;
inline void add(int id,int x){
    for(;x<=s[n];x+=x&-x) c[id][x]++;
}
inline int ask(int id,int x){
    int res=0;for(;x;x-=x&-x) res+=c[id][x];return res;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();for(int i=1;i<=n;++i) s[i]=s[i-1]+read();bin[0]=1;
    for(int i=1;i<=20;++i) bin[i]=bin[i-1]<<1;
    for(int k=0;k<=20;++k){
        if(bin[k]>s[n]) break;
        memset(c,0,sizeof(c));add(0,1);int cnt=0;
        for(int i=1;i<=n;++i){
            int x=s[i]&(bin[k]-1);++x;
            if(s[i]&bin[k]){
                int res=ask(0,x);
                res+=ask(1,s[n])-ask(1,x);
                if(res&1) cnt^=1;add(1,x);
            }else{
                int res=ask(0,s[n])-ask(0,x);
                res+=ask(1,x);
                if(res&1) cnt^=1;add(0,x);
            }
        }if(cnt) ans|=bin[k];
    }printf("%d\n",ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值