嘘~ 正在从服务器偷取页面 . . .

AT_dp_w Intervals 题解


AT_dp_w Intervals 题解

题目链接:AT_dp_w Intervals

题意

给定 $m$ 条规则形如 $(l_i,r_i,a_i)$,对于一个 01 串,其分数的定义是:

对于第 $i$ 条规则,若该串在 $[l_i,r_i]$ 中至少有一个 $1$,则该串的分数增加 $a_i$。

你需要求出长度为 $n$ 的 01 串中的最大分数。

输入格式

第一行 $n,m$ ,接下来 $m$ 行 $l_i,r_i,a_i$ 。

输出格式

输出字符串的最大得分。

数据范围

$1 \le n,m \le 2\times 10^5$ 且 $1\le l_i \le r_i\le n,~|a_i|\le 10^9$

这道题和 P9871 [NOIP2023] 天天爱打卡 思路挺像的。

设 $f(i,j)$ 为前 $i$ 个位置中,上一个 $1$ 放在 $j$ 的最优方案,则

后面这个东西可以把线段(即规则)按右端点排序,这样每次枚举 $r_k=i$ 的线段即可。

注意到对于线段 $(l_k,r_k)$ ,任何 $j \in [l_k, r_k]$ 的 $f(i-1,j)$ 都可以吃到 $a_k$ ,这就是一个区间修改

考虑用线段树维护每一轮的 $f(i,j)$ ,而不是用传统的滚动数组维护。

这样我们枚举 $i$ 的时候,就可以一轮区间加了,注意 $f(i,i)$ 要单独修改。

时间复杂度 $\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 N ((int)(2e5 + 15))
#define ls(at) ((at) << 1)
#define rs(at) ((at) << 1 | 1)

struct line { int l, r, v; } a[N];
struct node { int mx, tag; } tr[N * 4];
int n, m;
void push_up(int at) { up(tr[at].mx = tr[ls(at)].mx, tr[rs(at)].mx); }
void proc(int at, int v) { tr[at].mx += v; tr[at].tag += v; }
void push_down(int at)
{ 
    proc(ls(at), tr[at].tag);
    proc(rs(at), tr[at].tag); tr[at].tag = 0;
}
void build(int at, int l, int r)
{
    if(l == r) return ;
    push_down(at); 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)
{
    if(l == r) return tr[at].mx += v, void(0);
    push_down(at); int mid = (l + r) >> 1;
    if(x <= mid) Modify(x, v, ls(at), l, mid);
    else Modify(x, v, rs(at), mid + 1, r);
    push_up(at);
}
void update(int nl, int nr, int v, int at = 1, int l = 1, int r = n)
{
    if(nl <= l && r <= nr) return proc(at, v);
    push_down(at); int mid = (l + r) >> 1;
    if(nl <= mid) update(nl, nr, v, ls(at), l, mid);
    if(nr > mid) update(nl, nr, v, rs(at), mid + 1, r);
    push_up(at);
}
int query(int nl, int nr, int at = 1, int l = 1, int r = n)
{
    if(nl <= l && r <= nr) return tr[at].mx;
    push_down(at); int mid = (l + r) >> 1, mx = -INF;
    if(nl <= mid) up(mx, query(nl, nr, ls(at), l, mid));
    if(nr > mid) up(mx, query(nl, nr, rs(at), mid + 1, r));
    return mx;
}
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 >> m;
    for(int i = 1; i <= m; i++) cin >> a[i].l >> a[i].r >> a[i].v;
    sort(a + 1, a + 1 + m, [](line a, line b) { return a.r < b.r; });
    build(1, 1, n);
    for(int i = 1, p = 1; i <= n; i++)
    {
        Modify(i, max(0ll, tr[1].mx));
        while(p <= m && a[p].r == i) { update(a[p].l, i, a[p].v), ++p; }
    }
    cout << max(tr[1].mx, 0ll) << '\n';
    return 0;
}

参考文献

[1] https://www.luogu.com.cn/article/9gp9m3nn


文章作者: q779
版权声明: 本博客所有文章除特别声明外,均采用 CC BY-NC-ND 4.0 许可协议。转载请注明来源 q779 !
评论
  目录