LOJ6285 数列分块入门9


LOJ6285 数列分块入门 9


标签

  • 基础分块

前言


简明题意

  • 给一个序列,需要查询区间中的最小众数

思路

  • 首先要会写暴力。如果数据范围是1000但是有多组询问你会怎么写?先解决这个问题。我们可以开一个dp[][]数组,dp[i][j]表示区间[i,j]中的众数。具体过程是枚举所有的左端点,然后再枚举右端点,每次右端点右移就更新一下答案。
  • 暴力会写了,就考虑一下这题。实际上关于众数有这样一个性质:如果已知集合a的众数是x,那么集合a并上集合b的众数要么是x,要么是b中的一个数。那么我们就可以愉快的分块了,先根据上面所说的预处理出dp[i][j]第i块到第j块的众数,那么每次查询,答案要么是中间那一整块的答案,要么是左右两边不整块的那些数。
  • 所以现在问题就在求左右两边不整块的数的个数。这又是一个经典的二分问题了。直接开一个vector[]记录每一个数出现的所有位置,暴力枚举这些数,然后二分找到r和l的位置相减就可以了。
  • 最后就是数有些大,可以离散化一下。
  • 还有就是块的大小开成80才能AC

注意事项


总结

  • 如果已知集合a的众数是x,那么集合a并上集合b的众数要么是x,要么是b中的一个数
  • 求区间中数出现的次数,可以记录位置然后二分。

AC代码

#pragma GCC optimize(2)
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<vector>
#include<cstring>
#include<unordered_map>
using namespace std;

const int maxn = 1e5 + 10;

int read()
{
	int x = 0, f = 1; char ch = getchar();
	while (ch<'0' || ch>'9') { if (ch == '-')f = -1; ch = getchar(); }
	while (ch >= '0'&&ch <= '9') { x = x * 10 + ch - '0'; ch = getchar(); }
	return x * f;
}

int n, a[maxn];
int pos[maxn], len;
pair<int, int> dp[5000][5000];
vector<int> rec[maxn];//记录每个数出现的所有位置

vector<int> ls;
int mp[maxn];

int ask(int l, int r)
{
	int cur_max = dp[pos[l] + 1][pos[r] - 1].first, max_val = dp[pos[l] + 1][pos[r] - 1].second;

	for (int i = l; i <= min(pos[l] * len, r); i++)
	{
		auto& b = rec[a[i]];
		int cnt = upper_bound(b.begin(), b.end(), r) - lower_bound(b.begin(), b.end(), l);
		if (cnt > cur_max || (cnt == cur_max && a[i] < max_val))
			cur_max = cnt, max_val = a[i];
	}

	if (pos[l] != pos[r])
	{
		for (int i = pos[r] * len - len + 1; i <= r; i++)
		{
			auto& b = rec[a[i]];
			int cnt = upper_bound(b.begin(), b.end(), r) - lower_bound(b.begin(), b.end(), l);
			if (cnt > cur_max || (cnt == cur_max && a[i] < max_val))
				cur_max = cnt, max_val = a[i];
		}
	}

	return max_val;
}

void solve() 
{
	scanf("%d", &n);
	len = 80;
	for (int i = 1; i <= n; i++)
		a[i] = read(), pos[i] = (i - 1) / len + 1, ls.push_back(a[i]);

	sort(ls.begin(), ls.end());
	int ls_len = unique(ls.begin(), ls.end()) - ls.begin();

	for (int i = 1; i <= n; i++)
	{
		int t = lower_bound(ls.begin(), ls.begin() + ls_len, a[i]) - ls.begin() + 1;
		mp[t] = a[i];
		a[i] = t;
		rec[a[i]].push_back(i);
	}

	for (auto& it : rec)
		sort(it.begin(), it.end());

	//预处理dp
	int cnt[maxn];//记录一下每个数出现的次数
	for (int i = 1; i <= pos[n]; i++)//枚举起始块
	{
		int cur_max = -1e9, max_val;//分别是当前最多出现的次数、最多出现次数的那个数
		memset(cnt, 0, sizeof cnt);
		for (int j = i * len - len + 1; j <= n; j++)//枚举从第i块开始的每一个数
		{
			cnt[a[j]]++;
			if (cnt[a[j]] > cur_max || (cnt[a[j]] == cur_max && a[j] < max_val))//a[j]出现的次数多,或者a[j]出现的次数一样但更小,就更新
				cur_max = cnt[a[j]], max_val = a[j];
			dp[i][pos[j]].first = cur_max, dp[i][pos[j]].second = max_val;
		}
	}


	for (int i = 1; i <= n; i++) 
	{
		int l, r;
		l = read(), r = read();

		printf("%d\n", mp[ask(l, r)]);
	}
}

int main() {
	freopen("Testin.txt", "r", stdin);
	freopen("Testout.txt", "w", stdout);
	solve();
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值