[BZOJ1500][NOI2005]维修数列

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

题解

似乎这是一道裸的splay，但是空间卡的很紧，要优化。

调了4个多小时。。。

代码

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>

const int inf = +;
const int N = +;

inline int read(int &in) {
    in=;int f=;char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar()) if(ch=='-') f=-;
    for(;ch>='0'&&ch<='9';ch=getchar()) in=in*+ch-'0';
    return in*f;
}

using std::swap;
using std::max;

int n,m,rt,cnt;
int a[N],id[N],fa[N],c[N][];
int sum[N],size[N],v[N],mx[N],lx[N],rx[N];
bool tag[N],rev[N];
std::queue<int> q;

namespace Splay {
    void update(int x) {
        int l=c[x][],r=c[x][];
        sum[x]=sum[l]+sum[r]+v[x];
        size[x]=size[l]+size[r]+;
        mx[x]=max(mx[l],mx[r]);
        mx[x]=max(mx[x],rx[l]+v[x]+lx[r]);
        lx[x]=max(lx[l],sum[l]+v[x]+lx[r]);
        rx[x]=max(rx[r],sum[r]+v[x]+rx[l]);
    }
    void pushdown(int x)
    {
        int l=c[x][],r=c[x][];
        if(tag[x])
        {
            rev[x]=tag[x]=;
            if(l)tag[l]=,v[l]=v[x],sum[l]=v[x]*size[l];
            if(r)tag[r]=,v[r]=v[x],sum[r]=v[x]*size[r];
            if(v[x]>=)
            {
                if(l)lx[l]=rx[l]=mx[l]=sum[l];
                if(r)lx[r]=rx[r]=mx[r]=sum[r];
            }
            else 
            {
                if(l)lx[l]=rx[l]=,mx[l]=v[x];
                if(r)lx[r]=rx[r]=,mx[r]=v[x];
            }
        }
        if(rev[x])
        {
            rev[x]^=;rev[l]^=;rev[r]^=;
            swap(lx[l],rx[l]);swap(lx[r],rx[r]);
            swap(c[l][],c[l][]);swap(c[r][],c[r][]);
        }
    }
    void rotate(int x,int &k)
    {
        int y=fa[x],z=fa[y],l,r;
        l=(c[y][]==x);r=l^;
        if(y==k)k=x;
        else c[z][c[z][]==y]=x;
        fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
        c[y][l]=c[x][r];c[x][r]=y;
        update(y);update(x);
    }

    void splay(int x,int &k) {
        while(x!=k) {
            int y=fa[x],z=fa[y];
            if(y!=k) {
                if(c[y][]==x^c[z][]==y)rotate(x,k);
                else rotate(y,k);
            }
            rotate(x,k);
        }
    }

    int find(int x,int rk) {
        pushdown(x);
        int l=c[x][],r=c[x][];
        if(size[l]+==rk)return x;
        if(size[l]>=rk)return find(l,rk);
        return find(r,rk-size[l]-);
    }

    void rec(int x) {
        if(!x)return;
        int l=c[x][],r=c[x][];
        rec(l);rec(r);q.push(x);
        fa[x]=c[x][]=c[x][]=;
        tag[x]=rev[x]=;
    }

    int split(int k,int tot) {
        int x=find(rt,k),y=find(rt,k+tot+);
        splay(x,rt);splay(y,c[x][]);
        return c[y][];
    }

    void query(int k,int tot) {
        int x=split(k,tot);
        printf("%d\n",sum[x]);
    }

    void modify(int k,int tot,int val){
        int x=split(k,tot),y=fa[x];
        v[x]=val;tag[x]=;sum[x]=size[x]*val;
        if(val>=)lx[x]=rx[x]=mx[x]=sum[x];
        else lx[x]=rx[x]=,mx[x]=val;
        update(y);update(fa[y]);
    }

    void rever(int k,int tot) {
        int x=split(k,tot),y=fa[x];
        if(!tag[x]) {
            rev[x]^=;
            swap(c[x][],c[x][]);
            swap(lx[x],rx[x]);
            update(y);update(fa[y]);
        }
    }

    void erase(int k,int tot) {
        int x=split(k,tot),y=fa[x];
        rec(x);c[y][]=;
        update(y);update(fa[y]);
    }

    void build(int l,int r,int f) {
        if(l>r)return;
        int mid=(l+r)>>,now=id[mid],last=id[f];
        if(l==r) {
            sum[now]=a[l];size[now]=;
            tag[now]=rev[now]=;
            if(a[l]>=)lx[now]=rx[now]=mx[now]=a[l];
            else lx[now]=rx[now]=,mx[now]=a[l];
        }
        else build(l,mid-,mid),build(mid+,r,mid);
        v[now]=a[mid];fa[now]=last;update(now);
        c[last][mid>=f]=now;
    }

    void insert(int k,int tot) {
        for(int i=;i<=tot;i++)read(a[i]);
        for(int i=;i<=tot;i++)
            if(!q.empty())id[i]=q.front(),q.pop();
            else id[i]=++cnt;
        build(,tot,);int z=id[(+tot)>>];
        int x=find(rt,k+),y=find(rt,k+);
        splay(x,rt);splay(y,c[x][]);
        fa[z]=y;c[y][]=z;
        update(y);update(x);
    }

}


int main() {
    read(n);read(m);
    mx[]=a[]=a[n+]=-inf;
    for(int i=;i<=n;i++)read(a[i+]);
    for(int i=;i<=n+;i++)id[i]=i;
    Splay::build(,n+,);
    rt=(n+)>>;cnt=n+;
    int k,tot,val;
    char ch[];
    while(m--)
    {
        scanf("%s",ch);
        if(ch[]!='M'||ch[]!='X')read(k),read(tot);
        if(ch[]=='I')Splay::insert(k,tot);
        if(ch[]=='D')Splay::erase(k,tot);
        if(ch[]=='M') {
            if(ch[]=='X')printf("%d\n",mx[rt]);
            else read(val),Splay::modify(k,tot,val);
        }
        if(ch[]=='R')Splay::rever(k,tot);
        if(ch[]=='G')Splay::query(k,tot);
    }
    return ;
}