cf1136E. Nastya Hasn't Written a Legend(二分 线段树)

题意

题目链接

Sol

yy出了一个暴躁线段树的做法。

因为题目保证了

\(a_i + k_i <= a_{i+1}\)

那么我们每次修改时只需要考虑取max就行了。

显然从一个位置开始能影响到的位置是单调的,而且这些位置的每个改变量都是\((a_i + x) + \sum_{t=i}^{j-1} k_t\)

那么可以建两棵线段树分别维护这两部分的值

每次修改的时候二分出要修改的位置。

打cf一定要记得开数据结构题啊qwq

#include<bits/stdc++.h> 
#define Pair pair<int, int>
#define MP(x, y) make_pair(x, y)
#define fi first
#define se second
#define int long long 
#define LL long long 
#define ull unsigned long long 
#define Fin(x) {freopen(#x".in","r",stdin);}
#define Fout(x) {freopen(#x".out","w",stdout);}
using namespace std;
const int MAXN = 1e6 + 10, mod = 1e9 + 7, INF = 1e9 + 10;
const double eps = 1e-9;
template <typename A, typename B> inline bool chmin(A &a, B b){if(a > b) {a = b; return 1;} return 0;}
template <typename A, typename B> inline bool chmax(A &a, B b){if(a < b) {a = b; return 1;} return 0;}
template <typename A, typename B> inline LL add(A x, B y) {if(x + y < 0) return x + y + mod; return x + y >= mod ? x + y - mod : x + y;}
template <typename A, typename B> inline void add2(A &x, B y) {if(x + y < 0) x = x + y + mod; else x = (x + y >= mod ? x + y - mod : x + y);}
template <typename A, typename B> inline LL mul(A x, B y) {return 1ll * x * y % mod;}
template <typename A, typename B> inline void mul2(A &x, B y) {x = (1ll * x * y % mod + mod) % mod;}
template <typename A> inline void debug(A a){cout << a << '\n';}
template <typename A> inline LL sqr(A x){return 1ll * x * x;}
template <typename A, typename B> inline LL fp(A a, B p, int md = mod) {int b = 1;while(p) {if(p & 1) b = mul(b, a);a = mul(a, a); p >>= 1;}return b;}
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int N, a[MAXN], K[MAXN], S[MAXN], val[MAXN];
//第一棵线段树支持区间赋值 / 区间求和
//第二棵线段树支持区间赋值 / 区间求和 .....
struct Seg {
#define ls k << 1
#define rs k << 1 | 1
    int sum[MAXN], tag[MAXN], siz[MAXN];
    void update(int k) {
        sum[k] = sum[ls] + sum[rs];
    }
    void ps(int k, int f) {
        sum[k] = siz[k] * f;
        tag[k] = f;
    }
    void pushdown(int k) {
        if(tag[k] == -INF) return ;
        ps(ls, tag[k]); ps(rs, tag[k]);
        tag[k] = -INF;
    }
    void Build(int k, int l, int r) {
        siz[k] = r - l + 1; tag[k] = -INF;
        if(l == r) {sum[k] = val[l]; return ;}
        int mid = l + r >> 1;
        Build(ls, l, mid); Build(rs, mid + 1, r);
        update(k);
    }
    void IntMem(int k, int l, int r, int ql, int qr, int v) {
        if(ql <= l && r <= qr) {ps(k, v); return ;}
        int mid = l + r >> 1;
        pushdown(k);
        if(ql <= mid) IntMem(ls, l, mid, ql, qr, v);
        if(qr  > mid) IntMem(rs, mid + 1, r, ql, qr, v);
        update(k);
    }
    int IntQuery(int k, int l, int r, int ql, int qr) {
        if(ql <= l && r <= qr) return sum[k];
        int mid = l + r >> 1;
        pushdown(k);
        if(ql > mid) return IntQuery(rs, mid + 1, r, ql, qr);
        else if(qr <= mid) return IntQuery(ls, l, mid, ql, qr);
        else return IntQuery(ls, l, mid, ql, qr) + IntQuery(rs, mid + 1, r, ql, qr);
    }
}T[2];
int GetA(int x) {
    return T[0].IntQuery(1, 1, N, x, x) + K[x] - T[1].IntQuery(1, 1, N, x, x);
}
void Modify(int x, int v) {
    int l = x, r = N, ans = l, tmp = GetA(x);
    while(l <= r) {
        int mid = l + r >> 1;
        if(tmp + v + K[mid] - K[x] >= GetA(mid)) l = mid + 1, ans = mid;
        else r = mid - 1;   
    }
    T[0].IntMem(1, 1, N, x, ans, tmp + v);
    T[1].IntMem(1, 1, N, x, ans, K[x]);
}
signed main() {
    N = read();
    for(int i = 1; i <= N; i++) a[i] = read();
    for(int i = 2; i <= N; i++) 
        K[i] = read() + K[i - 1];
    memcpy(val, a, sizeof(a));
    T[0].Build(1, 1, N); 
    memcpy(val, K, sizeof(K));
    T[1].Build(1, 1, N);
    for(int i = 1; i <= N; i++) S[i] = S[i - 1] + K[i];
    int Q = read();
    while(Q--) {
        char c = 'g';
        while(c != 's' && c != '+') c = getchar();
        if(c == '+') {
            int x = read(), v = read(); 
            Modify(x, v);
        }
        else {
            int l = read(), r = read();
            int pre = T[0].IntQuery(1, 1, N, l, r);
            int nxt = S[r] - S[l - 1] - T[1].IntQuery(1, 1, N, l, r);
            cout << pre + nxt << '\n'; 
        }
        c = 'g';
    }
    return 0;
}

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏MasiMaro 的技术博文

C语言中处理结构体的原理

汇编中有几种寻址方式,分别是直接寻址:(ds:[idata])、寄存器间接寻址(ds:[bx])、寄存器相对寻址(ds:[bx + idata]、ds:[bx ...

15620
来自专栏一名合格java开发的自我修养

Linux之shell编程条件判断-if,while,for,case

sh if.sh intsmaze nihao 通过执行脚本时,在脚本命令的后面传入参数。

16220
来自专栏MasiMaro 的技术博文

向上取整算法

在进行内存分配的时候一般都需要在实际使用内存大小的基础上进行内存对齐,比如一般32位平台进行4字节对齐,而64位平台使用8字节对齐等等。 一般采用的算法是先利...

16320
来自专栏JAVA同学会

JAVA9模块化详解(一)——模块化的定义

java9已经出来有一段时间了,今天向大家介绍一下java9的一个重要特性——模块化。模块化系统的主要目的如下:

10310
来自专栏木溪知识加油站

笔记——安卓消息机制Handler(十六)

1、定义:Android的消息机制主要是指Handler的运行机制,Handler并不是专门用于更新UI的,它只是常被开发者用来更新UI,是同一个进程中线程间的...

12440
来自专栏游戏杂谈

从10W个数中随机抽走2个数,求出那两个数是多少

从1到10w(共10w个数)中随机抽走2个数,然后打乱剩下的数的顺序,问如果从这剩下的数中快速的找出抽走的是哪2个数?

11630
来自专栏一名合格java开发的自我修养

flink1.7自定义source实现

数据的来源是flink程序从中读取输入的地方。我们可以使用StreamExecutionEnvironment.addSource(sourceFunction...

15330
来自专栏游戏杂谈

JavaScript的类型错误:Illegal invocation

昨一看好像也没什么问题,在IE6~8下运行也没有任何的问题,可是在其它浏览器下就报错了,报了这样一个错:

1.5K10
来自专栏游戏杂谈

JavaScript求最大公共子串

对两个字符串的每一项都进行比较,若匹配则该项为1,不匹配则为0。然后求出对角线最长为1的那一段序列,即为最大公共子串。

11820
来自专栏木溪知识加油站

笔记——四大组件(十五)

1、Activity是一种展示型组件。Activity的启动由Intent触发,其中Intent可以分为显式Intent和隐式Intent,显式Intent可以...

10720

扫码关注云+社区

领取腾讯云代金券

年度创作总结 领取年终奖励