[bzoj1500 维修数列](NOI2005) (splay)
真的是太弱了TAT。。。光是把代码码出来就花了3h。。还调了快1h才弄完T_T
号称考你会不会splay(当然通过条件是1h内AC。。吓傻)。。。
黄学长的题解:http://hzwer.com/2841.html
当然了蒟蒻的splay模板全部都是借(抄)鉴(袭)黄学长的。。。。
主要是注意任何对子树的修改一定要记得pushup一定要记得pushup一定要记得pushup。。。计算最大子段和和bzoj1858一个姿势,就是维护子树的最大子段和,从左和从右开始的最大子段和(当然了平衡树和线段树在这里还是有点区别的);
对于修改和翻转可以打懒标记。。。插入的话可以先把要插入的数建出一颗splay树,然后把新树和原来的树连起来;
为了避免判断区间边界啊什么的可以多建两个边界节点。。。
据说要加上垃圾回收,这样的话还要新开一个id数组用来存节点的编号,那样的话数组开到50w就可以了。。
如果不要垃圾回收的话,删除时把对应区间的子树与原来的树断开连接就好了,垃圾回收还得把要砍掉的子树里面的节点遍历重置。。。oj上没写数据范围让人感觉有点虚。。
#include<stdio.h>
#include<math.h>
#include<queue>
#include<iostream>
#include<algorithm>
using namespace std;
const int inf=;
const int maxn=;//题目又没说数据范围= =//经测试加了垃圾优化开50w就好了= =
struct zs{
int c[],fa,mx,lmx,rmx,val,size,sum;
bool change,rev;
}tree[maxn];
queue<int>q;
int a[maxn],id[maxn];
char iid[];
int i,j,k,n,m,cnt,tot,x,c,root;
void update(int now){
int l=tree[now].c[],r=tree[now].c[];
tree[now].size=tree[l].size+tree[r].size+;
tree[now].sum=tree[l].sum+tree[r].sum+tree[now].val;
tree[now].mx=max(tree[l].mx,tree[r].mx);
tree[now].mx=max(tree[now].mx,tree[l].rmx+tree[now].val+tree[r].lmx);
tree[now].lmx=max(tree[l].lmx,tree[l].sum+tree[now].val+tree[r].lmx);
tree[now].rmx=max(tree[r].rmx,tree[r].sum+tree[now].val+tree[l].rmx);
}
void deal(int now,bool change,bool rev,int val){
if(change){
tree[now].change=;
tree[now].val=val;
tree[now].sum=tree[now].size*val;
if(val>)
tree[now].mx=tree[now].lmx=tree[now].rmx=tree[now].sum;
else tree[now].mx=val,tree[now].lmx=tree[now].rmx=;
}
else//如果整段覆盖了翻不翻转都一样
if(rev){
tree[now].rev^=;
// printf("! %d %d %d\n",now,tree[now].c[0],tree[now].c[1]);
swap(tree[now].c[],tree[now].c[]);
swap(tree[now].lmx,tree[now].rmx);
}
}
void pushdown(int now){
int l=tree[now].c[],r=tree[now].c[];
if(l)deal(l,tree[now].change,tree[now].rev,tree[now].val);
if(r)deal(r,tree[now].change,tree[now].rev,tree[now].val);
tree[now].change=tree[now].rev=;
}
void build(int l,int r,int premid){
int now,mid,fa;
if(l>r)return;
mid=(l+r)>>;
fa=id[premid];now=id[mid];//本来id[i]==i,但因为垃圾回收。。。
if(mid<premid)tree[fa].c[]=now;else tree[fa].c[]=now;
tree[now].fa=fa;
tree[now].val=a[mid];
if(l<r){
build(l,mid-,mid);build(mid+,r,mid);
update(now);
}else{
// printf(" !%d %d\n",now,tree[now].val);
tree[now].size=;
tree[now].mx=tree[now].sum=tree[now].val;
if(tree[now].val>)tree[now].lmx=tree[now].rmx=tree[now].mx;
else tree[now].lmx=tree[now].rmx=;
}
// printf("! %d %d %d %d %d %d %d %d\n",now,tree[now].c[0],tree[now].c[1],tree[now].mx,tree[now].lmx,tree[now].rmx,tree[now].change,tree[now].rev);
}
void rotate(int now,int&root){
int fa,gfa,l,r;
fa=tree[now].fa;gfa=tree[fa].fa;
if(fa==root)root=now;
else tree[gfa].c[tree[gfa].c[]!=fa]=now;
if(tree[fa].c[]==now)l=;else l=;
r=l^;
tree[now].fa=gfa;tree[fa].fa=now;
tree[fa].c[l]=tree[now].c[r];tree[now].c[r]=fa;
tree[tree[fa].c[l]].fa=fa;
update(fa);update(now);
}
void splay(int now,int&root){
int fa,gfa;
while(now!=root){
fa=tree[now].fa;gfa=tree[fa].fa;
if(fa!=root)
if((tree[fa].c[]==now)^(tree[gfa].c[]==fa))rotate(now,root);
else rotate(fa,root);
rotate(now,root);
}
}
int find(int now,int k){
if(tree[now].change||tree[now].rev)pushdown(now);
int l=tree[now].c[];
if(k<=tree[l].size)return find(l,k);
else if(k==tree[l].size+)return now;
else return find(tree[now].c[],k-tree[l].size-);
}
void ins(int pos,int tot){
// int now=find(root,pos+1);
for(int i=;i<=tot;i++)scanf("%d",&a[i]);
for(int i=;i<=tot;i++)
if(!q.empty())id[i]=q.front(),q.pop();
else id[i]=++cnt;
int troot=id[(tot+)>>];
build(,tot,);
int x=find(root,pos+),y=find(root,pos+);
splay(x,root);splay(y,tree[x].c[]);
tree[y].c[]=troot;tree[troot].fa=y;
update(y);update(x);
// printf(" !%d\n",root);
}
void recover(int &now){
tree[now].fa=;
if(tree[now].c[])recover(tree[now].c[]);
if(tree[now].c[])recover(tree[now].c[]);
q.push(now);
tree[now].change=tree[now].rev=;
now=;
}
void del(int pos,int tot){
// int now=find(root,pos+1);
int x=find(root,pos),y=find(root,pos+tot+);
splay(x,root);splay(y,tree[x].c[]);
// int z=tree[y].c[0];tree[y].c[0]=0;
recover(tree[y].c[]);
update(y);update(x);//要更新祖先TAT
}
void bilibili(int pos,int tot,int change,int rev,int val){//处理修改&&翻转
int x=find(root,pos),y=find(root,pos+tot+);
splay(x,root);splay(y,tree[x].c[]);
// printf(" ! %d %d %d %d\n",tree[y].c[0],change,rev,val);
deal(tree[y].c[],change,rev,val);
update(y);update(x);//也要更新祖先TAT
// printf("! %d %d\n",y,x);
}
void getsum(int pos,int tot){
// int now=find(root,pos+1);
int x=find(root,pos),y=find(root,pos+tot+);
splay(x,root);splay(y,tree[x].c[]);
int z=tree[y].c[];
printf("%d\n",tree[z].sum);
}
int main(){
scanf("%d%d",&n,&m);
for(i=;i<=n;i++)scanf("%d",&a[i+]);
a[]=a[n+]=-inf;
for(i=;i<=n+;i++)id[i]=i;
tree[].mx=-inf;
root=(n++)>>;cnt=n+;
build(,n+,);
// for(i=1;i<=cnt;i++)printf(" %d %d %d %d %d %d\n",i,tree[i].c[0],tree[i].c[1],tree[i].mx,tree[i].lmx,tree[i].rmx);
for(i=;i<=m;i++){
scanf("%s",iid);
if(iid[]=='S')scanf("%d%d",&x,&tot),ins(x,tot);
else if(iid[]=='L')scanf("%d%d",&x,&tot),del(x,tot);
else if(iid[]=='K')scanf("%d%d%d",&x,&tot,&c),bilibili(x,tot,,,c);
else if(iid[]=='V')scanf("%d%d",&x,&tot),bilibili(x,tot,,,);
else if(iid[]=='T')scanf("%d%d",&x,&tot),getsum(x,tot);
else if(iid[]=='X')printf("%d\n",tree[root].mx);
//if(i==4)
// for(int i=1;i<=cnt;i++)printf(" %d %d %d %d %d %d %d %d\n",i,tree[i].c[0],tree[i].c[1],tree[i].mx,tree[i].lmx,tree[i].rmx,tree[i].change,tree[i].rev);
// return 233;
}
return ;
}
//连续几天写到的题目都没有数据范围是什么鬼。。。
//至于读入的时候少打一个%d这种错显然只有傻叉才会犯的啊TAT
//变量名太长的悲剧。。。把几个过程缩成一个结果代码量还是别人的两倍T_T
感觉自己还不会splay。。。果然是太弱了。。也许过一阵子再来写这题效果会好点= =
2015.8.1:4h
2015.8.6:2h : 注意删点的时候那些新建节点时没有初始化的信息都要初始化,
2016.1.6:4h+。。TAT。。太久没玩懒标记结果卡了一晚上+一节电脑课。。。要好好想一下什么时候更新本节点的信息。。
[bzoj1500 维修数列](NOI2005) (splay)的更多相关文章
- BZOJ 1500/Luogu 2042 - 维修数列 - [NOI2005][Splay]
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1500 题目链接:https://www.luogu.org/problemnew/sho ...
- 【BZOJ1500】【NOI2005】维修数列(Splay)
[BZOJ1500][NOI2005]维修数列(Splay) 题面 不想再看见这种毒瘤题,自己去BZOJ看 题解 Splay良心模板题 真的很简单 我一言不发 #include<iostream ...
- [bzoj1269]文本编辑器editor [bzoj1500]维修数列
1269: [AHOI2006]文本编辑器editor Time Limit: 10 Sec Memory Limit: 162 MB Submit: 2540 Solved: 923 [Submit ...
- 【BZOJ1500】维修数列(splay)
题意: 输入的第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 ...
- NOI2005维修数列(splay)
题目描述: Description 请写一个程序,要求维护一个数列,支持以下 6 种操作: 请注意,格式栏 中的下划线‘ _ ’表示实际输入文件中的空格 Input 输入的第1 行包含两个数N 和M( ...
- BZOJ1500 维修数列
AC通道:http://www.lydsy.com/JudgeOnline/problem.php?id=1500 [前言] 据说没打这题就相当于没打过Splay,这题简直就是让你内心崩溃的... 这 ...
- BZOJ 1500 维修数列【Splay】
注意:1,内存限制,所以需要回收删除的点 2,当前节点的左连续区间和最大值=max(左子树的左连续区间和最大值,左子树的总和+当节点的值+max(右子树的左连续区间和最大值,0)):右连续区间和最大值 ...
- p2042 维修数列(SPLAY)
终于yy过去了 撒花 调了一天半,出了无数锅 以下是出锅列表 rotate的时候要判断父亲节点的父亲节点是否存在 rotate的时候记得修改节点的父亲信息 pushdown的时候注意特判有无左右子树 ...
随机推荐
- iOS屏幕适配(尺寸适配)
屏幕尺寸适配:一 在.pch中加入以下代码,在定义每个尺寸值的时候都调用下边的宏 //以iphone7为例 定义 view相关的宽高宏#define IPHONEHIGHT(b) [UIScreen ...
- [机器学习系列] k-近邻算法(K–nearest neighbors)
C++ with Machine Learning -K–nearest neighbors 我本想写C++与人工智能,但是转念一想,人工智能范围太大了,我根本介绍不完也没能力介绍完,所以还是取了他的 ...
- mysql-5.7.17.msi安装
mysql-5.7.17.msi安装,跟着截图摩擦,一步一步是爪牙,是魔鬼的步伐 开始: 可以创建其他用户 我自己改了日志名
- bzoj 4898: [Apio2017]商旅
Description 在广阔的澳大利亚内陆地区长途跋涉后,你孤身一人带着一个背包来到了科巴.你被这个城市发达而美丽的市场所 深深吸引,决定定居于此,做一个商人.科巴有个集市,集市用从1到N的整数编号 ...
- Spring之DAO一
前面博客把bean.aop简单了解了一下,今天主要是了解Spring中DAO层,如果使用传统的JDBC时需要创建连接.打开.执行sql.关闭连接这一系列的步骤,Spring框架对JDBC进行了封装,我 ...
- centos7 卸载home 扩大root空间
=============================================== 2017/11/1_第1次修改 ccb_warlock == ...
- NOI 2005维护数列
题目描述 请写一个程序,要求维护一个数列,支持以下 6 种操作:(请注意,格式栏 中的下划线‘ _ ’表示实际输入文件中的空格) 输入输出格式 输入格式: 输入文件的第 1 行包含两个数 N 和 M, ...
- Linux设置PHP环境变量
区分 环境变量从时间上可分为临时性和永久性,这里只说明永久性的设置 操作 PHP 安装目录 找到PHP的安装目录:我这里是/phpstudy/server/php 其bin目录为:/phpstudy/ ...
- Web开发入门学习笔记
公司web项目终于要启动了,本以为django学习可以在实战中进行,结果最终使用了Drupal框架,好吧,那我们就PHP走起,买了本<细说PHP>,先跟着过一遍Web开发入门. HTTP协 ...
- Java---hashCode()和equals()
1.hashCode()和equals() API hashCode()和equals()都来自上帝类Object, 所有的类都会拥有这两个方法,特定时,复写它们. 它们是用来在同一类中做比较用的,尤 ...