5301: [Cqoi2018]异或序列

本文介绍了一种使用莫队算法解决特定区间异或和等于k的问题的方法。通过对输入数据进行预处理,并利用前缀异或和的概念,实现了高效查询。代码示例展示了如何通过调整区间边界来更新当前异或和等于k的子区间的数量。

看完题,容易想到把求一段区间改成求 前r的异或和 ^ 前l-1的异或和 等于k,
那么这个东西发现可以直接用莫队统计,没了.
c++代码如下:

#include<bits/stdc++.h>
#define rep(i,x,y) for(register int i = x ; i <= y ; ++ i)
#define repd(i,x,y) for(register int i = x ; i >= y ; -- i)
using namespace std;
typedef long long ll;
template<typename T>inline void read(T&x)
{
    char c;int sign = 1;x = 0;
    do { c = getchar(); if(c == '-') sign = -1; }while(!isdigit(c));
    do { x = x * 10 + c - '0'; c = getchar(); }while(isdigit(c));
    x *= sign;
}

const int N = 1e5+850;
int n,m,k,a[N];
int num[N],ans[N],belong[N];
struct DATA
{
    int l,r,id;
}q[N];

const bool cmp(DATA a,DATA b) { return belong[a.l] < belong[b.l] || belong[a.l] == belong[b.l] && belong[a.r] < belong[b.r]; }

int main()
{
    read(n); read(m); read(k);
    rep(i,1,n) read(a[i]);
    rep(i,2,n) a[i] = a[i]^a[i-1];
    rep(i,1,m) read(q[i].l), read(q[i].r),q[i].id = i;

    int d = max(1,(int)pow(n,2.0/3));

    rep(i,1,n) belong[i] = i/d;

    sort(q + 1,q + 1 + m,cmp);

    int r = 0,l = 1,now = 0;
    rep(i,1,m)
    {
        --q[i].l;
        while(l < q[i].l)
        {
            --num[a[l]];
            now -= num[a[l++]^k];
        }

        while(l > q[i].l)
        {
            now += num[a[--l]^k];
            ++num[a[l]];
        }

        while(r < q[i].r)
        {
            now += num[a[++r]^k];
            ++num[a[r]];
        }

        while(r > q[i].r)
        {
            --num[a[r]];
            now -= num[a[r--]^k];
        }

        ans[q[i].id] = now;
    }

    rep(i,1,m) 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
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值