【bzoj5050】【bzoj九月月赛H】建造摩天楼

讲个笑话,这个题很休闲的。
大概是这样的,昨天看到这个题,第一眼星际把题目看反了然后感觉这是个傻逼题。
后来发现不对,这个修改一次的影响是很多的,可能导致一个数突然可以被改,也可能导致一个数不能被改。
大概就是一个不断拔高,最后拔得跟区间最大值一样高的过程。
后来开始想做法,感觉是不是可以维护一下最小的深度差,然后这个以内直接加,超过了就重构?
感觉似乎有点道理,写写写
自己想又感觉这个不太靠谱,可能会重构上天,然后就删了。
(flag++)
然后问了栋老师,栋老师一眼星际(跟我下午一样),然后忙着写集训队作业去了。
(天下神犇都神的相似,傻逼各有各的傻法)
后来问了lxl,挂张图:

他说的似乎有点道理,然而我没太弄懂他想表达什么
于是GG看其他题目去了。。。。。。
后来看了题解……

mdzz这不是我下午想的那个吗……?
然后想了下这个随机期望的复杂度,不太会,问了下栋老师,大概得到的结果是这样的:
首先分类讨论。

所以重构会在常数次,因此题解的复杂度是正确的~
然后谈一下维护:
由于变化会影响左右两边,所以维护一下左右两边的最小深度,每次在线段树上收集重构节点,大力计算就好。
#include<bits/stdc++.h>
#define lson (o<<1)
#define rson (o<<1|1)
const int N=;
const int inf=;
typedef long long ll;
inline int min(int x,int y){return x<y?x:y;}
using namespace std;
int n,m,a[N],exis[N],exl[N],exr[N],cnt,q[N];
ll sumv[N<<];int addv[N<<],minl[N<<],minr[N<<],size[N<<];
const int ch_top=4e7+;
typedef long long ll;
char ch[ch_top],*now_r=ch-,*now_w=ch-;
inline int read(){
while(*++now_r<'');
register int x=*now_r-'';
while(*++now_r>='')x=x*+*now_r-'';
return x;
}
inline void write(ll x){
static char st[];static int top;
while(st[++top]=''+x%,x/=);
while(*++now_w=st[top],--top);
*++now_w='\n';
}
inline void pushup(int o){
minr[o]=min(minr[lson],minr[rson]);
minl[o]=min(minl[lson],minl[rson]);
}
inline void puttag(int o,int v){
sumv[o]+=v*size[o];addv[o]+=v;minl[o]-=v;minr[o]-=v;
}
inline void pushdown(int o){
if(!addv[o])return;
puttag(lson,addv[o]);puttag(rson,addv[o]);
addv[o]=;
}
inline void build(int o,int l,int r){
if(l==r){
sumv[o]=a[l];size[o]=exis[l];minl[o]=exl[l];minr[o]=exr[l];
return;
}
int mid=(l+r)>>;
build(lson,l,mid);build(rson,mid+,r);
pushup(o);
sumv[o]=sumv[lson]+sumv[rson];
size[o]=size[lson]+size[rson];
}
inline void recalc(int o,int l,int r,int q){
if(l==r){
if(exis[l])a[l]+=addv[o];
addv[o]=;return;
}
pushdown(o);
int mid=(l+r)>>;
if(q<=mid)recalc(lson,l,mid,q);else recalc(rson,mid+,r,q);
}
inline ll querysum(int o,int l,int r,int ql,int qr){
if(ql<=l&&r<=qr)return sumv[o];
pushdown(o);
int mid=(l+r)>>;ll ans=;
if(ql<=mid)ans+=querysum(lson,l,mid,ql,qr);
if(qr>mid)ans+=querysum(rson,mid+,r,ql,qr);
return ans;
}
inline void change(int o,int l,int r,int ql,int qr){
if(ql<=l&&r<=qr)return void(puttag(o,));
pushdown(o);int mid=(l+r)>>;
if(ql<=mid)change(lson,l,mid,ql,qr);
if(qr>mid)change(rson,mid+,r,ql,qr);
pushup(o);
sumv[o]=sumv[lson]+sumv[rson];
}
inline void recexi(int o,int l,int r,int q){
if(l==r){
size[o]=exis[l];minl[o]=exl[l];minr[o]=exr[l];
return;
}
pushdown(o);int mid=(l+r)>>;
if(q<=mid)recexi(lson,l,mid,q);else recexi(rson,mid+,r,q);
pushup(o);
size[o]=size[lson]+size[rson];
}
inline void recml(int o,int l,int r,int q,int v){
if(l==r){minl[o]=v;return;}
pushdown(o);int mid=(l+r)>>;
if(q<=mid)recml(lson,l,mid,q,v);else recml(rson,mid+,r,q,v);
minl[o]=min(minl[lson],minl[rson]);
}
inline void recmr(int o,int l,int r,int q,int v){
if(l==r){minr[o]=v;return;}
pushdown(o);int mid=(l+r)>>;
if(q<=mid)recmr(lson,l,mid,q,v);else recmr(rson,mid+,r,q,v);
minr[o]=min(minr[lson],minr[rson]);
}
inline void dfs(int o,int l,int r){
if(minl[o]>&&minr[o]>=)return;
if(l==r){
if(minl[o]<=&&q[cnt]!=l)q[++cnt]=l;
if(minr[o]<)q[++cnt]=l+;
return;
}
pushdown(o);int mid=(l+r)>>;
dfs(lson,l,mid);dfs(rson,mid+,r);
}
inline bool flag(int i){return a[i-]>a[i];}
inline int calcl(int i){return exis[i]&&!exis[i-]?a[i-]-a[i]:inf;}
inline int calcr(int i){return exis[i]&&!exis[i+]?a[i+]-a[i]:inf;}
inline void calcexl(int i){if(i<||i>n)return;recml(,,n+,i,calcl(i));}
inline void calcexr(int i){if(i<||i>n)return;recmr(,,n+,i,calcr(i));}
inline void get(int i){if(i<||i>n)return;recalc(,,n+,i);}
inline void qwq(int i){
if(i<||i>n)return;get(i);get(i-);get(i+);exis[i]=flag(i);exl[i]=calcl(i);exr[i]=calcr(i);
recexi(,,n+,i);calcexl(i+);calcexr(i-);
}
int main(){
fread(ch,,ch_top,stdin);
n=read();m=read();
a[]=exl[]=exr[]=inf;a[n+]=exl[n+]=exr[n+]=inf;
for(int i=;i<=n;i++)a[i]=read(),exis[i]=flag(i);
for(int i=;i<=n;i++)exl[i]=calcl(i),exr[i]=calcr(i);
build(,,n+);
while(m--){
int opt=read(),x=read(),y=read();
if(opt==){
cnt=;change(,,n+,x,y);dfs(,,n+);
for(int i=;i<=cnt;i++)qwq(q[i]);
qwq(x);qwq(y+);
}
else{
write(querysum(,,n+,x,y));
}
}
fwrite(ch,,now_w-ch,stdout);
}
【bzoj5050】【bzoj九月月赛H】建造摩天楼的更多相关文章
- bzoj5050: 建造摩天楼
Description 属于小Q管辖的n座摩天楼从左往右排成一排,编号依次为1到n,第i座摩天楼的高度为h_i.小Q会进行m次以下两种 操作之一: 2 l r,询问h_l+h_{l+1}+...+h_ ...
- bzoj 1600 & Usaco 月赛 2008 建造栅栏 题解
[原题] 1600: [Usaco2008 Oct]建造栅栏 Time Limit: 5 Sec Memory Limit: 64 MB Submit: 785 Solved: 443 [Subm ...
- BZOJ 4070:[APIO2015]雅加达的摩天楼 最短路
4070: [Apio2015]雅加达的摩天楼 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 464 Solved: 164[Submit][Sta ...
- BZOJ 1600: [Usaco2008 Oct]建造栅栏( dp )
QAQ我没读过书...四边形都不会判定了 简单的dp.... --------------------------------------------------------------------- ...
- BZOJ 1600: [Usaco2008 Oct]建造栅栏
1600: [Usaco2008 Oct]建造栅栏 Time Limit: 5 Sec Memory Limit: 64 MB Description 勤奋的Farmer John想要建造一个四面的 ...
- 洛谷九月月赛II
题解:模拟 一旦不匹配就要break #include<iostream> #include<cstdio> #include<cstring> #include& ...
- bzoj 4070 [Apio2015]雅加达的摩天楼 Dijkstra+建图
[Apio2015]雅加达的摩天楼 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 644 Solved: 238[Submit][Status][D ...
- BZOJ 4070 [Apio2015]雅加达的摩天楼 ——分块 SPFA
挺有趣的分块的题目. 直接暴力建边SPFA貌似是$O(nm)$的. 然后考虑分块,$\sqrt n$一下用虚拟节点辅助连边, 以上的直接暴力连边即可. 然后卡卡时间,卡卡空间. 终于在UOJ上T掉辣. ...
- BZOJ 1600 [Usaco2008 Oct]建造栅栏:dp【前缀和优化】
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1600 题意: 给你一个长度为n的木板,让你把这个木板切割成四段(长度为整数),并且要求这四 ...
随机推荐
- bzoj3622-已经没有什么好害怕的的了
题意 给出两个长度为 \(n\) 的数列 \(a,b\) ,\(2n\) 个数都互不相同,求有多少种对应方式使得 \(a_i>b_i\) 的个数比 \(a_i<b_i\) 的个数恰好多 \ ...
- cdq分治学习
看了stdcall大佬的博客 传送门: http://www.cnblogs.com/mlystdcall/p/6219421.html 感觉cdq分治似乎很多时候都要用到归并的思想
- Unbuntu+nginx+mysql+php
1/准备 sudo su --切换到root 2/nginx安装 apt-get update apt-get install nginx 3/mysql 安装 apt-get install mys ...
- C++解析(15):二阶构造模式
0.目录 1.构造函数与半成品对象 2.二阶构造 3.小结 1.构造函数与半成品对象 关于构造函数: 类的构造函数用于对象的初始化 构造函数与类同名并且没有返回值 构造函数在对象定义时自动被调用 问题 ...
- CSS3制作3D水晶糖果按钮
本人仿照20个漂亮 CSS3 按钮效果及优秀的制作教程中的BonBon(Candy)Button实现了其棒棒糖果按钮,如下图所示: 在线演示地址见here. 使用完全使用CSS实现,无需JS.源码如下 ...
- mysql允许远程特定ip访问
1.登录 mysql -u root -p 之后输入密码进行登陆 2.权限设置及说明 2.1添加远程ip访问权限 GRANT ALL PRIVILEGES ON *.* TO 'root'@'10.1 ...
- linux 进程信号集合 sigset_t
sigset_t 号集及信号集操作函数:信号集被定义为一种数据类型: typedef struct { unsigned long sig[_NSIG_WORDS]: } sigset_t 信号集用来 ...
- PID控制算法的C语言实现八 变积分的PID控制算法C语言实现
变积分PID可以看成是积分分离的PID算法的更一般的形式.在普通的PID控制算法中,由于积分系数ki是常数,所以在整个控制过程中,积分增量是不变的.但是,系统对于积分项的要求是,系统偏差大时,积分作用 ...
- google的几道面试题
Question1:如何判断两个矩形是否有重叠区域? Answer(java): public boolean checkCross(Rectangle m, Rectangle n) { //求出m ...
- OpenCV---环境安装和初次使用
一:环境安装 pip3 install opencv-python #OpenCV模块,必须安装 pip3 install opencv-contrib-python #OpenCV扩展模块,选择安装 ...