讲个笑话,这个题很休闲的。

大概是这样的,昨天看到这个题,第一眼星际把题目看反了然后感觉这是个傻逼题。

后来发现不对,这个修改一次的影响是很多的,可能导致一个数突然可以被改,也可能导致一个数不能被改。

大概就是一个不断拔高,最后拔得跟区间最大值一样高的过程。

后来开始想做法,感觉是不是可以维护一下最小的深度差,然后这个以内直接加,超过了就重构?

感觉似乎有点道理,写写写

自己想又感觉这个不太靠谱,可能会重构上天,然后就删了。

(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】建造摩天楼的更多相关文章

  1. bzoj5050: 建造摩天楼

    Description 属于小Q管辖的n座摩天楼从左往右排成一排,编号依次为1到n,第i座摩天楼的高度为h_i.小Q会进行m次以下两种 操作之一: 2 l r,询问h_l+h_{l+1}+...+h_ ...

  2. bzoj 1600 &amp; Usaco 月赛 2008 建造栅栏 题解

    [原题] 1600: [Usaco2008 Oct]建造栅栏 Time Limit: 5 Sec  Memory Limit: 64 MB Submit: 785  Solved: 443 [Subm ...

  3. BZOJ 4070:[APIO2015]雅加达的摩天楼 最短路

    4070: [Apio2015]雅加达的摩天楼 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 464  Solved: 164[Submit][Sta ...

  4. BZOJ 1600: [Usaco2008 Oct]建造栅栏( dp )

    QAQ我没读过书...四边形都不会判定了 简单的dp.... --------------------------------------------------------------------- ...

  5. BZOJ 1600: [Usaco2008 Oct]建造栅栏

    1600: [Usaco2008 Oct]建造栅栏 Time Limit: 5 Sec  Memory Limit: 64 MB Description 勤奋的Farmer John想要建造一个四面的 ...

  6. 洛谷九月月赛II

    题解:模拟 一旦不匹配就要break #include<iostream> #include<cstdio> #include<cstring> #include& ...

  7. bzoj 4070 [Apio2015]雅加达的摩天楼 Dijkstra+建图

    [Apio2015]雅加达的摩天楼 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 644  Solved: 238[Submit][Status][D ...

  8. BZOJ 4070 [Apio2015]雅加达的摩天楼 ——分块 SPFA

    挺有趣的分块的题目. 直接暴力建边SPFA貌似是$O(nm)$的. 然后考虑分块,$\sqrt n$一下用虚拟节点辅助连边, 以上的直接暴力连边即可. 然后卡卡时间,卡卡空间. 终于在UOJ上T掉辣. ...

  9. BZOJ 1600 [Usaco2008 Oct]建造栅栏:dp【前缀和优化】

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1600 题意: 给你一个长度为n的木板,让你把这个木板切割成四段(长度为整数),并且要求这四 ...

随机推荐

  1. 高斯消元模板(pascal)

    洛谷P3389评测 program rrr(input,output); const eps=1e-8; var a:..,..]of double; n,i,j,k:longint; t:doubl ...

  2. Git Gerrit Repo User Manual

                      Git Repo Gerrit User Manual Revision History   Revision # Description Date Author ...

  3. Innobackupx工具命令简单解析

    --defaults-file 同xtrabackup的--defaults-file参数,指定mysql配置文件; --apply-log 对xtrabackup的--prepare参数的封装; - ...

  4. php-fpm: hundreds of seconds in the log

    favoriteI have nginx+php-fpm web serverSo I've noticed in php5-fpm.log many strange lines:[03-Sep-20 ...

  5. ping: unknown host 解决办法

    如果ping命令返回如下错误,那主要的可能性就是系统的DNS设置有误. [root@CentOS5 ~]# ping www.sina.com.cn ping: unknown host www.si ...

  6. 【bzoj4484】【jsoi2015】最小表示

    Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 432  Solved: 223[Submit][Status][Discuss] Descriptio ...

  7. 【bzoj4940】这是我自己的发明

    Portal --> bzoj4940 Solution (原题这题面到底是..怎么回事啊深深的套路qwq) 感觉自己对根号的算法还是很..没有感觉啊== 实际上这题和bzoj5016没有任何区 ...

  8. 驱动之SPI,UART,I2C的介绍与应用20170118

    这篇文章主要介绍基本的驱动也是用的最多的协议类驱动中的SPI,I2C和UART.首先从最简单的UART也就是串口讲起: 1.UART UART由两根线也就是TX,RX以及波特率产生器组成,操作比较简单 ...

  9. 图像处理之均值滤波介绍及C算法实现

    1 均值滤波介绍 滤波是滤波是将信号中特定波段频率滤除的操作,是从含有干扰的接收信号中提取有用信号的一种技术. 均值滤波是典型的线性滤波算法,它是指在图像上对目标像素给一个模板,该模板包括了其周围的临 ...

  10. vs下给生成的程序(exe)加入默认的申请管理员权限

    在vs下编程时,经常会用到一些特殊的权限,尤其是管理员权限,下面是在win7下笔者亲测通过的. 下面要分两种情况: 第一种情况就是在vs2003或者以下的编译器中运行 1.准备一个manifest文件 ...