P2894 [USACO08FEB]Hotel G

36 篇文章 0 订阅
20 篇文章 0 订阅

传送门

一个 N N N的数组

支持两种操作

  • 查询 1 1 1 ~ N N N 长度最少为 x x x 连续的 1 1 1
  • 将区间 [ l , r ] [l,r] [l,r] 的值置为 x , x ∈ { 0 , 1 } x,x \in \{0,1\} x,x{0,1}

分析

这道题和【Tunnel Warfare】有点相关

同样是线段树维护

这里有点像,区间连续最大值了。

对于每个位置,维护三个信息 { k , l , r } \{k, l, r\} {k,l,r}

  • k k k 为区间连续的 1 1 1长度
  • l l l 为区间左端点向右走,有多少连续的 1 1 1,也就是连续最后一个 1 1 1的位置 − - 左端点 l + 1 l+1 l+1
  • r r r 为区间右端点往左走,有多少连续的1

假设区间为大小为 1 1 1 { k = 1 , l = 1 , r = 1 } \{k=1, l=1,r=1\} {k=1,l=1,r=1}

推到区间长度为 2 2 2,对于三个信息讨论, l s ls ls为左儿子, r s rs rs为右儿子

  • k = max ⁡ ( t r [ l s ] . k , t r [ r s ] . k ) k = \max(tr[ls].k, tr[rs].k) k=max(tr[ls].k,tr[rs].k)​​, k k k肯定会 > = >= >= 左儿子和右儿子的最大值
  • l = t r [ l s ] . l l=tr[ls].l l=tr[ls].l l l l 肯定会 > = >= >=左儿子的 l l l
  • r = t r [ l s ] . r r=tr[ls].r r=tr[ls].r l l l 肯定会 > = >= >=左儿子的 l l l

但是,考虑最大值 k k k还可能是左右两边拼起来的中间,所以还有取个 max ⁡ ( k , t r [ l s ] . r + t r [ r s ] . l ) \max(k, tr[ls].r+tr[rs].l) max(k,tr[ls].r+tr[rs].l)

如果左边都是连续的 1 1 1,就能接上右边的 l l l了, l + = t r [ r s ] . l l+=tr[rs].l l+=tr[rs].l

同样如果右边都是连续的 1 1 1,接上左边为, r + = t r [ l s ] . r r+=tr[ls].r r+=tr[ls].r

综上

  • k = max ⁡ ( max ⁡ ( t r [ l s ] . k , t r [ r s ] . k ) , t r [ l s ] . r + t r [ r s ] . r ) k = \max(\max(tr[ls].k, tr[rs].k), tr[ls].r+tr[rs].r) k=max(max(tr[ls].k,tr[rs].k),tr[ls].r+tr[rs].r)​​
  • l = t r [ l s ] . l + ( t r [ l s ] . l e n = = t r [ l s ] . l ? t r [ r s ] . s : 0 ) l=tr[ls].l + (tr[ls].len == tr[ls].l?tr[rs].s:0) l=tr[ls].l+(tr[ls].len==tr[ls].l?tr[rs].s:0)
  • r = t r [ l s ] . r + ( t r [ r s ] . l e n = = t r [ r s ] . r ? t r [ l s ] . r : 0 ) r=tr[ls].r+(tr[rs].len==tr[rs].r?tr[ls].r:0) r=tr[ls].r+(tr[rs].len==tr[rs].r?tr[ls].r:0)

对于区间修改 s e t set set,再维护一个标记

查询的话

对于长度至少为 k k k的策略:

从区间 [ 1 , N ] [1, N] [1,N]开始查询,记左端点为 l l l,右端点为 r r r,中点为 m i d mid mid

  • 如果对于区间,左边的连续个数 > = >= >=k,往左边查询 k k k
  • 上面不满足后,如果右区间的 r + r+ r+右区间的 l > = k l>=k l>=k,直接返回 m i d − t r [ l s ] . r + 1 mid-tr[ls].r+1 midtr[ls].r+1
  • 上面不满足,往右区间查询 k k k
  • 如果发现区间的 k k k都小于查询的 k k k,直接返回 0 0 0,因为此区间最大值都小于 k k k,不存在解

以下是代码

//P2894
/*
  @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];

struct Tr {
	int k, l, r, len, lazy;
}tr[MAX_N];

void push_up(int rt) {
	int ls = rt << 1;
	int rs = ls|1;
	tr[rt].l = tr[ls].l;
	tr[rt].r = tr[rs].r;
	if (tr[rt].l == tr[ls].len) {
		tr[rt].l += tr[rs].l;
	}
	if (tr[rt].r == tr[rs].len) {
		tr[rt].r += tr[ls].r;
	}
	tr[rt].k = max(tr[ls].r + tr[rs].l, max(tr[ls].k, tr[rs].k));
}

void calc(int rt, int lazy) {
	if (lazy) {
		tr[rt].k =  tr[rt].l = tr[rt].r = tr[rt].len;
	} else {
		tr[rt].k = tr[rt].l = tr[rt].r = 0;
	}
	tr[rt].lazy = lazy;
}

void push_down(int rt) {
	if (tr[rt].lazy==-1) return;
	calc(rt<<1, tr[rt].lazy);
	calc(rt<<1|1, tr[rt].lazy);
	tr[rt].lazy = -1;
}

void build(int rt, int l, int r) {
	tr[rt].len = r - l + 1;
	tr[rt].lazy = -1;
	if (l == r) {
		tr[rt].k = 1;
		tr[rt].l = 1;
		tr[rt].r = 1;
		return;
	}
	int mid = l + ((r-l)>>1);
	build(rt<<1, l, mid);
	build(rt<<1|1, mid+1, r);
	push_up(rt);
}

void update(int rt, int l, int r, int x, int y, int k) {
	if (x <= l && r <= y) {
		calc(rt, k);
		return;
	}
	push_down(rt);
	int mid = l + ((r-l)>>1);
	if (x <= mid) update(rt<<1, l, mid, x, y, k);
	if (y  > mid) update(rt<<1|1, mid+1, r, x, y, k);
	push_up(rt);
}

int query(int rt, int l, int r, int k) {
	if (tr[rt].k < k) return 0;
	if (l == r) {
		return l;
	}
	push_down(rt);
	int mid = l + ((r-l)>>1);
	if (tr[rt<<1].k >= k) return query(rt<<1, l, mid, k);
	if (tr[rt<<1].r + tr[rt<<1|1].l >= k) return mid - tr[rt<<1].r + 1;
	return query(rt<<1|1, mid+1, r, k);
}

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

void solve(){
	sc("%lld%lld", &N, &M);
	int opt, x, y;
	int pos = 0;
	build(1, 1, N);
	for (int i = 1; i <= M; ++i) {
		sc("%lld", &opt);
		if (opt == 1) {
			sc("%lld", &x);
			pr("%lld\n", pos = query(1, 1, N, x));
			if (pos) update(1, 1, N, pos, pos+x-1, 0);
		} else if (opt == 2) {
			sc("%lld%lld", &x, &y);
			update(1, 1, N, x, x+y-1, 1);
		}
	}
}

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

	return AC;
}

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

Hexrt

客官,请不要给我小费!

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

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

打赏作者

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

抵扣说明:

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

余额充值