CF935F Fafa and Array 题解
题目链接:Fafa and Array
题意:
给出一个序列 $a$,定义函数 $f(a)=\sum_{i=1}^{n-1}|a_i-a_{i+1}|$。
给出两种操作:
$\texttt{1 l r x}$,表示在区间 $[l,r]$ 内找一个位置,
使得把这个位置的值加上 $x$ 后 $f(a)$ 最大,求这个最大值(修改是临时的)。
$\texttt{2 l r x}$,表示把区间 $[l,r]$ 加上 $x$。
输入格式
第一行一个正整数 $n~(3 \le n \le 10^5)$ ,表示数组的长度。
第二行 $n$ 个正整数 $a_{1},a_{2},…,a_{n}~(0 < a_i \le 10^9)$ ,表示数组的元素。
第三行一个正整数 $q~(1 \le q \le 10^5)$ ,表示询问的个数。
接下来 $q$ 行,第 $i$ 行有四个整数 $t_i,l_i,r_i,x_i~(t_i \in\{1, 2\},\,1<l_i \leq r_i<n,\, 0<x_i \leq 10^9)$ 。
数据保证至少有一个询问。
输出格式
对于每个询问,输出一行答案。
注意到在某个 $a_i$ 上修改造成的影响与 $a_{i-1},a_{i+1}$ 有关,因此考察 $a_i + x$ 与 $a_{i-1},a_{i+1}$ 间的大小关系
若 $a_i+x<a_{i-1}\, \land \, a_i+x<a_{i+1}$ ,则贡献为 $\Delta = -2x$ 。
若 $a_{i-1}<a_i+x<a_{i+1}$ ,则贡献为
若 $a_{i-1} \leq a_i+x\,\land\, a_{i+1} \leq a_i+x$ ,
- $a_{i-1} \leq a_i,~a_{i+1} \leq a_i$ ,那么 $\Delta = 2x$ 。
- $a_{i-1} \leq a_i<a_{i+1}$ ,那么 $\Delta = 2x - 2(a_{i+1} - a_i)$ 。
- $a_i<a_{i-1}\, \land \, a_i<a_{i+1}$ ,那么 $\Delta = 2x - 2((a_{i-1} - a_i) + (a_{i+1} - a_i))$ 。
显然只有第三种情况贡献为正,也就是对可以让 $f(a)$ 增大。
考虑用线段树维护 $\max \left(0, a_{i-1}-a_i\right)+\max \left(0, a_{i+1}-a_i\right)$ 的最小值。
这个差分一下就是 $\max \left(0,-a_{i-1}\right)+\max \left(0, a_i\right)$ 的最小值,需要特判 $l=1$ 和 $r=n$ 的情况。
那么贡献可能为负吗?显然当 $l=r$ 时我们只能强行修改,特判即可。
注意线段树单点修改了 $x$ 后 $x+1$ 的那个最小值会变更,需要额外跑一遍修改。
时间复杂度 $\mathcal{O}(n \log n)$
代码:
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define INF 0x3f3f3f3f3f3f3f3f
void up(int &x, int y) { x < y ? x = y : 0; }
void down(int &x, int y) { x > y ? x = y : 0; }
#define rep(i, a, b) for(int i = (a), i##END = (b); i <= i##END; i++)
#define Rep(i, a, b) for(int i = (a), i##END = (b); i >= i##END; i--)
#define N ((int)(1e5 + 15))
#define ls(at) ((at) << 1)
#define rs(at) ((at) << 1 | 1)
int n, ans, a[N], tr[N * 4];
void push_up(int at) { down(tr[at] = tr[ls(at)], tr[rs(at)]); }
void build(int at = 1, int l = 1, int r = n - 1)
{
if(l == r) {
if(l > 1) tr[at] = max(0ll, a[l]) + max(0ll, -a[l - 1]);
ans += abs(a[l]); return;
}
int mid = (l + r) >> 1;
build(ls(at), l, mid); build(rs(at), mid + 1, r); push_up(at);
}
void modify(int x, int v, int at = 1, int l = 1, int r = n - 1)
{
if(l == r)
{
ans += abs(a[l] + v) - abs(a[l]), a[l] += v;
if(x > 1) tr[at] = max(0ll, -a[l - 1]) + max(0ll, a[l]);
if(v && x < n - 1) modify(x + 1, 0); return;
}
int mid = (l + r) >> 1;
x <= mid ? modify(x, v, ls(at), l, mid) : modify(x, v, rs(at), mid + 1, r); push_up(at);
}
int qry(int nl, int nr, int at = 1, int l = 1, int r = n - 1)
{
if(nl <= l && r <= nr) return tr[at];
int mid = (l + r) >> 1;
if(nl > mid) return qry(nl, nr, rs(at), mid + 1, r);
if(nr <= mid) return qry(nl, nr, ls(at), l, mid);
return min(qry(nl, nr, ls(at), l, mid), qry(nl, nr, rs(at), mid + 1, r));
}
signed main()
{
ios::sync_with_stdio(0);
cin.tie(0); cout.tie(0);
// freopen("check.in","r",stdin);
// freopen("check.out","w",stdout);
cin >> n; rep(i, 1, n) cin >> a[i];
rep(i, 1, n - 1) a[i] = a[i + 1] - a[i];
build(); int q; cin >> q;
rep(i, 1, q)
{
static int op, l, r, x; cin >> op >> l >> r >> x;
if(op == 1)
{
if(l == r)
{
if(l == 1) cout << (ans - abs(a[l]) + abs(a[l] - x)) << '\n';
else if(l == n) cout << (ans - abs(a[l - 1]) + abs(a[l - 1] + x)) << '\n';
else cout << (ans - abs(a[l - 1]) - abs(a[l]) + abs(a[l - 1] + x) + abs(a[l] - x)) << '\n';
}else
{
int tmp = 2 * (x - qry(l > 1 ? l : 2, r));
if(l == 1) up(tmp, -abs(a[l]) + abs(a[l] - x));
if(r == n) up(tmp, -abs(a[r - 1]) + abs(a[r - 1] + x));
cout << (ans + max(0ll, tmp)) << '\n';
}
}else { if(l > 1) modify(l - 1, x); if(r < n) modify(r, -x); }
}
return 0;
}
参考文献:
[1] https://www.luogu.com.cn/article/xzmyel52
题外话:
Typora 怎么偷偷给我更新了?怪不得某个 bug 没了。