CF-617-E. XOR and Favorite Number【普通莫队】

给定一个长度为 N N N的数组,一个整数 K K K,求区间 l ∼ r l \sim r lr中有多少 p a i r { i , j } , i < = j pair\{i,j\},i<=j pair{i,j},i<=j,满足 a i ⊕ a i + 1 ⊕ . . . ⊕ a j = K a_i \oplus a_{i+1} \oplus ... \oplus a_j=K aiai+1...aj=K
传送门

分析

先考虑暴力
a i ⊕ a i + 1 ⊕ . . . ⊕ a j a_i \oplus a_{i+1} \oplus ... \oplus a_j aiai+1...aj 就是要求,异或前缀下的 o s u m i − 1 ⊕ o s u m j osum_{i-1} \oplus osum_j osumi1osumj
固定起点 l l l,暴力往后走
假如走到 i i i的位置,计算 o s u m l − 1 ⊕ o s u m i osum_{l-1} \oplus osum_i osuml1osumi,如果等于 K K K的话就加 1 1 1,依次计算,复杂度 M ∗ N 2 M*N^2 MN2
再想想优化
上面复杂度显然不行,考虑固定左端点,拓展右端点的时候直接把所有以右端点结尾的答案都计算出来
如果当前右端点走到了 i i i位置,因为我们要找到异或后等于 K K K的值,既然要让 i i i能做出贡献,那么此时一定在前面计算的某一个左端点 j j j,存在 o s u m j − 1 ⊕ o s u m i = K osum_{j-1} \oplus osum_i = K osumj1osumi=K,推出 o s u m j − 1 = o s u m i ⊕ K osum_{j-1} = osum_i \oplus K osumj1=osumiK,也就是说,所有以 i i i为结尾的区间异或等于 K K K的贡献,就是前面计算的所有,值等于 o s u m i ⊕ K osum_i \oplus K osumiK的数量

考虑添加一个数 x x x,计数器 c n t [ x ⊕ K ] + + cnt[x \oplus K]++ cnt[xK]++,答案 a n s + = c n t [ x ⊕ K ] ans+=cnt[x \oplus K] ans+=cnt[xK]
考虑删除一个数 x x x,计数器 c n t [ x ⊕ K ] − − cnt[x \oplus K]-- cnt[xK],答案 a n s − = c n t [ x ⊕ K ] ans -= cnt[x \oplus K] ans=cnt[xK]

注意细节,题目要求是 i = = j i==j i==j是合法的,所以得先计数器改变,再统计答案个数

那么计算一个区间内满足条件的个数,时间复杂度就是 N N N

接下来就是,如何求出多个区间的问题?
普通莫队就能够解决
添加,删除的操作上面应该已经说清楚了

代码

//CF617
/*
  @Author: YooQ
*/
#include <bits/stdc++.h>
using namespace std;
#define sc scanf
#define pr printf
#define ll long long
#define int long long
#define FILE_OUT freopen("out", "w", stdout);
#define FILE_IN freopen("in", "r", stdin);
#define debug(x) cout << #x << ": " << x << "\n";
#define AC 0
#define WA 1
#define INF 0x3f3f3f3f
const ll MAX_N = 1e6+5;
const ll MOD = 1e9+7;
int N, M, K;

int arr[MAX_N];
int cnt[MAX_N<<1];

int mo = 0;
int ans[MAX_N];

void add(int x) {
	mo += cnt[arr[x]^K];
	++cnt[arr[x]];
}

void del(int x) {
	--cnt[arr[x]];
	mo -= cnt[arr[x]^K];
}

int block;

int idb(int x) {
	return x/block;
}

struct Qr {
	int l, r, id;
	bool operator < (const Qr& B) const {
		return idb(l) ^ idb(B.l)? l < B.l : (idb(l)&1) ? r < B.r : r > B.r;
	}
}qr[MAX_N];

void solve(){
	sc("%lld%lld%lld", &N, &M, &K);
	block = max(sqrt(N), sqrt((N*N)/M));
	
	for (int i = 1; i <= N; ++i) {
		sc("%lld", &arr[i]);
		arr[i] ^= arr[i-1];
	}
	for (int i = 1; i <= M; ++i) {
		sc("%lld%lld", &qr[i].l, &qr[i].r);
		--qr[i].l;
		qr[i].id = i;
	}
	
	sort(qr+1, qr+1+M);
	
	int l = 1, r = 0;
	mo = 0;
	for (int i = 1; i <= M; ++i) {
		while (l > qr[i].l) add(--l);
		while (r < qr[i].r) add(++r);
		while (l < qr[i].l) del(l++);
		while (r > qr[i].r) del(r--);
		ans[qr[i].id] = mo;
	}
	
	for (int i = 1; i <= M; ++i) {
		pr("%lld\n", ans[i]);
	}
}

signed main()
{
	#ifndef ONLINE_JUDGE
	//FILE_IN
	FILE_OUT
	#endif
	int T = 1;//cin >> T;
	while (T--) solve();

	return AC;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Hexrt

客官,请不要给我小费!

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值