HDU 3308 LCIS (经典区间合并)【线段树】

时间:2023-03-09 16:36:37
HDU 3308 LCIS (经典区间合并)【线段树】

<题目链接>

题目大意:

给你一段序列,对其进行两种操作,一是修改某个序号的点的值;二是查询某个区间的LCIS(最长上升子序列)。

解题分析:

线段树区间合并的典型例题,用求某个区间的LCIS时,需要比较三个值,一是左区间的LCIS,二是右区间的LCIS,三是左右子区间合并的LCIS。最重要的是第三点如何实现,实现第三点需要维护一个最长后缀上升子序列和最长前缀上升子序列,总之,相对于一般的线段树,区间合并需要对Pushup()函数进行一些改动,query()的时候也要记得对三种情况进行讨论。

 #include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int siz=;
int a[siz],sum[siz<<],lsum[siz<<],rsum[siz<<]; void pushup(int rt,int l,int r) { //区间合并,更新该节点对应的最长上升子序列和最长前、后缀
lsum[rt]=lsum[rt<<]; //最长前缀
rsum[rt]=rsum[rt<<|]; //最长后缀
sum[rt]=max(sum[rt<<],sum[rt<<|]);
/*-- 先更新该节点需要维护的三个值 --*/ int m=(l+r)>>;
if(a[m]<a[m+]){ //如果左右子区间的LCIS能够合并,那么就进一步更新这三个点 if(lsum[rt<<]==(m-l+)) //如果左子区间的最长前缀==左子区间长度
lsum[rt]+=lsum[rt<<|]; //那么该节点的最长前缀除左子区间的最长前缀以外,还要加上右子区间的最长前缀
if(rsum[rt<<|]==(r-m)) //如果右子区间的最长后缀==右子区间长度
rsum[rt]+=rsum[rt<<]; //那么该节点的最长后缀除右子区间的最长后缀外,还要加上左子区间的最长后缀
/*-- 更新rt节点对应区间的最长前、后缀 --*/ sum[rt]=max(sum[rt],rsum[rt<<]+lsum[rt<<|]);
//更新rt节点对应区间的最长上升子序列
}
}
void build(int l,int r,int rt){
if(l==r){ //注意这里的初始化
sum[rt]=;
lsum[rt]=rsum[rt]=;
return;
}
int m=(l+r)>>;
build(l,m,rt<<);
build(m+,r,rt<<|);
pushup(rt,l,r);
}
void update(int loc,int val,int l,int r,int rt){ //单点更新
if(l==r){
a[l]=val;
return;
}
int m=(l+r)>>;
if(loc<=m)
update(loc,val,l,m,rt<<);
else
update(loc,val,m+,r,rt<<|);
pushup(rt,l,r);
}
int query(int L,int R,int l,int r,int rt){
if(L<=l&&r<=R)
return sum[rt];
int m=(l+r)>>;
int ans=; //注意这里ans初始化为1 /*-- 答案一共三种可能,在左区间或右区间或横跨两个区间 --*/
if(L<=m) //左、右区间的LCIS
ans=max(ans,query(L,R,l,m,rt<<));
if(R>m)
ans=max(ans,query(L,R,m+,r,rt<<|)); if(L<=m&&R>=m&&a[m]<a[m+]) //横跨左右两个子区间的情况
ans=max(ans,min(m-L+,rsum[rt<<])+min(R-m,lsum[rt<<|])); //rsum代表区间最长后缀、lsum为区间最长前缀
return ans;
} int main(){ //lsum为区间左端点开始长度
char c; //rsum为区间右端点开始长度
int t,n,m; //sum为区间最长长度
scanf("%d",&t);
while(t--){
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++)
scanf("%d",&a[i]);
build(,n,);
while(m--){
int u,v;
cin>>c>>u>>v;
if(c=='Q'){
u++,v++; //因为题目是从0开始编号
printf("%d\n",query(u,v,,n,));
}
else{
u++;
update(u,v,,n,);
}
}
}
return ;
}

2018-09-11