BZOJ4695:最假女选手
浅谈区间最值操作和历史最值问题: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:最假女选手的更多相关文章
- BZOJ4695 最假女选手(势能线段树)
		BZOJ题目传送门 终于体会到初步掌握势能分析思想的重要性了. 一开始看题,感觉套路还是很一般啊qwq.直接在线段树上维护最大值和最小值,每次递归更新的时候,如果不能完全覆盖就暴力递归下去.挺好写的欸 ... 
- 2018.07.27 bzoj4695: 最假女选手(线段树)
		传送门 线段树好题 支持区间加,区间取min" role="presentation" style="position: relative;"> ... 
- bzoj4695 最假女选手
		传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4695 [题解] SegmentTree beats!(见jiry_2论文/营员交流) 考虑只 ... 
- [BZOJ4695]最假女选手:segment tree beats!
		分析 segment tree beats!模板题. 看了gxz的博客突然发现自己写的mxbt和mnbt两个标记没用诶. 代码 #include <bits/stdc++.h> #defi ... 
- bzoj4695 最假女选手(势能线段树/吉司机线段树)题解
		题意: 已知\(n\)个数字,进行以下操作: \(1.\)给一个区间\([L,R]\) 加上一个数\(x\) \(2.\)把一个区间\([L,R]\) 里小于\(x\) 的数变成\(x\) \(3.\ ... 
- bzoj 4695  最假女选手 吉利线段树
		最假女选手 Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 480 Solved: 118[Submit][Status][Discuss] Desc ... 
- 【bzoj4695】最假女选手  线段树区间最值操作
		题目描述 给定一个长度为 N 序列,编号从 1 到 N .要求支持下面几种操作:1.给一个区间[L,R] 加上一个数x 2.把一个区间[L,R] 里小于x 的数变成x 3.把一个区间[L,R] 里大于 ... 
- 【bzoj4695】最假女选手
		zcy的励志故事.jpg 傻逼zcy突然想立一个flag,写一个segment-tree-beats的题娱乐一下 于是他就想起了这道题. 他打算今晚写完 然后光是写他就写的头昏脑涨,还犯了询问写反这种 ... 
- (WC2016模拟十一)【BZOJ4695】最假女选手
		ps:好久没更博啦……这几天连着有模拟赛,等初赛前后休息的时候来疯狂补坑吧……顺便补一下前面的数论啥的? 题解: mdzz我场上写了个15分暴力长度跟标算差不多... 线段树大法好啊!这题听说很多人做 ... 
随机推荐
- mysql 分组加行号
			SELECT end AS row_num_v_, case when @v <> _apply.PATIENT_ID then @v := _apply.PATIENT_ID else ... 
- Random Fourier Features
			K-DSN 深度堆叠网络 Random Features for Large-Scale Kernel Machines To accelerate the training of kernel ma ... 
- 洛谷 P2051 [AHOI2009]中国象棋
			题目描述 这次小可可想解决的难题和中国象棋有关,在一个N行M列的棋盘上,让你放若干个炮(可以是0个),使得没有一个炮可以攻击到另一个炮,请问有多少种放置方法.大家肯定很清楚,在中国象棋中炮的行走方式是 ... 
- java堆分析神器MAT
			Memory Analyzer(MAT) 基于Eclipse的软件 http://www.eclipse.org/mat/ 
- curl使用说明
			默认curl使用get请求,可以使用-d方式指定使用post方式传递数据 https://www.cnblogs.com/gbyukg/p/3326825.html 
- LeetCode:算法特辑——二分搜索
			LeetCode:算法特辑——二分搜索 算法模板——基础 int L =0; int R =arr.length; while(L<R) { int M = (R-L)/2+L; if(arr[ ... 
- UIImageView 获取图片的 宽 高
			该文章纯属这两天开发的经验之谈 并且也是平常没注意 这回发现的一个小方法 并且很实用 在开发中 提高了很大的效率 更加符合高保真的要求 通常 美术 切的一些图片 需要 :1还原的 现在 我们一般支持i ... 
- 【leetcode刷提笔记】Container With Most Water
			Given n non-negative integers a1, a2, ..., an, where each represents a point at coordinate (i, ai). ... 
- 剑指offer——和为s的两个数字VS和为s的连续正数序列
			两种方法都类似于快排的变形. #include <iostream> #include <string> using namespace std; bool FindNumbe ... 
- 算法(Algorithms)第4版 练习 1.5.3
			id数组和treesize数组变化情况: 0 1 2 3 4 5 6 7 8 9 1 1 1 1 1 1 1 1 1 1 10 components 9 0 1 2 3 4 5 6 7 8 9 1 1 ... 
