BZOJ3642 : [CEOI 2014] Cake
rank[i]表示第i美味的是哪块
left[i]表示在k左边美味度为i的是哪块
right[i]表示在k右边美味度为i的是哪块
用3棵线段树分别维护d序列的区间最大值、left序列的区间最大值、right序列的区间最小值
修改:
把第x块改成第y美味
把第y+1到第9美味的全部后移一位
然后把第x块美味度改成最大值+1
然后把第y-1到第1美味的美味度依次改成最大
查询:
设x到k这一段中美味度的最大值为y
求出k另一侧最靠近k的且美味度大于y的位置z
答案为|z-x|-1
时间复杂度为每次操作$O(\log n)$
一开始写了普通线段树被卡常数了TAT(最近怎么都被卡常数)
然后去学了zkw线段树后改成zkw线段树就过了,效果拔群,又快又短。
zkw线段树真是个好东西。
#include<cstdio>
const int N=250010,inf=~0U>>1;
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
inline int min(int a,int b){return a<b?a:b;}
inline int max(int a,int b){return a>b?a:b;}
inline void Max(int&a,int b){if(a<b)a=b;}
inline void Min(int&a,int b){if(a>b)a=b;}
int n,m,t,q,k,i,x,y,rank[N],d[N],left[5250010],right[5250010],M1,M2,ta_v[530000],tl_v[16777220],tr_v[16777220],que[500010][3];
char ch;
inline void ta_change(int x,int y){for(ta_v[x+=M1]=y,x>>=1;x;x>>=1)ta_v[x]=max(ta_v[x<<1],ta_v[x<<1|1]);}
inline int ta_ask(int x,int y){
int t=-inf;
for(x+=M1-1,y+=M1+1;x^y^1;x>>=1,y>>=1){
if(~x&1)Max(t,ta_v[x^1]);
if(y&1)Max(t,ta_v[y^1]);
}
return t;
}
inline void tl_change(int x,int y){for(tl_v[x+=M2]=y,x>>=1;x;x>>=1)tl_v[x]=max(tl_v[x<<1],tl_v[x<<1|1]);}
inline int tl_ask(int x,int y){
int t=-inf;
for(x+=M2-1,y+=M2+1;x^y^1;x>>=1,y>>=1){
if(~x&1)Max(t,tl_v[x^1]);
if(y&1)Max(t,tl_v[y^1]);
}
return t;
}
inline void tr_change(int x,int y){for(tr_v[x+=M2]=y,x>>=1;x;x>>=1)tr_v[x]=min(tr_v[x<<1],tr_v[x<<1|1]);}
inline int tr_ask(int x,int y){
int t=inf;
for(x+=M2-1,y+=M2+1;x^y^1;x>>=1,y>>=1){
if(~x&1)Min(t,tr_v[x^1]);
if(y&1)Min(t,tr_v[y^1]);
}
return t;
}
inline void modify(int x,int y){//d[x]改成y
if(x<k)tl_change(d[x],0);
if(x>k)tr_change(d[x],n+1);
ta_change(x,d[x]=y);
if(x<k)tl_change(y,x);
if(x>k)tr_change(y,x);
}
inline void change(int x,int y){//把第x块蛋糕设置成第y美味
int i,j=0;
for(i=1;i<=10;i++)if(rank[i]==x)j=i;
for(i=j?j:10;i>y;i--)rank[i]=rank[i-1];
for(modify(rank[y]=x,++t),i=y-1;i;i--)modify(rank[i],++t);
}
inline int ask(int x){
if(x==k)return 0;
return x<k?tr_ask(ta_ask(x,k-1)+1,m)-x-1:x-tl_ask(ta_ask(k+1,x)+1,m)-1;
}
int main(){
read(n),read(k);
for(i=1;i<=n;i++)read(d[i]),rank[n-d[i]+1]=i;
read(q);
m=(t=n)+1;
for(i=1;i<=q;i++){
while(!(((ch=getchar())=='E')||(ch=='F')));
if(ch=='E')read(que[i][1]),read(que[i][2]),m+=que[i][2];else que[i][0]=1,read(que[i][1]);
}
for(i=1;i<=m;i++)right[i]=n+1;
for(i=1;i<k;i++)left[d[i]]=i;
for(i=k+1;i<=n;i++)right[d[i]]=i;
for(M1=1;M1<n+2;M1<<=1);
for(i=1;i<=n;i++)ta_v[i+M1]=d[i];
for(i=M1-1;i;i--)ta_v[i]=max(ta_v[i<<1],ta_v[i<<1|1]);
for(M2=1;M2<m+2;M2<<=1);
for(i=1;i<=m;i++)tl_v[i+M2]=left[i],tr_v[i+M2]=right[i];
for(i=M2-1;i;i--)tl_v[i]=max(tl_v[i<<1],tl_v[i<<1|1]),tr_v[i]=min(tr_v[i<<1],tr_v[i<<1|1]);
for(i=1;i<=q;i++)if(!que[i][0])change(que[i][1],que[i][2]);else printf("%d\n",ask(que[i][1]));
return 0;
}
BZOJ3642 : [CEOI 2014] Cake的更多相关文章
- CEOI 2014 wall (最短路)
描述:给定一个网格图,每个区间可能会有城市,求在边上建墙使无法从外边到达所有城市切所有城市必须联通 n,m<=400 首先对于30%的数据,n,m<=10我们可以考虑用数位dp来解决这个问 ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
- 看完SQL Server 2014 Q/A答疑集锦:想不升级都难!
看完SQL Server 2014 Q/A答疑集锦:想不升级都难! 转载自:http://mp.weixin.qq.com/s/5rZCgnMKmJqeC7hbe4CZ_g 本期嘉宾为微软技术中心技术 ...
- Windows 7上执行Cake 报错原因是Powershell 版本问题
在Windows 7 SP1 电脑上执行Cake的的例子 http://cakebuild.net/docs/tutorials/getting-started ,运行./Build.ps1 报下面的 ...
- SQL Server 2014 新特性——内存数据库
SQL Server 2014 新特性——内存数据库 目录 SQL Server 2014 新特性——内存数据库 简介: 设计目的和原因: 专业名词 In-Memory OLTP不同之处 内存优化表 ...
- 2014年暑假c#学习笔记目录
2014年暑假c#学习笔记 一.C#编程基础 1. c#编程基础之枚举 2. c#编程基础之函数可变参数 3. c#编程基础之字符串基础 4. c#编程基础之字符串函数 5.c#编程基础之ref.ou ...
- SQL Server 2014聚集列存储索引
转发请注明引用和原文博客(http://www.cnblogs.com/wenBlog) 简介 之前已经写过两篇介绍列存储索引的文章,但是只有非聚集列存储索引,今天再来简单介绍一下聚集的列存储索引,也 ...
- Web 前端攻防(2014版)
在百度 FEX 团队时写的.有次让我写点前端开发中注意的安全规范,结果写着写着就跑题了,写了一堆纯前端实现的攻击方式...当然还有防范措施. 文章就懒得粘过来了~ 直接贴个地址算了:) http:// ...
- 谈谈我的微软特约稿:《SQL Server 2014 新特性:IO资源调控》
一.本文所涉及的内容(Contents) 本文所涉及的内容(Contents) 背景(Contexts) 撰写经历(Experience) 特约稿正文(Content-body) 第一部分:生活中资源 ...
随机推荐
- [BZOJ1220][POJ1091][HNOI2002]跳蚤
[BZOJ1220][POJ1091][HNOI2002]跳蚤 试题描述 Z城市居住着很多只跳蚤.在Z城市周六生活频道有一个娱乐节目.一只跳蚤将被请上一个高空钢丝的正中央.钢丝很长,可以看作是无限长. ...
- Bmob用户管理操作
注册用户 BmobUser bu = new BmobUser(); bu.setUsername("sendi"); bu.setPassword("123456&qu ...
- (0,eval)('this')与eval的区别
看doT源码的时候,看到了这么一句代码: global = (function(){ return this || (0,eval)('this'); }()); global.doT = doT; ...
- pro git 使用积累
http://www.zhihu.com/question/20070065 git相关问题的收集 Git 是 Linux 之父 Linus Trovalds,为管理 Linux 内核代码而建立的,被 ...
- 【leetcode】3Sum Closest
3Sum Closest Given an array S of n integers, find three integers in S such that the sum is closest t ...
- Windows命令行提取日期时间
参考: http://elicecn.blog.163.com/blog/static/174017473200931910320556/ SET str="%date:~0,4%%date ...
- XShell上传下载命令
参考:https://www.centos.bz/2012/12/xshell-securecrtrz-sz-upload-download/ 上传文件时,执行rz就会弹出文件选择对话框来选择文件.下 ...
- maven项目,导入的jar包,没有包含在pom文件中,install失败
[INFO] BUILD FAILURE[INFO] ------------------------------------------------------------------------[ ...
- 父页面刷新 保持iframe页面url不变
思路:点击父页面时写cookies-->刷新时从cookies中奖内容读取出来. 本文转自:http://blog.163.com/sdolove@126/blog/static/1146378 ...
- Java for LeetCode 033 Search in Rotated Sorted Array
Suppose a sorted array is rotated at some pivot unknown to you beforehand. (i.e., 0 1 2 4 5 6 7 migh ...