1500: [NOI2005]维修数列

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 11047  Solved: 3460
[Submit][Status][Discuss]

Description

Input

输入的第1 行包含两个数N 和M(M ≤20 000),N 表示初始时数列中数的个数,M表示要进行的操作数目。
第2行包含N个数字,描述初始时的数列。
以下M行,每行一条命令,格式参见问题描述中的表格。
任何时刻数列中最多含有500 000个数,数列中任何一个数字均在[-1 000, 1 000]内。
插入的数字总数不超过4 000 000个,输入文件大小不超过20MBytes。

Output

对于输入数据中的GET-SUM和MAX-SUM操作,向输出文件依次打印结果,每个答案(数字)占一行。

Sample Input

9 8
2 -6 3 5 1 -5 -3 6 3
GET-SUM 5 4
MAX-SUM
INSERT 8 3 -5 7 2
DELETE 12 1
MAKE-SAME 3 3 2
REVERSE 3 6
GET-SUM 5 4
MAX-SUM

Sample Output

-1
10
1
10

HINT

Source

Solution

splay模板题,exciting

插入和删除一个序列,有种方便快捷的方法

进行两次伸展,一次伸展到根,一次伸展到右子树上,这样会可以使得左子树空缺或需要删除是左子树

那么直接插入一整棵树到左子树,或者删除整个左子树即可

那么给定的序列,可以利用 线段树类似的build技巧,将其建成一个完美的二叉树,然后直接加入

BZOJ卡内存,所以删除的时候需要垃圾回收

具体的方法就是 使用一个队列,记录扔掉的点的编号,加入新东西时,先从队列中取编号,再建新编号即可

Code

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
using namespace std;
int read()
{
int x=,f=; char ch=getchar();
while (ch<'' || ch>'') {if (ch=='-') f=-; ch=getchar();}
while (ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}
return x*f;
}
#define maxn 1000100
queue<int>trashpool;
int n,m,cnt;
int fa[maxn],son[maxn][],key[maxn],size[maxn],sum[maxn],val[maxn],maxz[maxn],ll[maxn],rr[maxn],id[maxn];
bool rev[maxn],tag[maxn];
int sz,root;
void update(int now)
{
//if (!now) return;
size[now]=size[son[now][]]+size[son[now][]]+;
sum[now]=val[now]+sum[son[now][]]+sum[son[now][]];
maxz[now]=max(maxz[son[now][]],maxz[son[now][]]);
maxz[now]=max(maxz[now],rr[son[now][]]+val[now]+ll[son[now][]]);
ll[now]=max(ll[son[now][]],sum[son[now][]]+val[now]+ll[son[now][]]);
rr[now]=max(rr[son[now][]],sum[son[now][]]+val[now]+rr[son[now][]]);
}
void pushdown(int now)
{
int tg=tag[now],rv=rev[now];
if (tag[now])
{
tag[now]=; rev[now]=;
if (son[now][])
tag[son[now][]]=,val[son[now][]]=val[now],sum[son[now][]]=size[son[now][]]*val[now];
if (son[now][])
tag[son[now][]]=,val[son[now][]]=val[now],sum[son[now][]]=size[son[now][]]*val[now];
if(val[now]>=)
{
if(son[now][]) ll[son[now][]]=rr[son[now][]]=maxz[son[now][]]=sum[son[now][]];
if(son[now][]) ll[son[now][]]=rr[son[now][]]=maxz[son[now][]]=sum[son[now][]];
}
else
{
if(son[now][]) ll[son[now][]]=rr[son[now][]]=,maxz[son[now][]]=val[now];
if(son[now][]) ll[son[now][]]=rr[son[now][]]=,maxz[son[now][]]=val[now];
}
return;//
}
if (rev[now])
{
rev[now]^=; rev[son[now][]]^=; rev[son[now][]]^=;
swap(ll[son[now][]],rr[son[now][]]); swap(ll[son[now][]],rr[son[now][]]);
swap(son[son[now][]][],son[son[now][]][]); swap(son[son[now][]][],son[son[now][]][]);
}
}
void rotate(int x,int &k)
{
int y=fa[x],z=fa[y],l,r;
l=(son[y][]==x); r=l^;
if(y==k) k=x;
else son[z][son[z][]==y]=x;
fa[son[x][r]]=y; fa[y]=x;fa[x]=z;
son[y][l]=son[x][r]; son[x][r]=y;
update(y); update(x);
}
void splay(int x,int &k)
{
while(x!=k)
{
//printf("%d %d\n",x,k);
int y=fa[x],z=fa[y];
if(y!=k)
{
if(son[y][]==x^son[z][]==y)rotate(x,k);
else rotate(y,k);
}
rotate(x,k);
}
}
int findkth(int x,int rk)
{
pushdown(x);
int l=son[x][],r=son[x][];
if(size[l]+==rk) return x;
if(size[l]>=rk) return findkth(l,rk);
return findkth(r,rk-size[l]-);
}
void trash(int now)
{
if (!now) return;
trash(son[now][]),trash(son[now][]); trashpool.push(now);
fa[now]=son[now][]=son[now][]=tag[now]=rev[now]=;
}
int split(int pos,int tot)
{
int x=findkth(root,pos),y=findkth(root,pos+tot+);
splay(x,root);splay(y,son[x][]);
return son[y][];
}
void asksum(int pos,int tot)
{printf("%d\n",sum[split(pos,tot)]);}
void change(int pos,int tot,int va)
{
int x=split(pos,tot),y=fa[x];
val[x]=va; tag[x]=; sum[x]=size[x]*va;
if(va>=) ll[x]=rr[x]=maxz[x]=sum[x];
else ll[x]=rr[x]=,maxz[x]=va;
update(y); update(fa[y]);
}
void cover(int pos,int tot)
{
int x=split(pos,tot),y=fa[x];
if(!tag[x])
{
rev[x]^=;
swap(son[x][],son[x][]);swap(ll[x],rr[x]);
update(y); update(fa[y]);
}
}
void remove(int pos,int tot)
{
int x=split(pos,tot),y=fa[x];
trash(x); son[y][]=;
update(y); update(fa[y]);
}
void build(int l,int r,int father)
{
if(l>r) return;
int mid=(l+r)>>,now=id[mid],last=id[father];
if(l==r)
{
sum[now]=key[l]; size[now]=;
tag[now]=rev[now]=;
if(key[l]>=) ll[now]=rr[now]=maxz[now]=key[l];
else ll[now]=rr[now]=,maxz[now]=key[l];
// return;
}
else build(l,mid-,mid),build(mid+,r,mid);
val[now]=key[mid]; fa[now]=last; update(now);
son[last][mid>=father]=now;
}
void insert(int pos,int tot)
{
for (int i=; i<=tot; i++) key[i]=read();
for (int i=; i<=tot; i++)
if (!trashpool.empty()) id[i]=trashpool.front(),trashpool.pop(); else id[i]=++cnt;
build(,tot,); int mid=(+tot)>>;
int x=findkth(root,pos+),y=findkth(root,pos+),z=id[mid];
splay(x,root);splay(y,son[x][]);
fa[z]=y; son[y][]=z;
update(y); update(x);
} #define sb -100000000
int main()
{
n=read();m=read();
maxz[]=key[]=key[n+]=sb;
for (int i=; i<=n; i++) key[i+]=read();
for (int i=; i<=n+; i++) id[i]=i;
root=(n+)>>; cnt=n+;
build(,n+,);
char opt[];
while (m--)
{
scanf("%s",opt);int pos,tot,va;
switch (opt[])
{
case 'S': pos=read(),tot=read();insert(pos,tot); break;
case 'L': pos=read(),tot=read();remove(pos,tot); break;
case 'K': pos=read(),tot=read(),va=read(); change(pos,tot,va); break;
case 'V': pos=read(),tot=read();cover(pos,tot); break;
case 'T': pos=read(),tot=read();asksum(pos,tot); break;
case 'X': printf("%d\n",maxz[root]); break;
}
//for (int i=1; i<=n+1; i++) printf("%d ",val[i]);puts("");
//puts("---------------------------------------------------------------------");
}
return ;
}

真的是模板题,传说中NOI最难数据结构题??逗我?也就写了好几个小时罢了QAQ

一开始电脑爆炸,每次处理标记,n都会变化..自己瞎搞1h30min后,问char哥,char哥说,你的xx数组爆了,因为内存空间是连续的,所以波及了n...

然后30min后,电脑爆炸,于是强行考程序换YveH自带的笔记本,于是发现2h处理了电脑的智障问题....

然后发现很多问题...

such as:Splay死循啊,pushdown蛋疼啊,垃圾回收废物利用编号出毛病啊.....

以后有学弟找我要Splay模板题,不说啥了,就这道!!!

【BZOJ-1500】维修数列 Splay的更多相关文章

  1. [BZOJ 1500]维修数列 [Splay Tree从进阶到住院]

    历尽艰辛终于A掉了这题QwQ 贴COGS评论区几句话=.= 策爷:"splay/块状链表的自虐题.".深刻理解到如果没有M倾向就不要去写这题了.. -Chenyao2333 记得b ...

  2. BZOJ 1500 维修数列【Splay】

    注意:1,内存限制,所以需要回收删除的点 2,当前节点的左连续区间和最大值=max(左子树的左连续区间和最大值,左子树的总和+当节点的值+max(右子树的左连续区间和最大值,0)):右连续区间和最大值 ...

  3. bzoj 1500 维修数列

    splay乱搞. 调了两个多小时...这辈子再也不想写splay了... 维护左边最大连续和右边最大连续,维护两个标记,无脑push_down.push_up就行了. 注意最大连续和至少要包含一个数. ...

  4. bzoj 1500: [NOI2005]维修数列 splay

    1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 6556  Solved: 1963[Submit][Status ...

  5. BZOJ 1500: [NOI2005]维修数列 (splay tree)

    1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 4229  Solved: 1283[Submit][Status ...

  6. 【BZOJ1500】[NOI2005]维修数列 Splay

    [BZOJ1500][NOI2005]维修数列 Description Input 输入的第1 行包含两个数N 和M(M ≤20 000),N 表示初始时数列中数的个数,M表示要进行的操作数目.第2行 ...

  7. [NOI2005]维修数列 Splay tree 区间反转,修改,求和,求最值

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1500 Description Input 输入文件的第1行包含两个数N和M,N表示初始时数 ...

  8. BZOJ1500 [NOI2005]维修数列(Splay tree)

    [Submit][Status][Discuss] Description 请写一个程序,要求维护一个数列,支持以下 6 种操作: 请注意,格式栏 中的下划线‘ _ ’表示实际输入文件中的空格 Inp ...

  9. BZOJ1500: [NOI2005]维修数列[splay ***]

    1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 12278  Solved: 3880[Submit][Statu ...

随机推荐

  1. ToolkitScriptManager, ScriptManager AjaxControlToolkit 出错

    在使用ASP.NET AJAX应用时需添加ToolkitScriptManager或者ScriptManager.如果只是使用updatepanel,两者均可.如果有涉及到AJAX Control T ...

  2. windows Server 2008各版本区别详解

    Windows Server 2008 是专为强化下一代网络.应用程序和 Web 服务的功能而设计,是有史以来最先进的 Windows Server 操作系统.拥有 Windows Server 20 ...

  3. 产品经理技能之BRD的笔记之菜鸟入门

    链接:http://www.woshipm.com/pmd/178527.html?utm_source=tuicool 要学习MRD.PRD,先从BRD开始,才能做到知其然知其所以然. BRD是什么 ...

  4. 快速判断素数 --Rabin-Miller算法

    以前我在判断素数上一直只会 sqrt(n) 复杂度的方法和所谓的试除法(预处理出sqrt(n)以内的素数,再用它们来除). (当然筛选法对于判断一个数是否是素数复杂度太高) 现在我发现其实还有一种方法 ...

  5. struts2: 通过流输出实现exce导出

    参考下面代码,在Action中加一个方法: // 导出excel public String excel() throws Exception { StringBuffer excelBuf = ne ...

  6. python 播放 wav 文件

    未使用其他库, 只是使用 pywin32 调用系统底层 API 播放 wav 文件. # Our raison d'etre - playing sounds import pywintypes im ...

  7. TinyFrame再续篇:整合Spring AOP实现日志拦截

    上一篇中主要讲解了如何使用Spring IOC实现依赖注入的.但是操作的时候,有个很明显的问题没有解决,就是日志记录问题.如果手动添加,上百个上千个操作,每个操作都要写一遍WriteLog方法,工作量 ...

  8. 深入理解OOP(第一天):多态和继承(初期绑定和编译时多态)

    在本系列中,我们以CodeProject上比较火的OOP系列博客为主,进行OOP深入浅出展现. 无论作为软件设计的高手.或者菜鸟,对于架构设计而言,均需要多次重构.取舍,以有利于整个软件项目的健康构建 ...

  9. ASP.NET Web API 实现客户端Basic(基本)认证 之简单实现

    优点是逻辑简单明了.设置简单. 缺点显而易见,即使是BASE64后也是可见的明文,很容易被破解.非法利用,使用HTTPS是一个解决方案. 还有就是HTTP是无状态的,同一客户端每次都需要验证. 实现: ...

  10. redis连接数

    1.应用程序会发起多少个请求连接?1)对于php程序,以短连接为主.redis的连接数等于:所有web server接口并发请求数/redis分片的个数.2)对于java应用程序,一般使用JedisP ...