【BZOJ 1010】玩具装箱-CSDN博客

本文链接：https://blog.csdn.net/yhf_2015/article/details/53968770

本文介绍了一种结合动态规划与单调队列优化的方法，用于解决特定类型的最优化问题。通过状态定义、状态转移方程推导及代码实现，详细展示了如何利用单调队列优化动态规划中的计算过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

思路：

设 $f[i]$ 表示 $1$ 到 $i$ 号物品的最小代价， $s[i]$ 为 $c[i]$ 的前缀和。
状态转移方程： $f[i]=min\{f[j]+(s[i]-s[j]+i-j-1-L)^2\}$
考虑对方程进行优化。
令 $t[i]=s[i]+i$ ，则 $f[i]=min\{f[j]+(t[i]-t[j]-1-L)^2\}$
由于 $c[i]$ 非负，所以 $s[i]$ 在 $i$ 增长时单调不减。
令 $m=t[i]-1-L$ ，则 $f[i]=min\{f[j]+(m-t[j])^2\}$
设 $i$ 从 $p,q$ ，两个可能的 $j$ 的取值的地方转移，不妨令 $p<q$ 。
若方案 $q$ 优于 $p$ ，则有： $f[p]+(m-t[p])^2>f[q]+(m-t[q])^2$
展开： $f[p]-2m\times t[p]+t[p]^2>f[q]-2m\times t[q]+t[q]^2$
整理：

( f [ q ] + t [ q ] 2 ) - ( f [ p ] + t [ p ] 2 ) 2 \times ( t [ q ] - t [ p ] ) < m

$\frac{(f[q]+t[q]^2)-(f[p]+t[p]^2)}{2\times(t[q]-t[p])}<m$
因为

i $i$ 增加的的时候，

m $m$ 是单调不降的，所以可以用单调队列维护。

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
const int maxn = 50010;
int n, L;
ll c[maxn], f[maxn], s[maxn], q[maxn], t[maxn], l, r;
ll sqr(ll x){return x*x;}
double slope(ll p, ll q){
    return ((f[q]+sqr(t[q]))-(f[p]+sqr(t[p])))/(2.0*(t[q]-t[p]));
}
int main(){
    scanf("%d%d", &n, &L);
    for(int i = 1; i <= n; i ++){
        scanf("%lld", &c[i]);
        s[i] += s[i-1]+c[i];
        t[i] = s[i] + i;
    }
    l = 1, r = 1;
    for(int i = 1; i <= n; i ++){
        while(l<r && slope(q[l], q[l+1]) < t[i]-L-1) l ++;
        f[i] = f[q[l]] + sqr(t[i]-t[q[l]]-L-1);
        while(l<r && slope(q[r-1], q[r]) > slope(q[r], i)) r --;
        q[++ r] = i;
    }
    printf("%lld", f[n]);
    return 0;
}