……一道丧病线段树膜板题……

被常数卡的死去活来……QAQ

学到了些奇技淫巧:把取min标记 和 区间最小值 合并

可以快很多……

#include <bits/stdc++.h>
#define lc(t) ((t) << 1)
#define rc(t) (((t) << 1) | 1)
#define N 2000010
#define INF 1000000000
#define LL long long template <class T> inline T &read(T &x)
{
static int f;
static char c;
for (f = ; !isdigit(c = getchar()); ) {
if (c == '-')
f = -;
}
for (x = ; isdigit(c); c = getchar()) {
x = x * + c - ;
}
return x *= f;
} template <class T> inline void write(T x, const char p = '\n')
{
static int top;
static int s[];
if (x < ) {
x = -x;
putchar('-');
}
do s[++ top] = x % + ;
while (x /= );
while (top)
putchar(s[top --]);
putchar(p);
} using namespace std;
int n, m;
int mn[N], mx[N], cn[N], cx[N], sn[N], sx[N]; LL sum[N];
int t_a[N], t_n[N], t_x[N], i_n[N], i_x[N];
int lb[N], rb[N];
inline void upd(register int t)
{
sum[t] = sum[lc(t)] + sum[rc(t)];
if (mn[lc(t)] == mn[rc(t)])
mn[t] = mn[lc(t)], cn[t] = cn[lc(t)] + cn[rc(t)], sn[t] = min(sn[lc(t)], sn[rc(t)]);
else if (mn[lc(t)] < mn[rc(t)])
mn[t] = mn[lc(t)], cn[t] = cn[lc(t)], sn[t] = min(sn[lc(t)], mn[rc(t)]);
else
mn[t] = mn[rc(t)], cn[t] = cn[rc(t)], sn[t] = min(mn[lc(t)], sn[rc(t)]); if (mx[lc(t)] == mx[rc(t)])
mx[t] = mx[lc(t)], cx[t] = cx[lc(t)] + cx[rc(t)], sx[t] = max(sx[lc(t)], sx[rc(t)]);
else if (mx[lc(t)] > mx[rc(t)])
mx[t] = mx[lc(t)], cx[t] = cx[lc(t)], sx[t] = max(sx[lc(t)], mx[rc(t)]);
else
mx[t] = mx[rc(t)], cx[t] = cx[rc(t)], sx[t] = max(mx[lc(t)], sx[rc(t)]); }
void build(int t, int l, int r)
{
lb[t] = l; rb[t] = r;
if (l == r)
{
int x;
read(x);
mn[t] = mx[t] = sum[t] = x; cn[t] = cx[t] = ; sn[t] = INF; sx[t] = -INF;
}
else
{
build(lc(t), l, (l + r) / );
build(rc(t), (l + r) / + , r);
upd(t);
}
}
inline void add(register int t, int d)
{
sum[t] += (LL)d * (rb[t] - lb[t] + );
mn[t] += d; mx[t] += d; sn[t] += d; sx[t] += d;
t_a[t] += d;
} template <class T> inline void chkmin(T &x, T y) { x > y && (x = y); }
template <class T> inline void chkmax(T &x, T y) { x < y && (x = y); }
inline void m_x(register int t, int d)
{
sum[t] += (LL)cn[t] * (d - mn[t]);
mn[t] = d; chkmax(mx[t], d);
if (mx[t] == mn[t])
{
cx[t] = cn[t] = rb[t] - lb[t] + ;
sum[t] = (LL)mx[t] * (rb[t] - lb[t] + );
sn[t] = INF; sx[t] = -INF;
}
else chkmax(sx[t], d);
}
inline void m_n(register int t, int d)
{
sum[t] += (LL)cx[t] * (d - mx[t]);
mx[t] = d; chkmin(mn[t], d);
if (mx[t] == mn[t])
{
cx[t] = cn[t] = rb[t] - lb[t] + ;
sum[t] = (LL)mx[t] * (rb[t] - lb[t] + );
sn[t] = INF; sx[t] = -INF;
}
else chkmin(sn[t], d);
}
inline void psdw(register int t)
{
if (t_a[t])
{
add(lc(t), t_a[t]);
add(rc(t), t_a[t]);
t_a[t] = ;
}
if (mx[lc(t)] > mx[t] && sx[lc(t)] < mx[t]) m_n(lc(t), mx[t]);
if (mx[rc(t)] > mx[t] && sx[rc(t)] < mx[t]) m_n(rc(t), mx[t]);
if (mn[lc(t)] < mn[t] && sn[lc(t)] > mn[t]) m_x(lc(t), mn[t]);
if (mn[rc(t)] < mn[t] && sn[rc(t)] > mn[t]) m_x(rc(t), mn[t]);
}
namespace segment
{
int l, r, d;
void seg_add(int t)
{
if (l <= lb[t] && rb[t] <= r)
{
add(t, d);
return;
}
psdw(t);
if (l <= rb[lc(t)]) seg_add(lc(t));
if (r >= lb[rc(t)]) seg_add(rc(t));
upd(t);
}
void seg_m_x(int t)
{
if (mn[t] >= d) return;
if (l <= lb[t] && rb[t] <= r && sn[t] > d)
{
m_x(t, d);
return;
}
psdw(t);
if (l <= rb[lc(t)]) seg_m_x(lc(t));
if (r >= lb[rc(t)]) seg_m_x(rc(t));
upd(t);
}
void seg_m_n(int t)
{
if (mx[t] <= d) return;
if (l <= lb[t] && rb[t] <= r && sx[t] < d)
{
m_n(t, d);
return;
}
psdw(t);
if (l <= rb[lc(t)]) seg_m_n(lc(t));
if (r >= lb[rc(t)]) seg_m_n(rc(t));
upd(t);
}
LL get_sum(int t)
{
if (l <= lb[t] && rb[t] <= r)
{
return sum[t];
}
LL d = ;
psdw(t);
if (l <= rb[lc(t)]) d += get_sum(lc(t));
if (r >= lb[rc(t)]) d += get_sum(rc(t));
return d;
}
int get_max(int t)
{
if (l <= lb[t] && rb[t] <= r)
{
return mx[t];
}
int d = -INF;
psdw(t);
if (l <= rb[lc(t)]) d = max(d, get_max(lc(t)));
if (r >= lb[rc(t)]) d = max(d, get_max(rc(t)));
return d;
}
int get_min(int t)
{
if (l <= lb[t] && rb[t] <= r)
{
return mn[t];
}
int d = INF;
psdw(t);
if (l <= rb[lc(t)]) d = min(d, get_min(lc(t)));
if (r >= lb[rc(t)]) d = min(d, get_min(rc(t)));
return d;
}
} int main()
{
read(n);
build(, , n);
read(m);
while (m --)
{
using namespace segment;
int opt;
read(opt); read(l); read(r);
if (opt <= ) read(d);
switch (opt)
{
case : seg_add(); break;
case : seg_m_x(); break;
case : seg_m_n(); break;
case : write(get_sum()); break;
case : write(get_max()); break;
case : write(get_min()); break;
}
}
//write(c, ' '); write(d);
}

bzoj 4695: 最假女选手的更多相关文章

  1. bzoj 4695 最假女选手 吉利线段树

    最假女选手 Time Limit: 50 Sec  Memory Limit: 128 MBSubmit: 480  Solved: 118[Submit][Status][Discuss] Desc ...

  2. BZOJ.4695.最假女选手(线段树 Segment tree Beats!)

    题目链接 区间取\(\max,\ \min\)并维护区间和是普通线段树无法处理的. 对于操作二,维护区间最小值\(mn\).最小值个数\(t\).严格次小值\(se\). 当\(mn\geq x\)时 ...

  3. bzoj 4695: 最假女选手 && Gorgeous Sequence HDU - 5306 && (bzoj5312 冒险 || 小B的序列) && bzoj4355: Play with sequence

    算导: 核算法 给每种操作一个摊还代价(是手工定义的),给数据结构中某些东西一个“信用”值(不是手动定义的,是被动产生的),摊还代价等于实际代价+信用变化量. 当实际代价小于摊还代价时,增加等于差额的 ...

  4. BZOJ 4695 最假女选手 线段树

    题意: 给定一个长度为 N序列,编号从1 到 N.要求支持下面几种操作: 1.给一个区间[L,R] 加上一个数x  2.把一个区间[L,R] 里小于x 的数变成x  3.把一个区间[L,R] 里大于x ...

  5. BZOJ4695 最假女选手(势能线段树)

    BZOJ题目传送门 终于体会到初步掌握势能分析思想的重要性了. 一开始看题,感觉套路还是很一般啊qwq.直接在线段树上维护最大值和最小值,每次递归更新的时候,如果不能完全覆盖就暴力递归下去.挺好写的欸 ...

  6. bzoj4695 最假女选手

    传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4695 [题解] SegmentTree beats!(见jiry_2论文/营员交流) 考虑只 ...

  7. BZOJ4695:最假女选手

    浅谈区间最值操作和历史最值问题:https://www.cnblogs.com/AKMer/p/10225100.html 题目传送门:https://lydsy.com/JudgeOnline/pr ...

  8. 2018.07.27 bzoj4695: 最假女选手(线段树)

    传送门 线段树好题 支持区间加,区间取min" role="presentation" style="position: relative;"> ...

  9. 【bzoj4695】最假女选手 线段树区间最值操作

    题目描述 给定一个长度为 N 序列,编号从 1 到 N .要求支持下面几种操作:1.给一个区间[L,R] 加上一个数x 2.把一个区间[L,R] 里小于x 的数变成x 3.把一个区间[L,R] 里大于 ...

随机推荐

  1. 【codevs1106】 篝火晚会

    http://codevs.cn/problem/1106/ (题目链接) 题意 将1~n顺序排列的环改成另一个环,问n-不动点数. Solution 啊智障啦,不会做×_× 左转hzwer 代码 / ...

  2. [原]CentOS7.2最小安装环境部署Asp.NET Core笔记

    转载请注明原作者(think8848)和出处(http://think8848.cnblogs.com) 写在前面的话 不知不觉在cnblogs上注册已经10多年了,看我的园龄就直接暴露了我实际年龄, ...

  3. [LeetCode] Third Maximum Number 第三大的数

    Given a non-empty array of integers, return the third maximum number in this array. If it does not e ...

  4. Dapper学习笔记(一)

    https://github.com/StackExchange/dapper-dot-net Dapper是对IDbConnection的扩展,需要使用Dapper提供的扩展只需要把SqlMappe ...

  5. Android开发之延时执行

    new Handler().postDelayed(new Runnable() { @Override public void run() { //两秒后执行 startMainActivity() ...

  6. windows server 2008 各版本号下载地址(微软官网)

    前言: 微软官网上下载系统的镜像文件要远比百度网盘下载起来得更快. Windows Server 2008 32-bit Standard(标准版)

  7. C#图像处理笔记

    1.灰度拉伸 灰度拉伸又叫对比度拉伸,它是最基本的一种灰度变换,使用的是最简单的分段线性变换函数,它的主要思想是提高图像处理时灰度级的动态范围.

  8. OleDbDataReader快速数据读取方式

    查询得到OleDbDataReader后,有三种方式支持数据读取,如下: //方法一**速度中等 OleDbDataReader reader = command.ExecuteReader(); w ...

  9. R 查看函数源代码

    1:输入函数名. 2:methods( ),例如methods(quantile). 3:针对方法2得到的带*的,例stats:::quantile.default.

  10. ubuntu16.04文件形式安装mongodb

    下载文件:mongodb-linux-x86_64-ubuntu1604-3.4.1.tgz,解压到home目录. 在mongodb目录下新建data/db目录. 在桌面新建一个shell文件run- ...