前言


纪念一下我做的第二道斜率优化$dp$题,终于自己能把代码敲出来了,然而有很智障的$bug$,把$i$写成$q[i]$,找了半天QAQ。然后写$dp$公式并优化的能力稍微强了一点(自我感觉良好),对于斜率优化$dp$"去尾"的操作理解更深刻了

描述


$1∼N$号工厂,第$i$个工厂有$P_{i}$个成品,第$i$个工厂建立仓库需要$C_{i}$的费用,该工厂距离第一个工厂的距离为$X_{i}$,编号小的工厂只能往编号大的工厂搬用成品,每单位成品搬每单位距离需要花费1,问所有成品搬到工厂里面所需的最少费用是多少 [Link]

分析


设$f[i]$为第$i$个工厂建立仓库,前$i$个工厂的成品都搬到仓库中的最小花费,则容易得到动态转移方程:

         $f[i]=min(f[j]+P_{j+1}(X_{i}-X_{j+1})+P_{j+2}(X_{i}-X_{j+2})+\cdots +P_{i-1}(X_{i}-X_{i+1}))+Ci$

通式为

               $f[i]=min(f[j]+\sum_{k=j+1}^{i-1}P_{k}\cdot X_{i}-\sum_{k=j+1}^{i-1}P_{k}\cdot X_{k})+C_{i}$

令 $s[i]=\sum_{1}^{i}P[i],g[i]=\sum_{1}^{i}P_{i}\cdot X_{i}$
则方程变为
               $f[i]=min(f[j]+X_{i}\cdot(s[i−1]−s[j])-(g[i−1]−g[j]))+C_{i}$

则对于最优决策 $j$ ,有

                 $f[j]+g[j]=X_{i}\cdot s[j]+f[i]-X_{i}\cdot s[i−1]−Ci$

也就是要找 $y=kx+b$,$k$已知,找一对$x,y$使得截距最小

Code

#include <cstdio>
#define ll long long
#define empty (head>=tail)
const int maxn = 1e6+10;
ll n, head, tail, j;
ll x[maxn], p[maxn], c[maxn];
ll q[maxn], s[maxn], g[maxn], f[maxn];
inline long double X(ll i) {return s[i];}
inline long double Y(ll i) {return f[i]+g[i];}
inline long double rate(ll j,ll k) {return (Y(k)-Y(j))/(X(k)-X(j));}
int main()
{
scanf("%lld", &n);
for (int i = 1; i <= n; i++) {
scanf("%lld%lld%lld", &x[i], &p[i], &c[i]);
s[i] = s[i-1]+p[i], g[i] = g[i-1]+p[i]*x[i];
}
head = tail = 1;
for(int i = 1; i <= n; i++){
while(!empty&&rate(q[head],q[head+1])<x[i])head++;
j = q[head]; f[i] = f[j]+x[i]*(s[i-1]-s[j])-(g[i-1]-g[j])+c[i];
while(!empty&&rate(q[tail-1],q[tail])>rate(q[tail],i))tail--;
q[++tail] = i;
}
printf("%lld\n", f[n]);
}

思考


之前入门的题目说过,假如$g(c,b)<=g(b,a)$,那么这个$b$就没有任何用武之地了,就把这样的点去掉后再加入新的点。那么我会想到,能不能先不去掉这个$b$,直接加入新的点,就是在斜率小于$A[i]$的时候再去掉不也可以么,顶多时间长一点,然而现实给了我一发$WA$。我就去想为什么会是这个亚子?原因在于,必须保证通过前面的状态得到的$f[i]$为最优,而你的$b$可能会导致你通过原有的方式得到的$f[i]$并不是最优的,那么得到的$f[n]$自然不一定是最优解。

那我就去做了尝试,我的确没有去掉$b$,那么我只要从之前的状态里挑最优不也可以么,交上去之后发现的确可以(妥妥$WA$)。这就很头疼,也证明我的想法是不对的,因为在有$b$的前提下,你按照原来那种判断两点间斜率的方法去决定取哪个点的方法并不能确定哪个点是最优的,比如下图:

这个时候$C$是不满足要求的点,我没有去掉,之后我从$A,B,C,D$中挑选最优的点去更新$E$点,$E$对应的斜率$rate$是绿色的那条线,从最下面开始扫描,$AB$的斜率小于$rate$,扫到$B$这里停下,由于$BC$斜率大于$rate$,还是在$B$这里,之后$CD$的斜率比$rate$小,那就扫到$D$这里停下,然后通过$D$更新$E$点。但是你会发现其实$B$点比$D$点更优,这样选是不对的。但是如果我先把$C$点去掉,那么通过$AB,BD$的斜率去选择点,你会发现得到的$E$点是最优的

所以结果就是,正是因为满足要求的点之间的斜率是单调递增的才能保证每次得到的$f[i]$为最优,并且这样计算$f[i]$的确效率很高~

附上错误的代码(T_T):

#include <cstdio>
#define ll long long
#define empty (head>=tail)
const int maxn = 1e6+10;
ll n, head, tail, j;
ll x[maxn], p[maxn], c[maxn];
ll q[maxn], s[maxn], g[maxn], f[maxn];
inline long double X(ll i) {return s[i];}
inline long double Y(ll i) {return f[i]+g[i];}
inline long double rate(ll j,ll k) {return (Y(k)-Y(j))/(X(k)-X(j));}
int main()
{
scanf("%lld", &n);
for (int i = 1; i <= n; i++) {
scanf("%lld%lld%lld", &x[i], &p[i], &c[i]);
s[i] = s[i-1]+p[i], g[i] = g[i-1]+p[i]*x[i];
}
head = tail = 1;
for (int i = 1; i <= n; i++) {
int l = head, r = tail;
while(!empty&&l<r) {
if(rate(q[l],q[l+1])<x[i]) head = l+1;
l++;
}
j = q[head]; f[i] = f[j]+x[i]*(s[i-1]-s[j])-(g[i-1]-g[j])+c[i];
q[++tail] = i;
}
printf("%lld\n", f[n]);
}

参考文章:

https://www.cnblogs.com/1625--H/p/11267043.html

[ZJOI2007]仓库建设(斜率dp优化)的更多相关文章

  1. bzoj-1096 1096: [ZJOI2007]仓库建设(斜率优化dp)

    题目链接: 1096: [ZJOI2007]仓库建设 Description L公司有N个工厂,由高到底分布在一座山上.如图所示,工厂1在山顶,工厂N在山脚.由于这座山处于高原内陆地区(干燥少雨),L ...

  2. BZOJ 1096: [ZJOI2007]仓库建设 [斜率优化DP]

    1096: [ZJOI2007]仓库建设 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 4201  Solved: 1851[Submit][Stat ...

  3. BZOJ 1096: [ZJOI2007]仓库建设(DP+斜率优化)

    [ZJOI2007]仓库建设 Description L公司有N个工厂,由高到底分布在一座山上.如图所示,工厂1在山顶,工厂N在山脚.由于这座山处于高原内陆地区(干燥少雨),L公司一般把产品直接堆放在 ...

  4. bzoj1096[ZJOI2007]仓库建设 斜率优化dp

    1096: [ZJOI2007]仓库建设 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 5482  Solved: 2448[Submit][Stat ...

  5. 【BZOJ1096】[ZJOI2007]仓库建设 斜率优化

    [BZOJ1096][ZJOI2007]仓库建设 Description L公司有N个工厂,由高到底分布在一座山上.如图所示,工厂1在山顶,工厂N在山脚.由于这座山处于高原内陆地区(干燥少雨),L公司 ...

  6. bzoj 1096: [ZJOI2007]仓库建设 斜率優化

    1096: [ZJOI2007]仓库建设 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 2242  Solved: 925[Submit][Statu ...

  7. 【BZOJ】1096: [ZJOI2007]仓库建设(dp+斜率优化)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1096 首先得到dp方程(我竟然自己都每推出了QAQ)$$d[i]=min\{d[j]+cost(j+ ...

  8. 【bzoj1096】[ZJOI2007]仓库建设 斜率优化dp

    题目描述 L公司有N个工厂,由高到底分布在一座山上.如图所示,工厂1在山顶,工厂N在山脚.由于这座山处于高原内陆地区(干燥少雨),L公司一般把产品直接堆放在露天,以节省费用.突然有一天,L公司的总裁L ...

  9. P2120 [ZJOI2007]仓库建设 斜率优化dp

    好题,这题是我理解的第一道斜率优化dp,自然要写一发题解.首先我们要写出普通的表达式,然后先用前缀和优化.然后呢?我们观察发现,x[i]是递增,而我们发现的斜率也是需要是递增的,然后就维护一个单调递增 ...

  10. [ZJOI2007]仓库建设(斜率优化)

    L公司有N个工厂,由高到底分布在一座山上. 工厂1在山顶,工厂N在山脚. 由于这座山处于高原内陆地区(干燥少雨),L公司一般把产品直接堆放在露天,以节省费用. 突然有一天,L公司的总裁L先生接到气象部 ...

随机推荐

  1. Syn_Flood攻击&防御手段

    Syn_Flood攻击原理 攻击者首先伪造地址对服务器发起SYN请求(我可以建立连接吗?),服务器就会回应一个ACK+SYN(可以+请确认).而真实的IP会认为,我没有发送请求,不作回应.服务器没有收 ...

  2. SpringBoot同时接收单个对象和List<object>参数

    最近做项目的有个需求,是把多个文件移动到另一个文件夹下,这需要把 新的文件夹id -- Long类型 多个文件的信息 -- List< Object > 类型 这两个参数传给后台,我的后台 ...

  3. 【JS学习】var let const声明变量的异同点

    [JS学习]var let const声明变量的异同点 前言: 本博客系列为学习后盾人js教程过程中的记录与产出,如果对你有帮助,欢迎关注,点赞,分享.不足之处也欢迎指正,作者会积极思考与改正. 总述 ...

  4. python3 最基本且简单的实现组合设计模式

    组合模式是把一个类别归为一个整体,并且组织多个整体之间的关系,使用通过树形结构来描述所有整体. 一般的编写方式为一个类别对象下添加多个该对象的元素,该对象又为其它对象下的元素.简单来说就是一个学校有一 ...

  5. Redis中哈希分布不均匀该怎么办

    前言 Redis 是一个键值对数据库,其键是通过哈希进行存储的.整个 Redis 可以认为是一个外层哈希,之所以称为外层哈希,是因为 Redis 内部也提供了一种哈希类型,这个可以称之为内部哈希.当我 ...

  6. linux 文件目录权限

    文件目录权限: 什么是文件权限: 在Linux中,每个文件都有所属的所有者,和所有组,并且规定了文件的所有者,所有组以及其他人对文件的,可读,可写,可执行等权限. 对于目录的权限来说,可读是读取目录文 ...

  7. Kubernetes之GlusterFS集群文件系统高可用安装,提供动态卷存储

    GlusterFS高可用安装 一. 准备工作 安装好的k8s集群,提供其中三个节点给GFS,这三个节点都至少有一个可用的裸块设备 在k8s所有节点安装所需要的组件 # ubuntu16.04 add- ...

  8. html新特性笔记

    HTML5知识总结 l  文档类型声明:<!DOCTYPE HTML> l  新绘制元素: Canvas:标签定义图形,比如图表和其他图像.该标签基于 JavaScript 的绘图 API ...

  9. 路由协议-RIP协议

    一.路由协议的发展历程和分类 距离矢量路由协议--听信"谣言",使用跳数作为度量值,最大16(0-15)跳:RIP 链路状态路由协议--"地图"路由协议:OSP ...

  10. jQuery 移入显示div,移出当前div,移入到另一个div还是显示。

    jQuery 移入移出 操作div 1 <style type="text/css"> 2 .box{ 3 position: relative; 4 } 5 .box ...