2019HDU多校训练第二场 Longest Subarray

题意:给你一个串,问满足以下条件的子串中最长的是多长:对于每个数字,要么在这个子串没出现过,要么出现次数超过k次。

思路:对于这类问题,常常转化为数据结构的询问问题。我们考虑枚举右端点,对于当前右端点,我们单独考虑每一种数的合法区间。假设当前枚举的右端点是i,考虑的数字是c,在右端点左边离i最近的数字c的位置是p1,离i第k远的数字c的位置是p2, 容易发现,数字c的合法区间为[1, p2]和[p1 + 1, i],对应的情况是选择这个数至少k个和不选这个数。那么,如果我们用线段树来维护覆盖的区间,对于每一种数的合法区间在线段树上+1,这样我们只要找到在i前面值为c的最小的位置就是右端点为i的最优解。由于每次右端点只移动1,所以可以在O(logn)时间内维护一个数合法区间的变化。最小位置的找法可以通过维护区间最大值然后在线段树上二分即可。

代码:

#include <bits/stdc++.h>
#define ls (o << 1)
#define rs (o << 1 | 1) 
using namespace std;
const int maxn = 100010;
vector<int> pos[maxn];
int now[maxn];
int n, c, k;
int a[maxn];
struct SegmenTree {
	int mx, lz;
};
SegmenTree tr[maxn * 4];
void pushup(int o) {
	tr[o].mx = max(tr[ls].mx, tr[rs].mx);
}
void pushdown(int o) {
	if(tr[o].lz != 0) {
		tr[ls].mx += tr[o].lz;
		tr[rs].mx += tr[o].lz;
		tr[ls].lz += tr[o].lz;
		tr[rs].lz += tr[o].lz;
		tr[o].lz = 0;
	}
}
void build(int o, int l, int r) {
	tr[o].mx = 0;
	tr[o].lz = 0;
	if(l == r) {
		return;
	}
	int mid = (l + r) >> 1;
	build(ls, l, mid);
	build(rs, mid + 1, r);
	pushup(o);
}
void update(int o, int l, int r, int ql, int qr, int val) {
	if(ql > qr) return;
	if(l >= ql && r <= qr) {
		tr[o].lz += val;
		tr[o].mx += val;
		return;
	}
	pushdown(o);
	int mid = (l + r) >> 1;
	if(ql <= mid) update(ls, l, mid, ql, qr, val);
	if(qr > mid) update(rs, mid + 1, r, ql, qr, val);
	pushup(o); 
}
int query(int o, int l, int r) {
	if(l == r) return l;
	int ans = -1, mid = (l + r) >> 1;
	pushdown(o);
	if(tr[ls].mx == c) ans = query(ls, l, mid);
	else if(tr[rs].mx == c) ans = query(rs, mid + 1, r);
	return ans;
}
int main() {
	while(~scanf("%d%d%d", &n, &c, &k)) {
		for (int i = 1; i <= c; i++) {
			pos[i].clear();
			pos[i].push_back(0);
		}
		for (int i = 1; i <= n; i++) {
			scanf("%d", &a[i]);
			pos[a[i]].push_back(i);
		}
		build(1, 1, n);
		for (int i = 1; i <= c; i++) {
			pos[i].push_back(n + 1);
			update(1, 1, n, pos[i][0] + 1, pos[i][1] - 1, 1);
			now[i] = 0;
		}
		int ans = 0;
		for (int i = 1; i <= n; i++) {
			int t = a[i];
			update(1, 1, n, pos[t][now[t]] + 1, pos[t][now[t] + 1] - 1, -1);
			if(now[t] >= k) update(1, 1, n, 1, pos[t][now[t] - k + 1], -1);
			now[t]++;
			update(1, 1, n, pos[t][now[t]] + 1, pos[t][now[t] + 1] - 1, 1);
			if(now[t] >= k) update(1, 1, n, 1, pos[t][now[t] - k + 1], 1);
			int tmp = query(1, 1, n);
			if(tmp != -1)
				ans = max(ans, i - tmp + 1);
		}
		printf("%d\n", ans);
	}
} 

  

转载于:https://www.cnblogs.com/pkgunboat/p/11242904.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值