浅谈区间最值操作和历史最值问题:https://www.cnblogs.com/AKMer/p/10225100.html

题目传送门:https://lydsy.com/JudgeOnline/problem.php?id=4695

吉司机线段树板子大集合。所有信息都封装在一个结构体里会比开多个数组快14秒。

注意暴力\(dfs\)子树时要\(pushdown\)。

时间复杂度:\(O(nlog^2n)\)

空间复杂度:\(O(n)\)

代码如下:

#include <cstdio>
#include <algorithm>
using namespace std;
typedef long long ll; const int maxn=5e5+5,inf=1e9; int n,m;
int a[maxn]; inline int read() {
int x=0,f=1;char ch=getchar();
for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+ch-'0';
return x*f;
} struct segment_tree {
struct tree_node {
ll sum;
int cntA,cntZ,tag,A,B,Y,Z;//A最大值,B严格次大值,Y严格次小值,Z最小值
}tree[maxn<<2]; inline void update(int p) {
tree[p].A=max(tree[p<<1].A,tree[p<<1|1].A);
tree[p].Z=min(tree[p<<1].Z,tree[p<<1|1].Z);
tree[p].sum=tree[p<<1].sum+tree[p<<1|1].sum;
tree[p].cntA=(tree[p<<1].A==tree[p].A)*tree[p<<1].cntA;
tree[p].cntA+=(tree[p<<1|1].A==tree[p].A)*tree[p<<1|1].cntA;
tree[p].cntZ=(tree[p<<1].Z==tree[p].Z)*tree[p<<1].cntZ;
tree[p].cntZ+=(tree[p<<1|1].Z==tree[p].Z)*tree[p<<1|1].cntZ;
tree[p].B=tree[p<<1|1].A!=tree[p].A?tree[p<<1|1].A:tree[p<<1|1].B;
tree[p].B=max(tree[p].B,tree[p<<1].A!=tree[p].A?tree[p<<1].A:tree[p<<1].B);
tree[p].Y=tree[p<<1|1].Z!=tree[p].Z?tree[p<<1|1].Z:tree[p<<1|1].Y;
tree[p].Y=min(tree[p].Y,tree[p<<1].Z!=tree[p].Z?tree[p<<1].Z:tree[p<<1].Y);
} inline void build(int p,int l,int r) {
if(l==r) {
tree[p].cntA=tree[p].cntZ=1;
tree[p].sum=tree[p].A=tree[p].Z=a[l];
tree[p].B=-inf,tree[p].Y=inf;
return;
}
int mid=(l+r)>>1;
build(p<<1,l,mid);
build(p<<1|1,mid+1,r);
update(p);
} inline void add_tag(int p,int l,int r,int v) {
tree[p].A+=v,tree[p].Z+=v,tree[p].tag+=v;
if(tree[p].Y!=inf)tree[p].Y+=v;
if(tree[p].B!=-inf)tree[p].B+=v;
tree[p].sum+=1ll*(r-l+1)*v;
} inline void max_tag(int p,int v) {
tree[p].sum+=1ll*tree[p].cntZ*(v-tree[p].Z),tree[p].Z=v;
if(tree[p].Y==inf)tree[p].A=v;
else tree[p].B=max(tree[p].B,v);
} inline void min_tag(int p,int v) {
tree[p].sum-=1ll*tree[p].cntA*(tree[p].A-v),tree[p].A=v;
if(tree[p].B==-inf)tree[p].Z=v;
else tree[p].Y=min(tree[p].Y,v);
} inline void solveMax(int p,int l,int r,int limit) {
if(limit<=tree[p].Z)return;
if(limit>tree[p].Z&&limit<tree[p].Y) {
max_tag(p,limit);
return;
}
int mid=(l+r)>>1;push_down(p,l,r);
solveMax(p<<1,l,mid,limit);
solveMax(p<<1|1,mid+1,r,limit);
update(p);
} inline void solveMin(int p,int l,int r,int limit) {
if(limit>=tree[p].A)return;
if(limit<tree[p].A&&limit>tree[p].B) {
min_tag(p,limit);
return;
}
int mid=(l+r)>>1;push_down(p,l,r);
solveMin(p<<1,l,mid,limit);
solveMin(p<<1|1,mid+1,r,limit);
update(p);
} inline void push_down(int p,int l,int r) {
int mid=(l+r)>>1;
if(tree[p].tag) {
add_tag(p<<1,l,mid,tree[p].tag);
add_tag(p<<1|1,mid+1,r,tree[p].tag);
tree[p].tag=0;
}
solveMin(p<<1,l,mid,tree[p].A);
solveMin(p<<1|1,mid+1,r,tree[p].A);
solveMax(p<<1,l,mid,tree[p].Z);
solveMax(p<<1|1,mid+1,r,tree[p].Z);
} inline void ADD(int p,int l,int r,int L,int R,int v) {
if(L<=l&&r<=R) {
add_tag(p,l,r,v);
return;
}
int mid=(l+r)>>1;push_down(p,l,r);
if(L<=mid)ADD(p<<1,l,mid,L,R,v);
if(R>mid)ADD(p<<1|1,mid+1,r,L,R,v);
update(p);
} inline void MAX(int p,int l,int r,int L,int R,int v) {
if(tree[p].Z>=v)return;
if(L<=l&&r<=R) {
solveMax(p,l,r,v);
return;
}
int mid=(l+r)>>1;push_down(p,l,r);
if(L<=mid)MAX(p<<1,l,mid,L,R,v);
if(R>mid)MAX(p<<1|1,mid+1,r,L,R,v);
update(p);
} inline void MIN(int p,int l,int r,int L,int R,int v) {
if(tree[p].A<=v)return;
if(L<=l&&r<=R) {
solveMin(p,l,r,v);
return;
}
int mid=(l+r)>>1;push_down(p,l,r);
if(L<=mid)MIN(p<<1,l,mid,L,R,v);
if(R>mid)MIN(p<<1|1,mid+1,r,L,R,v);
update(p);
} inline ll querySum(int p,int l,int r,int L,int R) {
if(L<=l&&r<=R)return tree[p].sum;
int mid=(l+r)>>1;ll res=0;push_down(p,l,r);
if(L<=mid)res=querySum(p<<1,l,mid,L,R);
if(R>mid)res+=querySum(p<<1|1,mid+1,r,L,R);
return res;
} inline int queryMax(int p,int l,int r,int L,int R) {
if(L<=l&&r<=R)return tree[p].A;
int mid=(l+r)>>1,res=-inf;push_down(p,l,r);
if(L<=mid)res=max(res,queryMax(p<<1,l,mid,L,R));
if(R>mid)res=max(res,queryMax(p<<1|1,mid+1,r,L,R));
return res;
} inline int queryMin(int p,int l,int r,int L,int R) {
if(L<=l&&r<=R)return tree[p].Z;
int mid=(l+r)>>1,res=inf;push_down(p,l,r);
if(L<=mid)res=min(res,queryMin(p<<1,l,mid,L,R));
if(R>mid)res=min(res,queryMin(p<<1|1,mid+1,r,L,R));
return res;
}
}T; int main() {
n=read();
for(int i=1;i<=n;i++)
a[i]=read();
T.build(1,1,n);
m=read();
for(int i=1;i<=m;i++) {
int opt=read(),l=read(),r=read(),v=(opt<4?read():0);
if(opt==1)T.ADD(1,1,n,l,r,v);
if(opt==2)T.MAX(1,1,n,l,r,v);
if(opt==3)T.MIN(1,1,n,l,r,v);
if(opt==4)printf("%lld\n",T.querySum(1,1,n,l,r));
if(opt==5)printf("%d\n",T.queryMax(1,1,n,l,r));
if(opt==6)printf("%d\n",T.queryMin(1,1,n,l,r));
}
return 0;
}

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

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

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

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

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

  3. bzoj4695 最假女选手

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

  4. [BZOJ4695]最假女选手:segment tree beats!

    分析 segment tree beats!模板题. 看了gxz的博客突然发现自己写的mxbt和mnbt两个标记没用诶. 代码 #include <bits/stdc++.h> #defi ...

  5. bzoj4695 最假女选手(势能线段树/吉司机线段树)题解

    题意: 已知\(n\)个数字,进行以下操作: \(1.\)给一个区间\([L,R]\) 加上一个数\(x\) \(2.\)把一个区间\([L,R]\) 里小于\(x\) 的数变成\(x\) \(3.\ ...

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

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

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

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

  8. 【bzoj4695】最假女选手

    zcy的励志故事.jpg 傻逼zcy突然想立一个flag,写一个segment-tree-beats的题娱乐一下 于是他就想起了这道题. 他打算今晚写完 然后光是写他就写的头昏脑涨,还犯了询问写反这种 ...

  9. (WC2016模拟十一)【BZOJ4695】最假女选手

    ps:好久没更博啦……这几天连着有模拟赛,等初赛前后休息的时候来疯狂补坑吧……顺便补一下前面的数论啥的? 题解: mdzz我场上写了个15分暴力长度跟标算差不多... 线段树大法好啊!这题听说很多人做 ...

随机推荐

  1. One Person Game(概率+数学)

    There is a very simple and interesting one-person game. You have 3 dice, namelyDie1, Die2 and Die3.  ...

  2. 大数据学习系列(7)-- hadoop集群搭建

    1.配置ssh免登陆 #进入到我的home目录 cd ~/.ssh ssh-keygen -t rsa 执行完这个命令后,会生成两个文件id_rsa(私钥).id_rsa.pub(公钥) 将公钥拷贝到 ...

  3. Jaxb2 实现JavaBean与xml互转

    一.简介 JAXB(Java Architecture for XML Binding) 是一个业界的标准,是一项可以根据XML Schema产生Java类的技术.该过程中,JAXB也提供了将XML实 ...

  4. zookeepeer ID生成器 (一)

    目录 写在前面 1.1. ZK 的分布式命名服务 1.1.1. 分布式 ID 生成器的类型 UUID方案 1.1.2. ZK生成分布式ID 写在最后 疯狂创客圈 亿级流量 高并发IM 实战 系列 疯狂 ...

  5. Maven项目启动报错

    错误信息如下: 六月 , :: 下午 org.apache.tomcat.util.digester.SetPropertiesRule begin 警告: [SetPropertiesRule]{S ...

  6. 003-unity3d 物理引擎简介以及示例

    一.概述 物理引擎就是模拟真实世界中物体碰撞.跌落等反应的引擎,通过ballence.愤怒的小鸟等理解.Unity3D的物理引擎使用的是Nvidia的PhysX. 物理引擎是一个计算机程序模拟牛顿力学 ...

  7. 变分推断(Variational Inference)

    变分 对于普通的函数f(x),我们可以认为f是一个关于x的一个实数算子,其作用是将实数x映射到实数f(x).那么类比这种模式,假设存在函数算子F,它是关于f(x)的函数算子,可以将f(x)映射成实数F ...

  8. eclipse svn 分支合并到主干

    最近公司产品上线,整个系统架构包含有七八个子系统,并且子系统都是集群部署.所以每次升级维护都要确保尽可能不出问题.因为整个系统刚上线不久,意味着新系统不定期有BUG需修复,但新功能模块也在持续的开发中 ...

  9. VHDL基础 学习笔记

    最近一直忙着学校里的活动,所以没怎么更新,上周活动忙完了,正好也借着数电实验的机会,重新学习一下VHDL的编程.以下是转自360doc的教程: ———————————————————————————— ...

  10. vim 操作符命令和位移(如删除和修改)

    一.字符删除命令:x,d,D x 删除光标下的字符,前面可以加入命令计数,如:5x 代表删除从当前光标到后面的5个字符,包括空格: X 删除光标前面的一个字符: dw 删除光标下到word后的数据: ...