并不对劲的bzoj1500: [NOI2005]维修数列
这题没什么好说的……小清新数据结构题……并不对劲的人太菜了,之前照着标程逐行比对才过了这道题,前几天刚刚把这题一遍写对……
其实这题应该口胡很容易。操作1,2,3,4,5就是普通的splay支持的那些操作。操作6要转一个小弯,对于每个点记lmax,rmax,xmax,分别表示贴左边界的最大连续子段和(可能为空)、贴右边界的最大连续子段和(可能为空)、最大子段和(不能为空)。现在考虑将lson[u]和rson[u]子树对应的序列的lmax,rmax,xmax经过一些神奇的操作合并,得到u的子树对应的这些值。lmax[u]:跨越u时是左边整段+key[u]+右边lmax;不跨越时是左边的lmax。rmax[u]:同上。xmax[u]:跨越u时是左边贴右边界+右边贴左边界+key[u];不跨越时左边、右边的xmax二选一。
听上去这题已经口胡完了,但是调它的路还长这呢。要注意几点小小的问题:
1.一开始建树时不要一个个插入,可以每次取中点作为根,再递归处理左右两边,这样会形成更加平衡的树。会发现操作1是插入一个序列,可以类似地先将这个序列建成平衡的树。(不这么做可能TLE)
2.发现数据略坑,无脑新建点会需要4*10^6个点。但是只会有5*10^5个点被用到,所以对于那些删掉的点,可以将它们的编号存在【随便什么简单数据结构】里,用的时候再拿出来。(不这么做会MLE或RE)
3.在做2时,假设将编号x存到了【随便什么简单数据结构】里,发现它被删除前的父亲、儿子、lmax、rmax等一系列信息还在,注意要记得初始化。(不这么做当然会WA)
4.在pushdown时,发现当有一段被改成同一个数后就不用管翻转标记了。(不这么做会常数极大)
5.pushup/down时懒得判某个儿子为0怎么办?在pushup/down前后强行将编号为0的点的各种值全部设置为初始值!(不这么做会写烦)
6.在做5时,xmax[0]=-inf,其余都为0。(不这么做会WA)
7.对于某点打上修改标记时,pushdown前它的左、右儿子都没有改变,这时pushup就会得到修改前的答案。不能盲目地pushup。(不这么做会WA)
8.pushdown翻转标记时,记得交换lmax,rmax。(不这么做会WA)
#include<algorithm> #include<cmath> #include<cstdio> #include<cstdlib> #include<cstring> #include<iomanip> #include<iostream> #include<map> #include<queue> #include<stack> #include<vector> #define rep(i,x,y) for(register LL i=(x);i<=(y);i++) #define dwn(i,x,y) for(register LL i=(x);i>=(y);i--) #define re register #define Int re LL #define maxn 500010 #define mi (l+r>>1) #define ls son[u][0] #define rs son[u][1] #define s0 fa[0]=son[0][0]=son[0][1]=rev[0]=sum[0]=siz[0]=key[0]=lm[0]=rm[0]=yes[0]=0,xm[0]=-inf[0] #define rks son[rk][0] #define LL long long using namespace std; inline LL read() { LL x=0,f=1; char ch=getchar(); while(isdigit(ch)==0 && ch!='-')ch=getchar(); if(ch=='-')f=-1,ch=getchar(); while(isdigit(ch))x=x*10+ch-'0',ch=getchar(); return x*f; } inline void write(LL x) { LL f=0;char ch[20]; if(!x){puts("0");return;} if(x<0){putchar('-');x=-x;} while(x)ch[++f]=x%10+'0',x/=10; while(f)putchar(ch[f--]); putchar('\n'); } LL n,m; LL fa[maxn],son[maxn][2],rev[maxn],sum[maxn],xm[maxn],siz[maxn],inf[5],x,y,cnt; LL st[maxn],top,tmp[maxn],key[maxn],rt,lm[maxn],rm[maxn],yes[maxn]; void mark(LL u,LL k){if(u)yes[u]=1ll,key[u]=k,sum[u]=k*siz[u],xm[u]=max(k*siz[u],k),lm[u]=rm[u]=max(k*siz[u],0ll);} void rever(LL u){if(u)rev[u]^=1ll,swap(ls,rs),swap(lm[u],rm[u]); } void pu(LL u) { s0; siz[u]=siz[ls]+siz[rs]+1ll; sum[u]=sum[ls]+sum[rs]+key[u]; lm[u]=max(lm[ls],sum[ls]+key[u]+lm[rs]); rm[u]=max(rm[rs],sum[rs]+key[u]+rm[ls]); xm[u]=max(max(xm[ls],xm[rs]),lm[rs]+rm[ls]+key[u]); s0; } void pd(LL u) { s0; if(yes[u])mark(ls,key[u]),mark(rs,key[u]),rev[u]=yes[u]=0; else if(rev[u])rever(ls),rever(rs),rev[u]=0; s0; } LL res(LL k,LL f) { LL u; if(top>0ll)u=st[top--]; else u=++cnt; fa[u]=f,ls=rs=rev[u]=yes[u]=0;siz[u]=1; sum[u]=xm[u]=k;key[u]=k; lm[u]=rm[u]=max(0ll,k); return u; } LL build(LL l,LL r) { if(l>r)return 0; LL u=res(tmp[mi],0ll); ls=build(l,mi-1ll),rs=build(mi+1ll,r); fa[ls]=fa[rs]=u,pu(u); return u; } LL kth(LL k) { LL u=rt; while(u) { pd(u); LL rk=siz[ls]+1ll; if(rk==k)break; if(rk>k)u=ls; else k-=rk,u=rs; } return u; } void getst(LL u) { if(!u)return ; pd(u),getst(ls),st[++top]=u,getst(rs); } inline LL getso(LL u) {return son[fa[u]][0]!=u;} void rot(LL u) { LL fu=fa[u],ffu=fa[fu],l=getso(u),fl=getso(fu),r=l^1ll,rson=son[u][r]; fa[fu]=u,fa[u]=ffu,fa[rson]=fu,son[ffu][fl]=u,son[fu][l]=rson,son[u][r]=fu; pu(fu),pu(u); } void splay(LL u,LL k) { while(fa[u]!=k){LL fu=fa[u];if(fa[fu]!=k)rot(getso(u)^getso(fu)?u:fu);rot(u);} if(!k)rt=u; } void ins(LL pos,LL num) { LL lk=kth(pos),rk=kth(pos+1ll); splay(lk,0),splay(rk,lk); rep(i,1,num)tmp[i]=read(); rks=build(1,num),fa[rks]=rk,pu(rk),pu(lk); } void del(LL pos,LL num) { LL lk=kth(pos-1ll),rk=kth(pos+num); splay(lk,0),splay(rk,lk),getst(rks),rks=0,pu(rk),pu(lk); } void mksa(LL pos,LL num) { LL c=read(),lk=kth(pos-1ll),rk=kth(pos+num); splay(lk,0),splay(rk,lk),mark(rks,c),pd(rks),pu(rk),pu(lk); } void reverse(LL pos,LL num) { LL lk=kth(pos-1ll),rk=kth(pos+num); splay(lk,0),splay(rk,lk),rever(rks),pd(rks),pu(rk),pu(lk); } LL getsum(LL pos,LL num) { LL lk=kth(pos-1ll),rk=kth(pos+num); splay(lk,0),splay(rk,lk); return sum[rks]; } int main() { char s[20]; memset(inf,63,sizeof(inf)); n=read(),m=read(); rep(i,1,n)tmp[i]=read();tmp[0]=tmp[n+1]=-inf[0]; rt=build(0,n+1); while(m--) { scanf("%s",s); if(s[0]=='I')x=read()+1,y=read(),ins(x,y); if(s[0]=='D')x=read()+1,y=read(),del(x,y); if(s[0]=='M'&&s[2]=='K')x=read()+1,y=read(),mksa(x,y); if(s[0]=='R')x=read()+1,y=read(),reverse(x,y); if(s[0]=='G')x=read()+1,y=read(),write(getsum(x,y)); if(s[0]=='M'&&s[2]=='X'){if(siz[rt]>2)write(xm[rt]);else write(0);} } return 0; }
最后祝您身体健康,再见。