玩具装箱题解 - 洛谷
玩具装箱题解 - cnblogs
斜率优化 - OIWiki
玩具装箱（HAOI2008）
P 教授要去看奥运，但是他舍不下他的玩具，于是他决定把所有的玩具运到北京。他使用自己的压缩器进行压缩，其可以将任意物品变成一堆，再放到一种特殊的一维容器中。
P 教授有编号为 \(1 \cdots n\) 的 \(n\) 件玩具，第 \(i\) 件玩具经过压缩后的一维长度为 \(C_i\)。
为了方便整理，P教授要求：
- 在一个一维容器中的玩具编号是连续的。
- 同时如果一个一维容器中有多个玩具，那么两件玩具之间要加入一个单位长度的填充物。形式地说，如果将第 \(i\) 件玩具到第 \(j\) 个玩具放到一个容器中，那么容器的长度将为 \(x=j-i+\sum\limits_{k=i}^{j}C_k\)。
制作容器的费用与容器的长度有关，根据教授研究，如果容器长度为 \(x\)，其制作费用为 \((x-L)^2\)。其中 \(L\) 是一个常量。P 教授不关心容器的数目，他可以制作出任意长度的容器，甚至超过 \(L\)。但他希望所有容器的总费用最小。
设 \(f_i\) 表示前 \(i\) 个玩具装箱的最小费用，\(s_i\) 表示 \(c_i\) 的前缀和，则有
\[f_i=\min_{j<i}\{f_j+(s_i-s_j+i-j-1-L)^2\} \]如一个一个遍历，复杂度为 \(O(n^2)\)。展开后会出现 \(-2s_is_j\) 这一项，不满足单调队列优化判定式 \(\displaystyle f_i=\min_{j<i}\{a_i+b_j\}\)。即需要最小化的多项式和 \(i,j\) 均有关。
将 \(()^2\) 括号内与 \(i\) 有关的设为 \(A=s_i+i\)，与 \(j\) 有关，或与 \(i,j\) 均无关的设为 \(B=s_j+j+1+L\)，先不考虑 \(\min\)，则有
\[f_i=f_j+(A-B)^2 \]展开后移项得到
\[f_j+B^2=2AB+f_i-A^2 \]上式只有 \(f_i-A^2\) 未知。通过考虑其的几何意义，在均摊 \(O(1)\) 的时间内转移求得最小值，便是「斜率优化」。
将 \((B,f_j+B^2)\) 视作平面中的点，记为 \(P_j\)。对于 \(j\)，当等式成立时，相当于有一条斜率为 \(2A\) 的直线经过了点 \(P_j\)。而这条直线在 \(y\) 轴上的截距便是 \(f_i-A^2\)。
要求出最小的截距，可以把斜率为 \(2A\) 的直线从下往上平移，直到碰到第一个点，此时的 \(y\) 轴截距即为最小。
考虑凸包的几何意义（点集的「边界」），可以发现使得截距最小的点一定在点集的下凸包上。
继续观察，还可以发现一个事实。若设构成下凸包的点集为 \(S\)，并从左往右标号为 \(S_1,S_2,\dots, S_n\)，设两点 \(P_u,P_v\) 连成的直线的斜率为 \(K(P_u,P_v)\)。则使得截距最小的点 \(S_k\) 必满足 \(K(S_{k-1},S_k)<2A\)，\(K(S_k,S_{k+1})> 2A\)。
凸包本身可以用队列维护。维护方法是在加入一个点 \(P_i\)，判断队尾的点 \(P_r\) 是否满足 \(K(P_r,P_{r-1})>K(P_r,P_i)\)。如果满足，则弹出 \(P_r\)。重复此过程直到 \(K(P_r,P_{r-1})<K(P_r,P_i)\) 或队列中元素不多于一个。
而此题的斜率 \(2A\) 单调递增，则 \(S_k\) 左边点的个数也单调递增，所以可以不断 pop 使得 \(K(P_l,P_{l+1})<2A\) 成立的队头。pop 完之后，队头即为 \(S_k\)。
此时来归纳一下本题的做法，对于每个点 \(P_i\)：
- 若 \(K(P_l,P_{l+1})<2A\)，则将队头弹出，直到队列中元素数量不多于一个或条件不成立。
- 取出队头，计算 \(f_i\)。
- 若 \(K(P_r,P_{r-1})>K(P_r,P_i)\)，则弹出队尾，直到队列中元素不多于一个或条件不成立。
- 将 \(P_i\) 加入队尾。
Code：
- 队列中元素多于一个的「代码意义」为
head < tail。 - 将简化后的式子中的变量用函数写出来。
换成注释里的写法就莫名其妙是错的，可能是因为精度问题，我暂且谔谔。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 5e4 + 5;
int n, L, q[N], hd, tl;
ll s[N], f[N];
inline double A(int i) { return s[i] + i; }
inline double B(int i) { return s[i] + i + L + 1; }
inline double Y(int i) { return B(i) * B(i) + f[i]; }
inline double X(int i) { return B(i); }
inline double K(int i, int j) { return (Y(j) - Y(i))/(X(j) - X(i)); }
int main() {
ios::sync_with_stdio(false); cin.tie(nullptr);
cin >> n >> L;
for(int i = 1; i <= n; i++)
cin >> s[i], s[i] += s[i - 1];
hd = tl = 1; // 相当于 q[1] = f[0] = s[0] = 0
for(int i = 1; i <= n; i++) {
while(hd < tl && K(q[hd], q[hd + 1]) < 2 * A(i)) ++hd;
f[i] = f[q[hd]] + (ll)(pow((A(i) - B(q[hd])), 2) + 0.1); // 避免精度误差
// f[i] = f[q[hd]] + B(q[hd]) * B(q[hd]) - 2 * A(i) * B(q[hd]) + A(i) * A(i);
while(hd < tl && K(q[tl - 1], q[tl]) > K(q[tl], i)) --tl;
q[++tl] = i;
}
cout << f[n] << '\n';
return 0;
}
标签：容器,截距,int,玩具,斜率,tl,笔记,DP 来源： https://www.cnblogs.com/huaruoji/p/15466805.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。