[HNOI2008]玩具装箱TOY

题面描述

传送门

思路

理解题意之后，发现

斜率优化真是一个大坑

$F_i$为最近一个箱子装到第$i$件物品的花费

首先我们来看看状态转移方程
($sum$为前缀和)

$$F_i=\min (F_i,F_j+(su{m}_i-su{m}_j+(i-(j+1))-L{)}^2)$$

$$=\min (F_i,F_j+((su{m}_i+i)-(su{m}_j+j)-(L+1){)}^2)$$

设$s_i=su{m}_i+i,s_j=su{m}_j+j,L=L+1$，

则$$F_i=\min (F_i,F_j+(s_i-s_j-L{)}^2)$$

设有$F_j+(s_i-s_j-L{)}^2\ge F_k+(s_i-s_k-L{)}^2(j<k<i)①$

说明在$i$时，$k$优于$j$。

那么对于后面状态$t$呢，是否有$F_j+(s_t-s_j-L{)}^2\ge F_k+(s_t-s_k-L{)}^2$

证明:

显而易见的，$s_t=s_i+val(val>0)$，则

$$\begin{gathered} F_j+(s_t-s_j-L{)}^2=F_j+(s_i-s_j-L+val{)}^2② \\ {F}_k+(s_t-s_k-L{)}^2=F_k+(s_i-s_j-L+val{)}^2③ \end{gathered}$$

由$②$得:

$$F_j+(s_i-s_j-L{)}^2+2\ast val\ast (s_i-s_j-L)+va{l}^2$$

由$③$得

$$F_k+(s_i-s_j-L{)}^2+2\ast val\ast (s_i-s_k-L)+va{l}^2$$

比较上述二式，并结合$①$,可知

只要证明$s_i-s_j-L>s_i-s_k-L$即可

$-s_j>-s_k$

由于$s_j<s_k$，则不等式成立。

所以在$i$时，$k$优于$j$，那么在后面状态时，$k$永远优于$j$。

如何淘汰$j$(队头)

$$\begin{gathered} F_j+(s_i-s_j-L{)}^2\ge F_k+(s_i-s_k-L{)}^2 \\ {F}_j+(s_i-L{)}^2-2\ast s_j\ast (s_i-L)+{s_j}^2\ge F_k+(s_i-L{)}^2-2\ast s_k\ast (s_i-L)+{s_k}^2 \\ {F}_j-F_k+{s_j}^2-{s_k}^2\ge 2\ast (s_j-s_k)\ast (s_i-L) \end{gathered}$$

由$s_k>s_j$可知，

$$\frac{F_j-F_k+{s_j}^2-{s_k}^2}{s_j-s_k}\le 2\ast (s_i-L)$$

设$x_1=F_j+{s_j}^2,x_2=F_k+{s_k}^2,y_1=s_j,y_2=s_k$

则$$calc(j,k)=\frac{x_1-x_2}{y_1-y_2}\le 2\ast (s_i-L)$$

这句话的意思就是$k$比$j$更优，所以$j$(先进队列)就没有存在的必要，要踢去。

由于$s_i$随$i$增大而增大，那么$calc(q_{head},q_{head+1})$随$head$增大而增大，才符合$head$为最优解，因此斜率是不断递增的。

踢队尾

我们首先一定要确定$i$这个状态，要优于$q_{tail}$

要维护的队列由于$q_{tail-1}$一定优于$q_{tail}$

现有一个标准值$2\ast (s_t-L)(tail<t)$

$calc(q_{tail-1},q_{tail)}\ge 2\ast (s_t-L)$

如果$q_{tail}$要优于$i$，

则$calc(q_{tail},i)\ge calc(q_{tail-1},q_{tail})$

才能确保$q_{tail}$优于$i$

反之，当$calc(q_{tail},i)\le calc(q_{tail-1},q_{tail})$时，需要删去队尾。

貌似就没了哦。

最后几句话

本题维护的是一个下凸壳(我的做法是这样的)，所以斜率优化真是一个神奇的东西。

AC code

#pragma GCC optimize("Ofast")
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
#define gc getchar()
#define ll long long
using namespace std;
const int N=5e4+10;
inline void qr(ll &x)
{
	x=0;int f=1;char c=gc;
	while(c<'0'||c>'9'){if(c=='-')f=-1;c=gc;}
	while(c>='0'&&c<='9'){x=x*10+(c^48);c=gc;}
	x*=f;
}
void qw(ll x)
{
	if(x<0)x=-x,putchar('-');
	if(x/10)qw(x/10);
	putchar(x%10+48);
}
ll s[N],f[N];
int q[N],l,r;
double calc(int j,int k)
{
	return (f[j]-f[k]+s[j]*s[j]-s[k]*s[k])/(double)(s[j]-s[k]);
}
int n,L;
int main()
{
	scanf("%d%d",&n,&L);++L;
	for(int i=1;i<=n;i++)qr(s[i]),s[i]+=s[i-1]+1;
	l=1;r=1;q[1]=0;
	for(int i=1;i<=n;i++)
	{
		while(l<r&&calc(q[l],q[l+1])<=2.0*(s[i]-L))++l;
		f[i]=f[q[l]]+(s[i]-s[q[l]]-L)*(s[i]-s[q[l]]-L);
		while(l<r&&calc(q[r],i)<=calc(q[r-1],q[r]))--r;
		q[++r]=i;
	}
	qw(f[n]);puts("");
	return 0;
}