Bzoj1500 [NOI2005]维修数列

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 12469  Solved: 3955

Description

Input

输入的第1 行包含两个数N 和M(M ≤20 000)，N 表示初始时数列中数的个数，M表示要进行的操作数目。
第2行包含N个数字，描述初始时的数列。
以下M行，每行一条命令，格式参见问题描述中的表格。
任何时刻数列中最多含有500 000个数，数列中任何一个数字均在[-1 000, 1 000]内。
插入的数字总数不超过4 000 000个，输入文件大小不超过20MBytes。

Output

对于输入数据中的GET-SUM和MAX-SUM操作，向输出文件依次打印结果，每个答案（数字）占一行。

Sample Input

9 8
2 -6 3 5 1 -5 -3 6 3
GET-SUM 5 4
MAX-SUM
INSERT 8 3 -5 7 2
DELETE 12 1
MAKE-SAME 3 3 2
REVERSE 3 6
GET-SUM 5 4
MAX-SUM

Sample Output

-1
10
1
10

HINT

Source

Splay模板题

真·一道题写一天

每个指令都是基本操作，放在一起就迷之难写……

由于空间限制，必须要重复使用del掉的旧结点

几个WA点：

　　申请新结点时，没有初始化ch[]

　　序列st，ed编号记录错误，导致MAX-SUM操作转错

　　pushdown只处理当前层，并把标记传给下层←会WA，非要把下一层也处理掉才行。

　　pushdown没有判“当前结点为空时直接跳过”

/*by SilverN*/
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<vector>
using namespace std;
const int INF=1e9;
const int mxn=1000100;
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int stk[mxn],top;//回收结点 
struct node{
    int ch[2];
    int fa,w,mk;
    int rx,lx,max;
    int size,sum;
    bool rev;
}t[mxn];
int a[mxn];//存储数值 
int cnt=0;
int root=0;
int n,m,st,ed;
void pd_cov(int rt){//区间修改
    if(!rt)return;
    t[rt].w=t[rt].mk;
    t[rt].sum=t[rt].size*t[rt].w;
    t[rt].lx=t[rt].rx=t[rt].max=max(t[rt].w,t[rt].sum);
    return;
}
void pd_rev(int rt){//区间翻转 
    if(!rt)return;
    t[rt].rev^=1;
    swap(t[rt].ch[0],t[rt].ch[1]);
    swap(t[rt].lx,t[rt].rx);//翻转左右连续和 
    return;
}
void pushdown(int x){
    if(!x)return;
    int lc=t[x].ch[0],rc=t[x].ch[1];
    if(t[x].mk!=-INF){
        t[lc].mk=t[rc].mk=t[x].mk;
        pd_cov(lc);pd_cov(rc);
        t[x].mk=-INF;
    }
    if(t[x].rev){
        pd_rev(lc);pd_rev(rc);
        t[x].rev=0;
    }
    return;
}
inline void pushup(int x){
    if(!x)return;
    int lc=t[x].ch[0],rc=t[x].ch[1];
    t[x].lx=max(t[lc].lx,t[lc].sum+t[x].w+max(0,t[rc].lx));
    t[x].rx=max(t[rc].rx,t[rc].sum+t[x].w+max(0,t[lc].rx));
    t[x].sum=t[lc].sum+t[x].w+t[rc].sum;
    t[x].max=max(max(t[lc].max,t[rc].max),max(t[rc].lx,0)+t[x].w+max(t[lc].rx,0));
    t[x].size=t[lc].size+t[rc].size+1;
    return;
}
void rotate(int x,int &k){
    int y=t[x].fa,z=t[y].fa,lc,rc;
    if(t[y].ch[0]==x)lc=0;else lc=1;
    rc=lc^1;
    if(y==k)k=x;
    else t[z].ch[t[z].ch[1]==y]=x;
    t[x].fa=z;t[y].fa=x;  t[t[x].ch[rc]].fa=y;
    t[y].ch[lc]=t[x].ch[rc];t[x].ch[rc]=y;
    pushup(y);
}
void relax(int x,int k){if(x!=k)relax(t[x].fa,k);pushdown(x);return;}
void Splay(int x,int &k){
    relax(x,k);
    while(x!=k){
        int y=t[x].fa,z=t[y].fa;
        if(y!=k){
            if((t[z].ch[0]==y)^(t[y].ch[0]==x))rotate(x,k);
            else rotate(y,k);
        }
        rotate(x,k);
    }
    pushup(x);
    return;
}
inline void split(int x,int y){
    Splay(x,root);
    Splay(y,t[x].ch[1]);
    return;
}
inline int find(int x,int num){//查询序列中第num个数 
    pushdown(x);
    int lc=t[x].ch[0],rc=t[x].ch[1];
    if(t[lc].size>=num)return find(lc,num);
    if(t[lc].size+1==num)return x;
    return find(rc,num-t[lc].size-1);
}
inline int newnode(int x){//申请新结点 
    int rt;
    if(top)rt=stk[top--];else rt=++cnt;
    t[rt].w=a[x];
    t[rt].lx=t[rt].rx=t[rt].max=-INF;
    t[rt].rev=0;t[rt].mk=-INF;
    return rt;
}
int Build(int l,int r,int fa){
    if(l>r)return 0;
    int mid=(l+r)>>1,lc,rc;
    int x=newnode(mid);
    if(l<mid) t[x].ch[0]=Build(l,mid-1,x);
        else t[x].ch[0]=0;
    if(r>mid) t[x].ch[1]=Build(mid+1,r,x);
        else t[x].ch[1]=0;
    t[x].fa=fa;
    pushup(x);
    return x;
}
void del(int &x){
    if(!x)return;
    stk[++top]=x;
    t[x].fa=0;
    del(t[x].ch[0]),del(t[x].ch[1]);
    x=0;
    return;
}
void insert(int pos,int tot){
    for(int i=1;i<=tot;i++)a[i]=read();
    int tmp=Build(1,tot,0);
    int x=find(root,pos+1),y=find(root,pos+2);
    split(x,y);
    t[y].ch[0]=tmp;t[tmp].fa=y;
    pushup(y);pushup(x);
    return;
}
void msame(int pos,int tot,int c){
    int x=find(root,pos),y=find(root,pos+tot+1);
    split(x,y);
    t[t[y].ch[0]].mk=c;
    pd_cov(t[y].ch[0]);
    pushup(y);pushup(x);
    return;
}
void reverse(int pos,int tot){
    int x=find(root,pos),y=find(root,pos+tot+1);
    split(x,y);
    pd_rev(t[y].ch[0]);
    pushup(y);pushup(x);
    return;
}
void solvedel(int pos,int tot){
    int x=find(root,pos),y=find(root,pos+tot+1);
    split(x,y);
    del(t[y].ch[0]);
    pushup(y);pushup(x);
    return;
}
inline int getmax(int pos,int tot){
    split(st,ed);
    return t[t[ed].ch[0]].max;
}
inline int getsum(int pos,int tot){
    int x=find(root,pos),y=find(root,pos+tot+1);
    split(x,y);
    return t[t[y].ch[0]].sum;
}
void DEBUG(){
    printf("root:%d
",root);
    for(int i=1;i<=cnt;i++){
        printf("id:%d  fa:%d lc:%d rc:%d w:%d
",i,t[i].fa,t[i].ch[0],t[i].ch[1],t[i].w);
        printf("  lx:%d rx:%d mx:%d
",t[i].lx,t[i].rx,t[i].max);
        printf("  size:%d
",t[i].size);
    }
    return;
}
int main(){
    int i,j,x,y,c;
    n=read();m=read();
    for(i=1;i<=n;i++)a[i+1]=read();
    st=1;ed=n+2;
    a[st]=a[ed]=0;
    t[0].lx=t[0].rx=t[0].max=-INF;
    root=Build(1,n+2,0);
    st=find(root,st);ed=find(root,ed);//记录头尾结点的编号 
    char op[15];
    while(m--){
        scanf("%s",op);
        switch(op[0]){
            case 'I':x=read();y=read();insert(x,y);break;
            case 'D':{x=read();y=read();solvedel(x,y);break;}
            case 'M':{
                if(op[2]=='K'){x=read();y=read();c=read();msame(x,y,c);}
                else{printf("%d
",getmax(st,ed));}
                break;
            }
            case 'R':{x=read();y=read();reverse(x,y);break;}
            case 'G':{x=read();y=read();printf("%d
",getsum(x,y));break;}
            case 'B':DEBUG();break;
        }
    }
    return 0;
}