SP3267 DQUERY - D-query【莫队】或【主席树】

20 篇文章 0 订阅
7 篇文章 0 订阅

传送门
给定一个长度为 N N N的数组,求区间 l ∼ r l \sim r lr 不同的数字

分析

对于此类后面出现会影响(覆盖)前面出现效果的问题
需要往当前位置左边最后一次出现位置思考
比如, m e x mex mex 也是如此
这道题求 区间mex 也是同样处理

这种位置会覆盖前一个位置影响的,都考虑左边最近一个的位置
静态区间在线查询 M E X MEX MEX
主席树,维护最近值域出现的位置
r r r 树上二分找值域的出现位置小于 l l l 的最小值
证明:
假如一个数 x x x 在区间 l ∼ r l \sim r lr 出现过,那么它在 r r r数上对应值域的值,一定是大于等于 l l l
推得,一个数在 l ∼ r l \sim r lr区间出现过的条件是,在 r r r树上值域的值一定大于 l l l
反之,则不出现。 m e x mex mex则是不出现的最小数,对应在树上就是值域最小且出现时间小于 l l l

所以,此题有主席树解法

主席树

主席树第 r r r个版本的信息是, 1 ∼ r 1 \sim r 1r区间内,所有出现的数,最后一次出现的位置
假如,这里建立一个权值线段树,要在权值线段树上找到所有在 l ∼ r l \sim r lr 区间的个数(由于是第 r r r棵树,所以不存在大于 r r r的位置,找大于等于 l l l的即可),查询复杂度为 n l o g n nlogn nlogn,因为要遍历到每一个叶子节点进行判断,但是我们这里不需要管权值的大小,只需要知道在大于 l l l 的个数,因此我们可以在此基础上,忽略权值大小,直接对位置进行维护,因为 一个下标唯一确定一个权值,所以更新的时候对上一次出现的下标位置进行删除 − 1 -1 1,更新当前位置 + 1 +1 +1,查询的时候直接查区间 l ∼ r l \sim r lr的和即可

莫队

莫队写法就比较正常了
由于是不带修改的,分块正常分配为 N \sqrt{N} N 即可,对不同块的 l l l进行从小到大排序,同一块的根据 r r r从小到大奇偶排序
假如说区间中,增加一个数的时候,发现是第一次出现,计数器 + 1 +1 +1
删除一个数的时候,发现目前只有一个了,删了就没了,计数器 − 1 -1 1
答案易知

代码

主席树

//SP3267 
/*
  @Author: YooQ
*/
#include <bits/stdc++.h>
using namespace std;
#define sc scanf
#define pr printf
#define ll 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];

struct Tr {
	int k, l, r;
}tr[MAX_N<<1];
int indx = 0;
int root[MAX_N];

int clone(int x) {
	tr[++indx] = tr[x];
	return indx;
}

void push_up(int rt) {
	tr[rt].k = tr[tr[rt].l].k + tr[tr[rt].r].k;
}

#define mid l + ((r-l)>>1)
void update(int& rt, int l, int r, int x, int k) {
	rt = clone(rt);
	if (l == r) {
		tr[rt].k += k;
		return;
	}
	if (x <= mid) update(tr[rt].l, l, mid, x, k);
	if (x  > mid) update(tr[rt].r, mid+1, r, x, k);
	push_up(rt);
}

int query(int rt, int l, int r, int x, int y) {
	if (!rt) return 0;
	if (x <= l && r <= y) {
		return tr[rt].k;
	}
	if (y <= mid) return query(tr[rt].l, l, mid, x, y);
	if (x  > mid) return query(tr[rt].r, mid+1, r, x, y);
	return query(tr[rt].l, l, mid, x, y) + query(tr[rt].r, mid+1, r, x, y);
}

int pre[MAX_N];

void solve(){
	sc("%d", &N);
	for (int i = 1; i <= N; ++i) {
		sc("%d", &arr[i]);
		root[i] = root[i-1];
		if (pre[arr[i]]) {
			update(root[i], 1, N, pre[arr[i]], -1);
		}
		update(root[i], 1, N, i, 1);
		pre[arr[i]] = i;
	}
	
	sc("%d", &M);
	int l, r;
	for (int i = 1; i <= M; ++i) {
		sc("%d%d", &l, &r);
		pr("%d\n", query(root[r], 1, N, l, r));
	}
}

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

	return AC;
}

莫队

较基础的简单莫队

pending...
  • 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、付费专栏及课程。

余额充值