BZOJ 1500 维修数列

Problem

BZOJ请戳我

Solution

又是一道splay的恶心题，大模拟(真相了)，细节多。

首先这道题卡空间，而且节点多，你需要开个队列回收被删除的节点，写个get函数。
然后按照套路，你需要两个哨兵节点。
对于插入操作，由于插入节点较多，用insert操作会导致复杂度变成$O(nlogn)$，不如直接直接调用build函数，旋出需要插入的相邻节点，然后对pos+1的左儿子建树，复杂度$O(n)$。
对于删除操作，旋出相应区间，然后再dfs回收。
对于修改操作和翻转操作，打标记。在下放标记时，先下放修改标记，因为在全修改后，翻转就没有意义了。
对于求和操作，维护一个sum数组表示答案。

对于求最大子区间，考虑维护两个辅助数组，lx[x]表示从左端开始的最大连续子区间，rx[x]表示从右段开始的最大连续子区间。那么我们就有$mx[x]=max(lx[l],rx[r],key[x]+rx[l]+lx[r])$。对于lx和rx的转移我们有
$lx[x]=max(lx[l],sum[l]+key[x]+lx[r])$，$rx[x]=max(rx[r],sum[r]+key[x]+rx[l])$。
而在翻转操作时，需要swap(lx[x],rx[x])。
然而mx操作要求至少选取一个数字。因此当修改操作修改成负值时，特殊处理一下，lx[x]=rx[x]=0,mx[x]=key[x]。而且为了使得哨兵节点和空节点不影响答案，我们设它的权值为-INF。

具体的实现可以看代码，感觉各种长啊。本来写完的时候只有3.5k，结果调试了很久，结果就变成4.6k了……
调不出来，建议写拍子……或者上Codevs？

Code

#include <cstdio>
#include <queue>
using namespace std;
const int maxn=500100,INF=0x3f3f3f3f;
queue<int> q;
int n,m,rt,size,a[maxn],key[maxn],f[maxn],sz[maxn],ch[maxn][2];
int lazy[maxn],rev[maxn],sum[maxn],lx[maxn],rx[maxn],mx[maxn];
char op[10];
template <typename Tp> inline void read(Tp &x)
{
    x=0;int f=0;char ch=getchar();
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') f=1,ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    if(f) x=-x;
}
inline int max(int x,int y){return x>y?x:y;}
inline int get()
{
    if(size<500010) return ++size;
    int x=q.front();
    q.pop();
    return x;
}
inline void clear(int x){key[x]=sz[x]=ch[x][0]=ch[x][1]=lazy[x]=rev[x]=sum[x]=lx[x]=rx[x]=mx[x]=0;}
inline void pushup(int x)
{
    int l=ch[x][0],r=ch[x][1];
    sz[x]=sz[l]+sz[r]+1;
    sum[x]=sum[l]+sum[r]+key[x];
    mx[x]=max(mx[l],max(mx[r],key[x]+rx[l]+lx[r]));
    lx[x]=max(lx[l],sum[l]+key[x]+lx[r]);
    rx[x]=max(rx[r],sum[r]+key[x]+rx[l]);
}
inline void pushdown(int x)
{
    int l=ch[x][0],r=ch[x][1];
    if(lazy[x])
    {
        rev[x]=lazy[x]=0;
        if(l) lazy[l]=1,key[l]=key[x],sum[l]=key[x]*sz[l];
        if(r) lazy[r]=1,key[r]=key[x],sum[r]=key[x]*sz[r];
        if(key[x]>=0)
        {
            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]=0,mx[l]=key[x];
            if(r) lx[r]=rx[r]=0,mx[r]=key[x];
        }
    }
    if(rev[x])
    {
        rev[x]=0;rev[l]^=1;rev[r]^=1;
        swap(lx[l],rx[l]);swap(lx[r],rx[r]);
        swap(ch[l][0],ch[l][1]);swap(ch[r][0],ch[r][1]);
    }
}
void rotate(int x,int &k)
{
    int fa=f[x],ff=f[fa],l,r;
    if(ch[fa][0]==x)l=0;else l=1;r=l^1;
    if(fa==k)k=x;
    else {if(ch[ff][0]==fa)ch[ff][0]=x;else ch[ff][1]=x;}
    f[x]=ff;f[fa]=x;f[ch[x][r]]=fa;
    ch[fa][l]=ch[x][r];ch[x][r]=fa;
    pushup(fa);pushup(x);
}
void splay(int x,int &k)
{
    while(x!=k)
    {
        int fa=f[x],ff=f[fa];
        if(fa!=k)
        {
            if(ch[ff][0]==fa^ch[fa][0]==x) rotate(x,k);
            else rotate(fa,k);
        }
        rotate(x,k);
    }
}
int build(int l,int r)
{
    if(l>r) return 0;
    if(l==r)
    {
        int now=get();
        lazy[now]=rev[now]=0;lx[now]=rx[now]=max(a[l],0);
        key[now]=mx[now]=sum[now]=a[l];sz[now]=1;
        return now;
    }
    int m=(l+r)>>1,now=get();
    f[ch[now][0]=build(l,m-1)]=now;
    f[ch[now][1]=build(m+1,r)]=now;
    key[now]=a[m];pushup(now);
    return now;
}
void dfs(int x)
{
    if(!x) return ;
    q.push(x);
    if(ch[x][0]) dfs(ch[x][0]);
    if(ch[x][1]) dfs(ch[x][1]);
    clear(x);
}
int find_th(int x)
{
    int now=rt;x++;
    while(1)
    {
        if(lazy[now]||rev[now]) pushdown(now);
        if(sz[ch[now][0]]+1==x) return now;
        if(sz[ch[now][0]]+1<x) x-=sz[ch[now][0]]+1,now=ch[now][1];
        else now=ch[now][0];
    }
}
void insert()
{
    int pos,tot;
    read(pos);read(tot);
    for(int i=1;i<=tot;i++) read(a[i]);
    splay(find_th(pos),rt);splay(find_th(pos+1),ch[rt][1]);
    f[ch[ch[rt][1]][0]=build(1,tot)]=ch[rt][1];
    pushup(ch[rt][1]);pushup(rt);
}
void recycle()
{
    int pos,tot;
    read(pos);read(tot);
    splay(find_th(pos-1),rt);splay(find_th(pos+tot),ch[rt][1]);
    dfs(ch[ch[rt][1]][0]);ch[ch[rt][1]][0]=0;
    pushup(ch[rt][1]);pushup(rt);
}
void update()
{
    int pos,tot,c,x;
    read(pos);read(tot);read(c);
    splay(find_th(pos-1),rt);splay(find_th(pos+tot),ch[rt][1]);
    key[x=ch[ch[rt][1]][0]]=c;lazy[ch[ch[rt][1]][0]]=1;sum[x]=sz[x]*c;
    if(c>=0) lx[x]=rx[x]=mx[x]=sum[x];
    else lx[x]=rx[x]=0,mx[x]=c;
    pushup(ch[rt][1]);pushup(rt);
}
void reverse()
{
    int pos,tot,x;
    read(pos);read(tot);
    splay(find_th(pos-1),rt);splay(find_th(pos+tot),ch[rt][1]);
    rev[x=ch[ch[rt][1]][0]]^=1;swap(ch[x][0],ch[x][1]);
    swap(lx[x],rx[x]);pushup(ch[rt][1]);pushup(rt);
}
void query()
{
    int pos,tot;
    read(pos);read(tot);
    splay(find_th(pos-1),rt);splay(find_th(pos+tot),ch[rt][1]);
    printf("%d\n",sum[ch[ch[rt][1]][0]]);
}
int main()
{
    #ifndef ONLINE_JUDGE
    freopen("in","r",stdin);
    #endif
    read(n);read(m);
    for(int i=1;i<=n;i++) read(a[i]);
    mx[0]=a[0]=a[n+1]=-INF;
    rt=build(0,n+1);
    while(m--)
    {
        scanf("%s",op);
        switch(op[2])
        {
            case 'S':insert();break;
            case 'L':recycle();break;
            case 'K':update();break;
            case 'V':reverse();break;
            case 'T':query();break;
            case 'X':printf("%d\n",mx[rt]);break;
        }
    }
    return 0;
}