Hdu3308--LCIS(线段树区间合并)

题意:

给你n个数,m个操作。操作有两种:
1.U x y 将数组第x位变为y (下标从0开始)
2. Q x y 问数组第x位到第y位连续最长子序列的长度。
对于每次询问,输出一个答案

思路:

线段树单点修改区间合并
对于每个区间节点维护3个值:
ls[root]:以区间左端点开始的最长连续上升序列长度(前缀序列)
rs[root]:以区间右端点结束的最长连续上升序列长度(后缀序列)
Tree[root]:该区间最终的LCIS
区间[l, r]的LCIS序列可能存在于左半区间[l, mid],也可能在右半区间[mid+1, r],还可能跨左右区间,在中间。

如果原数组A[mid] >= A[mid+1],说明LCIS序列不可能在中间。
Tree[root] = max(Tree[root], Tree[root<<1|1);

如果A[mid] < A[mid+1],说明LCIS序列可能在中间,
Tree[root] = max(前缀序列, 后缀序列,中间序列),此时还要注意:如果当前区间中所有的值都符合上升序列,那么更新它的父节点时,它可以加上其他区间的边界值而进行扩充,比如:[l, mid]区间全部符合LCIS,如果此时A[mid] < A[mid+1],那么[l, r]的前缀序列就不是[l, mid]的前缀序列了,而是再要加上右区间[mid+1, r]的前缀。

#include <cstdio>
#include <cstring>
#include <algorithm>
#define lson root<<1, l, mid
#define rson root<<1|1, mid+1, r
using namespace std;
const int INF = 0x3f3f3f3f;
const int maxn = 1e5+5;
int A[maxn], Tree[maxn<<2], ls[maxn<<2], rs[maxn<<2];	// Tree[]:总的LCIS;ls[]:前缀的LCIS;rs[]:后缀的LCIS 
void push_up(int root, int l, int r){
	ls[root] = ls[root<<1]; rs[root] = rs[root<<1|1];
	Tree[root] = max(Tree[root<<1], Tree[root<<1|1]);
	int mid = (l+r) >> 1;
	if(A[mid] < A[mid+1]){	// LCIS可能在中间 
		if(ls[root] == mid-l+1) ls[root] += ls[root<<1|1];
		if(rs[root] == r-mid) rs[root] += rs[root<<1];
		Tree[root] = max(Tree[root], rs[root<<1]+ls[root<<1|1]);
	}
}
void Stree_build(int root, int l, int r){
	if(l == r){
		Tree[root] = ls[root] = rs[root] = 1;
		return;
	}
	int mid = (l+r) >> 1;
	Stree_build(lson);
	Stree_build(rson);
	push_up(root, l, r); 
}

void update(int k, int l, int r, int root){
	if(l == r) return;
	int mid = (l+r) >> 1;
	if(k <= mid)
		update(k, l, mid, root<<1);
	else
		update(k, mid+1, r, root<<1|1);
	push_up(root, l, r); 
}

int Query(int la, int rb, int l, int r, int root){
	if(la <= l&&rb >= r) return Tree[root];
	if(l == r) return 0;
	int mid = (l+r) >> 1;
	if(rb <= mid) return Query(la, rb, l, mid, root<<1);
	if(la > mid) return Query(la, rb, mid+1, r, root<<1|1);
	int lmax = Query(la, rb, l, mid, root<<1);
	int rmax = Query(la, rb, mid+1, r, root<<1|1);
	int ans = max(lmax, rmax);
	if(A[mid] < A[mid+1]){
		ans = max(ans, min(mid-la+1, rs[root<<1]) + min(rb-mid, ls[root<<1|1]));
	}
	return ans;
}


int main()
{
	//freopen("in.txt","r",stdin);
	int n, m, a,b;
	char op[5];
	int T; scanf("%d", &T); 
	while(T--){
		scanf("%d%d", &n,&m);
		for(int i = 1; i <= n; ++i) scanf("%d", &A[i]);
		Stree_build(1, 1, n);
		while(m--){
			scanf("%s", op);
			scanf("%d%d",&a,&b);
			++a;
			if(op[0] == 'U'){
				A[a] = b;
				update(a, 1, n, 1);
			}
			else{
				++b;
				printf("%d\n", Query(a, b, 1, n, 1));
			}
		}
	}

	return 0;
}


  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值