[斜率优化]BZOJ 1010——玩具装箱toy

1010: [HNOI2008]玩具装箱toy

题目描述

　　P教授要去看奥运，但是他舍不下他的玩具，于是他决定把所有的玩具运到北京。他使用自己的压缩器进行压缩，其可以将任意物品变成一堆，再放到一种特殊的一维容器中。P教授有编号为1…N的N件玩具，第i件玩具经过压缩后变成一维长度为Ci.为了方便整理，P教授要求在一个一维容器中的玩具编号是连续的。同时如果一个一维容器中有多个玩具，那么两件玩具之间要加入一个单位长度的填充物，形式地说如果将第i件玩具到第j个玩具放到一个容器中，那么容器的长度将为 x=j-i+Sigma(Ck) i<=K<=j 制作容器的费用与容器的长度有关，根据教授研究，如果容器长度为x,其制作费用为(X-L)^2.其中L是一个常量。P教授不关心容器的数目，他可以制作出任意长度的容器，甚至超过L。但他希望费用最小.

解题思路

根据题目描述不难想到是DP，很快就会推出方程

$$f[i]=min(f[j]+(i-j-1-L+s[i]-s[j])^2)$$
进一步化简 $$f[i]=min(f[j]+((i-1-L+s[i])+(-j-s[j]))^2)$$
把无关项移出来 $$f[i]=min(f[j]-2*(i-1-L+s[i])(j+s[j])+(j+s[j])^2)+(i-1-L+s[i])^2$$

设$k=(i-1-L+s[i]),x=(j+s[j]),y=(j+s[j])^2+f[j]$
通过$b=-kx+y,y=kx+b$

接下来套用斜率优化就可以了。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
using namespace std;
const int maxn=50005;
LL s[maxn],f[maxn];
int n,L,now,til;
struct jz{
    LL x,y;
    jz(LL a=0,LL b=0):x(a),y(b){}
    jz operator-(const jz &b){return jz(x-b.x,y-b.y);}
}a[maxn];
inline int _read(){
    int num=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') num=num*10+ch-48,ch=getchar();
    return num;
}
LL getb(jz x,LL k){return -k*x.x+x.y;}
LL check(jz x,jz y){return x.x*y.y-x.y*y.x;}
int main(){
    freopen("exam.in","r",stdin);
    freopen("exam.out","w",stdout);
    n=_read(),L=_read();
    for (int i=1;i<=n;i++) s[i]=s[i-1]+_read();
    memset(f,63,sizeof(63));
    f[0]=0;now=til=1;a[1]=jz(0,0);
    for (int i=1;i<=n;i++){
        LL w=s[i]+i-L-1,k=2*w;
        while (now<til&&getb(a[now+1],k)<getb(a[now],k)) now++;
        f[i]=getb(a[now],k)+w*w;
        jz x(s[i]+i,f[i]+(s[i]+i)*(s[i]+i));
        while (til>1&&check(a[til]-a[til-1],x-a[til-1])<=0) til--;
        now=min(now,til);a[++til]=x;
    }
    printf("%lld\n",f[n]);
    return 0;
}