• 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);
        }
    }
  • 相关阅读:
    MSMQ实现自定义序列化存储
    流程部署的查询、删除、流程
    使用trello管理你的项目
    CentOS 6.4 编译安装 gcc 4.8.1
    jquery实现无限滚动瀑布流实现原理
    系统分析员备考之经济管理篇(二)
    架构、架构师和架构设计
    ASP.net Web API综合示例
    c中函数参数传递
    准备抽象NHibernate和EntityFramework
  • 原文地址:https://www.cnblogs.com/nbwzyzngyl/p/8043687.html
Copyright © 2020-2023  润新知