NC200195 区区区间

题目链接

题目

题目描述

Keven 特别喜欢线段树,他给你一个长度为 n 的序列,对序列进行 m 次操作。

操作有两种:

1 1 l r k :表示将下标在 [l,r] 区间内的数字替换成 [k,k+1,,k+rl]

2 l r :表示查询区间 [l,r] 的区间和

输入描述

第一行两个整数 nm ,表示序列的长度和操作次数 (1n,m2×105) 第二行 n 个整数,表示序列的初始值 a1,a2,an(1ai2×105) 接下来 m 行,每行三或四个数字,若第一个数字是 1 ,则表示操作 1 ,反之则表示操作 2(1lrn1k2e5)

输出描述

对于每个操作 2 ,输出一行一个整数表示区间和。

示例1

输入

5 5
1 1 1 1 1
2 1 5
1 1 5 1
2 1 5
1 1 3 3
2 1 3

输出

5
15
12

说明

第一次1操作后,序列是1 2 3 4 5
第二次1操作后,序列是3 4 5 4 5

题解

知识点:线段树,数学。

本题的区间修改不具有区间无关性,因为在一次修改的不同的子区间,修改维护的值是不同的,因此需要在 update 函数内做特殊处理,是非标准的线段树。

区间信息需要维护区间长度 len 、区间和 sum ,区间合并直接加即可。

区间修改需要维护数值起始点 k ,新区间和可以直接公式求出 (2k+len1)len2

区间修改可以设置懒标记,只需要记录最新一次的修改即可。

注意,要考虑标记的单位元值,即表示未修改过的标记值,这里取 1 。在 push_down 函数中需先判断标记是否为 1 ,再进行下传。

区间更新函数有两种写法:

  1. 不更改修改值,在更新时根据子区间位置直接求出实际修改值。
void update(int rt, int l, int r, int x, int y, F f) {
if (r < x || y < l) return;
if (x <= l && r <= y) return node[rt] = F{ f.k + l - x }(node[rt]), lazy[rt] = F{ f.k + l - x }(lazy[rt]), void();
push_down(rt);
int mid = l + r >> 1;
update(rt << 1, l, mid, x, y, f);
update(rt << 1 | 1, mid + 1, r, x, y, f);
node[rt] = node[rt << 1] + node[rt << 1 | 1];
}
  1. 每次更新子区间时,同时更改相应的修改值。注意,这种方法需要控制更新的子区间,否则无法计算相应的修改值。
void update(int rt, int l, int r, int x, int y, F f) {
if (x == l && r == y) return node[rt] = f(node[rt]), lazy[rt] = f(lazy[rt]), void();
push_down(rt);
int mid = l + r >> 1;
if (y <= mid) update(rt << 1, l, mid, x, y, f);
else if (x >= mid + 1) update(rt << 1 | 1, mid + 1, r, x, y, f);
else update(rt << 1, l, mid, x, mid, f), update(rt << 1 | 1, mid + 1, r, mid + 1, y, { f.k + mid + 1 - x });
node[rt] = node[rt << 1] + node[rt << 1 | 1];
}

实际上,后者精确控制更新区间的方法,在遇到其他因为修改与区间位置有关等问题,需要魔改 update 函数时,更好用。

时间复杂度 O((n+m)logn)

空间复杂度 O(n)

代码

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
struct T {
int len;
ll sum;
static T e() { return { 0,0 }; }
friend T operator+(const T &a, const T &b) { return { a.len + b.len,a.sum + b.sum }; }
};
struct F {
int k;
static F e() { return{ -1 }; }
T operator()(const T &x) { return { x.len, (2LL * k + x.len - 1) * x.len / 2 }; }
F operator()(const F &g) { return { k }; }
};
class SegmentTreeLazy {
int n;
vector<T> node;
vector<F> lazy;
void push_down(int rt) {
if (!~lazy[rt].k) return;
node[rt << 1] = lazy[rt](node[rt << 1]);
lazy[rt << 1] = lazy[rt](lazy[rt << 1]);
lazy[rt].k += node[rt << 1].len;
node[rt << 1 | 1] = lazy[rt](node[rt << 1 | 1]);
lazy[rt << 1 | 1] = lazy[rt](lazy[rt << 1 | 1]);
lazy[rt] = F::e();
}
/* void update(int rt, int l, int r, int x, int y, F f) {
if (r < x || y < l) return;
if (x <= l && r <= y) return node[rt] = F{ f.k + l - x }(node[rt]), lazy[rt] = F{ f.k + l - x }(lazy[rt]), void();
push_down(rt);
int mid = l + r >> 1;
update(rt << 1, l, mid, x, y, f);
update(rt << 1 | 1, mid + 1, r, x, y, f);
node[rt] = node[rt << 1] + node[rt << 1 | 1];
} */
void update(int rt, int l, int r, int x, int y, F f) {
if (x == l && r == y) return node[rt] = f(node[rt]), lazy[rt] = f(lazy[rt]), void();
push_down(rt);
int mid = l + r >> 1;
if (y <= mid) update(rt << 1, l, mid, x, y, f);
else if (x >= mid + 1) update(rt << 1 | 1, mid + 1, r, x, y, f);
else update(rt << 1, l, mid, x, mid, f), update(rt << 1 | 1, mid + 1, r, mid + 1, y, { f.k + mid + 1 - x });
node[rt] = node[rt << 1] + node[rt << 1 | 1];
}
T query(int rt, int l, int r, int x, int y) {
if (r < x || y < l) return T::e();
if (x <= l && r <= y) return node[rt];
push_down(rt);
int mid = l + r >> 1;
return query(rt << 1, l, mid, x, y) + query(rt << 1 | 1, mid + 1, r, x, y);
}
public:
SegmentTreeLazy(int _n = 0) { init(_n); }
SegmentTreeLazy(const vector<T> &src) { init(src); }
void init(int _n) {
n = _n;
node.assign(n << 2, T::e());
lazy.assign(n << 2, F::e());
}
void init(const vector<T> &src) {
assert(src.size());
init(src.size() - 1);
function<void(int, int, int)> build = [&](int rt, int l, int r) {
if (l == r) return node[rt] = src[l], void();
int mid = l + r >> 1;
build(rt << 1, l, mid);
build(rt << 1 | 1, mid + 1, r);
node[rt] = node[rt << 1] + node[rt << 1 | 1];
};
build(1, 1, n);
}
void update(int x, int y, F f) { update(1, 1, n, x, y, f); }
T query(int x, int y) { return query(1, 1, n, x, y); }
};
int main() {
std::ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
int n, m;
cin >> n >> m;
vector<T> a(n + 1);
for (int i = 1;i <= n;i++) {
int x;
cin >> x;
a[i] = { 1,x };
}
SegmentTreeLazy sgt(a);
while (m--) {
int op, l, r;
cin >> op >> l >> r;
if (op == 1) {
int k;
cin >> k;
sgt.update(l, r, { k });
}
else cout << sgt.query(l, r).sum << '\n';
}
return 0;
}
posted @   空白菌  阅读(21)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 开源Multi-agent AI智能体框架aevatar.ai,欢迎大家贡献代码
· Manus重磅发布:全球首款通用AI代理技术深度解析与实战指南
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· 没有Manus邀请码?试试免邀请码的MGX或者开源的OpenManus吧
点击右上角即可分享
微信分享提示