• NOI2005 维护数列(splay)


    学了半天平衡树,选择了一道题来写一写,发现题目是裸的splay模板,但是还是写不好,这个的精髓之处在于在数列的某一个位置加入一个数列,类似于treap里面的merge,然后还学到了题解里面的的回收空间做法,利用时间换取空间也是一种非常不错的做法,由于题目中有区间翻转的操作,所以只能用splay来写,然后对于一些区间最大值和区间最大子序列的询问我们可以用类似于线段树的合并操作来维护,总的来说这是一道对于树形数据结构非常好的运用的题目。——by VANE

    #include<bits/stdc++.h>
    using namespace std;
    const int inf=0x3f3f3f3f;
    const int N=1e6+7;
    int n,m,rt,cnt;
    int a[N],id[N],fa[N],c[N][2];
    int sum[N],sz[N],v[N],mx[N],lx[N],rx[N];
    bool tag[N],rev[N];
    queue<int> q;
    inline int read()
    {
        int x=0,f=1;char ch=getchar();
        while(ch<'0'||ch>'9') {if(ch=='-') f=-f; ch=getchar();}
        while(ch>='0'&&ch<='9') {x=x*10+ch-'0';ch=getchar();}
        return x*f;
    }
    inline void pushup(int x)
    {
        int l=c[x][0],r=c[x][1];
        sum[x]=sum[l]+sum[r]+v[x];
        sz[x]=sz[l]+sz[r]+1;
        mx[x]=max(mx[l],max(mx[r],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]);
    }
    inline void pushdown(int x)
    {
        int l=c[x][0],r=c[x][1];
        if(tag[x])
        {
            rev[x]=tag[x]=0;
            if(l) tag[l]=1,v[l]=v[x],sum[l]=v[x]*sz[l];
            if(r) tag[r]=1,v[r]=v[x],sum[r]=v[x]*sz[r];
            if(v[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]=v[x];
                if(r) lx[r]=rx[r]=0,mx[r]=v[x];
            }
        }
        if(rev[x])
        {
            rev[x]=0;rev[l]^=1;rev[r]^=1;
            swap(lx[l],rx[l]);swap(lx[r],rx[r]);
            swap(c[l][0],c[l][1]);swap(c[r][0],c[r][1]);
        }
    }
    inline void rotate(int x,int &k)
    {
        int y=fa[x],z=fa[y],l=(c[y][1]==x),r=l^1;
        if(y==k) k=x;else c[z][c[z][1]==y]=x;
        fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
        c[y][l]=c[x][r];c[x][r]=y;
        pushup(y);pushup(x);
    }
    inline void splay(int x,int &k)
    {
        while(x!=k)
        {
            int y=fa[x],z=fa[y];
            if(y!=k)
            {
                if(c[z][0]==y ^ c[y][0]==x) rotate(x,k);
                else rotate(y,k);
            }
            rotate(x,k);
        }
    }
    inline int find(int x,int rk)
    {
        pushdown(x);
        int l=c[x][0],r=c[x][1];
        if(sz[l]+1==rk) return x;
        if(sz[l]>=rk) return find(l,rk);
        return find(r,rk-sz[l]-1);
    }
    inline void recycle(int x)
    {
        int &l=c[x][0],&r=c[x][1];
        if(l) recycle(l);
        if(r) recycle(r);
        q.push(x);
        fa[x]=l=r=tag[x]=rev[x]=0;
    }
    inline int split(int k,int tot)
    {
        int x=find(rt,k),y=find(rt,k+tot+1);
        splay(x,rt);splay(y,c[x][1]);
        return c[y][0];
    }
    inline void query(int k,int tot)
    {
        int x=split(k,tot);
        printf("%d
    ",sum[x]);
    }
    inline void modify(int k,int tot,int val)
    {
        int x=split(k,tot),y=fa[x];
        v[x]=val;tag[x]=1;sum[x]=sz[x]*val;
        if(val>=0) lx[x]=rx[x]=mx[x]=sum[x];
        else lx[x]=rx[x]=0,mx[x]=val;
        pushup(y);pushup(fa[y]);
    }
    inline void rever (int k,int tot)
    {
        int x=split(k,tot),y=fa[x];
        if(!tag[x])
        {
            rev[x]^=1;
            swap(c[x][0],c[x][1]);
            swap(lx[x],rx[x]);
            pushup(y);pushup(fa[y]);
        }
    }
    inline void erase(int k,int tot)
    {
        int x=split(k,tot),y=fa[x];
        recycle(x);c[y][0]=0;
        pushup(y);pushup(fa[y]);
    }
    inline void build(int l,int r,int f)
    {
        int mid=l+r>>1,now=id[mid],pre=id[f];
        if(l==r)
        {
            mx[now]=sum[now]=a[l];
            lx[now]=rx[now]=max(a[l],0);
            sz[now]=1;
        }
        if(l<mid) build(l,mid-1,mid);
        if(mid<r) build(mid+1,r,mid);
        v[now]=a[mid];fa[now]=pre;
        pushup(now);
        c[pre][mid>=f]=now;
    }
    inline void insert(int k,int tot)
    {
        for(int i=1;i<=tot;++i) a[i]=read();
        for(int i=1;i<=tot;++i)
        {
            if(!q.empty()) id[i]=q.front(),q.pop();
            else id[i]=++cnt;
        }
        build(1,tot,0);
        int z=id[1+tot>>1];
        int x=find(rt,k+1),y=find(rt,k+2);
        splay(x,rt);splay(y,c[x][1]);
        fa[z]=y;c[y][0]=z;
        pushup(y);pushup(x);
    }
    int main()
    {
        n=read();m=read();
        mx[0]=a[1]=a[n+2]=-inf;
        for(int i=1;i<=n;++i) a[i+1]=read();
        for(int i=1;i<=n+2;++i) id[i]=i;
        build(1,n+2,0);
        rt=n+3>>1;cnt=n+2;
        int k,tot,val;char ch[10];
        while(m--)
        {
            scanf("%s",ch);
            if(ch[0]!='M'||ch[2]!='X') k=read(),tot=read();
            if(ch[0]=='I') insert(k,tot);
            if(ch[0]=='D') erase(k,tot);
            if(ch[0]=='M')
            {
                if(ch[2]=='X') printf("%d
    ",mx[rt]);
                else val=read(),modify(k,tot,val);
            }
            if(ch[0]=='R') rever(k,tot);
            if(ch[0]=='G') query(k,tot);
        }
    }
  • 相关阅读:
    hutool 解析 Excel
    上传文件
    Cannot construct instance of `com.**` (although at least one Creator exists)
    Java8之Optional
    java8之Stream
    java8之Lambda
    springboot+mybatis事务管理
    queryWrapper in like
    Java 组装 Tree
    JWT
  • 原文地址:https://www.cnblogs.com/nbwzyzngyl/p/8043687.html
Copyright © 2020-2023  润新知