公园遛狗 / 小白逛公园【ybt高效进阶4-4-3】【luogu P4513】【题解未完成】
题目大意:
给你一个序列,要维护两个操作。
单点修改和在一个区间中找它的最大子段和。
思路:
这题是一个比较经典的线段树做法。
在线段树中我们维护四个值:
v
a
l
(
x
)
=
该
区
间
的
区
间
和
val(x)=该区间的区间和
val(x)=该区间的区间和
l
s
(
x
)
=
该
区
间
包
含
左
端
点
的
最
大
子
段
和
ls(x)=该区间包含左端点的最大子段和
ls(x)=该区间包含左端点的最大子段和
r
s
(
x
)
=
该
区
间
包
含
左
端
点
的
最
大
子
段
和
rs(x)=该区间包含左端点的最大子段和
rs(x)=该区间包含左端点的最大子段和
l
e
n
(
x
)
=
区
间
的
最
大
子
段
和
len(x)=区间的最大子段和
len(x)=区间的最大子段和
那么我们有:
{
v
a
l
(
x
)
=
v
a
l
(
x
<
<
1
)
+
v
a
l
(
x
<
<
1
∣
)
l
s
(
x
)
=
m
a
x
(
l
s
(
x
<
<
1
)
,
v
a
l
(
x
<
<
1
)
+
l
s
(
x
<
<
1
∣
1
)
)
r
s
(
x
)
=
m
a
x
(
r
s
(
x
<
<
1
∣
1
)
,
v
a
l
(
x
<
<
1
∣
1
)
+
r
s
(
x
<
<
1
)
)
l
e
n
(
x
)
=
m
a
x
(
l
e
n
(
x
<
<
1
)
,
l
e
n
(
x
<
<
1
∣
1
)
,
r
s
(
x
<
<
1
)
+
l
s
(
x
<
<
1
∣
1
)
)
\left\{\begin{matrix} val(x)=val(x<<1)+val(x<<1|)& \\ ls(x)=max(ls(x<<1),val(x<<1)+ls(x<<1|1))& \\ rs(x)=max(rs(x<<1|1),val(x<<1|1)+rs(x<<1))& \\ len(x)=max(len(x<<1),len(x<<1|1),rs(x<<1)+ls(x<<1|1))& \end{matrix}\right.
⎩⎪⎪⎨⎪⎪⎧val(x)=val(x<<1)+val(x<<1∣)ls(x)=max(ls(x<<1),val(x<<1)+ls(x<<1∣1))rs(x)=max(rs(x<<1∣1),val(x<<1∣1)+rs(x<<1))len(x)=max(len(x<<1),len(x<<1∣1),rs(x<<1)+ls(x<<1∣1))
该段代码:
val(x)=val(x<<1)+val(x<<1|1); //区间和=左区间和+右区间和
ls(x)=max(ls(x<<1),val(x<<1)+ls(x<<1|1)); //=max(左区间左最大子段和,左区间和+右区间左最大子段和)
rs(x)=max(rs(x<<1|1),val(x<<1|1)+rs(x<<1));//同上
len(x)=max(max(len(x<<1),len(x<<1|1)),ls(x<<1|1)+rs(x<<1));
//=max(左区间最大子段和,右区间最大子段和,左区间右最大子段和+右区间左最大子段和)
在查询答案的时候,我们要像更新答案一样去查询,代码与修改神似,具体注释见代码:
代码:
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#define r register
#define rep(i,x,y) for(r ll i=x;i<=y;++i)
#define per(i,x,y) for(r ll i=x;i>=y;--i)
using namespace std;
typedef long long ll;
const ll V=5*1e5+10;
ll a[V],n,m,k,x,y;
struct node
{
ll ls,rs,val,len;
#define ls(x) tree[x].ls //为减少码量的无奈之举
#define rs(x) tree[x].rs
#define val(x) tree[x].val
#define len(x) tree[x].len
}tree[V<<3];
void update(ll x) //更新答案的函数
{
val(x)=val(x<<1)+val(x<<1|1);
ls(x)=max(ls(x<<1),val(x<<1)+ls(x<<1|1));
rs(x)=max(rs(x<<1|1),val(x<<1|1)+rs(x<<1));
len(x)=max(max(len(x<<1),len(x<<1|1)),ls(x<<1|1)+rs(x<<1));
}
void build(ll now,ll le,ll ri) //初始建树
{
if(le==ri)
{
ls(now)=rs(now)=val(now)=len(now)=a[le];
return ;
}
ll mid=le+ri>>1;
build(now<<1,le,mid);
build(now<<1|1,mid+1,ri);
update(now);
}
void add(ll now,ll le,ll ri,ll p,ll v) //修改,思路见上面
{
if(le==ri)
{
ls(now)=rs(now)=val(now)=len(now)=v;
return ;
}
ll mid=le+ri>>1;
if(p<=mid) add(now<<1,le,mid,p,v);
else add(now<<1|1,mid+1,ri,p,v);
update(now);
}
node ask(ll now,ll le,ll ri,ll x,ll y)
{
if(le>=x&&ri<=y) return tree[now];
ll mid=le+ri>>1;
node res;
if(mid>=y) return ask(now<<1,le,mid,x,y); //都在左区间
else if(mid+1<=x) return ask(now<<1|1,mid+1,ri,x,y); //都在右区间
//在两个区间之间
node r1=ask(now<<1,le,mid,x,y); //递归求解
node r2=ask(now<<1|1,mid+1,ri,x,y);
res.val=r1.val+r2.val;
res.ls=max(r1.ls,r2.ls+r1.val);//就像修改一样更新数据,比较大小
res.rs=max(r2.rs,r1.rs+r2.val);
res.len=max(max(r1.len,r2.len),r1.rs+r2.ls);
return res;
}
int main()
{
scanf("%lld%lld",&n,&m);
rep(i,1,n)
scanf("%lld",&a[i]);
build(1,1,n);
rep(i,1,m)
{
scanf("%lld%lld%lld",&k,&x,&y);
if(k==2) add(1,1,n,x,y);
else
{
if(x>y) swap(x,y);
node ans=ask(1,1,n,x,y);
printf("%lld\n",ans.len);
}
}
return 0;
}