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;
}