BZOJ 1500 [NOI2005]维修数列 FHQ Treap
终于A了这题。。。这题还是很好。。。但是我太菜。。。重构了三遍qwq
FHQ Treap大法好!qwq。。。~~
Ins:直接拿输入造一棵树,把原来的树split成[1,pos],[pos+1,n],然后merge三棵树;
Del:把要删的区间split出来,merge他两边的树,记着要回收内存;
Mk-Same:把要改的区间split出来,打上标记,更新这棵树根的信息(见cover(x,v)),再merge回去;
Rev:把要翻转的区间split出来,打上标记,更新这棵树根的信息(见reverse(x)),再merge回去;
Get:把要求的区间split出来,输出sum[根],然后再merge回去;
Mx:输出mx[根]
维护最大子段和的思路跟线段树是类似的(链接)
顺便说一下自己当初犯的错误:
1.没有及时下放rev标记;重构x1
2.由于读取字符串用的getchar(),但是操作中有些串有负号。。。导致我读进来一堆负数。。。以后乖乖用字符串吧。。;重构x2,第三次重构才发现
3.stk内存池开小了。。。以为自己哪里递归写错了MLE+RE,后来看到网上有用queue的,就水了一波结果水过了。。后来又换成了数组stk。。;重构x3
代码:
#include<cstdio>
#include<iostream>
#include<queue>
#include<cstdlib>
#define R register int
#define ls(x) (ch[x][0])
#define rs(x) (ch[x][1])
using namespace std;
const int N=,Inf=;
inline int g() {
R ret=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
}
int n,m,tot,rt,top,sz[N],ch[N][],vl[N],mxl[N],mxr[N],mx[N],sum[N],tg[N],cov[N],dat[N],a[N>>],stk[N<<];
inline int max(int a,int b) {return a>b?a:b;}
inline int cre(int v) {R x; if(top) x=stk[top],--top; else x=++tot;
sz[x]=,ls(x)=rs(x)=tg[x]=; cov[x]=Inf,dat[x]=rand();
vl[x]=sum[x]=mx[x]=mxl[x]=mxr[x]=v; return x;
}
inline void upd(int x) {
if(!x) return ; sz[x]=sz[ls(x)]+sz[rs(x)]+; sum[x]=sum[ls(x)]+sum[rs(x)]+vl[x];
mx[x]=max(max(mxl[rs(x)],)+vl[x]+max(mxr[ls(x)],),max(mx[ls(x)],mx[rs(x)]));
mxl[x]=max(mxl[ls(x)],sum[ls(x)]+vl[x]+max(mxl[rs(x)],));
mxr[x]=max(mxr[rs(x)],sum[rs(x)]+vl[x]+max(mxr[ls(x)],));
}
inline void cover(int x,int v) {vl[x]=cov[x]=v,sum[x]=sz[x]*v,mxl[x]=mxr[x]=max(sum[x],),mx[x]=max(vl[x],sum[x]);}
inline void reverse(int x) {swap(ls(x),rs(x)),swap(mxl[x],mxr[x]),tg[x]^=;}
inline void spread(int x) {
if(tg[x]) {if(ls(x)) reverse(ls(x)); if(rs(x)) reverse(rs(x));}
if(cov[x]!=Inf) {if(ls(x)) cover(ls(x),cov[x]); if(rs(x)) cover(rs(x),cov[x]);} tg[x]=,cov[x]=Inf;
}
inline int build(int l,int r) {
if(l>r) return ; R md=l+r>>,v=a[md]; R x=cre(v);
ls(x)=build(l,md-),rs(x)=build(md+,r); upd(x); return x;
}
inline void split(int o,int rk,int& x,int& y) {
if(!o) {x=y=; return ;} spread(o);
if(rk>sz[ls(o)]) x=o,split(rs(x),rk-sz[ls(x)]-,rs(x),y);
else y=o,split(ls(y),rk,x,ls(y)); upd(o);
}
inline int merge(int x,int y) {
if(!x||!y) return x|y; spread(x),spread(y);
if(dat[x]<dat[y]) {rs(x)=merge(rs(x),y); upd(x); return x;}
else {ls(y)=merge(x,ls(y)); upd(y); return y;}
}
inline void del(int x) {if(!x) return; del(ls(x)),stk[++top]=x,del(rs(x));}
signed main() { R w,x,y,z; srand(); vl[]=mx[]=-Inf,sum[]=sz[]=tg[]=cov[]=;
n=g(),m=g(); for(R i=,x;i<=n;++i) a[i]=g(); rt=build(,n); while(m--) { register char ch;
while(!isalpha(ch=getchar())); if(ch=='M') {
getchar(),ch=getchar(); if(ch=='K') {
while(!isspace(ch=getchar())); R pos=g(),tot=g(),c=g(); split(rt,pos-,x,y),split(y,tot,y,z);
cover(y,c); rt=merge(x,merge(y,z));
} else if(ch=='X') {printf("%d\n",mx[rt]); while(!isspace(ch=getchar()));}
} else if(ch=='I') {R pos=g(),tot=g(); for(R i=,x;i<=tot;++i) a[i]=g(); z=build(,tot);split(rt,pos,x,y),rt=merge(merge(x,z),y);}
else if(ch=='D') {while(!isspace(ch=getchar())); R pos=g(),tot=g(); split(rt,pos-,x,y),split(y,tot,y,z),rt=merge(x,z); del(y);}
else if(ch=='R') {while(!isspace(ch=getchar())); R pos=g(),tot=g(); split(rt,pos-,x,y),split(y,tot,y,z); reverse(y); rt=merge(x,merge(y,z));}
else if(ch=='G') {while(!isspace(ch=getchar())); R pos=g(),tot=g(); split(rt,pos-,x,y),split(y,tot,y,z); printf("%d\n",sum[y]); rt=merge(x,merge(y,z));}
}
}
2019.05.11
BZOJ 1500 [NOI2005]维修数列 FHQ Treap的更多相关文章
- bzoj 1500: [NOI2005]维修数列 splay
1500: [NOI2005]维修数列 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 6556 Solved: 1963[Submit][Status ...
- BZOJ 1500: [NOI2005]维修数列 (splay tree)
1500: [NOI2005]维修数列 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 4229 Solved: 1283[Submit][Status ...
- [BZOJ 1500] [NOI2005] 维修数列
题目链接:BZOJ - 1500 题目分析 我要先说一下,这道题我写了一晚上,然后Debug了一整个白天..........再一次被自己的蒟蒻程度震惊= = 这道题是传说中的Splay维护数列的Bos ...
- 【BZOJ】1500: [NOI2005]维修数列
[算法]splay [题解]数据结构 感谢Occult的模板>_<:HYSBZ 1500 维修数列 #include<cstdio> #include<cctype> ...
- 1500: [NOI2005]维修数列
Description Input 输入的第1 行包含两个数N 和M(M ≤20 000),N 表示初始时数列中数的个数,M表示要进行的操作数目.第2行包含N个数字,描述初始时的数列.以下M行,每行一 ...
- 【BZOJ】1500: [NOI2005]维修数列(splay+变态题)
http://www.lydsy.com/JudgeOnline/problem.php?id=1500 模板不打熟你确定考场上调试得出来? 首先有非常多的坑点...我遇到的第一个就是,如何pushu ...
- bzoj千题计划221:bzoj1500: [NOI2005]维修数列(fhq treap)
http://www.lydsy.com/JudgeOnline/problem.php?id=1500 1.覆盖标记用INF表示无覆盖标记,要求可能用0覆盖 2.代表空节点的0号节点和首尾的两个虚拟 ...
- [NOI2005] 维修数列
1500: [NOI2005]维修数列 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 8397 Solved: 2530 Description In ...
- [BZOJ1500][NOI2005]维修数列---解题报告
Portal Gun:[BZOJ1500][NOI2005]维修数列 有一段时间没写博客了,最近在刚数据结构......各种板子背得简直要起飞,题目也是一大堆做不完,这里就挑一道平衡树的题来写写好了 ...
随机推荐
- unity渲染层级关系小结
http://blog.csdn.net/meegomeego/article/details/42060389 最近连续遇到了几个绘制图像之间相互遮挡关系不正确的问题,网上查找的信息比较凌乱,所以这 ...
- poj 2390 Bank Interest(计算本利和)
一.Description Farmer John made a profit last year! He would like to invest it well but wonders how m ...
- NBIOT回答
该部分分享的是物联网各垂直应用领域里,NB-IoT技术的部署,看看适合NB-IoT技术的垂直应用场景有哪些?垂直应用服务商又该如何部署? 1 NB-IoT适合的垂直应用场景有哪些? 2 NB-IoT垂 ...
- Code:目录
ylbtech-Code:目录 1.返回顶部 1. https://github.com/ 2.返回顶部 1. https://gitee.com 2. 3.返回顶部 4.返回顶部 5.返 ...
- mouseenter与mouseover的区别
mouseover 事件:只有在鼠标指针穿过被选元素时,才会触发. mouseover 事件:鼠标指针穿过任何子元素,都会触发. 请看例子的演示.
- Java探索之旅(3)——选择与循环
1.选择结构与输出 ❶Switch语句: Switch表达式必须算出 char,byte,short,int类型数值之一,总是括号括住:Value1----ValueN,对应有相同数据类型且为常量或者 ...
- new/delete 和malloc/free 的区别
new/delete 和malloc/free 的区别 一.基本概念malloc/free:1.函数原型及说明: void *malloc(long NumBytes):该函数分配了NumB ...
- [poj2104]可持久化线段树入门题(主席树)
解题关键:离线求区间第k小,主席树的经典裸题: 对主席树的理解:主席树维护的是一段序列中某个数字出现的次数,所以需要预先离散化,最好使用vector的erase和unique函数,很方便:如果求整段序 ...
- nodejs assert 模块
assert模块是Node的内置模块,主要用于断定.如果表达式不符合预期,就抛出一个错误.该模块提供11个方法,但只有少数几个是常用的. 1.assert() assert(value[, messa ...
- p4503&bzoj3555 企鹅QQ
传送门(洛谷) 传送门(bzoj) 题目 PenguinQQ是中国最大.最具影响力的SNS(Social Networking Services)网站,以实名制为基础,为用户提供日志.群.即时通讯.相 ...