[CQOI2018] 异或序列

题目描述:

雾.

题目分析:

看上去像是个莫队啊.
维护一个前缀异或和 val[]
a[l]^a[l+1]^…^a[r]=val[r]^val[l-1]
如果a^b=c 那么 a^c=b
这样我们开一个桶记录下当前前缀异或和为X的数量,就可以O(1)转移惹.

题目链接:

BZOJ 5301
Luogu 4462

Ac 代码:

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
const int maxm=1e5+100;
int block_len;
int val[maxm],cnt[maxm],pos[maxm],Ans[maxm];
int n,m,ans,k;
struct node{
    int l,r,id;
}q[maxm];
inline bool comp(node x,node y)
{
    return pos[x.l]==pos[y.l]?x.r<y.r:x.l<y.l;
}
inline void update(int x,int f)
{
    if(f==1)
    {
        ans+=cnt[k^x];
        cnt[x]++;
    }
    else
    {
        cnt[x]--;
        ans-=cnt[k^x];
    }
}
int main()
{
    scanf("%d%d%d",&n,&m,&k);
    block_len=std::sqrt(n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&val[i]),val[i]^=val[i-1];
        pos[i]=(i-1)/block_len+1;
    }
    for(int i=1;i<=m;i++)
     scanf("%d%d",&q[i].l,&q[i].r),q[i].id=i;
    std::sort(q+1,q+m+1,comp);
    int l=1,r=0;
    for(int i=1;i<=m;i++)
    {
        q[i].l--;
        while(l<q[i].l) update(val[l++],-1);
        while(l>q[i].l) update(val[--l],1);
        while(r>q[i].r) update(val[r--],-1);
        while(r<q[i].r) update(val[++r],1);
        Ans[q[i].id]=ans;
    }
    for(int i=1;i<=m;i++)
     printf("%d\n",Ans[i]);
    return 0;
}
# P4462 [CQOI2018] 异或序列 ## 题目描述 已知一个长度为 $n$ 的整数数列 $a_1,a_2,\dots,a_n$,给定查询参数 $l,r$,问在 $a_l,a_{l+1},\dots,a_r$ 区间内,有多少子区间满足异或等于 $k$。也就是说,对于所有的 $x,y (l \leq x \leq y \leq r)$,能够满足 $a_x \oplus a_{x+1} \oplus \dots \oplus a_y = k$ 的 $x,y$ 有多少组。 ## 输入格式 输入文件第一行,为 $3$ 个整数 $n,m,k$。 第二行为空格分开的 $n$ 个整数,即 $a_1,a_2,..a_n$。 接下来 $m$ 行,每行两个整数 $l_j,r_j$,表示一查询。 ## 输出格式 输出文件共 $m$ 行,对应每个查询的计算结果。 ## 输入输出样例 #1 ### 输入 #1 ``` 4 5 1 1 2 3 1 1 4 1 3 2 3 2 4 4 4 ``` ### 输出 #1 ``` 4 2 1 2 1 ``` ## 说明/提示 对于 $30\%$ 的数据,$1 \leq n, m \leq 1000$。 对于 $100\%$ 的数据,$1 \leq n, m \leq 10^5$,$0 \leq k, a_i \leq 10^5$,$1 \leq l_j \leq r_j \leq n$。 #include <bits/stdc++.h> #define ll long long using namespace std; const ll N=2e5; struct node{ ll l,r,f,ans; }b[N+5]; ll n,m,k; ll blk[N+5],lenb; ll a[N+5],pre[N+5]; ll cnt[N+5]; ll lp,rp,s; bool cmp(node l1,node l2){ if(blk[l1.l]==blk[l2.l]) return l1.r<l2.r; return blk[l1.l]<blk[l2.l]; } bool cmpf(node l1,node l2){ return l1.f<l2.f; } void solve(ll l,ll r,ll &ans){ while(lp<l){ s-=cnt[k^pre[lp-1]]; cnt[pre[lp]]--; lp++; } while(lp>l){ lp--; cnt[pre[lp]]++; s+=cnt[k^pre[lp-1]]; } while(rp<r){ rp++; s+=cnt[k^pre[rp]]; cnt[pre[rp]]++; if((pre[lp-1]^pre[rp])==k) s++; } while(rp>r){ s-=cnt[k^pre[rp]]; cnt[pre[rp]]--; if((pre[lp-1]^pre[rp])==k) s--; rp--; } ans=s; return; } int main(){ scanf("%lld%lld%lld",&n,&m,&k); ll sx=1; lenb=(ll)sqrt(n); for(ll i=1;i<=n;i++){ blk[i]=sx; if(sx%lenb==0) sx++; } for(ll i=1;i<=n;i++){ scanf("%lld",&a[i]); pre[i]=(pre[i-1]^a[i]); } for(ll i=1;i<=m;i++){ scanf("%lld%lld",&b[i].l,&b[i].r); } sort(b+1,b+m+1,cmp); lp=rp=cnt[pre[1]]=1; if(k==0) s=1; for(ll i=1;i<=m;i++){ solve(b[i].l,b[i].r,b[i].ans); } sort(b+1,b+m+1,cmpf); for(ll i=1;i<=m;i++) printf("%lld\n",b[i].ans); return 0; } 不改变原来代码写法情况下标出并修正错误
最新发布
07-31
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值