题目描述





分析

  • 对于 \(Subtask\ 1\),可以写一个 \(n^3\) 的 \(DP\),\(f[i][j]\) 代表第 \(i\) 个建筑高度为 \(j\) 时的最小花费,随便转移即可

    时间复杂度 \(O(n \times h^2)\)
  • 对于 \(Subtask\ 2\),我们沿用 \(Subtask\ 1\)的思路,记录前缀后缀 \(min\),将复杂度优化至 \(O(n \times h)\)

    但是显然两维的定义无法继续进行优化,我们可以考虑改变一下定义的方式

    设 \(f[i]\) 表示考虑前 \(i\) 个建筑,并且第 \(i\) 个建筑高度不变的最优答案

    可以发现,枚举两个不变的边界,那么中间的建筑必定被提高成相同的小于等于边界的高度

    也就是说我们需要把一些坑填平

    因为增加峰的高度既花人力又不能提高观赏度

    增加坡的高度花人力但不能提高观赏度

    只有增加坑的高度才会有贡献,而且增加的值不能超过边界的高度

    因为超过边界的高度又变成了峰

    因此我们可以枚举边界然后再枚举填平的高度

    转移方程为 \(f[i]=f[j]+(sum2[i-1]-sum2[j])+(i-j-1)*h*h-(sum1[i-1]-sum1[j])*2*h+c*abs(a[i]+a[j]-2*h)\)

    其中 \(sum1[i]\) 为 \(a\) 数组的前缀和,\(sum2[i]\) 为 \(a\) 数组平方的前缀和

    上面的式子是展开后的式子,原式子并不难推

    注意这种做法我们需要把 \(a[0]\) 和 \(a[n+1]\) 置为无穷大,因为我们有可能提高第一个和最后一个建筑的高度

    要特判 \(i,j\) 等于 \(0\)或 \(n+1\) 的情况

    最后的答案为 \(f[n+1]\)

    时间复杂度 \(O(n^2 \times h)\),不知道为什么能过这个子任务
  • 对于 \(Subtask\ 3\),我们把一维 \(DP\) 的状态转移方程化简得到

    \(f[i]=f[j]+(i-j-1)*h*h-2*(sum1[i-1]-sum1[j]+c)*h+(sum2[i-1]-sum2[j])+c*(a[i]+a[j])\)

    我们发现前半部分是一个关于高度 \(h\) 的二次函数

    可以直接由对称轴求出最小值

    时间复杂度 \(O(n^2)\)
  • 对于 \(Subtask\ 4\) 和 \(Subtask\ 5\),我们会发现只有两个高度较大的建筑夹着一堆高度较小的建筑才有贡献

    因此可以用单调队列(栈)维护

    时间复杂度 \(O(n logn)\)

    复杂度的瓶颈在\(ST\) 表查询最值上

代码

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<cstring>
#define rg register
inline int read(){
rg int x=0,fh=1;
rg char ch=getchar();
while(ch<'0' || ch>'9'){
if(ch=='-') fh=-1;
ch=getchar();
}
while(ch>='0' && ch<='9'){
x=(x<<1)+(x<<3)+(ch^48);
ch=getchar();
}
return x*fh;
}
typedef long long ll;
const int maxn=1e6+5;
int n,c,a[maxn],maxh,lg[maxn],st[maxn][22],head,tail,q[maxn];
ll f[maxn],sum1[maxn],sum2[maxn];
int zhao(ll i,ll j){
if(j==0 && i==n+1) return ((double)(sum1[i-1]-sum1[j])/(double)(i-j-1)+0.5);
else if(j==0 || i==n+1) return ((double)(2*sum1[i-1]-2*sum1[j]+c)/(double)(2.0*(i-j-1))+0.5);
return ((double)(c+sum1[i-1]-sum1[j])/(double)(i-j-1)+0.5);
}
int cx(int l,int r){
rg int k=lg[r-l+1];
return std::max(st[l][k],st[r-(1<<k)+1][k]);
}
int main(){
memset(f,0x3f,sizeof(f));
n=read(),c=read();
for(rg int i=1;i<=n;i++){
a[i]=read();
maxh=std::max(maxh,a[i]);
sum1[i]=sum1[i-1]+a[i];
sum2[i]=sum2[i-1]+1LL*a[i]*a[i];
st[i][0]=a[i];
}
for(rg int i=2;i<=n;i++){
lg[i]=lg[i/2]+1;
}
for(rg int j=1;j<=20;j++){
for(rg int i=1;i+(1<<j)-1<=n;i++){
st[i][j]=std::max(st[i][j-1],st[i+(1<<(j-1))][j-1]);
}
}
a[0]=a[n+1]=0x3f3f3f3f;
f[0]=f[1]=0;
head=1,tail=1;
rg int mmax,mmin,now;
for(rg int i=1;i<=n+1;i++){
while(head<=tail){
if(q[tail]==i-1){
if(i<=n) f[i]=std::min(f[i],f[q[tail]]+1LL*c*std::abs(a[i]-a[i-1]));
else f[i]=f[i-1];
} else {
mmax=cx(q[tail]+1,i-1);
mmin=std::min(a[i],a[q[tail]]);
if(mmin>=mmax){
now=zhao(i,q[tail]);
if(now<mmax) now=mmax;
if(now>mmin) now=mmin;
if(q[tail]==0 && i==n+1) f[i]=std::min(f[i],f[q[tail]]+1LL*(i-q[tail]-1)*now*now-2LL*(sum1[i-1]-sum1[q[tail]])*now+1LL*(sum2[i-1]-sum2[q[tail]]));
else if(q[tail]==0) f[i]=std::min(f[i],f[q[tail]]+1LL*(i-q[tail]-1)*now*now-1LL*(2*sum1[i-1]-2*sum1[q[tail]]+c)*now+1LL*(sum2[i-1]-sum2[q[tail]])+1LL*c*a[i]);
else if(i==n+1) f[i]=std::min(f[i],f[q[tail]]+1LL*(i-q[tail]-1)*now*now-1LL*(2*sum1[i-1]-2*sum1[q[tail]]+c)*now+1LL*(sum2[i-1]-sum2[q[tail]])+1LL*c*a[q[tail]]);
else f[i]=std::min(f[i],f[q[tail]]+1LL*(i-q[tail]-1)*now*now-2LL*(sum1[i-1]-sum1[q[tail]]+c)*now+1LL*(sum2[i-1]-sum2[q[tail]])+1LL*c*(a[i]+a[q[tail]]));
}
}
if(a[i]>=a[q[tail]])tail--;
else break;
}
q[++tail]=i;
}
printf("%lld\n",f[n+1]);
return 0;
}

联赛模拟测试18 A. 施工 单调队列(栈)优化DP的更多相关文章

  1. 单调队列 &amp;&amp; 斜率优化dp 专题

    首先得讲一下单调队列,顾名思义,单调队列就是队列中的每个元素具有单调性,如果是单调递增队列,那么每个元素都是单调递增的,反正,亦然. 那么如何对单调队列进行操作呢? 是这样的:对于单调队列而言,队首和 ...

  2. HDU 3507 单调队列 斜率优化

    斜率优化的模板题 给出n个数以及M,你可以将这些数划分成几个区间,每个区间的值是里面数的和的平方+M,问所有区间值总和最小是多少. 如果不考虑平方,那么我们显然可以使用队列维护单调性,优化DP的线性方 ...

  3. 【P2422】良好的感觉(单调栈优化DP//奇怪的暴力)

    话说正解是单调栈优化DP,然而貌似根据某种玄学的推算,这个题暴力出解貌似也是可以的.首先,我们枚举所有的点作为最小点,然后横向展开,遇到更小的就停止...然后再操作一下,看上去时间O(N^2),然而由 ...

  4. 【BZOJ-2892&amp;1171】强袭作战&amp;大sz的游戏 权值线段树+单调队列+标记永久化+DP

    2892: 强袭作战 Time Limit: 50 Sec  Memory Limit: 512 MBSubmit: 45  Solved: 30[Submit][Status][Discuss] D ...

  5. luogu 2216 理想的正方形 单调队列(其实没有DP)

    #include<bits/stdc++.h> using namespace std; ; ; int a,b,n; int g[A][A],q[A][N],Q[A][N]; int h ...

  6. 洛谷 P2254 [NOI2005]瑰丽华尔兹(单调栈优化DP)

    题目描述 不妨认为舞厅是一个N行M列的矩阵,矩阵中的某些方格上堆放了一些家具,其他的则是空地.钢琴可以在空地上滑动,但不能撞上家具或滑出舞厅,否则会损坏钢琴和家具,引来难缠的船长.每个时刻,钢琴都会随 ...

  7. 联赛模拟测试20 C. Weed

    题目描述 \(duyege\) 的电脑上面已经长草了,经过辨认上面有金坷垃的痕迹. 为了查出真相,\(duyege\) 准备修好电脑之后再进行一次金坷垃的模拟实验. 电脑上面有若干层金坷垃,每次只能在 ...

  8. [考试反思]0812NOIP模拟测试18:稀释

    []200 [4]190[5]180 [6]170[7]150 [11]140[16]130[19]120[24]100 有些事情,看淡了,也就简单了. 连续爆炸之后,没什么感觉了. 把态度放正,把结 ...

  9. [CSP-S模拟测试]:小P的单调数列(树状数组+DP)

    题目描述 小$P$最近喜欢上了单调数列,他觉得单调的数列具有非常多优美的性质.经过小$P$复杂的数学推导,他计算出了一个单调增数列的艺术价值等于该数列中所有书的总和.并且以这个为基础,小$P$还可以求 ...

随机推荐

  1. css应对已有class和特殊class的冲突

    类之间是没有优先级的,当需要应用特殊样式时,可以先删除通用class后增加特殊class

  2. Hibernate 异常 —— Unable to instantiate default tuplize

    出现这个异常 —— Unable to instantiate default tuplizer ,是 Hibernate 的映射文件(*.hbm.xml)导致的.仔细检查一下工程里的映射文件吧. 笔 ...

  3. windows下一个erlang包装镜像启动

    于linux环境,erlang经systools:make_script("",[])和systools:make_tar()命令生成图像包,安装镜像包,图片包的安装过程,通过替换 ...

  4. c++学籍管理系统

    程序在编译时出错(vc++ 6.0) 求哪位大神帮忙改改 #include<iostream> #include <string> #include<conio.h> ...

  5. [转]js 取得 Unix时间戳(Unix timestamp)

    本文转自:https://blog.csdn.net/o0snow/article/details/6858829 js 取得 Unix时间戳 Unix时间戳(Unix timestamp),或称Un ...

  6. 【转载】论文笔记系列-Tree-CNN: A Deep Convolutional Neural Network for Lifelong Learning

    一. 引出主题¶ 深度学习领域一直存在一个比较严重的问题——“灾难性遗忘”,即一旦使用新的数据集去训练已有的模型,该模型将会失去对原数据集识别的能力.为解决这一问题,本文提出了树卷积神经网络,通过先将 ...

  7. SUID、SGID、粘滞位

    粘滞位(Stikybit) +t,只有用户自己可以删除自己创建文件,其他用户只能查看,不能删除.        1:创建两个用户  useradd oo                         ...

  8. CONTINUOUS MIGRATION

    转自:https://pgloader.io/blog/continuous-migration/ After having been involved in many migration proje ...

  9. 【转】学习使用:before和:after伪元素

    如果你一直密切关注着各种网页设计的博客,你可能已经注意到了:before和:after伪元素已经在前端开发中获得了相当多的关注.特别是在Nicolas Gallagher的博客中,后期运用了很多伪类元 ...

  10. 青蛙的约会(poj1061+欧几里德)

    青蛙的约会 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 94174   Accepted: 17412 Descripti ...