维修数列 bzoj-1500 NOI-2005

题目大意:给定n个数,m个操作,支持:在指定位置插入一段数;删除一个数;区间修改;区间翻转。查询:区间和;全局最大子序列。

注释:$1\le n_{max} \le 5\cdot 10^5$,$1\le m \le 2\cdot 10^4$。

想法:据说是... ...最GB的平衡树裸题,如果把这题切了,话说所有的平衡树都不怕了。

具体的,

插入和删除对于非旋转Treap来讲都是基本操作;

区间修改的话我们就先删除然后插入即可;

区间反转就打标记;

区间和就是单点维护子树和;

全局最大子序列的话像小白逛公园一样维护一下从左右开始的最长子序列。

最后,附上丑陋的代码... ...

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 502333
#define MAXN 4002333
using namespace std;
char c;int v;
inline void read(int &x)
{
x=0;c=getchar();v=1;
while(c<'0'||c>'9')
{
if(c=='-')v=-1;
c=getchar();
}
while(c>='0'&&c<='9')x=x*10+c-'0',c=getchar();
x*=v;
}
int n,m,root;
struct Node
{
int ls,rs,size,key,val,sum,maxx,lmax,rmax,f;
bool turn;
}a[N];
struct par{int x,y;};
int q[MAXN],l=1,r,cnt,temp[N];
inline void clear(int x)
{
a[x].ls=a[x].rs=a[x].size=a[x].key=a[x].val=0;
a[x].maxx=a[x].lmax=a[x].rmax=a[x].turn=0;
a[x].f=2333;
}
inline void update(int x)
{
int ls=a[x].ls,rs=a[x].rs;
a[x].size=1;a[x].maxx=a[x].sum=a[x].val;
if(ls)a[x].size+=a[ls].size,a[x].sum+=a[ls].sum,a[x].maxx=max(a[x].maxx,a[ls].maxx);
if(rs)a[x].size+=a[rs].size,a[x].sum+=a[rs].sum,a[x].maxx=max(a[x].maxx,a[rs].maxx);
a[x].maxx=max(a[x].maxx,a[ls].rmax+a[x].val+a[rs].lmax);
a[x].lmax=max(a[ls].lmax,a[ls].sum+a[x].val+a[rs].lmax);
a[x].rmax=max(a[rs].rmax,a[rs].sum+a[x].val+a[ls].rmax);
}
inline void pushdown(int x)
{
if(!x)return;
int ls=a[x].ls,rs=a[x].rs;
if(a[x].turn)
{
if(ls)
{
a[ls].turn^=1;
swap(a[ls].lmax,a[ls].rmax);
}
if(rs)
{
a[rs].turn^=1;
swap(a[rs].lmax,a[rs].rmax);
}
swap(a[x].ls,a[x].rs);
a[x].turn=0;
}
ls=a[x].ls,rs=a[x].rs;
if(a[x].f!=2333)
{
if(ls)
{
a[ls].f=a[ls].val=a[x].f;a[ls].sum=a[ls].size*a[ls].f;
if(a[ls].f>0)a[ls].maxx=a[ls].lmax=a[ls].rmax=a[ls].sum;
else a[ls].maxx=a[ls].f,a[ls].lmax=a[ls].rmax=0;
}
if(rs)
{
a[rs].f=a[rs].val=a[x].f;a[rs].sum=a[rs].size*a[rs].f;
if(a[rs].f>0)a[rs].maxx=a[rs].lmax=a[rs].rmax=a[rs].sum;
else a[rs].maxx=a[rs].f,a[rs].lmax=a[rs].rmax=0;
}
a[x].f=2333;
}
}
int merge(int x,int y)
{
if(!x||!y)return x|y;
pushdown(x),pushdown(y);
if(a[x].key>a[y].key)
{
a[x].rs=merge(a[x].rs,y);update(x);
return x;
}
a[y].ls=merge(x,a[y].ls);update(y);
return y;
}
par split(int x,int k)
{
if(!k)return (par){0,x};
pushdown(x);
int ls=a[x].ls,rs=a[x].rs;
if(k==a[ls].size)
{
a[x].ls=0;update(x);
return (par){ls,x};
}
if(k==a[ls].size+1)
{
a[x].rs=0;update(x);
return (par){x,rs};
}
if(k<a[ls].size)
{
par t=split(ls,k);
a[x].ls=t.y;update(x);
return (par){t.x,x};
}
par t=split(rs,k-a[ls].size-1);
a[x].rs=t.x;update(x);
return (par){x,t.y};
}
inline int makenew(int x)
{
int p;
if(l<=r)p=q[l++];
else p=++cnt;
a[p].size=1,a[p].maxx=a[p].sum=a[p].val=x;
a[p].lmax=a[p].rmax=max(x,0);
a[p].turn=0,a[p].f=2333;
a[p].key=rand();
return p;
}
int build(int L,int R)
{
if(L==R)return makenew(temp[L]);
int mid=(L+R)>>1;
return merge(build(L,mid),build(mid+1,R));
}
void del(int x)
{
if(!x)return;
del(a[x].ls);
del(a[x].rs);
clear(x);
q[++r]=x;
}
int main()
{
read(n),read(m);
for(int i=1;i<=n;i++)read(temp[i]);
root=build(1,n);
char flag[10];
for(int x,num,val,j,i=1;i<=m;i++)
{
scanf("%s",flag);
if(flag[0]=='I')
{
read(x);read(num);
for(j=1;j<=num;j++)read(temp[j]);
par t=split(root,x);
root=merge(merge(t.x,build(1,num)),t.y);
}
else if(flag[0]=='D')
{
read(x);read(num);x--;
par t1=split(root,x),t2=split(t1.y,num);
del(t2.x);
root=merge(t1.x,t2.y);
}
else if(flag[0]=='R')
{
read(x);read(num);x--;
par t1=split(root,x),t2=split(t1.y,num);
a[t2.x].turn^=1;swap(a[t2.x].lmax,a[t2.x].rmax);
root=merge(t1.x,merge(t2.x,t2.y));
}
else if(flag[0]=='G')
{
read(x);read(num);x--;
par t1=split(root,x),t2=split(t1.y,num);
printf("%d\n",a[t2.x].sum);
root=merge(t1.x,merge(t2.x,t2.y));
}
else if(flag[2]=='X')printf("%d\n",a[root].maxx);
else
{
read(x);read(num);read(val);x--;
par t1=split(root,x),t2=split(t1.y,num);
a[t2.x].f=a[t2.x].val=val;a[t2.x].sum=val*a[t2.x].size;
if(val>0)a[t2.x].maxx=a[t2.x].lmax=a[t2.x].rmax=a[t2.x].sum;
else a[t2.x].maxx=val,a[t2.x].lmax=a[t2.x].rmax=0;
root=merge(t1.x,merge(t2.x,t2.y));
}
}
return 0;
}

小结:裸题就是开心... ...

P.S.:非旋转Treap就是比splay优越!

[bzoj1500][NOI2005]维修数列_非旋转Treap的更多相关文章

  1. 2018.08.06 bzoj1500: [NOI2005]维修数列(非旋treap)

    传送门 平衡树好题. 我仍然是用的fhqtreap,感觉速度还行. 维护也比线段树splay什么的写起来简单. %%%非旋treap大法好. 代码: #include<bits/stdc++.h ...

  2. [BZOJ1500][NOI2005]维修数列---解题报告

    Portal Gun:[BZOJ1500][NOI2005]维修数列 有一段时间没写博客了,最近在刚数据结构......各种板子背得简直要起飞,题目也是一大堆做不完,这里就挑一道平衡树的题来写写好了 ...

  3. [BZOJ1500][NOI2005]维修数列 解题报告 Splay

    Portal Gun:[BZOJ1500][NOI2005]维修数列 有一段时间没写博客了,最近在刚数据结构......各种板子背得简直要起飞,题目也是一大堆做不完,这里就挑一道平衡树的题来写写好了 ...

  4. BZOJ1500[NOI2005]维修数列——非旋转treap

    题目描述 请写一个程序,要求维护一个数列,支持以下 6 种操作: 请注意,格式栏 中的下划线‘ _ ’表示实际输入文件中的空格 输入 输入的第1 行包含两个数N 和M(M ≤20 000),N 表示初 ...

  5. bzoj千题计划221:bzoj1500: [NOI2005]维修数列(fhq treap)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1500 1.覆盖标记用INF表示无覆盖标记,要求可能用0覆盖 2.代表空节点的0号节点和首尾的两个虚拟 ...

  6. [bzoj3173]最长上升子序列_非旋转Treap

    最长上升子序列 bzoj-3173 题目大意:有1-n,n个数,第i次操作是将i加入到原有序列中制定的位置,后查询当前序列中最长上升子序列长度. 注释:1<=n<=10,000,开始序列为 ...

  7. [bzoj4864][BeiJing2017Wc]神秘物质_非旋转Treap

    神秘物质 bzoj-4864 BeiJing-2017-Wc 题目大意:给定一个长度为n的序列,支持插入,将相邻两个元素合并并在该位置生成一个指定权值的元素:查询:区间内的任意一段子区间的最大值减最小 ...

  8. [bzoj1552\bzoj2506][Cqoi2014]robotic sort 排序机械臂_非旋转Treap

    robotic sort 排序机械臂 bzoj-1552 bzoj-2506 Cqoi-2014 题目大意:给定一个序列,让你从1到n,每次将[1,p[i]]这段区间反转,p[i]表示整个物品权值第i ...

  9. [bzoj1861][Zjoi2006]Book 书架_非旋转Treap

    Book 书架 bzoj-1861 Zjoi-2006 题目大意:给你一个序列,支持:将指定编号的元素抽出,放到序列顶(底):将指定编号元素左右篡位:查询指定编号元素位置:查询指定数量位置元素编号. ...

随机推荐

  1. JAVA8与lambda表达式

    一.lambda表达式的来源 Lambda 表达式(拉姆达表达式)来源于lambda演算是Alonzo Church给出的一套图灵机等价的形式计算系统.lambda演算系统以函数和变量为基础,也可以进 ...

  2. Quartz实现执行任务记录数据库,方便计算任务的执行次数以及成功次数

    任务执行实体 /** * 任务执行情况详情 */ public class JobExecuteDetail implements Serializable{ /** * */ private sta ...

  3. 强迫症!一行代码拿到url特定query的值

    简单的说一下背景,看到小伙伴给我发的项目中有一段获取当前url特定query值的代码,本着能写1行代码就不写5行代码的原则,我把这个获取方法给改了一下 之前的代码如下: const queryArr ...

  4. POJ 2286 The Rotation Game IDA*

    (再一次感谢学长幻灯片) ID A* 随便自己yy了一下. 额嗯 思路什么的都没有问题 就是改不对.. 无奈地删代码...边删边交. 删啊删 哎呦 AC了 ... ... ... 找删的那一段 . o ...

  5. Django html页面 'ascii' codec can't encode characters in position 8-10: ordinal not

    用Django开发的页面,之前用的是python3.X,后来又换成python2.X后各种报错,编码问题,于是在所有python文件开头加了编码:#coding=utf-8 但是后来发现,有些文件加了 ...

  6. shell 实用命令学习

    查找文件 -iname 大小写不敏感 “*.log” .log后缀的文件 -type d 文件类型为目录的 find ./ -name 'index.html' 查找当前目录,及其子目录下文件

  7. ios -使用NSLayoutConstraint实现多个view等宽等高等间距

    @interface ViewController () { UIView *firstView; UIView *secondView; UIView *thirdView; } @end @imp ...

  8. 2013款MacBook Air装Windows7单系统

    经过两天的摸索,查找无数资料终于把2013款的MacBook Air装上了WIN 7,虽然网上有很多的资料但是都不是我想要的,第一个我的是2013款的MacBook Air,跟原来2012 11款Ma ...

  9. (1)dotnet开源电商系统-brnshop&brnMall 和老外开发的nopCommerce(dotnet两套电商来PK--第一篇)

    一直想做电商软件,但是实在不想学PHP了,所以前后关注了这两个开源电商系统.一个是国人出品的,一个据说是俄罗斯人写得(不知道对不对).目前两个开源软件都在学习了解中,以下的博文可能会涉及到这两套系统, ...

  10. OpenWRT 常用软件安装

    root@Jack:/tmp/opkg-lists# opkg--help opkg:unrecognized option `--help' opkgmust have one sub-comman ...