zkw线段树
zkw线段树
————循环写法
zkw线段树不能处理有运算优先级的问题(加法乘法混合处理)
我们观察一下递归式线段树的代码,很容易就会发现:无论是建树、修改还是查询,都是自顶向下的。
zkw线段树则正好反过来,即自底向上
https://blog.csdn.net/keshuqi/article/details/52205884
https://www.cnblogs.com/Judge/p/9514862.html#_label0
#include<cstdio>
#include<iostream>
#define ll long long
using namespace std;
const int M=1e5+111;
char buf[1<<21],*p1=buf,*p2=buf;
inline ll read(){
ll x=0,f=1; char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;
}
char sr[1<<21],z[20];int C=-1,Z;
inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
inline void print(ll x){
if(C>1<<20)Ot();if(x<0)sr[++C]=45,x=-x;
while(z[++Z]=x%10+48,x/=10);
while(sr[++C]=z[Z],--Z);sr[++C]='\n';
}
ll n,m,q;
ll sum[M<<2],add[M<<2];
inline void build(){
for(m=1;m<=n;m<<=1);
for(int i=m+1;i<=m+n;++i) sum[i]=read();
for(int i=m-1;i;--i) sum[i]=sum[i<<1]+sum[i<<1|1];
}
inline void update_part(int s,int t,ll v){
ll A=0,lc=0,rc=0,len=1;
for(s+=m-1,t+=m+1;s^t^1;s>>=1,t>>=1,len<<=1){
if(s&1^1) add[s^1]+=v,lc+=len;
if(t&1) add[t^1]+=v,rc+=len;
sum[s>>1]+=v*lc,sum[t>>1]+=v*rc;
} for(lc+=rc,s>>=1;s;s>>=1) sum[s]+=v*lc;
}
inline ll query_sum(int s,int t){
ll lc=0,rc=0,len=1,ans=0;
for(s+=m-1,t+=m+1;s^t^1;s>>=1,t>>=1,len<<=1){
if(s&1^1) ans+=sum[s^1]+len*add[s^1],lc+=len;
if(t&1) ans+=sum[t^1]+len*add[t^1],rc+=len;
if(add[s>>1]) ans+=add[s>>1]*lc;
if(add[t>>1]) ans+=add[t>>1]*rc;
} for(lc+=rc,s>>=1;s;s>>=1) if(add[s]) ans+=add[s]*lc;
return ans;
}
signed main(){
n=read(),q=read(),build();
int opt,x,y; ll k;
while(q--){
opt=read(),x=read(),y=read();
if(opt&1) k=read(),update_part(x,y,k);
else print(query_sum(x,y));
} Ot(); return 0;
}
建树
inline void build(int n){
for(m=1;m<=n;m<<=1);
for(int i=m+1;i<=m+n;++i)
sum[i]=mn[i]=mx[i]=read();
//建树
for(int i=m-1;i;--i){
sum[i]=sum[i<<1]+sum[i<<1|1];
mn[i]=min(mn[i<<1],mn[i<<1|1]),
mn[i<<1]-=mn[i],mn[i<<1|1]-=mn[i];
mx[i]=max(mx[i<<1],mx[i<<1|1]),
mx[i<<1]-=mx[i],mx[i<<1|1]-=mx[i];
}
//维护区间和,最大、最小值
}
单点修改
inline void upd(int x,int v,int A=0){
x+=m,mx[x]+=v,mn[x]+=v;
for(;x>1;x>>=1){
sum[x]+=v;
A=min(mn[x],mn[x^1]);
mn[x]-=A,mn[x^1]-=A,mn[x>>1]+=A;
A=max(mx[x],mx[x^1]),
mx[x]-=A,mx[x^1]-=A,mx[x>>1]+=A;
}
}