[BZOJ3064][Tyvj1518] CPU监控

题目：[BZOJ3064][Tyvj1518] CPU监控

思路：

线段树专题讲的。以下为讲课时的课件：

给出序列，要求查询一些区间的最大值、历史最大值，支持区间加、区间修改。序列长度和操作数<=1e5。
维护标记的线段树裸题。出这道题的目的就是为了让大家加强对于线段树标记下传的理解，以及线段树维护信息时分类讨论一定要明确。
在确定了要写线段树后，先考虑要维护哪些标记。然后思考这些标记应当如何维护，何时下传，何时修改。特别注意，当区间加、区间覆盖操作同时存在时，标记之间的先后处理顺序非常重要。
这道题细节非常多，难点集中于标记下传。
具体做法：对于既有set又有add的线段树，在任何时刻，同一个节点上必然不能同时出现set和add两种标记，因为set和add会由顺序不同导致冲突。
• 我们维护六个信息：
• mx[o]表示o点当前最大值；
• gmx[o]表示o点全局最大值；
• add[o]表示o点当前增加的值；
• gadd[o]表示o点从上一次pushdown后直到现在，出现的最大add值；
• set[o]表示o点当前赋的值；
• gset[o]表示o点从上一次pushdown后直到现在，出现的最大set值。
• ※对于不存在set值的节点，均赋为-inf。
接着，先不管查询历史操作，只考虑标记的合并与传递。既然两种标记不能并存，那在传递的时候一定要讨论当前节点和子节点的标记情况，分情况传递。
• 1.父节点有add，子节点有add：直接将父节点的add累加到子节点上。
• 2.父节点有add，子节点有set：将父节点的add值累加到子节点的set值上。
• 3.父节点有set：无论子节点的情况，直接覆盖子节点的set值，同时清掉子节点的add值。
这样，每次访问到一个非叶子节点时，首先进行标记下传更新子节点信息，清空当前节点标记，而每次修改时，如果区间被完全包含，则直接给节点打上标记，更新当前节点信息即可。
最后来考虑历史信息的维护。我们需要在进行操作时，不遗漏没有进行操作的节点出现的最大值，所以维护了“从上次pushdown该节点直到现在，区间中出现的最大set与add标记”，这样，区间的历史最值（答案、set、add）就可能来自于父节点没有pushdown的最大set或者add，在pushdown父节点的时候考虑就好了。另外，在完成对子节点set、add标记更新之后还要再考虑对子节点历史标记的影响。这样，即可保证访问到的节点的信息
都是最新的。
• 你会写到自闭的。
• 就算是看标程也要自己写一遍。
• 做维护信息比较复杂的线段树题目时，一定要找对维护标记的时机，清晰而彻底的认真思考会不会有遗漏，可以避免长时间的盲目调试。

方法一：
果然写自闭了...
最后差不多是对着题解看上几遍然后抄的。
最初容易把问题想简单，以为只要随便用个gmax表示历史最大值，然后就能直接维护了。
实际上，考虑这种情况：我们对一段区间曾做过一些操作，标记停留在某层还没有被下放。之后又进行一项新的操作，把原有的标记覆盖掉了，于是在维护子节点的gmax值时就会出错。
那么该怎么做呢？
为了防止信息丢失，必须保证每次pushdown之前，还没有被下放的标记能够更新子节点的gmax，然后才能把这个标记覆盖掉。
如果每次pushdown都先把原来标记下放一遍，复杂度会退化->TLE.
我们可以维护“从上次pushdown到现在为止没有被下放的最大add值和set值”，问题就解决了。
以上是基本思想，只涉及线段树最基本的lazytag，但具体实现需要对lazy标记的下放有深刻的理解。
题目中有两个操作：add和set，我们考虑下放顺序，共四种情况：

1. add - add
2. set - set
3. set - add
4. add - set
   前两种直接维护，第三种可以把add累加到set上。
   于是能做了，有的题解就是直接分两类讨论：add-add，add-set。
   这里我们也可以按讲课时的方法，用set盖掉add，于是任意时刻线段树节点上只会出现一种标记。可以在提交的代码中加入断言进行验证。
   注意区分：
   mx、add、set标记只会在Add操作、Set操作、pushdown时父节点的add和set标记时更新，gmax、gadd、gset在各种操作时都可能被更新。
   把不同的更新分别写成函数，可以帮助整理思维。

---

Code：

#include <bits/stdc++.h>
using namespace std;
const int N=1e6,inf=0x3f3f3f3f;
int n,T,ans_max,ans_gmax;
struct Tree{
	int l,r,mx,gmx,add,gadd,set,gset;
#define l(p) (node[p].l)
#define r(p) (node[p].r)
#define mx(p) (node[p].mx)
#define gmx(p) (node[p].gmx)
#define add(p) (node[p].add)
#define gadd(p) (node[p].gadd)
#define set(p) (node[p].set)
#define gset(p) (node[p].gset)
#define ls(p) (p<<1)
#define rs(p) (p<<1|1)
#define mid ((l(p)+r(p))>>1)
}node[N<<2];
inline void cmax(int &a,int b){ a=((b>a)?b:a); }
void pup(int p){
	mx(p)=max(mx(ls(p)),mx(rs(p)));
	cmax(gmx(p),max(gmx(ls(p)),gmx(rs(p))));
}
void build(int p,int l,int r){
	l(p)=l;r(p)=r;
	set(p)=gset(p)=-inf;
	if(l==r){
		scanf("%lld",&mx(p));
		gmx(p)=mx(p);
		return;
	}
	build(ls(p),l,mid);
	build(rs(p),mid+1,r);
	pup(p);
}
void m_add(int p,int val){/*p节点实际增加val*/
	cmax(gmx(p),mx(p)+=val);
	if(set(p)!=-inf) cmax(gset(p),set(p)+=val);
	else cmax(gadd(p),add(p)+=val);
}
void m_set(int p,int val){/*p节点实际被val覆盖*/
	cmax(gmx(p),mx(p)=val);
	cmax(gset(p),set(p)=val);
	add(p)=0;
}
void g_add(int p,int val){/*用父亲的gadd更新p节点*/
	cmax(gmx(p),mx(p)+val);
	if(set(p)!=-inf) cmax(gset(p),set(p)+val);
	else cmax(gadd(p),add(p)+val);
}
void g_set(int p,int val){/*用父亲的gset更新p节点*/
	cmax(gmx(p),val);
	cmax(gset(p),val);
}
void pdown(int p){
	if(gadd(p)){
		g_add(ls(p),gadd(p));
		g_add(rs(p),gadd(p));
		gadd(p)=0;
	}
	if(gset(p)!=-inf){
		g_set(ls(p),gset(p));
		g_set(rs(p),gset(p));
		gset(p)=-inf;
	}
	if(add(p)){
		m_add(ls(p),add(p));
		m_add(rs(p),add(p));
		add(p)=0;
	}
	if(set(p)!=-inf){
		m_set(ls(p),set(p));
		m_set(rs(p),set(p));
		set(p)=-inf;
	}
}
void Add(int p,int L,int R,int val){
	if(l(p)>R||r(p)<L) return;
	if(L<=l(p)&&r(p)<=R) return m_add(p,val); 
	pdown(p);
	Add(ls(p),L,R,val);
	Add(rs(p),L,R,val);
	pup(p);
}
void Set(int p,int L,int R,int val){
	if(l(p)>R||r(p)<L) return;
	if(L<=l(p)&&r(p)<=R) return m_set(p,val);
	pdown(p);
	Set(ls(p),L,R,val);
	Set(rs(p),L,R,val);
	pup(p);
}
void query(int p,int L,int R){
	if(l(p)>R||r(p)<L) return;
	if(L<=l(p)&&r(p)<=R) {
		ans_max=max(ans_max,mx(p));
		ans_gmax=max(ans_gmax,gmx(p));
		return;
	}
	pdown(p);
	query(ls(p),L,R);
	query(rs(p),L,R);
}
int main(){
	scanf("%d",&n);
	build(1,1,n);
	scanf("%d",&T);
	char op[5];
	int x,y,z;
	while(T--){
		ans_max=ans_gmax=-inf;
		scanf("%s",op);
		if(op[0]=='Q') {
			scanf("%d%d",&x,&y);
			query(1,x,y);
			printf("%d
",ans_max);
		} else if(op[0]=='A') {
			scanf("%d%d",&x,&y);
			query(1,x,y);
			printf("%d
",ans_gmax);
		} else if(op[0]=='P') {
			scanf("%d%d%d",&x,&y,&z);
			Add(1,x,y,z);
		} else if(op[0]=='C') {
			scanf("%d%d%d",&x,&y,&z);
			Set(1,x,y,z);
		}
	}
	return 0;
}

---

方法二：
康题解的时候发现还有一种巧妙的标记。
以下转载某dalao的博客。

定义标记 ((a,b)) 表示先给区间内所有数加上 (a) ，再与 (b) 取 (max) ，即 (Ai=max(Ai+a,b))。
那么区间加操作就相当于打标记 ((z,−infty))，区间赋值操作就相当于打标记((-infty ,z))
考虑两个标记合并：将标记 ((c,d))合并到 ((a,b)) 上（即先有 ((a,b)) 后有 ((c,d))），可以得到结果 ((a+c,max(b+c,d)))。
考虑两个标记取最大值：将标记作用后的值看作原值的函数，那么显然是一个分段函数，第一段是 (y=b)，第二段是 (y=x+a) 。那么标记取最大值就是取两个函数最上方的部分，显然是两段分别取最大值，即 ((a,b)) 与 ((c,d)) 取最大值后为 ((max(a,c),max(b,d))) 。

这种标记很妙啊，好像又是出自集训队论文，抽空再研究。
先对着题解抄一份。
注意细节：两个标记合并时，add与-inf取max，防止赋值操作add值连加几个-inf爆了int。

---

Code

#include <bits/stdc++.h>
using namespace std;
const int N=1e6,inf=0x3f3f3f3f;
int n,T,ans_max,ans_gmax;
struct tag{
	int add,set;
	inline tag(int add=0,int set=-inf):add(add),set(set){}
	inline tag operator + (const tag &b){
		return tag(max(-inf,add+b.add),max(set+b.add,b.set));
	}
	inline tag operator & (const tag &b){
		return tag(max(add,b.add),max(set,b.set));
	} 
};
struct Tree{
	int l,r,mx,gmx;tag now,pre;
#define l(p) (node[p].l)
#define r(p) (node[p].r)
#define mx(p) (node[p].mx)
#define gmx(p) (node[p].gmx)
#define now(p) (node[p].now)
#define pre(p) (node[p].pre)
#define ls(p) (p<<1)
#define rs(p) (p<<1|1)
#define mid ((l(p)+r(p))>>1)
}node[N<<2];
inline void cmax(int &a,int b){ a=((b>a)?b:a); }
void pup(int p){
	mx(p)=max(mx(ls(p)),mx(rs(p)));
	cmax(gmx(p),max(gmx(ls(p)),gmx(rs(p))));
}
void build(int p,int l,int r){
	l(p)=l;r(p)=r;
	now(p)=pre(p)=tag();
	if(l==r){
		scanf("%lld",&mx(p));
		gmx(p)=mx(p);
		return;
	}
	build(ls(p),l,mid);
	build(rs(p),mid+1,r);
	pup(p);
}
void calc(int p,tag n,tag h){
	pre(p)=pre(p)&(now(p)+h);
	now(p)=now(p)+n;
	cmax(gmx(p),max(mx(p)+h.add,h.set));
	mx(p)=max(mx(p)+n.add,n.set);
}
void pdown(int p){
	calc(ls(p),now(p),pre(p));
	calc(rs(p),now(p),pre(p));
	now(p)=pre(p)=tag();
}
void change(int p,int L,int R,tag key){
	if(l(p)>R||r(p)<L) return;
	if(L<=l(p)&&r(p)<=R) return calc(p,key,key); 
	pdown(p);
	change(ls(p),L,R,key);
	change(rs(p),L,R,key);
	pup(p);
}
void query(int p,int L,int R){
	if(l(p)>R||r(p)<L) return;
	if(L<=l(p)&&r(p)<=R) {
		ans_max=max(ans_max,mx(p));
		ans_gmax=max(ans_gmax,gmx(p));
		return;
	}
	pdown(p);
	query(ls(p),L,R);
	query(rs(p),L,R);
}
int main(){
	scanf("%d",&n);
	build(1,1,n);
	scanf("%d",&T);
	char op[5];
	int x,y,z;
	while(T--){
		ans_max=ans_gmax=-inf;
		scanf("%s",op);
		if(op[0]=='Q') {
			scanf("%d%d",&x,&y);
			query(1,x,y);
			printf("%d
",ans_max);
		} else if(op[0]=='A') {
			scanf("%d%d",&x,&y);
			query(1,x,y);
			printf("%d
",ans_gmax);
		} else if(op[0]=='P') {
			scanf("%d%d%d",&x,&y,&z);
			change(1,x,y,tag(z,-inf));
		} else if(op[0]=='C') {
			scanf("%d%d%d",&x,&y,&z);
			change(1,x,y,tag(-inf,z));
		}
	}
	return 0;
}