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. mysql新建用户的方法

    新增 insert into mysql.user(Host,User,Password,ssl_cipher,x509_issuer,x509_subject) values("local ...

  2. BZOJ 1066 【SCOI2007】 蜥蜴

    Description 在一个r行c列的网格地图中有一些高度不同的石柱,一些石柱上站着一些蜥蜴,你的任务是让尽量多的蜥蜴逃到边界外. 每行每列中相邻石柱的距离为$1$,蜥蜴的跳跃距离是d,即蜥蜴可以跳 ...

  3. Exchange WebSerivce Usage

    //ExchangeService版本为2007SP1 ExchangeService service = new ExchangeService(ExchangeVersion.Exchange20 ...

  4. 转:用WCAT进行IIS压力测试

    Microsoft的Web容量分析工具(WCAT) 是测试你的客户-服务器网络配置的必备工具.这个工具在你的网络上对多种工作量的场景进行仿真,允许你确定你的网络和服务器的最佳配置.WCAT是专门为 评 ...

  5. ssh生成key不交互

    ssh-keygen -t rsa -f ~/.ssh/id_rsa -P ""   首次执行不交互 第二次再次执行会让输入y

  6. codevs1842 递归第一次

    难度等级:白银 1842 递归第一次 题目描述 Description 同学们在做题时常遇到这种函数 f(x)=5 (x>=0) f(x)=f(x+1)+f(x+2)+1 (x<0) 下面 ...

  7. QTableView 添加进度条

    记录一下QTableView添加进度条 例子很小,仅供学习 使用QItemDelegate做的实现 有自动更新进度 要在.pro文件里添加 CONFIG += c++ ProgressBarDeleg ...

  8. Windows 10 后台音频

    UWP版本的网易云音乐已经上架,虽然还不支持Windows Phone但是整体而言功能已经比较齐全了! 那么如何在Windows 10 UWP实现后台播放呢? 我之前是一直在做Windows Phon ...

  9. vue.js+boostrap最佳实践

    一.为什么要写这篇文章 最近忙里偷闲学了一下vue.js,同时也复习了一下boostrap,发现这两种东西如果同时运用到一起,可以发挥很强大的作用,boostrap优雅的样式和丰富的组件使得页面开发变 ...

  10. HDInsight - 1,简介

    最近工作需要,要看HDInsight部分,这里要做笔记.自然是官网资料最权威,所以内容都从这里搬过来:https://azure.microsoft.com/en-us/documentation/a ...