JZOJ 4933. 【NOIP2017提高组模拟12.24】C

时间:2022-12-17 16:49:42

题目

JZOJ 4933. 【NOIP2017提高组模拟12.24】C

题解

这是一道想都能想到方法的码农题
然而,我们要注意几个要点。
在维护值的时候,我们只用打2个LAZY标志。一个更新最新值,另一个判断是修改还是覆盖。
还有,以后在打线段树的时候,不要开二维数组,尽量用一维。这样,各数组的名称可以直观看出维护哪个值。
重点讲一下第7个操作。
我们维护这样几个值:最左边的数ln,最右边的数rn,左边的lwz,右边的rwz,这段区间的答案lg。
JZOJ 4933. 【NOIP2017提高组模拟12.24】C
当合并的时候,比一下两端区间的lg,然后新的区间的ln,ls等于做区间的ln,ls,新的区间的rn,rs等于做区间的rn,rs。如果左区间的rn=右区间的ln,那么尝试用ls_左区间+rs _右区间 更新区间的lg。
在查询的时候,递归要变成函数,比较一下左边lg大还是右边lg大,如果此区间能够合并,那么看看合并后能不能更新lg(注意越界)。
我说的越界如下图:
JZOJ 4933. 【NOIP2017提高组模拟12.24】C
还有其他很多地方要注意的,十分恶心,希望大家打代码时注意一下。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#define N 100010
#define LL long long
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
LL a[N];
LL mx[N*15],mn[N*15],sum[N*15],sum1[N*15],la[N*15],rp[N*15];
LL ln[N*15],rn[N*15],ls[N*15],rs[N*15],lg[N*15];
LL i,j,K,L,R,n,q,op,ans,cnt;
inline LL read()
{
LL p,data=0;
char ch=0;
while ((ch!='-') && ch<'0' || ch>'9') ch=getchar();
if (ch=='-')
{
p=-1;
ch=getchar();
} else p=1;
while (ch>='0' && ch<='9') data=data*10+ch-'0',ch=getchar();
return data*p;
}
void down(LL s,LL l,LL r)
{
if (rp[s]==1)
{
sum[s]+=la[s]*(r-l+1);
mx[s]+=la[s];
mn[s]+=la[s];
ln[s]+=la[s];
rn[s]+=la[s];
if (!rp[s*2]) rp[s*2]=1;
if (!rp[s*2+1]) rp[s*2+1]=1;
la[s*2]+=la[s];
la[s*2+1]+=la[s];
rp[s]=la[s]=0;
} else
if (rp[s]==2)
{
sum[s]=la[s]*(r-l+1);
mx[s]=mn[s]=ln[s]=rn[s]=la[s];
lg[s]=ls[s]=rs[s]=r-l+1;
rp[s*2]=rp[s*2+1]=2;
la[s*2]=la[s*2+1]=la[s];
rp[s]=la[s]=0;
}
}
void up(LL s,LL l,LL r)
{
sum[s]=sum[s*2]+sum[s*2+1];
mx[s]=max(mx[s*2],mx[s*2+1]);
mn[s]=min(mn[s*2],mn[s*2+1]);
ls[s]=ls[s*2];
rs[s]=rs[s*2+1];
ln[s]=ln[s*2];
rn[s]=rn[s*2+1];
lg[s]=max(lg[s*2],lg[s*2+1]);
if (rn[s*2]==ln[s*2+1])
{
lg[s]=max(lg[s],rs[s*2]+ls[s*2+1]);
if (ls[s*2]==(l+r)/2-l+1) ls[s]+=ls[s*2+1];
if (rs[s*2+1]==r-(l+r)/2) rs[s]+=rs[s*2];
}
}
void build(LL s,LL l,LL r)
{
if (l==r)
{
mx[s]=mn[s]=sum[s]=ln[s]=rn[s]=a[l];
ls[s]=rs[s]=lg[s]=1;
return;
}
LL mid=(l+r)/2;
build(s*2,l,mid);
build(s*2+1,mid+1,r);
up(s,l,r);
}
void inspls(LL s,LL l,LL r,LL x,LL y,LL z)
{
LL wz=(l+r)/2;
down(s*2,l,wz);down(s*2+1,wz+1,r);
if (l==x && r==y)
{
la[s]=z;rp[s]=1;
down(s,l,r);
return;
}
if (y<=wz) inspls(s*2,l,wz,x,y,z);
else if (x>wz) inspls(s*2+1,wz+1,r,x,y,z);
else
{
inspls(s*2,l,wz,x,wz,z);
inspls(s*2+1,wz+1,r,wz+1,y,z);
}
up(s,l,r);
}
void replace(LL s,LL l,LL r,LL x,LL y,LL z)
{
LL wz=(l+r)/2;
down(s*2,l,wz);down(s*2+1,wz+1,r);
if (l==x && r==y)
{
la[s]=z;rp[s]=2;
down(s,l,r);
return;
}
if (y<=wz) replace(s*2,l,wz,x,y,z);
else if (x>wz) replace(s*2+1,wz+1,r,x,y,z);
else
{
replace(s*2,l,wz,x,wz,z);
replace(s*2+1,wz+1,r,wz+1,y,z);
}
up(s,l,r);
}
void find(LL s,LL l,LL r,LL x,LL y,LL op)
{
LL wz=(l+r)/2;
down(s*2,l,wz);down(s*2+1,wz+1,r);
if (l==x && r==y)
{
if (op==4) ans+=sum[s];
if (op==5) ans=min(ans,mn[s]);
if (op==6) ans=max(ans,mx[s]);
return;
}
if (y<=wz) find(s*2,l,wz,x,y,op);
else if (x>wz) find(s*2+1,wz+1,r,x,y,op);
else
{
find(s*2,l,wz,x,wz,op);
find(s*2+1,wz+1,r,wz+1,y,op);
}
up(s,l,r);
}
LL findlg(LL s,LL l,LL r,LL x,LL y)
{
LL wz=(l+r)/2;
down(s*2,l,wz);down(s*2+1,wz+1,r);
if (l==x && r==y) return lg[s];
if (y<=wz) return findlg(s*2,l,wz,x,y);
else if (x>wz) return findlg(s*2+1,wz+1,r,x,y);
else
{
LL ans;
ans=max(findlg(s*2,l,wz,x,wz),findlg(s*2+1,wz+1,r,wz+1,y));
if (rn[s*2]==ln[s*2+1])
ans=max(ans,min(rs[s*2],wz-x+1)+min(ls[s*2+1],y-wz));
return ans;
}
up(s,l,r);
}
int main()
{
n=read();q=read();
fo(i,1,n) a[i]=read();
build(1,1,n);
fo(i,1,q)
{
op=read();
if (op<4)
{
L=read();R=read();K=read();
if (op==1) inspls(1,1,n,L,R,K);
if (op==2) inspls(1,1,n,L,R,-K);
if (op==3) replace(1,1,n,L,R,K);
} else
{
L=read();R=read();
if (op==4) ans=0;
if (op==5) ans=2147483647;
if (op==6) ans=-2147483647;
if (op==7) ans=findlg(1,1,n,L,R);
else find(1,1,n,L,R,op);
printf("%lld\n",ans);
}
}
}