区间众数(分块)

问题

n个数构成一个序列,查询 [ L , R ] [L, R] [L,R] 区间内出现次数最多的那个数的出现次数。时间、空间复杂度要求均为: O ( n ⋅ n ) O(n\cdot \sqrt{n}) O(nn )

分析

  • 基本思想:分块
  • 离散化,使得序列元素 a n a_n an 转换成 a n ′ ∈ [ 1 , n ] a_n^{'} \in [1, n] an[1,n]
  • 计算块 i 、 j i、j ij 间的最大次数,用 f [ i ] [ j ] f[i][j] f[i][j] 记录
  • 计算元素 a n ′ a_n^{'} an 在前 i i i 块中出现的次数(前缀和),用 g [ a ] [ i ] g[a][i] g[a][i] 记录
  • [ L , R ] [L, R] [L,R] 区间包含两头的不完整块及中的完整块;中间完整块的最大次数可查 f [ i ] [ j ] f[i][j] f[i][j],不完整块部分暴力枚举,然后查询 g [ a ] [ i ] g[a][i] g[a][i],得出相应数据的最大次数;最后,取两最大值中的最大值。

代码

#include<bits/stdc++.h>
using namespace std;
#define MAXN 60010
#define MAXS 256
const int K = 9; // 块长2^K
int st[MAXS], en[MAXS]; // 块首、尾位置
int n, m, a[MAXN], c[MAXN], f[MAXS][MAXS], g[MAXN][MAXS];
void lsh(){ // 离散化
	for(int i = 1; i <= n; ++i) c[i] = i;
	sort(c+1, c+n+1, [](int x, int y){return a[x] < a[y];});
	int tmp = a[c[1]];
	a[c[1]] = 1;
	for(int i = 2; i <= n; ++i){
		if(a[c[i]] == tmp) a[c[i]] = a[c[i-1]];
		else tmp = a[c[i]], a[c[i]] = a[c[i-1]]+1;
	}
}
void pre_f(){
	int b = n >> K, tmp;
	for(int i = 0; i <= b; ++i){		
		tmp = 0;
		for(int j = i; j <= b; ++j){
			for(int k = st[j]; k <= en[j]; ++k)
				if(tmp < ++c[a[k]]) tmp = c[a[k]];
			f[i][j] = tmp;
		}
		memset(c, 0, sizeof c);
	}
}
void pre_g(){
	memset(g, 0, sizeof g);
	for(int i = 1; i <= n; ++i) g[a[i]][i>>K]++;
	int b = n >> K;
	for(int i = 1; i <= n; ++i){
		for(int j = 1; j <= b; ++j) g[i][j] += g[i][j-1];
	}
}
int ask(int l, int r){
	int L, R, ans = 0;
	L = l >> K, R = r >> K;
	if(L+1 >= R){
		for(int i = l; i <= r; ++i) if(ans < ++c[a[i]]) ans = c[a[i]];
		for(int i = l; i <= r; ++i) c[a[i]] = 0;
		return ans;
	}
	for(int i = l; i <= en[L]; ++i) ++c[a[i]];
	for(int i = st[R]; i <= r; ++i) ++c[a[i]];
	for(int i = l; i <= en[L]; ++i){
		if(c[a[i]] == 0) continue;
		int tmp = c[a[i]]+g[a[i]][R-1]-g[a[i]][L];
		if(ans < tmp) ans = tmp;
		c[a[i]] = 0;
	}
	for(int i = st[R]; i <= r; ++i){
		if(c[a[i]] == 0) continue;
		int tmp = c[a[i]]+g[a[i]][R-1]-g[a[i]][L];
		if(ans < tmp) ans = tmp;
		c[a[i]] = 0;
	}
	ans = max(ans, f[L+1][R-1]);
	return ans;
}
int main(){
	int t;
	scanf("%d", &t);
	while(t--){
		scanf("%d%d", &n, &m);
		for(int i = 1; i <= n; ++i) en[i>>K] = i;
		for(int i = n; i >= 1; --i) st[i>>K] = i;
		for(int i = 1; i <= n; ++i) scanf("%d", a+i);
		lsh(), pre_f(), pre_g();
		int l, r, ans = 0;
		memset(c, 0, sizeof c);
		while(m--) {
			scanf("%d%d", &l, &r);
			ans = ask(l^ans, r^ans);
			printf("%d\n", ans);
		}
	}
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

jpphy0

你的鼓励将是我创作的最大动力

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

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

打赏作者

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

抵扣说明:

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

余额充值