BZOJ1500: [NOI2005]维修数列（平衡树）

最新推荐文章于 2019-09-10 20:14:22 发布
最新推荐文章于 2019-09-10 20:14:22 发布
阅读量300 收藏
点赞数
分类专栏： 平衡树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_35649707/article/details/78128300
版权

传送门
题意：
平衡树的一系列操作。

题解：
Splay/无旋Treap
Splay:

#include<bits/stdc++.h>
using namespace std;
struct IO{
    streambuf *ib,*ob;
    inline void init(){
        ios::sync_with_stdio(false);
        cin.tie(NULL);cout.tie(NULL);
        ib=cin.rdbuf();ob=cout.rdbuf();
    }
    inline int read(){
        char ch=ib->sbumpc();int i=0,f=1;
        while(!isdigit(ch)){if(ch=='-')f=-1;ch=ib->sbumpc();}
        while(isdigit(ch)){i=(i<<1)+(i<<3)+ch-'0';ch=ib->sbumpc();}
        return i*f;
    }
    inline void W(long long x){
        static int buf[50];
        if(!x){ob->sputc('0');ob->sputc('\n');return;}
        if(x<0){ob->sputc('-');x=-x;}
        while(x)buf[++buf[0]]=x%10,x/=10;
        while(buf[0])ob->sputc(buf[buf[0]--]+'0');
        ob->sputc('\n');
    }
}io;

const int Maxn=5e5+50,INF=0x3f3f3f3f;

int n,m,nowtot;

struct node{
    node *lc,*rc,*fa;
    int val,sum,revtag,addtag,maxx,lmax,rmax,sze,isadd;
    node();
    inline void upt();
    inline void add(int v){
        val=v;sum=v*sze;
        int t=(v>0)?sze:1;
        maxx=lmax=rmax=t*v;
        addtag=v;isadd=1;
    }
    inline void rev(){
        swap(lc,rc);revtag^=1;
        swap(lmax,rmax);
    }
    inline void pushdown();
}Pool[Maxn],*pool=Pool,*null=Pool;
node::node():lc(null),rc(null),val(-INF),maxx(-INF),lmax(-INF),rmax(-INF),isadd(0){}
void node::upt() {
    sum=lc->sum+rc->sum+val;
    lmax=max(lc->lmax,max(lc->sum+val,lc->sum+val+rc->lmax));
    rmax=max(rc->rmax,max(rc->sum+val,rc->sum+val+lc->rmax));
    sze=lc->sze+rc->sze+1;
    maxx=max(val+max(lc->rmax,0)+max(rc->lmax,0),max(lc->maxx,rc->maxx));
}
void node::pushdown(){
    if(revtag){
        revtag=0;
        if(lc!=null)lc->rev();
        if(rc!=null)rc->rev();
    }
    if(isadd){
        if(lc!=null)lc->add(addtag);
        if(rc!=null)rc->add(addtag);
        isadd=0;
    }
}
struct RecyclePool{
    node *que[Maxn];
    int tail;
    inline node* newnode(int val){
        static node* tmp;
        tmp=(tail?que[tail--]:++pool);
        tmp->lc=tmp->rc=tmp->fa=null;tmp->isadd=0;
        tmp->revtag=tmp->addtag=0;tmp->sze=1;
        tmp->val=tmp->sum=tmp->maxx=tmp->lmax=tmp->rmax=val;
        return tmp;
    }
    inline void recycle(node *u){
        if(u==null)return;
        que[++tail]=u;
        recycle(u->lc);recycle(u->rc);
    }
}recyclepool;
struct Splay{
    node *rt;
    inline bool which(node *x){return x->fa->lc==x;}
    inline void rotate(node *x){
        node *y=x->fa,*z=y->fa;
        if(z!=null)(z->lc==y?z->lc:z->rc)=x;
        x->fa=z;y->fa=x;node *b;
        if(y->lc==x){b=x->rc;x->rc=y;y->lc=b;}
        else{b=x->lc;x->lc=y;y->rc=b;}
        if(b)b->fa=y;
        y->upt();x->upt();
    }
    inline void splay(node *x,node *tar){
        static node* que[Maxn];static int tail;
        que[tail=1]=x;
        for(node* y=x;y->fa!=null;y=y->fa)que[++tail]=y->fa;
        for(int i=tail;i>=1;i--)que[i]->pushdown();
        while(x->fa!=tar){
            node *y=x->fa,*z=y->fa;
            if(z!=tar)(which(x)^which(y))?(rotate(x)):(rotate(y));
            rotate(x);
        }
        if(tar==null)rt=x;
    }
    inline void insert(node *&now,node *fa,int pos,int val){
        now->pushdown();
        if(now==null){
            now=recyclepool.newnode(val);now->fa=fa;
            splay(now,null);
            return;
        }
        if(now->lc->sze+1>pos)insert(now->lc,now,pos,val);
        else insert(now->rc,now,pos-now->lc->sze-1,val);
    }
    inline node* findpos(node *now,int pos){
        now->pushdown();
        if(now->lc->sze+1==pos)return now;
        else if(now->lc->sze+1<pos)return findpos(now->rc,pos-now->lc->sze-1);
        else return findpos(now->lc,pos);
    }
    inline int querysum(int l,int r){
        splay(findpos(rt,l),null);splay(findpos(rt,r+2),rt);
        return rt->rc->lc->sum;
    }
    inline void del(int l,int r){
        splay(findpos(rt,l),null);splay(findpos(rt,r+2),rt);
        recyclepool.recycle(rt->rc->lc);rt->rc->lc->fa=null;rt->rc->lc=null;
        rt->rc->upt();rt->upt();
    }
    inline void rev(int l,int r){
        splay(findpos(rt,l),null);splay(findpos(rt,r+2),rt);
        rt->rc->lc->rev();rt->rc->upt();rt->upt();
    }
    inline void add(int l,int r,int v){
        splay(findpos(rt,l),null);splay(findpos(rt,r+2),rt);
        rt->rc->lc->add(v);rt->rc->upt();rt->upt();
    }
    inline node* build(node *fa,int *a,int l,int r){
        int mid=(l+r)>>1;
        node *tmp=recyclepool.newnode(a[mid]);tmp->fa=fa;
        if(l==r)return tmp;
        if(l<mid)tmp->lc=build(tmp,a,l,mid-1);
        if(r>mid)tmp->rc=build(tmp,a,mid+1,r);
        tmp->upt();
        return tmp;
    }
    inline void insert(int pos,int *a,int tot){
        node *tmp=build(null,a,1,tot);
        splay(findpos(rt,pos+1),null);splay(findpos(rt,pos+2),rt);
        tmp->fa=rt->rc;rt->rc->lc=tmp;rt->rc->upt();rt->upt();
    }
}splay;
int a[Maxn];
inline void print(node *now){
    if(now->lc!=null)print(now->lc);
    cout<<now->val<<" ";
    if(now->rc!=null)print(now->rc);
}
int main(){
io.init();n=io.read(),m=io.read();splay.rt=null;nowtot=n;
    splay.insert(splay.rt,null,0,-INF);
    splay.insert(splay.rt,null,1,-INF);
    for(int i=1;i<=n;i++){splay.insert(splay.rt,null,i,io.read());}
    for(int i=1;i<=m;i++){
        static char ch[15];
        cin>>(ch+1);
        if(ch[1]=='G'){
            int pos=io.read(),tot=io.read();
            if(!tot){io.W(0);continue;}
            io.W(splay.querysum(pos,pos+tot-1));
        }
        else if(ch[1]=='M'&&ch[3]=='X'){
            if(!nowtot){io.W(0);continue;}
            io.W(splay.rt->maxx);
        }
        else if(ch[1]=='I'){
            int pos=io.read(),tot=io.read();nowtot+=tot;
            for(int j=1;j<=tot;j++){
                a[j]=io.read();
            }
            splay.insert(pos,a,tot);
        }
        else if(ch[1]=='D'){
            int pos=io.read(),tot=io.read();nowtot-=tot;
            splay.del(pos,pos+tot-1);
        }
        else if(ch[1]=='R'){
            int pos=io.read(),tot=io.read();
            splay.rev(pos,pos+tot-1);
        }
        else{
            int pos=io.read(),tot=io.read(),c=io.read();
            if(!tot)continue;
            splay.add(pos,pos+tot-1,c);        
        }
    }
}

无旋Treap:

#include<bits/stdc++.h>
using namespace std;
typedef unsigned int uint;
struct IO{
    streambuf *ib,*ob;
    inline void init(){
        ios::sync_with_stdio(false);
        cin.tie(NULL);cout.tie(NULL);
        ib=cin.rdbuf();ob=cout.rdbuf();
    }
    inline int read(){
        char ch=ib->sbumpc();int i=0,f=1;
        while(!isdigit(ch)){if(ch=='-')f=-1;ch=ib->sbumpc();}
        while(isdigit(ch)){i=(i<<1)+(i<<3)+ch-'0';ch=ib->sbumpc();}
        return i*f;
    }
    inline void W(long long x){
        static int buf[50];
        if(!x){ob->sputc('0');ob->sputc('\n');return;}
        if(x<0){ob->sputc('-');x=-x;}
        while(x)buf[++buf[0]]=x%10,x/=10;
        while(buf[0])ob->sputc(buf[buf[0]--]+'0');
        ob->sputc('\n');
    }
    inline uint unit(){
        static uint status0=23333;
        status0^=status0<<13;
        status0^=status0>>17;
        status0^=status0<<5;
        return status0;
    }
}io;

const int Maxn=5e5+50,INF=0x3f3f3f3f;

int n,m,nowtot;

struct node{
    node *lc,*rc;
    int val,sum,revtag,addtag,maxx,lmax,rmax,sze,isadd;
    uint pri;
    node();
    inline void upt();
    inline void add(int v){
        val=v;sum=v*sze;
        int t=(v>0)?sze:1;
        maxx=lmax=rmax=t*v;
        addtag=v;isadd=1;
    }
    inline void rev(){
        swap(lc,rc);revtag^=1;
        swap(lmax,rmax);
    }
    inline void pushdown();
}Pool[Maxn],*pool=Pool,*null=Pool;
typedef pair<node*,node*> pii;
node::node():lc(null),rc(null),val(-INF),maxx(-INF),lmax(-INF),rmax(-INF),isadd(0), revtag(0){}
void node::upt() {
    sum=lc->sum+rc->sum+val;
    lmax=max(lc->lmax,max(lc->sum+val,lc->sum+val+rc->lmax));
    rmax=max(rc->rmax,max(rc->sum+val,rc->sum+val+lc->rmax));
    sze=lc->sze+rc->sze+1;
    maxx=max(val+max(lc->rmax,0)+max(rc->lmax,0),max(lc->maxx,rc->maxx));
}
void node::pushdown(){
    if(revtag){
        revtag=0;
        if(lc!=null)lc->rev();
        if(rc!=null)rc->rev();
    }
    if(isadd){
        if(lc!=null)lc->add(addtag);
        if(rc!=null)rc->add(addtag);
        isadd=0;
    }
}
struct RecyclePool{
    node *que[Maxn];
    int tail;
    inline node* newnode(int val){
        static node* tmp;
        tmp=(tail?que[tail--]:++pool);
        tmp->lc=tmp->rc=null;tmp->isadd=0;tmp->pri=io.unit();
        tmp->revtag=tmp->addtag=0, tmp->sze=1;
        tmp->val=tmp->maxx=tmp->lmax=tmp->rmax=val;
        return tmp;
    }
    inline void recycle(node *u){
        if(u==null)return;
        que[++tail]=u;
        recycle(u->lc);recycle(u->rc);
    }
}recyclepool;
struct Treap{
    node *rt;
    inline node* merge(node *x,node *y){
        if(x==null)return y;
        if(y==null)return x;
        x->pushdown(), y->pushdown();
        if(x->pri < y->pri){
            x->rc = merge(x->rc, y), x->upt();
            return x;
        }
        else{
            y->lc = merge(x, y->lc), y->upt();
            return y;
        }
    }
    inline pii split(node *u,int k){
        if(u==null) return pii(null, null);
        u->pushdown();
        if (u->lc->sze + 1 <= k){
            pii tr = split(u->rc, k - u->lc->sze - 1);
            u->rc = tr.first, u->upt();
            return pii(u,tr.second);
        }
        else {
            pii tr = split(u->lc,k);
            u->lc = tr.second;u->upt();
            return pii(tr.first,u);
        }
    }
    inline int querysum(int l,int tot){
        pii tr=split(rt,l-1);pii tr2=split(tr.second,tot);
        int ans = tr2.first->sum;
        rt=merge(tr.first,merge(tr2.first,tr2.second));
        return ans;
    }
    inline void del(int l,int tot){
        pii tr=split(rt,l-1);pii tr2=split(tr.second,tot);
        recyclepool.recycle(tr2.first);
        rt=merge(tr.first,tr2.second);
    }
    inline void rev(int l,int tot){
        pii tr=split(rt,l-1);pii tr2=split(tr.second,tot);
        tr2.first->rev();
        rt=merge(tr.first,merge(tr2.first,tr2.second));
    }
    inline void add(int l,int tot,int val){
        pii tr=split(rt,l-1);pii tr2=split(tr.second,tot);
        tr2.first->add(val);
        rt=merge(tr.first,merge(tr2.first,tr2.second));
    }
    inline node* build(int *a,int n){
        static node* st[Maxn],*pre;static int tail;
        tail=0;
        for(int i=0;i<n;i++){
            node *u=recyclepool.newnode(a[i]);pre=null;
            while(tail&&st[tail]->pri>u->pri)
                st[tail]->upt(),pre=st[tail--];
            if(tail)st[tail]->rc=u;
            u->lc=pre;st[++tail]=u;
        }
        while(tail)st[tail--]->upt();
        return st[1];
    }
    inline void insert(int *a, int n, int pos) {
        node *u = build(a, n);
        pii t = split(rt, pos);
        rt = merge(merge(t.first, u), t.second);
    }
}treap;
int a[Maxn];

int main(){
    io.init();n=io.read(),m=io.read();treap.rt=null;nowtot=n;
    for(int i=0;i<n;i++){a[i] = io.read();}
    treap.rt=treap.build(a,n);
    for(int i=1;i<=m;i++){
        static char ch[15];
        cin>>(ch+1);
        if(ch[1]=='G'){
            int pos=io.read(),tot=io.read();
            if(!tot){io.W(0);continue;}
            io.W(treap.querysum(pos,tot));
        }
        else if(ch[1]=='M'&&ch[3]=='X'){
            if(!nowtot){io.W(0);continue;}
            io.W(treap.rt->maxx);
        }
        else if(ch[1]=='I'){
            int pos=io.read(),tot=io.read();nowtot+=tot;
            for(int j=0;j<tot;j++){
                a[j] = io.read();
            }
            treap.insert(a, tot, pos);
        }
        else if(ch[1]=='D'){
            int pos=io.read(),tot=io.read();nowtot-=tot;
            if(!tot)continue;
            treap.del(pos,tot);
        }
        else if(ch[1]=='R'){
            int pos=io.read(),tot=io.read();
            if(!tot)continue;
            treap.rev(pos,tot);
        }
        else{
            int pos=io.read(),tot=io.read(),c=io.read();
            if(!tot)continue;
            treap.add(pos,tot,c);        
        }
    }
}
DZYO
关注
专栏目录
【BZOJ】1500: [NOI2005]维修数列-splay
远行客
12-05 144
splay模板
BZOJ1500: [NOI2005]维修数列（洛谷P2042）
forezxl的博客
12-09 348
Splay
[BZOJ1500][NOI2005]维修数列（平衡树）
薇小薇
03-15 269
题目： 我是超链接 题解： md这题调的我哭出来，最后发现是all数组del的时候没清零Orz 反正是平衡树的基本操作吧 因为害怕转着转着出现奇怪的错误，我们把1和n+2当做空节点，权值为-INF，但是sum,maxx都是0，避免影响结果，注意这里0节点也要处理一下，因为我们用ch[now][0/1]的时候节点也是0，容易GG Insert：只需把pos点splay到根，pos+1sp...
bzoj 1500
diyan5166的博客
03-08 125
再写一次毒瘤题，改变代码风格改用了指针，以为会不太适应，事实上不会撒。。 对拍数据要等20分钟才出来= = 1 //#include<bits/stdc++.h> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #inc...
bzoj1500
lpf_as_a_oier的博客
12-27 972
这题是裸的平衡树。不过我却调了很久，果然还是太渣。 最大和子序列比较难维护，我是先维护LMAX表示从左边起区间的连续最大和，RMAX表示从右边起区间的连续最大和，然后SUM1(最大和子序列)就是等于MAX（SUM1[LC[X]],SUM1[RC[X]],LMAX[RC[X]]+KEY[X]+RMAX[LC[X]]） 其他的好像都不难，记得插入的时候不要一个一个插 很明显我们要设定2个标记，翻
BZOJ 1500 [NOI2005]维修数列 （无旋treap良心题解）
shadiao.的博客
09-10 263
题目传送门 上面链接沉了幸亏洛谷还有这道题 无旋treap解题的原理： 1.无旋treap可以在维护平衡的同时保持中序遍历不变（treap也可以，应该都可以），因此我们用下标代表插入二叉树的值，那么我们可以按照中序遍历获得原来的数组。 2.无旋treap的分裂操作可以分裂出包含前k个节点的树和剩下元素的树，就相当于可以把数组切开。 你可以切取想要的区间进行操作，操作完再merge回去即可...
NOI2005 维护数列
Holylight_Knight
10-10 634
NOI2005 维护数列tags :【Splay】,【BZOJ1500】,【洛谷2042】在线评测地址：洛谷，大视野题目描述1500: [NOI2005]维修数列Time Limit: 10 Sec Memory Limit: 64 MBDescriptionInput输入的第1 行包含两个数N 和M(M ≤20 000)，N 表示初始时数列中数的个数，M表示要进行的操作数目。第2行包含N个数字，
splay伸展树 指针型 平衡树基本操作 序列维护 详细讲解+总结
anantheparty的博客
12-21 3486
转载请保留本博客源地址：http://blog.csdn.net/u011327397/article/details/53783700 本来是要去学lct，然后看到要用splay我又不会，就开始看splay了，然后发现splay艹起来特别爽，但是序列操作网上没看到有详讲的，理解的时候很多奇怪的问题（无限RE）再加上指针调试复杂度高，花了很多时间，就觉得来写一个讲解吧。
最强平衡树——Treap[以我的最弱击败你的最强]
lemonoil的博客
05-13 671
以我的最弱战胜你的最强！！——Treap本人蒟蒻，在平衡树坑中深陷数年。为了早日逃离此天坑，特作此文。什么是平衡树?度娘传送门 什么是treap？ACdreamers%%% 注：本篇所有代码都在片尾！！（醒目）CMP那么了解了这些，我们先列出一个list NAME 优势 劣势 splay LCT，序列之王 常数大，代码量稍大 RBT 自适应深度平衡树，速度在同类BS
【NOI 2005】【JZOJ 2413】维护数列 平衡树操作模板题
LYD729
01-24 528
平衡树模板题 splay操作大全
BZOJ 1500([NOI2005]维修数列-Splay的数列维护)
nike0good |Oier&ACMer | 熟能生巧
02-13 2020
1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 3087  Solved: 920 [Submit][Status][Discuss] Description Input 输入文件的第1行包含两个数N和M，N表示初始时数列中数的个数，M表示要进行的操作数目。 第2行包含N个数字，描述初始
poj3274 找平衡数列（哈希加一点数学思维）
weixin_30872157的博客
05-09 110
题目传送门题目大意：有n只牛，每只牛有k个属性，接下来n个数字，每个数字的二进制位上的1和0分别表示某种属性的有或者无，然后一个特殊数列就是，一个区间内所有牛的各种属性的总和相等（有e种1属性 e种2属性andsoon），问你这排牛的最长的特殊数列长度是多少。思路：看上去像dp，但思路走不通，然后看网上大佬的思路，仿佛推开新世界的大门。数组sum[i][j]...
【BZOJ1500】[NOI2005]维修数列【Splay】
BraketBN
03-16 426
http://www.lydsy.com/JudgeOnline/problem.php?id=1500 这神题终于过了... 主要是： find()那里少写了一个pushdown()。 辣鸡回收写跪了。 /* Footprints In The Blood Soaked Snow */ #include #include using namespace s
[BZOJ]1500 维修数列
Le Petit Prince
07-25 399
1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 14433  Solved: 4701 [Submit][Status][Discuss] Description Input 输入的第1 行包含两个数N 和M(M ≤20 000)，N 表示初始时数列中数的个数，M表示要进行的操作数目
bzoj1500（妥妥的splay模板题）
cq_phqg
09-04 1680
1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 6366  Solved: 1910 [Submit][Status] Description Input 输入文件的第1行包含两个数N和M，N表示初始时数列中数的个数，M表示要进行的操作数目。第2行包含N个数字，描述初始时的数列。以下
[BZOJ1500][NOI2005]维护数列(无旋Treap)
Ephemeral
03-06 453
1500: [NOI2005]维修数列Time Limit: 10 Sec&nbsp;&nbsp;Memory Limit: 64 MBSubmit: 16199&nbsp;&nbsp;Solved: 5391 Description请写一个程序，要求维护一个数列，支持以下 6 种操作： 请注意，格式栏 中的下划线‘ _ ’表示实际输入文件中的空格 Input输入的第1 行包...
BZOJ2839:集合计数(容斥,组合数学)
最新发布
06-08
这个问题是一个比较经典的容斥原理的应用题目。我们可以考虑使用容斥原理来求解这个问题。具体来说，我们可以考虑求出所有集合的个数，然后减去至少包含一个特定元素的集合的个数，再加上至少包含两个特定元素的集合...
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值