首页 > > 程序设计 > C/C++ >

BZOJ4355: Play with sequence(吉司机线段树)

2018-09-19 02:44:35来源：博客园阅读 ()

题意

题目链接

Sol

传说中的吉司机线段树？？感觉和BZOJ冒险那题差不多，就是强行剪枝。。。

这题最坑的地方在于对于操作1，$C >= 0$, 操作2中需要对0取max，$a[i] >= 0$，这不就是统计最小值出现的次数么？？

按照套路

维护好区间赋值标记 / 区间加法标记 / 区间max标记 / 区间最小值 / 区间最小值出现的次数 / 区间次小值

对于第二个操作就拆成区间加和区间max

区间max是一个很神奇的操作

设当前加入的数为val

若val>=mn，那该操作对该区间无影响

若se < val < mn，该操作只会对次小值产生影响，因为对其他的标记均不会产生影响，因此打一个额外的标记即可

否则暴力递归

时间复杂度：$O(n log^2n)$??

另外这东西可以做区间加 / 查询历史版本，前者维护一下标记就行，，后者嘛，，，等长大了再研究吧qwq

#include<cstdio>
#include<algorithm>
#define LL long long 
//#define int long long 
using namespace std;
const int MAXN = 3 * 1e5 + 10;
const LL INF = 1e10 +10;
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;
}
#define ls k << 1
#define rs k << 1 | 1
int N, Q;
int a[MAXN];
struct Node {
    int l, r, siz, cnt;
    LL Mx, add, mn, se, cov;
}T[MAXN << 2];
void update(int k) {
    //T[k].mn = min(T[ls].mn, T[rs].mn);
    if(T[ls].mn < T[rs].mn)      T[k].cnt = T[ls].cnt, T[k].mn = T[ls].mn, T[k].se = min(T[rs].mn, T[ls].se);
    else if(T[ls].mn > T[rs].mn) T[k].cnt = T[rs].cnt, T[k].mn = T[rs].mn, T[k].se = min(T[ls].mn, T[rs].se);
    else              T[k].cnt = T[ls].cnt + T[rs].cnt, T[k].mn = T[ls].mn, T[k].se = min(T[ls].se, T[rs].se);
}
void MemP(int k, LL val) {
    T[k].cov = val;
    T[k].cnt = T[k].siz;
    T[k].se = INF;
    T[k].Mx = -INF;
    T[k].add = 0;
    T[k].mn = val;
}
void AddP(int k, LL val) {
    if(T[k].se != INF)  T[k].se += val;
    if(T[k].Mx != -INF) T[k].Mx += val;
    T[k].mn += val;
    T[k].add += val;
}
void MaxP(int k, LL val) {
    T[k].mn = max(T[k].mn, val);//
    T[k].Mx = max(T[k].Mx, val);
}
void pushdown(int k) {
    if(T[k].cov != INF) MemP(ls, T[k].cov), MemP(rs, T[k].cov), T[k].cov = INF;
    if(T[k].add) AddP(ls, T[k].add), AddP(rs, T[k].add), T[k].add = 0;
    if(T[k].Mx != -INF) MaxP(ls, T[k].Mx), MaxP(rs, T[k].Mx), T[k].Mx = -INF;
}
void Build(int k, int ll, int rr) {
    T[k] = (Node) {ll, rr, rr - ll + 1, 0, -INF, 0, 0, INF, INF};
    if(ll == rr) {
        T[k].mn = a[ll]; T[k].cnt = 1;
        return ;
    }
    int mid = T[k].l + T[k].r >> 1;
    Build(ls, ll, mid); Build(rs, mid + 1, rr);
    update(k);
}
void IntMem(int k, int ll, int rr, LL val) {
    if(ll <= T[k].l && T[k].r <= rr) {MemP(k, val); return ;}
    pushdown(k);
    int mid = T[k].l + T[k].r >> 1;
    if(ll <= mid) IntMem(ls, ll, rr, val); 
    if(rr >  mid) IntMem(rs, ll, rr, val);
    update(k);
}
void IntAdd(int k, int ll, int rr, LL val) {
    if(ll <= T[k].l && T[k].r <= rr) {AddP(k, val); return ;}
    pushdown(k);
    int mid = T[k].l + T[k].r >> 1;
    if(ll <= mid) IntAdd(ls, ll, rr, val);
    if(rr >  mid) IntAdd(rs, ll, rr, val);
    update(k);
}
void IntMax(int k, int ll, int rr, LL val) {
    if(val <= T[k].mn) return ;
    if(ll <= T[k].l && T[k].r <= rr && T[k].se > val) {//tag
        MaxP(k, val); 
        return ;
    }
    int mid = T[k].l + T[k].r >> 1;
    pushdown(k);
    if(ll <= mid) IntMax(ls, ll, rr, val); 
    if(rr >  mid) IntMax(rs, ll, rr, val);
    update(k);
}
LL Query(int k, int ll, int rr) {
    // int ans = 0;
    if(ll <= T[k].l && T[k].r <= rr) 
        return (T[k].mn == 0 ? T[k].cnt : 0);
    pushdown(k);
    int mid = T[k].l + T[k].r >> 1;
    if(ll > mid) return Query(rs, ll, rr);
    else if(rr <= mid) return Query(ls, ll, rr);
    else return Query(ls, ll, rr) + Query(rs, ll, rr);
}
main() {
    // freopen("4355.in", "r", stdin);
    // freopen("4355.out", "w", stdout);
    N = read(); Q = read();
    for(int i = 1; i <= N; i++) a[i] = read();
    Build(1, 1, N);
    while(Q--) {
        int opt = read(), l = read(), r = read(), val;
        if(opt == 3) printf("%d\n", Query(1, l, r));
        else {
            val = read();
            if(opt == 1) IntMem(1, l, r, val);
            else {
                IntAdd(1, l, r, val);
                IntMax(1, l, r, 0);
            }
        }
    }
    return 0;
}