NOI2005 维修数列 Splay

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

题目大意：
给出一个初始序列，要求实现下面的操作：
1.在当前序列的某个位置插入连续的若干个元素；
2.删除当前序列的某一段元素；
3.将当前序列中某一段元素全部修改为某个值；
4.翻转某个区间；
5.求当前序列某一段元素的总和；
6.求出当前序列中连续和的最大值。

一看就是Splay裸题，一看就知道操作很多很恶心。

各种操作当然要写准确了，主要注意下面几点：

1.显然有lazy操作，那么如何下放？何时下放？

下放时可以参见线段树的处理方法，打上标记时，当前节点已经更新完毕。
由于几乎每个操作都要首先提取区间，而提取区间总要用到查找第k大数的函数。这个函数是自根向下进行的，这时遇到lazy标记显然必须下放。同时，找到编号后就会立即Splay,而这个点上跳时经过的点一定恰好是执行查找第k大函数时访问到的点，也就是说，这个时候该下放的一定都下放完毕了。
当然，为了保险，也可以在Splay操作的旋转前把节点到根路径上的点放在栈里面，先下放再旋转。

2.如何更新某个节点的各种值？

尤其注意Splay维护区间时更新与线段树有不同。线段树中，叶子节点的值是容易操作的，非叶子节点的的各种值仅由它的左右儿子的值决定，而且一定既有左儿子，又有右儿子。
在Splay中，节点本身的各种值不仅受儿子影响，也有本身的影响。而且，非叶子节点有可能只有一个儿子。这个时候可以写四种情况逐个判断，也可以把编号0的节点（代表没有儿子）的值适当处理，使得不会对答案产生影响即可。

3.关于虚拟节点。

为了每次都能把想要的区间提取到根的右儿子的左子树上，分别在序列的首位添加一个占位元素。占位元素的Size必须是1，而其他值可以任取。处理操作6时，虽然已知讯问区间，还是要把整个区间提取出来，否则占位元素的值会对答案产生影响。（我就是这个操作没有处理对，调了2h+）

此外，插入操作和初始序列的建树可以采用与线段树类似的方法构建，这样不难实现，而且插入的是一棵严格的平衡树。

代码：

#include<stdio.h>
#include<algorithm>
#include<cstring>
#define MAXT 5000005
#define MAXN 500005
using namespace std;

int N,M;

int Tmax(int x,int y,int z)
{
    if(x>=y&&x>=z)return x;
    if(y>=z)return y;
    return z;
}

int rt,tot,fa[MAXT],ls[MAXT],rs[MAXT],Size[MAXT];
int val[MAXT],Sum[MAXT],MaxL[MAXT],MaxR[MAXT],MaxV[MAXT];
int lazy[MAXT],rev[MAXT];bool cl[MAXT];

void Update(int p)
{
    int l=ls[p],r=rs[p];

    Size[p]=Size[l]+Size[r]+1;
    Sum[p]=Sum[l]+Sum[r]+val[p];

    MaxL[p]=Tmax(MaxL[l],Sum[l]+val[p],Sum[l]+val[p]+MaxL[r]);
    MaxR[p]=Tmax(MaxR[r],Sum[r]+val[p],Sum[r]+val[p]+MaxR[l]);
    MaxV[p]=Tmax(MaxV[l],MaxV[r],Tmax(MaxR[l]+val[p],MaxL[r]+val[p],MaxR[l]+val[p]+MaxL[r]));
    if(MaxV[p]<val[p])MaxV[p]=val[p];
}

void Putdown(int p)
{
    if(cl[p]==false&&rev[p]==0)return;
    int l=ls[p],r=rs[p];
    if(rev[p])
    {
        rev[p]=0;rev[l]^=1;rev[r]^=1;
        swap(MaxL[l],MaxR[l]);swap(ls[l],rs[l]);
        swap(MaxL[r],MaxR[r]);swap(ls[r],rs[r]);
    }
    if(cl[p])
    {
        cl[p]=false;cl[l]=true;cl[r]=true;
        lazy[l]=lazy[r]=lazy[p];
        val[l]=val[r]=lazy[p];
        Sum[l]=Size[l]*lazy[p];
        Sum[r]=Size[r]*lazy[p];
        if(lazy[p]>0)
        {
            MaxL[l]=MaxR[l]=MaxV[l]=Size[l]*lazy[p];
            MaxL[r]=MaxR[r]=MaxV[r]=Size[r]*lazy[p];
        }
        else
        {
            MaxL[l]=MaxR[l]=MaxV[l]=lazy[p];
            MaxL[r]=MaxR[r]=MaxV[r]=lazy[p];
        }
    }
}

void Zig(int x)
{
    int y=fa[x],z=fa[y];
    if(z)
    {
        if(ls[z]==y)ls[z]=x;
        else rs[z]=x;
    }
    fa[x]=z;fa[y]=x;fa[rs[x]]=y;
    ls[y]=rs[x];rs[x]=y;
    Update(y);
    Update(x);
}

void Zag(int x)
{
    int y=fa[x],z=fa[y];
    if(z)
    {
        if(ls[z]==y)ls[z]=x;
        else rs[z]=x;
    }
    fa[x]=z;fa[y]=x;fa[ls[x]]=y;
    rs[y]=ls[x];ls[x]=y;
    Update(y);
    Update(x);
}

void Splay(int x,int t)
{
    int y,z;
    while(fa[x]!=t)
    {
        y=fa[x];z=fa[y];
        if(z==t)
        {
            if(ls[y]==x)Zig(x);
            else Zag(x);
        }
        else
        {
            if(ls[z]==y){
                if(ls[y]==x)Zig(y),Zig(x);
                else Zag(x),Zig(x);
            }
            else
            {
                if(rs[y]==x)Zag(y),Zag(x);
                else Zig(x),Zag(x);
            }
        }
    }
    if(!t)rt=x;
}

int A[MAXN];
int Build(int x,int y,int f)
{
    int p=++tot,mid=x+y>>1;
    fa[p]=f;val[p]=A[mid];
    if(x==y)
    {
        Size[p]=1;
        Sum[p]=MaxL[p]=MaxR[p]=MaxV[p]=A[mid];
        return p;
    }
    if(x<=mid-1)ls[p]=Build(x,mid-1,p);
    if(mid+1<=y)rs[p]=Build(mid+1,y,p);
    Update(p);
    return p;
}

int GetKth(int k)
{
    int p=rt;
    while(p)
    {
        Putdown(p);
        if(Size[ls[p]]+1==k)return p;
        if(Size[ls[p]]>=k)p=ls[p];
        else k-=Size[ls[p]]+1,p=rs[p];
    }
}

void Extract(int l,int r)//将区间[l,r]提取到根的右儿子的左子树上
{
    int p,q;
    p=GetKth(l);
    Splay(p,0);
    q=GetKth(r+2);
    Splay(q,rt);
}

void Ins()
{
    int pos,sz,i,p,q;

    scanf("%d%d",&pos,&sz);
    for(i=1;i<=sz;i++)scanf("%d",&A[i]);
    p=GetKth(pos+1);Splay(p,0);
    q=GetKth(pos+2);Splay(q,rt);
    ls[rs[rt]]=Build(1,sz,rs[rt]);
    Update(rs[rt]);
    Update(rt);
}

void Del(int l,int r)
{
    Extract(l,r);
    int p=ls[rs[rt]];
    fa[p]=0;ls[rs[rt]]=0;
    Update(rs[rt]);
    Update(rt);
}

void Modify(int l,int r,int k)
{
    Extract(l,r);
    int p=ls[rs[rt]];

    lazy[p]=k;cl[p]=true;
    val[p]=k;
    Sum[p]=Size[p]*k;
    if(k>0)MaxL[p]=MaxR[p]=MaxV[p]=Size[p]*k;
    else MaxL[p]=MaxR[p]=MaxV[p]=k;

    Update(rs[rt]);
    Update(rt);
}

void Reverse(int l,int r)
{
    Extract(l,r);
    int p=ls[rs[rt]];
    rev[p]^=1;
    swap(MaxL[p],MaxR[p]);
    swap(ls[p],rs[p]);
}

int GetSum(int l,int r)
{
    Extract(l,r);
    return Sum[ls[rs[rt]]];
}


int main()
{
    Sum[0]=0;
    MaxL[0]=MaxR[0]=MaxV[0]=val[0]=-1e9;
    Size[0]=0;
    //上面三行就是“对0号点的适当处理”

    int i,pos,sz,c;
    char op[20];

    scanf("%d%d",&N,&M);
    for(i=1;i<=N;i++)scanf("%d",&A[i+1]);

    rt=Build(1,N+2,0);
    for(i=1;i<=M;i++)
    {
        scanf("%s",op);
        if(op[0]=='I')Ins();
        else if(op[0]=='D')
        {
            scanf("%d%d",&pos,&sz);
            Del(pos,pos+sz-1);
        }
        else if(op[0]=='R')
        {
            scanf("%d%d",&pos,&sz);
            Reverse(pos,pos+sz-1);
        }
        else if(op[0]=='G')
        {
            scanf("%d%d",&pos,&sz);
            printf("%d\n",GetSum(pos,pos+sz-1));
        }
        else if(op[2]=='K') 
        {
            scanf("%d%d%d",&pos,&sz,&c);
            Modify(pos,pos+sz-1,c);
        }
        else
        {
            Extract(1,Size[rt]-2);
            printf("%d\n",MaxV[ls[rs[rt]]]);
        }
    }
}

此外，Update函数也可以用分类讨论写，只不过略恶心：

void Update(int p)
{
    int l=ls[p],r=rs[p];

    if(l)
    {
        if(r)
        {
            Size[p]=Size[l]+Size[r]+1;
            Sum[p]=Sum[l]+Sum[r]+val[p];
            MaxL[p]=Tmax(MaxL[l],Sum[l]+val[p],Sum[l]+val[p]+MaxL[r]);
            MaxR[p]=Tmax(MaxR[r],Sum[r]+val[p],Sum[r]+val[p]+MaxR[l]);
            MaxV[p]=Tmax(MaxV[l],MaxV[r],Tmax(val[p]+MaxR[l],val[p]+MaxL[r],MaxR[l]+val[p]+MaxL[r]));
            MaxV[p]=max(MaxV[p],val[p]);
        }
        else
        {
            Size[p]=Size[l]+1;
            Sum[p]=Sum[l]+val[p];
            MaxL[p]=max(MaxL[l],Sum[l]+val[p]);
            MaxR[p]=max(val[p],MaxR[l]+val[p]);
            MaxV[p]=Tmax(val[p],MaxV[l],MaxR[l]+val[p]);
        }
    }
    else
    {
        if(r)
        {
            Size[p]=Size[r]+1;
            Sum[p]=Sum[r]+val[p];
            MaxL[p]=max(val[p],val[p]+MaxL[r]);
            MaxR[p]=max(MaxR[r],Sum[r]+val[p]);
            MaxV[p]=Tmax(val[p],MaxV[r],MaxL[r]+val[p]);
        }
        else
        {
            Size[p]=1;
            Sum[p]=val[p];
            MaxL[p]=MaxR[p]=MaxV[p]=val[p];
        }
    }
}
rgnoH
关注
专栏目录
[NOI2005] 维护数列（Splay）
hipamp的博客
10-11 375
题意 要求维护一个数列，支持以下操作： 其中，插入总数字个数小于 4e64e64e6，每个数在 [−1000,1000][-1000,1000][−1000,1000] 范围，保证每次操作合法。 分析 这道题是一道 SplaySplaySplay 码农题，细节颇多，让我们一个一个操作分析。 插入操作 我们可以先对要插入的这 tottottot 个数建树，得到树根 trttrttrt，再将 trttrttrt 插入原序列。 删除操作 就是普通的删除，将 r+1r+1r+1 旋到根，将 l−1l-1l−1
【BZOJ】1500: [NOI2005]维修数列-splay
远行客
12-05 142
splay模板
【平衡树】【NOI 2005】维护数列
morestep的专栏
07-27 965
339. [NOI2005] 维护数列★★★★☆ 输入文件：seq2005.in 输出文件：seq2005.out 简单对比 时间限制：3 s 内存限制：256 MB 【问题描述】请写一个程序，要求维护一个数列，支持以下 6 种操作：（请注意，格式栏 中的下划线‘ _ ’表示实际输入文件中的空格） 【输入格式】输入文件的第 1 行包含两个数 N 和 M，N 表示初始时数列中数的个数
[NOI2005]维护数列——Splay
weixin_30852419的博客
07-22 71
题面 　　洛谷P2042 解析 　　　　这道题集齐了大部分Splay的经典操作，上课时wys老师讲了一下，这个题还是很值得一做的 　　基操： 　　　　Splay上的所有区间操作都需要提取区间，如果要提取一个[l, r]的区间，只需要把l-1的点转到根，再把r+1的点转到根的右儿子，那么根的右儿子的左儿子的整颗子树就是我们要提取的区间 　　操作1：插入 　　　　这个题的插入并非插入一个...
[NOI2005]维修数列
a5163273的博客
08-15 123
Description 请写一个程序，要求维护一个数列，支持以下 6 种操作： 请注意，格式栏 中的下划线‘ _ ’表示实际输入文件中的空格 Input 输入的第1 行包含两个数N 和M(M ≤20 000)，N 表示初始时数列中数的个数，M表示要进行的操作数目。 第2行包含N个数字，描述初始时的数列。 以下M行，每行一条命令，格式参见问题描述中的表格。 任何时刻数列中最多含有500 ...
1500: [NOI2005]维修数列（Splay）
qq_43544481的博客
12-25 142
Description 请写一个程序，要求维护一个数列，支持以下 6 种操作： 请注意，格式栏 中的下划线‘ _ ’表示实际输入文件中的空格 Input 输入的第1 行包含两个数N 和M(M ≤20 000)，N 表示初始时数列中数的个数，M表示要进行的操作数目。 第2行包含N个数字，描述初始时的数列。 以下M行，每行一条命令，格式参见问题描述中的表格。 任何时刻数列中最多含有500 000个数...
BZOJ1500: [NOI2005]维修数列（洛谷P2042）
forezxl的博客
12-09 343
Splay
大视野在线测评 1500 [NOI2005]维修数列（splay）
madaidao的专栏
05-27 1387
题目链接 1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 8379  Solved: 2521 [Submit][Status][Discuss] Description Input 输入文件的第1行包含两个数N和M，N表示初始时数列中数的个数，M表示要进行的操作数目。
NOI2005 维护数列
Holylight_Knight
10-10 628
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个数字，
noi2005测试数据
09-13
noi2005测试数据
[NOI 2005] 维护数列 splay
九野的博客
05-30 1695
题目链接：http://www.lydsy.com:808/JudgeOnline/problem.php?id=1500 写的比较小心，还好没有经历虐心的wa。。 数据结构的经典神题。。敲完才发现搞了250行。。 思路还是很简单的，和线段树的求最大子段和一样，左连续右连续什么的，， 对于区间操作都把这个区间转到根的右子树的左子树下。。 1、这题需要搞个内存池，开始t
NOI2005 维护数列(splay)
weixin_30662539的博客
12-15 96
学了半天平衡树，选择了一道题来写一写，发现题目是裸的splay模板，但是还是写不好，这个的精髓之处在于在数列的某一个位置加入一个数列，类似于treap里面的merge，然后还学到了题解里面的的回收空间做法，利用时间换取空间也是一种非常不错的做法，由于题目中有区间翻转的操作，所以只能用splay来写，然后对于一些区间最大值和区间最大子序列的询问我们可以用类似于线段树的合并操作来维护，总的来说这是一道...
noi2005维护数列 splay
黑历史
03-09 283
用splay来维护数列的一道很好的模版练习题 有些要注意的地方： 不用指针来实现splay的话要注意内存回收，不然就会无限RE+TLE凡是有对序列进行修改的操作都要up一下reverse的时候要注意交换左起最大连续和跟右起最大连续和 代码： #include #include #include #include #include #include #include #include
P2042 [NOI2005]维护数列 splay
haiku的博客
11-09 194
学块状链表的时候找到的这个题，看了大佬的思路发现splay真香，块状链表虽然好理解，但写起来真的有点烦。 P2042 [NOI2005]维护数列 对一个数列进行插入，删除，修改（统一置为某数），区间翻转，区间求和，整体最大子列和。 用splay维护序列的基本操作：先通过旋转把要操作的区间单独拎出来，比如要操作从k开始的tot个数，那就先把下标k-1对应的结点转到整棵树的根，然后把k+tot转到根的...
[NOI2005] 维修数列 (Splay)
xyyxyyx的博客
03-22 171
题意 思路 注意 代码 ver. 1 #include<bits/stdc++.h> using namespace std; const int inf = 10000; const int N = 5e5+10; const int M = 4e6+10; #define rdi read<int> #define rdl read<long long> t...
【NOI2005】维护数列
WorldWide_D的博客
05-30 1504
Description请写一个程序，要求维护一个数列，支持以下6种操作：（请注意，格式栏中的下划线‘ _ ’表示实际输入文件中的空格） 1. 插入 INSERT_posi_tot_c1_c2_…_ctot 在当前数列的第posi个数字后插入tot个数字：c1, c2, …, ctot；若在数列首插入，则posi为0 2. 删除 DELETE_posi_tot 从当前数列的第posi个数字开始连续
[BZOJ1500][NOI2005]维修数列（splay）
zyf2000
12-04 915
题目描述传送门题解splay模板题（笑 这题是我今年1月左右看到的。记得当时只调过了样例然后交上去过了一个点就用了我一天的时间。果然我还是太傻逼。 填上了好大一个坑。 前5个都是splay的基础操作，就是第6个比较麻烦。因为要维护区间最大连续子序列和，可以维护3个值：Max子树最大连续子序列和，maxl子树自左端起最大连续子序列和，maxr子树自右端起最大连续子序列和。注意这里的Max是不允许
BZOJ 1500 [NOI2005]维修数列(splay)
weixin_33753003的博客
01-08 77
题目链接：http://61.187.179.132/JudgeOnline/problem.php?id=1500 题意：一个数列，维护以下操作。 思路：以下A代表根节点，B在根节点右子树，C代表根节点右子树的左子树即B的左子树 （1）插入时，将pos-1和pos+cnt分别调至A和B，将插入节点建立一棵子树插在C即可； （2）删除时，将pos-1和pos+cnt分别调至A和B，将C...
NOI比赛专用平台：NOI Linux详解
"NOI系列比赛平台NOI Linux是一款专为NOI（全国青少年信息学奥林匹克竞赛）设计的定制化操作系统，旨在提供一个统一、标准化的竞赛环境，消除因平台差异导致的评分不公。该系统由NOI科学委员会授权，由评测小组基于...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值