【笔记篇】斜率优化dp（三） APIO特别行动队

最新推荐文章于 2022-02-13 15:31:38 发布

Enzymii

最新推荐文章于 2022-02-13 15:31:38 发布

阅读量502

点赞数

分类专栏：ＤＰ之神虐杀众生文章标签：斜率优化 dp及优化

本文链接：https://blog.csdn.net/enzymii/article/details/79167201

版权

ＤＰ之神虐杀众生专栏收录该内容

11 篇文章 0 订阅

订阅专栏

旁听了一波给舒老师和学弟的pkuwc面试讲座…

这里有一段隐身的吐槽, 请反白观看. 不想看的跳过这一段看似空白的东西就好了…
刚开始ATP学姐给我们讲了自己面试的时候的事情..描绘了一下当时面试的场面和当时问的问题…~~ATP学姐太可爱了ＯvO~~可能准备的也不是很充足用了好多的”然后” “就是”之类的, 自己也在吐槽, 老师后来也槽来着, ~~不过还是超可爱n(≧▽≦)n~~旁边老师在记录, 尤其记下了一个什么”你认为学习信息男生和女生有什么区别?”的问题..~~可是这种问题怎么可能问男生啊讲道理~~, 然后学姐讲完之后就是shallwe大爷来讲, shallwe大爷还准备了稿子, 然后就讲了一些面基技巧.. 但是也是比较老生常谈的东西, 之前参加各种各样的风采大赛回答问题的时候也就是那一套之类的小时候还能做到, 现在反而越来越不行了.. 然后讲了各种各样的问题, 说是问的问题都在blog里面.. 老师还说要blog的地址来给学弟学妹什么的(那也就是包括我们咯), 但是听完回来一看这blog里的东西能给老师看?! 各种吐槽甚至还有吐槽面试培训的.. 然后他们就回去学习了= =
老师就开始种讲评, 还让两个人分别回答了上面说的男女区别问题然后两个都做出了比较”片面”的回答, 老师就说要往什么三个方面想之类的说了一堆, 反正就是不怼人的同时优雅地装逼的技巧, ~~什么在北大不说清华比北大强之类的, 说了确定不是石乐志么= =~~然后再就是”你有什么缺点”, 听说年年问..~~大家准备好的答案都能倒着背了好么~~就说要找个”不是缺点的缺点”, 然后两个人都说自己比较执着啊… 就不能找点有创意的么←←, 感觉全世界都有执着得有点傻这个缺点… 要我的话可能会说自己好奇心过强? 反正就这么着吧. 然后还有什么进门关门鞠躬瞅瞅有没有废纸倒拖把什么的… 反正就是一些鸡汤文里描写的面试的细节问题, 还有着装问题, 我觉得只要不是明显违和的穿女装应该都不是什么大问题吧←←~~对没错哪怕是可爱的男孩纸穿女装其实也是可以哒~(≧▽≦)/~~~两个人又练了一波关门鞠躬之类的终于煎熬地熬过去了.. 然后赶紧去了个厕所毕竟半下午+半晚上没有去了.. 回来之后就调这个题不过很快就调A了这还比较和善..~~嗯主要原因是这次的luogu没有出锅...~~

===我===是===一===条===分===割===线===我===也===是===一===扇===传===送===门===

继续吐槽, 不过这次跟题目有关了.. 听完回来就继续写这道题啊…
这道题显然支持 $O(n)$ 的做法, (感觉1e6的话套个log就不是很靠谱了)
那就是要优化, 就是要化式子嘛..然后刚开始以为这个题是

\sum i = 1 n (a x 2 i + b x i + c) = a \sum i = 1 n x 2 i + b s n + n * c

$\sum_{i=1}^n(ax_i^2+bx_i+c)=a\sum_{i=1}^nx_i^2+bs_n+n*c$
然后最大化

∑ni=1x2i $\sum_{i=1}^nx_i^2$ 之后直接输出

a∗f[n]+b∗s[n]+n∗c $a*f[n]+b*s[n]+n*c$ 的大水题,
但是发现好像样例差了好多好多…结果发现自己快读没读负数…
然后改完之后输出负数了… 手玩了一下发现好像这个

n∗c $n*c$ 的

n $n$ 好像不一定是题目里给的

n $n$ …
那就不能这么化式子了, 就开始重新化.. 还是得从状态转移方程入手…

f [i] = m a x {f [j] + a (s [i] - s [j]) 2 + b (s [i] - s [j]) + c} (j \in [1, i))

$f[i]=max\{f[j]+a(s[i]-s[j])^2+b(s[i]-s[j])+c\}(j\in[1,i))$
然后去括号

f [i] = f [j] + a s [i] 2 - 2 a s [i] s [j] + a s [j] 2 + b s [i] - b s [j] + c

$f[i]=f[j]+as[i]^2-2as[i]s[j]+as[j]^2+bs[i]-bs[j]+c$
移项能得到

f[j]+as[j]2−bs[j] $f[j]+as[j]^2-bs[j]$ =

2as[i] $2as[i]$

s[j] $s[j]$ +

f[i]−(as[i]2+bs[i]+c) $f[i]-(as[i]^2+bs[i]+c)$
由于题目中已经给了

−5<a<−1 $-5<a<-1$ , 所以斜率一定是<0的, 而

f[j]+as[j]2−bs[j] $f[j]+as[j]^2-bs[j]$ 直观上很明显是递减的. 我们考虑凸性对结果的影响.
这里写图片描述

我们考虑中间的点在位于两个斜率中间(大于大斜率小于小斜率显然不优就不考虑了)的情况.
发现如果是位于下凸壳的点是不优的, 而位于上凸壳的点是优的, 我们就要维护一个上凸壳.

而这里的斜率 $2as[i]$ 显然是递减的(绝对值递增), 所以
这里写图片描述
如果一个斜率一旦已经小于(绝对值大)于队首的边, 那么队首的点在以后也不会优了我们让队首出队即可.
依然是维护一个单调队列.
然后就没有了. 代码: 1A了哈哈哈哈 (放肆地笑ing…)

#include <cstdio>
const int N=16+6;
typedef long long LL;
LL f[N],s[N]; int q[N],h,t,a,b,c,n;
inline int gn(int a=0,char c=0,int f=1){
    for(;(c<48||c>57)&&c!='-';c=getchar());if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar()) a=a*10+c-'0'; return a*f;
}
inline double slope(int x,int y){
    return 1.0*(f[x]+a*s[x]*s[x]-b*s[x]-f[y]-a*s[y]*s[y]+b*s[y])/(s[x]-s[y]);
}
int main(){ n=gn(); a=gn(); b=gn(); c=gn();
    for(int i=1;i<=n;++i) s[i]=s[i-1]+gn();
    for(int i=1;i<=n;++i){
        while(h<t&&slope(q[h],q[h+1])>=2*a*s[i]) ++h;
        f[i]=f[q[h]]+a*(s[i]-s[q[h]])*(s[i]-s[q[h]])+b*(s[i]-s[q[h]])+c;
        while(h<t&&slope(q[t],q[t-1])<=slope(q[t],i)) --t;
        q[++t]=i;
    } printf("%lld",f[n]);
}