更多文章可以在本人的个人小站:https://kaiserwilheim.github.io 查看。
转载请注明出处。
什么是线段树?
线段树是一种数据结构,其可以以
O
(
log
n
)
O( \log n)
O(logn) 的时间复杂度处理维护区间信息的问题。
其维护的信息满足可以快速合并两个区间的信息(例如区间和、区间最大(小)值。
线段树同时可以支持区间修改,但是区间修改时必须能够快速求出每个区间维护的信息。
线段树的结构
线段树在建立的时候,会将每一个区间长度不为1的区间进行二分,并递归对这两个区间进行判定。
这是一棵区间长度为12的线段树:
我们可以明显看到,线段树是一棵二叉树。
线段树的高度一般是
log
n
\log n
logn 的,这就需要我们建立
2
⌈
log
n
⌉
+
1
2^{\lceil \log n \rceil + 1}
2⌈logn⌉+1 个节点,这将会耗费我们大量的空间。
一般情况下,我们为了访问方便,也为了节省空间,我们一般将一个节点(设其编号为 p p p)的左儿子和右儿子分别设为 p < < 1 p << 1 p<<1 和 p < < 1 ∣ 1 p << 1 | 1 p<<1∣1,就相当于是在一棵满二叉树上面访问。
于是我们就给出编号:
下面给出示例代码:
假设我们维护的是区间和,并使用类似
struct SegTree
{
int l, r;
int sum;
}tr;
的结构体来存储数据,那么建立线段树的函数就是像下面这个样子的:
void build(int p, int l, int r)
{
tr[p].l = l, tr[p].r = r;
if(l == r)
{
tr[p].sum = a[l];
return;
}
int mid = (l + r) >> 1;
build(p << 1, l, mid);
build(p << 1 | 1, mid + 1, r);
tr[p].sum = tr[p << 1].sum + tr[p << 1 | 1].sum;
return;
}
线段树的区间询问
我们这里以区间求和为例。
我们先手动模拟一下区间求和的过程:
假设我们对这样一个长度为12的序列 { 1 , 1 , 4 , 5 , 1 , 4 , 1 , 9 , 1 , 9 , 8 , 1 } \lbrace 1,1,4,5,1,4,1,9,1,9,8,1 \rbrace {1,1,4,5,1,4,1,9,1,9,8,1} 建立线段树,那么我们每一个区间的数据是这个样子的:
我们假设需要对 KaTeX parse error: Undefined control sequence: \[ at position 1: \̲[̲ 2,8 \] 这个区间内的数据进行求和,那么我们覆盖的区间就是这个样子的:
我们想象一下将区间向上推:
最后我们将 KaTeX parse error: Undefined control sequence: \[ at position 1: \̲[̲ 2,8 \] 这个大区间分为 KaTeX parse error: Undefined control sequence: \[ at position 1: \̲[̲ 2,2 \] \cup \[… 这几个小区间。
我们每一次二分递归的时候,如果当前区间被我们询问的区间完全包含就返回当前区间的值,如果没有完全包含就递归询问左右区间,并输出两个区间的返回值之和。
在递归询问左右两区间的时候需要注意一下子区间是否与询问的区间有重合部分,如果两者交集为空,那就不需要询问了。
下面给出示例代码:
int segadd(int p, int l, int r)
{
if(tr[p].l >= l && tr[p].r <= r)return tr[p].sum;
int mid = (tr[p].l + tr[p].r) >> 1;
int res = 0;
if(l <= mid)res += segsum(p << 1, l, r);
if(r > mid)res += segsum(p << 1 | 1, l, r);
return res;
}
线段树的区间修改
区间修改和区间求和有着相同的思路。
我们这里以区间加为例。
还是之前的例子,如果我们想要给区间 KaTeX parse error: Undefined control sequence: \[ at position 1: \̲[̲ 2,8 \] 同时加上一个数会怎样?
让我们再来看一下这个图:
其中白色的区间是完全和我们询问没有任何关系的区间,黄色的区间是与我们询问的区间有交集,但是没有被新闻区间所包含的区间。
红色的区间和橙色的区间都被询问区间所包含,其中红色的区间是我们在区间求和的时候返回了值的区间,橙色的区间虽然被询问区间包含,但是我们没有访问到他们。
区间修改也是一样的思路。我们只需要将需要的区间修改即可,剩下的被询问区间包含的区间只需要标记一下就可以了。
这里需要用到懒标记的思想。
懒标记
懒标记的意思就是很懒的标记。我们用到它的时候才会将它兑现,否则他就在那里静静地待着就可以了。
我们考虑给结构体里面加上一个懒标记,就像这样:
struct SegTree
{
int l, r;
int sum;
int tag;//<- 就是这个
}tr;
当某个区间有懒标记的时候,就意味着它的子树里面的所有区间都应该被修改,但是目前没有应用上。
但是拥有懒标记的区间的值在修改的时候已经被修改好了。
每一次我们递归询问子区间的时候,我们就需要将这个懒标记下传(如果有的话)。
我们需要一个pushdown
函数来完成这件事。
我们这里仍然以区间加为例:
void pushdown(int p)
{
auto &root = tr[p], &left = tr[p << 1], &rght = tr[p << 1 | 1];
if(root.tag)
{
left.sum += root.tag * (left.r - left.l + 1);
rght.sum += root.tag * (rght.r - rght.l + 1);
//更新子区间值
left.tag += root.tag;
rght.tag += root.tag;
//更新子区间懒标记
root.tag = 0;//懒标记清空
}
return;
}
例题
例题1:
#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int N = 100010;
struct SegTree
{
int l, r;
ll sum, tag;
}tr[N << 3];
ll a[N];
void pushup(int p)
{
tr[p].sum = tr[p << 1].sum + tr[p << 1 | 1].sum;
}
void pushdown(int p)
{
auto &root = tr[p], &left = tr[p << 1], &rght = tr[p << 1 | 1];
if(root.tag)
{
left.sum += root.tag * (left.r - left.l + 1);
rght.sum += root.tag * (rght.r - rght.l + 1);
left.tag += root.tag;
rght.tag += root.tag;
root.tag = 0;
}
return;
}
void build(int p, int l, int r)
{
tr[p].l = l, tr[p].r = r;
if(l == r)
{
tr[p].sum = a[l];
return;
}
int mid = (l + r) >> 1;
build(p << 1, l, mid);
build(p << 1 | 1, mid + 1, r);
pushup(p);
}
void segadd(int p, int l, int r, int k)
{
if(tr[p].l >= l && tr[p].r <= r)
{
tr[p].tag += k;
tr[p].sum += k * (tr[p].r - tr[p].l + 1);
return;
}
pushdown(p);
int mid = (tr[p].l + tr[p].r) >> 1;
if(l <= mid)segadd(p << 1, l, r, k);
if(r > mid)segadd(p << 1 | 1, l, r, k);
pushup(p);
}
ll segsum(int p, int l, int r)
{
if(tr[p].l >= l && tr[p].r <= r)return tr[p].sum;
pushdown(p);
ll res = 0;
int mid = (tr[p].l + tr[p].r) >> 1;
if(l <= mid)res += segsum(p << 1, l, r);
if(r > mid)res += segsum(p << 1 | 1, l, r);
return res;
}
int main()
{
int n, m;
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i++)scanf("%lld", &a[i]);
build(1, 1, n);
while(m--)
{
int op, x, y, k;
scanf("%d%d%d", &op, &x, &y);
if(op == 1)
{
scanf("%d", &k);
segadd(1, x, y, k);
}
else if(op == 2)
{
printf("%lld\n", segsum(1, x, y));
}
}
return 0;
}
例题2:
#include <bits/stdc++.h>
using namespace std;
int n, m, d;
const int N = 100010;
int a[N];
struct segtree
{
int l, r;
long long v, ta, tm;
}node[N << 2];
void build(int p, int l, int r)
{
node[p].l = l, node[p].r = r, node[p].ta = 0, node[p].tm = 1;
if(l == r)
{
node[p].v = a[r];
return;
}
int mid = (l + r) >> 1;
build(p << 1, l, mid);
build(p << 1 | 1, mid + 1, r);
node[p].v = (node[p << 1].v + node[p << 1 | 1].v) % d;
return;
}
void tag(int p)
{
auto &root = node[p], &left = node[p << 1], &rght = node[p << 1 | 1];
if((root.ta != 0) || (root.tm != 1))
{
left.tm = (left.tm * root.tm) % d, left.ta = (left.ta * root.tm + root.ta) % d;
rght.tm = (rght.tm * root.tm) % d, rght.ta = (rght.ta * root.tm + root.ta) % d;
left.v = (left.v * root.tm + root.ta * (left.r - left.l + 1)) % d;
rght.v = (rght.v * root.tm + root.ta * (rght.r - rght.l + 1)) % d;
root.ta = 0, root.tm = 1;
}
return;
}
void segadd(int p, int l, int r, int k)
{
if((node[p].r <= r) && (node[p].l >= l))
{
node[p].ta = (node[p].ta + k) % d;
node[p].v = (node[p].v + k * (node[p].r - node[p].l + 1)) % d;
return;
}
if((node[p].r < l) || (node[p].l > r))return;
tag(p);
int mid = (node[p].r + node[p].l) >> 1;
if(mid >= l)segadd(p << 1, l, r, k);
if(mid < r)segadd(p << 1 | 1, l, r, k);
node[p].v = (node[p << 1].v + node[p << 1 | 1].v) % d;
return;
}
void segmul(int p, int l, int r, int k)
{
if((node[p].r <= r) && (node[p].l >= l))
{
node[p].tm = (node[p].tm * k) % d;
node[p].ta = (node[p].ta * k) % d;
node[p].v = (node[p].v * k) % d;
return;
}
if((node[p].r < l) || (node[p].l > r))return;
tag(p);
int mid = (node[p].r + node[p].l) >> 1;
if(mid >= l)segmul(p << 1, l, r, k);
if(mid < r)segmul(p << 1 | 1, l, r, k);
node[p].v = (node[p << 1].v + node[p << 1 | 1].v) % d;
return;
}
long long segsum(int p, int l, int r)
{
if((node[p].r <= r) && (node[p].l >= l))return node[p].v;
if((node[p].r < l) || (node[p].l > r))return 0;
tag(p);
int mid = (node[p].r + node[p].l) >> 1;
long long ans = 0;
if(mid >= l)ans = ans + segsum(p << 1, l, r);
if(mid < r)ans = ans + segsum(p << 1 | 1, l, r);
return ans % d;
}
int main()
{
scanf("%d%d%d", &n, &m, &d);
for(int i = 1; i <= n; i++)scanf("%d", &a[i]);
build(1, 1, n);
for(int i = 1; i <= m; i++)
{
int op;
scanf("%d", &op);
if(op == 1)
{
int x, y, k;
scanf("%d%d%d", &x, &y, &k);
segmul(1, x, y, k);
}
else if(op == 2)
{
int x, y, k;
scanf("%d%d%d", &x, &y, &k);
segadd(1, x, y, k);
}
else if(op == 3)
{
int x, y;
scanf("%d%d", &x, &y);
printf("%lld\n", segsum(1, x, y));
}
}
}
扩展
标记永久化
我们可以尝试一个思路:不下传标记。
我们重新思考一下懒标记的意义。
懒标记代表的是,这个点的所有子区间都需要被这个标记更新他们的值。
于是我们就可以将标记放进函数里面,随着函数下传,到时候也可以起到更新值的作用。
这样的方法适用于不能pushdown
的时候。
具体情况我也没有遇见过,据说是需要配合可持久化权值线段树(aka主席树)使用。
示例代码:
ll segsum(int p, int l, int r, int tag)
{
if(tr[p].l >= l && tr[p].r <= r)return tr[p].sum + tag * (tr[p].r - tr[p].l + 1);
ll res = 0;
int mid = (tr[p].l + tr[p].r) >> 1;
if(l <= mid)res += segsum(p << 1, l, r, tag + tr[p].tag);
if(r > mid)res += segsum(p << 1 | 1, l, r, tag + tr[p].tag);
return res;
}
线段树上二分
想不到吧,线段树还能优化二分!
我们考虑这样一个问题:我们需要找到一个位置 p o s pos pos,使得 KaTeX parse error: Undefined control sequence: \[ at position 1: \̲[̲ 1,pos \] 这个区间内数字的和不小于 v a l val val。保证所有权值非负。
那么我们可以直接进行二分,这样每做一次的复杂度是 O ( log 2 n ) O(\log^2 n) O(log2n) 的。
我们可以用线段树来优化这样的过程。
我们带着
v
a
l
val
val 在线段树上搜索。
如果当前节点的左儿子的权值
≥
v
a
l
\geq val
≥val,就带着
v
a
l
val
val 继续跳到左儿子搜索;
如果当前节点的左儿子的权值
<
v
a
l
< val
<val,就带着
v
a
l
−
t
r
[
t
r
[
p
]
.
l
s
]
.
v
a
l
val-tr[tr[p].ls].val
val−tr[tr[p].ls].val 继续跳到右儿子搜索。
动态开点线段树
有时候我们需要维护的下标区域非常广(比如说 1 0 9 10^9 109),但是真正有值的地方其实很少。
我们当然不可能为这个序列直接创建一棵完整的线段树,这样子会耗费我们大量的空间维护许多无用的数据。
比如说维护一个序列里面有多少个数字大于某个数,但是这个数列的值域是 1 0 9 10^9 109,但是数列长度最多只有 1 0 5 10^5 105。
我们当然不可能直接创建一个覆盖整个值域的线段树,这样子会需要我们创建 2 3 1 2^31 231 个节点,即使我们不考虑区间修改,只维护一个区间和的话也需要高达24GiB的空间,十分的浪费。
于是我们可以动态开点,只建立用到的区间,其他的就不会被声明,大大减少了我们需要用的空间。
或者说,我们什么时候需要用到某个节点了,我们就建立它。
但是这样做的代价就是我们无法再使用简单的 p<<1
和 p<<1|1
来访问当前节点的左右子节点了。
我们需要记录下当前节点的左右节点的编号。
我们以维护一个序列里面有多少个数字大于某个数为例。
那么我们的结构体就会变成这个样子:
struct SegTree
{
int l, r;
int ls, rs;
int sum;
}tr;
其中 ls
和 rs
分别代表了当前节点的左儿子和右儿子的编号。
那么我们就不需要build
了,只需要一开始将第一个节点初始化为根节点即可。
然后我们的区间修改函数就是这个样子的。
我们可以有两种风格:
- 提前开点
我们在访问子节点之前,首先看一下我们将要访问的这个子节点到底是否存在,如果不存在的话就新开一个。
这样子的话支持把当前区间的左右端点放在结构体里面存着,也支持将区间左右端点当做参数放到函数里面传。
void segadd(int p, int l, int r, int k)
{
if(tr[p].l >= l && tr[p].r <= r)
{
tr[p].sum += k * (tr[p].r - tr[p].l + 1);
tr[p].add += k;
return;
}
int mid = (tr[p].l + tr[p].r) >> 1;
if(l <= mid)
{
if(!tr[p].ls)
{
tr[p].ls = ++idx;
tr[tr[p].ls].l = tr[p].l;
tr[tr[p].ls].r = mid;
}
segadd(tr[p].ls, l, r, k);
}
if(r > mid)
{
if(!tr[p].rs)
{
tr[p].rs = ++idx;
tr[tr[p].rs].l = mid + 1;
tr[tr[p].rs].r = tr[p].r;
}
segadd(tr[p].rs, l, r, k);
}
tr[p].sum = tr[tr[p].ls].sum + tr[tr[p].rs].sum;
return;
}
- 用到再开
我们在访问子节点的时候才看一下,我们的这个子节点是否存在。如果不存在的话就创建一个。
这样的话我们就在函数里面传节点编号的引用,并且同时把当前区间的左右端点放进函数里面当参数传。
void segadd(int &p, int l, int r, int al, int ar, int k)
{
if(!p)p = ++idx;
if(al >= l && ar <= r)
{
tr[p].sum += k * (ar - al + 1);
tr[p].add += k;
return;
}
int mid = (al + ar) >> 1;
if(l <= mid)segadd(tr[p].ls, l, r, al, mid, k);
if(r > mid)segadd(tr[p].rs, l, r, mid + 1, ar, k);
tr[p].sum = tr[tr[p].ls].sum + tr[tr[p].rs].sum;
return;
}
而查询函数不需要做太大的变动,只需要看一下当前访问的节点是否存在就可以了。
int segsum(int p, int l, int r)
{
if(!p)return 0;
if(tr[p].l >= l && tr[p].r <= r)return tr[p].sum;
int mid = (tr[p].l + tr[p].r) >> 1;
int res = 0;
if(l <= mid)res += segsum(tr[p].ls, l, r);
if(r > mid)res += segsum(tr[p].rs, l, r);
return res;
}
可持久化权值线段树
面对眼前的区间第 k k k 小问题,我们应该怎么办?
我们的一个思路就是使用主席树。
主席树,全称为可持久化权值线段树,通过保存每次插入操作时的历史版本,以便查询区间第 k k k 小。
怎么保存呢?
一种暴力的方法就是可以每一次插入的时候都开一棵线段树。
但是这样对空间的使用不是很友好。
那么我们分析一下,每一次修改的时候涉及的都只有一条从根节点到我们想要修改的节点的路径。
那么我们每一次只需要新建从根节点一直到我们想要修改的节点的一整条路径就可以了。
我们每一次就只新建了 O ( log n ) O(\log n) O(logn) 个节点,大大节省了我们需要用的空间。
这种情况下,我们就不能采用堆式存储法(就是那个把左儿子和右儿子分别设为 p < < 1 p << 1 p<<1 和 p < < 1 ∣ 1 p << 1 | 1 p<<1∣1 的)了,必须使用动态开点。
对于预计需要的节点个数,我们可以采用这种办法:
我们每一次操作最多需要新开 2 ⌈ log n ⌉ + 1 2^{\lceil \log n \rceil +1} 2⌈logn⌉+1 个节点,我们只需要根据 n n n 和 m m m 的大小来声明就可以了。
怎么维护区间第 k k k 小?
一般情况下,我们可以用对区间建立一棵权值线段树的方法来维护当前区间的区间第
k
k
k 小。
原理就是,权值线段树维护的是区间内每一个数有多少个。
我们只需要用类似线段树上二分的思想在线段树上搜索即可。
但是我们不能每一次都根据当前区间建立一棵线段树,这样又麻烦又浪费空间。
所以我们就像利用前缀和减去另一个前缀和得到区间和一样,将给当前区间建立一棵线段树的问题转化成了给分别建立在 KaTeX parse error: Undefined control sequence: \[ at position 1: \̲[̲ 1,l-1 \] 和 KaTeX parse error: Undefined control sequence: \[ at position 1: \̲[̲ 1,r \] 上的两棵线段树做差。
而这样的两棵线段树可以通过每一次向线段树维护的区间最后追加一个数字来得到。
而这也就是主席树维护区间第 k k k 小的主要思想。
例题:
Luogu P3834 【模板】可持久化线段树 2:https://www.luogu.com.cn/problem/P3834
#include<bits/stdc++.h>
using namespace std;
const int N = 200010, M = 200010;
int n, m; int a[N];
vector<int>nums;
struct Node
{
int l, r;
int cnt;
}tr[N * 4 + N * 17]; int root[N], idx;
int find(int x) { return lower_bound(nums.begin(), nums.end(), x) - nums.begin(); }
int build(int l, int r)
{
int p = ++idx;
if(l == r)return p;
int mid = l + r >> 1;
tr[p].l = build(l, mid);
tr[p].r = build(mid + 1, r);
return p;
}
int insert(int p, int l, int r, int x)
{
int q = ++idx;
tr[q] = tr[p];
if(l == r)
{
tr[q].cnt++;
return q;
}
int mid = l + r >> 1;
if(x <= mid)tr[q].l = insert(tr[p].l, l, mid, x);
else tr[q].r = insert(tr[p].r, mid + 1, r, x);
tr[q].cnt = tr[tr[q].l].cnt + tr[tr[q].r].cnt;
return q;
}
int query(int q, int p, int l, int r, int k)
{
if(l == r)return r;
int cnt = tr[tr[q].l].cnt - tr[tr[p].l].cnt;
int mid = l + r >> 1;
if(k <= cnt)return query(tr[q].l, tr[p].l, l, mid, k);
else return query(tr[q].r, tr[p].r, mid + 1, r, k - cnt);
}
int main()
{
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i++)
{
scanf("%d", &a[i]);
nums.push_back(a[i]);
}
sort(nums.begin(), nums.end());
nums.erase(unique(nums.begin(), nums.end()), nums.end());
root[0] = build(0, nums.size() - 1);
for(int i = 1; i <= n; i++)
root[i] = insert(root[i - 1], 0, nums.size() - 1, find(a[i]));
while(m--)
{
int l, r, k;
scanf("%d%d%d", &l, &r, &k);
printf("%d\n", nums[query(root[r], root[l - 1], 0, nums.size() - 1, k)]);
}
return 0;
}
可持久化的普通线段树例题:
Luogu P3919 【模板】可持久化线段树 1(可持久化数组):https://www.luogu.com.cn/problem/P3919
#include<bits/stdc++.h>
using namespace std;
const int N = 1000010, M = N;
int n, m;
struct node
{
int l, r;
int v;
}tr[(N << 2) + (M * 20)];
int root[M], idx = 0;
int a[N];
int build(int l, int r)
{
int p = ++idx;
if(l == r)
{
tr[p].v = a[r];
return p;
}
int mid = (l + r) >> 1;
tr[p].l = build(l, mid);
tr[p].r = build(mid + 1, r);
//tr[p].v = tr[tr[p].l].v + tr[tr[p].r].v;
return p;
}
int change(int p, int l, int r, int loc, int val)
{
int q = ++idx;
tr[q] = tr[p];
if(l == r)
{
tr[q].v = val;
return q;
}
int mid = (l + r) >> 1;
if(loc <= mid)tr[q].l = change(tr[p].l, l, mid, loc, val);
else tr[q].r = change(tr[p].r, mid + 1, r, loc, val);
//tr[q].v = tr[tr[q].l].v + tr[tr[q].r].v;
return q;
}
int search(int p, int l, int r, int loc)
{
if(l ==r)return tr[p].v;
int mid = (l + r) >> 1;
if(loc <= mid)return search(tr[p].l, l, mid, loc);
else return search(tr[p].r, mid + 1, r, loc);
}
int main()
{
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i++)scanf("%d", &a[i]);
root[1] = build(1, n);
root[0]++;
int op, ver, loc, val;
while(m--)
{
scanf("%d%d%d", &ver, &op, &loc);
ver++;
if(op == 1)
{
scanf("%d", &val);
root[++root[0]] = change(root[ver], 1, n, loc, val);
}
else if(op == 2)
{
printf("%d\n", search(root[ver], 1, n, loc));
root[++root[0]] = root[ver];
}
else
{
puts("Youwike AK IOI!");
}
}
return 0;
}
区间最大子段和
区间最大子段和也可以使用线段树来维护。
我们考虑记录一下每一个区间内的最大前缀和、最大后缀和、最大子段和与区间和。
区间合并的时候这样合并:
- 最大前缀和:左儿子的最大前缀和与(右儿子的最大前缀和与左儿子的区间和之和)取最大值;
- 最大后缀和:右儿子的最大后缀和与(左儿子的最大前缀和与右儿子的区间和之和)取最大值;
- 最大子段和:左儿子的最大子段和、右儿子的最大子段和与(左儿子的最大后缀和与右儿子的最大前缀和之和)取最大值;
- 区间和:直接加。
写了个重载加号:
struct SegTree
{
int l, r;
int lms, rms;
int ms, sum;
friend SegTree operator + (const SegTree &lhs, const SegTree &rhs)
{
return
{
lhs.l,
rhs.r,
max(lhs.lms, rhs.lms + lhs.sum),
max(rhs.rms, lhs.rms + rhs.sum),
max(lhs.rms + rhs.lms, max(lhs.ms, rhs.ms)),
lhs.sum + rhs.sum
};
}
}tr[N << 3];
然后查询函数就是这个样子的:
SegTree query(int p, int l, int r)
{
if(tr[p].l >= l && tr[p].r <= r)return tr[p].ms;
int mid = (tr[p].l + tr[p].r) >> 1;
if(r <= mid)return query(p << 1, l, r);
else if(l > mid)return query(p << 1 | 1, l, r);
else return query(p << 1, l, r) + query(p << 1 | 1, l, r);
}
如果用线段树维护区间最大子段和的话还可以支持单点修改。
思路跟正常的单点修改一样。
void modify(int p, int pos, int k)
{
if(tr[p].l == tr[p].r)
{
tr[p] = { tr[p].l,tr[p].r,max(k,0),max(k,0),max(k,0),k };
return;
}
int mid = (tr[p].l + tr[p].r) >> 1;
if(pos <= mid)modify(p << 1, pos, k);
if(pos > mid)modify(p << 1 | 1, pos, k);
tr[p] = tr[p << 1] + tr[p << 1 | 1];
return;
}
暴力修改,不是正解。
能过官方数据,但是不能过Hack数据。
#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N = 100010;
struct SegTree
{
int l, r;
ll lms, rms;
ll ms, sum;
friend SegTree operator + (const SegTree &lhs, const SegTree &rhs)
{
return
{
lhs.l,
rhs.r,
max(lhs.lms, rhs.lms + lhs.sum),
max(rhs.rms, lhs.rms + rhs.sum),
max(lhs.rms + rhs.lms, max(lhs.ms, rhs.ms)),
lhs.sum + rhs.sum
};
}
}tr[N << 3];
ll a[N];
void build(int p, int l, int r)
{
tr[p].l = l, tr[p].r = r;
if(l == r)
{
tr[p] = { l,r,max(0ll,a[l]),max(0ll,a[l]),max(0ll,a[l]),a[l] };
return;
}
int mid = (l + r) >> 1;
build(p << 1, l, mid);
build(p << 1 | 1, mid + 1, r);
tr[p] = tr[p << 1] + tr[p << 1 | 1];
return;
}
void modify(int p, int pos, ll k)
{
if(tr[p].l == tr[p].r)
{
tr[p] = { tr[p].l,tr[p].r,max(k,0ll),max(k,0ll),max(k,0ll),k };
return;
}
int mid = (tr[p].l + tr[p].r) >> 1;
if(pos <= mid)modify(p << 1, pos, k);
if(pos > mid)modify(p << 1 | 1, pos, k);
tr[p] = tr[p << 1] + tr[p << 1 | 1];
return;
}
SegTree query(int p, int l, int r)
{
if(tr[p].l >= l && tr[p].r <= r)return tr[p];
int mid = (tr[p].l + tr[p].r) >> 1;
if(r <= mid)return query(p << 1, l, r);
else if(l > mid)return query(p << 1 | 1, l, r);
else return query(p << 1, l, r) + query(p << 1 | 1, l, r);
}
int main()
{
int m, n;
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i++)
scanf("%lld", &a[i]);
build(1, 1, n);
while(m--)
{
int op, l, r;
scanf("%d%d%d", &op, &l, &r);
if(op == 0)
{
ll k;
scanf("%lld", &k);
for(int i = l; i <= r; i++)
{
if(k > a[i])
{
a[i] = k;
modify(1, i, k);
}
}
}
else if(op == 1)
{
SegTree res = query(1, l, r);
printf("%lld\n", res.ms);
}
}
return 0;
}