1500: [NOI2005]维修数列

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 12544  Solved: 3970
[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

打掉大boss了,一边抄一边写。

因为插入的数总共为4000000个,所以我们可以暴力地每个点直接插入,插入时,我们把新插入的元素建成一棵小树,再挂到大树上。

删除一样,转到那个节点,然后删除,但是这里要一个个删除,因为内存要循环利用。直接暴力删除,把节点放进一个栈中,用的时候拿出来,因为一共只有4000000个元素,绝对不会超时。

修改是打一个懒标记,tag[x]=c,然后pushdown即可。

翻转同上。

求和也是利用标记,sum[x]=sum[child[x][0]]+sum[child[x][1]]+key[x]

最后一个我们需要维护三个值,x节点表示的区间,lm[x]表示从x节点的区间最左端开始,最大的和是多少,rm[x]则是右边,mx[x]则表示这个节点的最大和

mx[x]=max(mx[child[x][0]],mx[child[x][1]],max(rm[child[x][0]],0)+key[x]+max(0,lm[child[x][1]])) 自己拿手画一下,想一下即可。

注意:pushdown要更新自己的sum[x],key[x],然后也要把自己的儿子节点child[x][0],child[x][1]也更新掉,翻转也是,我不知道为什么,这里查了两个小时才发现。

makesame的初始值要设成-1001,因为更新的数在[-1000,1000]之间。

#include<cstdio>
#include<cstring>
#include<stack>
#include<vector>
using namespace std;
#define N 1000010
#define inf -1101
stack<int> s;
int n,m,root,cnt;
int size[N],child[N][],lm[N],rm[N],key[N],mx[N],sum[N];
int tag1[N],tag2[N],fa[N],a[N];
void update(int x)
{
size[x]=size[child[x][]]+size[child[x][]]+;
sum[x]=sum[child[x][]]+sum[child[x][]]+key[x];
mx[x]=max(max(mx[child[x][]],mx[child[x][]]),
max(rm[child[x][]],)+key[x]+max(lm[child[x][]],));
lm[x]=max(lm[child[x][]],sum[child[x][]]+key[x]+max(,lm[child[x][]]));
rm[x]=max(rm[child[x][]],sum[child[x][]]+key[x]+max(,rm[child[x][]]));
}
void paint(int x)
{
tag2[x]^=;
swap(child[x][],child[x][]);
swap(lm[x],rm[x]);
}
void paint(int x,int c)
{
if(!x) return;
tag1[x]=key[x]=c; sum[x]=size[x]*c;
lm[x]=rm[x]=mx[x]=max(sum[x],key[x]);
}
void pushdown(int x)
{
if(tag1[x]!=inf)
{
paint(child[x][],tag1[x]);//更新儿子节点
paint(child[x][],tag1[x]);//更新儿子节点
tag1[x]=inf;
}
if(tag2[x])
{
tag2[x]^=;
if(child[x][]) paint(child[x][]);//更新儿子节点
if(child[x][]) paint(child[x][]);//更新儿子节点
}
}
void zig(int x)
{
int y=fa[x];
fa[x]=fa[y]; child[fa[x]][child[fa[x]][]==y]=x;
child[y][]=child[x][];
fa[child[x][]]=y;
child[x][]=y; fa[y]=x;
update(y); update(x);
}
void zag(int x)
{
int y=fa[x];
fa[x]=fa[y]; child[fa[x]][child[fa[x]][]==y]=x;
child[y][]=child[x][];
fa[child[x][]]=y;
child[x][]=y; fa[y]=x;
update(y); update(x);
}
void splay(int x,int t)
{
while(fa[x]!=t)
{
pushdown(x);
int y=fa[x],z=fa[y];
if(z==t)
{
x==child[y][]?zig(x):zag(x); break;
}
x==child[y][]?zig(x):zag(x);
x==child[z][]?zig(x):zag(x);
}
if(!t) root=x;
update(root);
}
int find(int x,int rank)
{
pushdown(x);
if(size[child[x][]]+==rank) return x;
if(size[child[x][]]>=rank) return find(child[x][],rank);
else return find(child[x][],rank-size[child[x][]]-);
}
int newnode()
{
int ret=;
if(!s.empty())
{
ret=s.top();
s.pop();
} else ret=++cnt;
return ret;
}
void build(int l,int r,int&x,int last)
{
int mid=(l+r)>>;
x=newnode(); key[x]=a[mid]; fa[x]=last;
if(l<=mid-) build(l,mid-,child[x][],x);
if(mid+<=r) build(mid+,r,child[x][],x);
update(x);
}
void insert(int l,int r)
{
int x=find(root,l),y=find(root,l+);//在l之后在l+1之前
splay(x,); splay(y,root);
for(int i=;i<=r-l+;i++) scanf("%d",&a[i]);
build(,r-l+,child[child[root][]][],child[root][]);
update(child[root][]); update(root);
}
void erase(int&x)
{
if(!x) return;
s.push(x);
fa[x]=; key[x]=size[x]=sum[x]=tag2[x]=;
tag1[x]=inf;
lm[x]=rm[x]=mx[x]=-(<<);
erase(child[x][]);
erase(child[x][]);
x=;
}
void del(int l,int r)
{
int x=find(root,l-),y=find(root,r+);
splay(x,); splay(y,root);
erase(child[child[root][]][]);
update(child[root][]); update(root);
}
void makesame(int l,int r,int c)
{
int x=find(root,l-),y=find(root,r+);
splay(x,); splay(y,root);
paint(child[child[root][]][],c);//更新儿子节点
update(child[root][]); update(root);
splay(child[child[root][]][],);
}
void reverse(int l,int r)
{
int x=find(root,l-),y=find(root,r+);
splay(x,); splay(y,root);
int b=child[child[root][]][];
if(b) paint(b);//更新儿子节点
splay(child[child[root][]][],);
}
void getsum(int l,int r)
{
int x=find(root,l-),y=find(root,r+);
splay(x,); splay(y,root);
printf("%d\n",sum[child[child[root][]][]]);
}
void maxsum()
{
printf("%d\n",mx[root]);
}
int main()
{
for(int i=;i<=N;i++) tag1[i]=inf;
scanf("%d%d",&n,&m);
for(int i=;i<=n+;i++)
{
scanf("%d",&a[i]);
}
mx[]=mx[n+]=lm[]=lm[n+]=rm[]=rm[n+]=mx[]=
lm[]=rm[]=a[]=a[n+]=-(<<);
int x; build(,n+,x,);
root=;
while(m--)
{
char s[]; int pos,tot,c; scanf("%s",s);
if(s[]=='S')
{
scanf("%d%d",&pos,&tot);
insert(pos+,pos+tot);
}
if(s[]=='L')
{
scanf("%d%d",&pos,&tot);
del(pos+,pos+tot);
}
if(s[]=='K')
{
scanf("%d%d%d",&pos,&tot,&c);
if(tot>) makesame(pos+,pos+tot,c);
}
if(s[]=='V')
{
scanf("%d%d",&pos,&tot);
reverse(pos+,pos+tot);
}
if(s[]=='T')
{
scanf("%d%d",&pos,&tot);
getsum(pos+,pos+tot);
}
if(s[]=='X')
{
maxsum();
}
}
return ;
}

bzoj1500的更多相关文章

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

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

  2. 【BZOJ1500】【NOI2005】维修数列(Splay)

    [BZOJ1500][NOI2005]维修数列(Splay) 题面 不想再看见这种毒瘤题,自己去BZOJ看 题解 Splay良心模板题 真的很简单 我一言不发 #include<iostream ...

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

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

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

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

  5. [您有新的未分配科技点] 无旋treap:从单点到区间(例题 BZOJ1500&NOI2005 维护数列 )

    1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MB Description Input 输入的第1 行包含两个数N 和M(M ≤20 ...

  6. 【BZOJ1500】维修数列(splay)

    题意: 输入的第1 行包含两个数N 和M(M ≤20 000),N 表示初始时数列中数的个数,M表示要进行的操作数目.第2行包含N个数字,描述初始时的数列.以下M行,每行一条命令,格式参见问题描述中的 ...

  7. [bzoj1500][NOI2005]维修数列_非旋转Treap

    维修数列 bzoj-1500 NOI-2005 题目大意:给定n个数,m个操作,支持:在指定位置插入一段数:删除一个数:区间修改:区间翻转.查询:区间和:全局最大子序列. 注释:$1\le n_{ma ...

  8. BZOJ1500[NOI2005]维修数列

    Description Input 输入的第1 行包含两个数N 和M(M ≤20 000),N 表示初始时数列中数的个数,M表示要进行的操作数目.第2行包含N个数字,描述初始时的数列.以下M行,每行一 ...

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

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

随机推荐

  1. doT js 模板引擎【初探】要优雅不要污

    js中拼接html,总是感觉不够优雅,本着要优雅不要污,决定尝试js模板引擎. JavaScript 模板引擎 JavaScript 模板引擎作为数据与界面分离工作中最重要一环,越来越受开发者关注. ...

  2. php函数获取真实客户端IP地址

    function getIPaddress(){ $IPaddress=''; if (isset($_SERVER)){ if (isset($_SERVER["HTTP_X_FORWAR ...

  3. .NET Core 2.0版本预计于2017年春季发布

    英文原文: NET Core 2.0 Planned for Spring 2017 微软项目经理 Immo Landwerth 公布了即将推出的 .NET Core 2.0 版本的细节,该版本预计于 ...

  4. Node.js之sails框架

    先开一坑,有空更新,记录最近钉钉项目上对node及sails框架的学习记录和理解

  5. sencha ext js 6 入门

    Sencha Ext JS号称是目前世界上最先进和最强大的.支持多平台多设备的JavaScript应用程序开发框架.首先看一下Ext JS的发展简史. 1 Ext JS发展简史 YUI-Ext的作者J ...

  6. C#Excel文件加密实现,支持xlsx、docx、pptx(C#\Net\Asp.Net)

    从此刻开始,我已封闭!概不接客! 像风一样的男人,像风一样的性格,无拘无束,不拘一格.那么问题来了,当风遇到沙,不一定你是风儿,我是沙儿的缠缠绵绵,.也许是漫天黄沙,飞粒走石.如果我们期望擒住这漫天的 ...

  7. Web Mercator Non-Conformal, Non-Mercator

    public static void XYtoGL(Coordinate coordinate) { double R = 6378137; coordinate.x = coordinate.x / ...

  8. iOS中的各种id

     网卡Mac地址会随机化 

  9. ListView之性能优化

    listview加载的核心是其adapter,本文通过减少adapter中创建.处理view的次数来提高listview加载的性能,总共分四个层次: 0.最原始的加载 1.利用convertView ...

  10. ASP.NET MVC 3 网站优化总结(三)Specify Vary: Accept-Encoding header

    继续进行 ASP.NET MVC 3 网站优化工作,使用 Google Page 检测发现提示 You should Specify Vary: Accept-Encoding header,The ...