1022 - 线段树 - Rikka with Sequence（HDU 5828）

传送门

题意

维护区间开方，区间加，区间求和。

分析

区间开方：对于一个区间，显然在开几次方之后区间的极差会很小，并且对整个区间进行区间加操作对区间的极差并没有什么影响，这个极差随着修改数的增多仍然是单调不增的，也就是说，在对这个区间进行了一个近似于常数的开方操作之后整个区间的极差就不会超过1了。
那么可以维护区间最值amax,amin，每次区间开方的时候，先看当前的amax-amin是否不超过1，如果超过就递归修改。
否则，我们令当前的amax开方之后变成bmax,amin开方之后变成bmin
接下来判断，如果amax-bmax=amin-bmin，相当于给整个区间打一个区间减标记。
否则，就只有一种可能了，amax=amin+1&&bmax=bmin，这样就相当于一个区间覆盖操作。
这样实际上我们将原问题转化成了区间加，区间覆盖，询问区间和的简单情况。
摘自ldx讲稿

下面结合自己的理解再粗浅的谈一谈
一般情况下，我们遇到区间开方的问题都只涉及这一种修改操作，那么由于我们知道一个数在$log(log(n))$次的开方后就会变为1，我们就维护一个区间最大值，当这个最大值大于1的时候递归修改，否则不予以操作
但这个题有点点……嗯……讨厌，它还要求区间加，显然上面的方法此刻就有问题了
再仔细想一下，维护区间的最小值和最大值，就按照ldx大佬的方法操作咯

犯了两个sb错误

代码

#include<bits/stdc++.h>
#define in read()
#define N 100009
#define lc (k<<1)
#define rc (k<<1)|1
#define ll long long
using namespace std;
inline int read(){
	char ch;int f=1,res=0;
	while((ch=getchar())<'0'||ch>'9') if(ch=='-') f=-1;
	while(ch>='0'&&ch<='9'){
		res=(res<<3)+(res<<1)+ch-'0';
		ch=getchar();
	}
	return f==1?res:-res;
}
int t,n,m,a[N];
ll sum[N<<2],maxn[N<<2],minn[N<<2],lzy[N<<2],tag[N<<2]; 
void pushup(int k){
	sum[k]=sum[lc]+sum[rc];
	maxn[k]=max(maxn[lc],maxn[rc]);
	minn[k]=min(minn[lc],minn[rc]);
	return ;
}
void Cover(int k,int l,int r,ll v){
	lzy[k]=0;
	sum[k]=(r-l+1)*v;
	tag[k]=maxn[k]=minn[k]=v;
}
void Add(int k,int l,int r,ll v){
	sum[k]+=(r-l+1)*v;
	maxn[k]+=v;minn[k]+=v;
	lzy[k]+=v;
}
void pushdown(int k,int l,int r){
	int mid=l+r>>1;
	if(tag[k]!=-1) Cover(lc,l,mid,tag[k]),Cover(rc,mid+1,r,tag[k]);
	if(lzy[k]) Add(lc,l,mid,lzy[k]),Add(rc,mid+1,r,lzy[k]);
	tag[k]=-1;lzy[k]=0;
}
void build(int k,int l,int r){
	if(l==r){
		sum[k]=maxn[k]=minn[k]=a[l];
		return;
	}
	int mid=l+r>>1;
	build(lc,l,mid);build(rc,mid+1,r);
	pushup(k);
}
void modify(int k,int l,int r,int x,int y,ll v){//qujian jia
	if(x<=l&&r<=y){
		Add(k,l,r,v);return;
	}
	pushdown(k,l,r);
	int mid=l+r>>1;
	if(x<=mid) modify(lc,l,mid,x,y,v);
	if(y>mid) modify(rc,mid+1,r,x,y,v);
	pushup(k);
}
void Sqrt(int k,int l,int r,int x,int y){
	if(x>r||y<l) return;
	if(x<=l&&r<=y&&maxn[k]-minn[k]<=1){
		ll bmax=sqrt(maxn[k]),bmin=sqrt(minn[k]);
		if(bmax!=bmin)	Add(k,l,r,bmax-maxn[k]);
		else	Cover(k,l,r,bmax);
		return;///
	}
	pushdown(k,l,r);
	int mid=l+r>>1;
	if(x<=mid) Sqrt(lc,l,mid,x,y);
	if(y>mid) Sqrt(rc,mid+1,r,x,y);
	pushup(k);
}
ll query(int k,int l,int r,int x,int y){
	if(x<=l&&r<=y) return sum[k];
	pushdown(k,l,r);
	int mid=l+r>>1;ll res=0;
	if(x<=mid) res+=query(lc,l,mid,x,y);
	if(y>mid) res+=query(rc,mid+1,r,x,y);
	return res;
}
int main(){
	t=in;
	while(t--){
		memset(tag,-1,sizeof(tag));
		memset(lzy,0,sizeof(lzy));
		n=in;m=in;
		int i,j,op,l,r;
		for(i=1;i<=n;++i) a[i]=in;
		build(1,1,n);
		for(i=1;i<=m;++i){
			op=in;l=in;r=in;
			if(op==1){
				ll x=in;
				modify(1,1,n,l,r,x);//qujian jia
			}
			else if(op==2) Sqrt(1,1,n,l,r);
			else{ printf("%lld\n",query(1,1,n,l,r));	}
		}
	}
	
	return 0;
}