【动态规划18】bzoj1010 [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。但他希望费用最小.

输出输出格式

第一行输入两个整数N，L.接下来N行输入Ci.1<=N<=50000,1<=L,Ci<=10^7
输出最小费用

O(n^2)的动规方程依然显然。
$f[i]=min(f[j]+(i-j-1+pre[i]-pre[j]-L)^2)$$(j<i)$
设$k<j<i$，且从j转移比从k转移情况更优
$f[j]+(i-j-1+pre[i]-pre[j]-L)^2<f[k]+(i-k-1+pre[i]-pre[k]-L)^2$
好像可以优化？不过蒟蒻不是很会，强行刚。
然而似乎刚不动，我们还是考虑一下优化..
我们可以发现原式子可以变为$f[j]+(pre[i]+i)-(pre[j]+j)-(L+1)$
那我们处理pre的时候再将pre[x]的值变为pre[x]+x，并将L++
则原式变为了$f[j]+(pre[i]-pre[j]-L)^2<f[k]+(pre[i]-pre[k]-L)^2$
那么这个时候强刚还是可以做到的。
$f[j]+pre[i]^2-2*pre[i]*(pre[j]+L)+(pre[j]+L)^2<=f[k]+pre[i]^2-2*pre[i]*(pre[k]+L)+(pre[k]+L)^2$
移项，乱搞之后
$pre[i]>=$${f[j]+(pre[j]+L)^2-f[k]-(pre[k]+L)^2)}\over{2*(pre[j]-pre[k])}$
由于Ci>=1，所以不用考虑pre[j]与pre[k]相同的情况

#include<bits/stdc++.h>
#define fer(i,j,n) for(int i=j;i<=n;i++)
#define far(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
const int maxn=100010;
const int INF=1e9+7;
using namespace std;
/*----------------------------------------------------------------------------*/
inline ll read()
{
    char ls;ll x=0,sng=1;
    for(;ls<'0'||ls>'9';ls=getchar())if(ls=='-')sng=-1;
    for(;ls>='0'&&ls<='9';ls=getchar())x=x*10+ls-'0';
    return x*sng;
}
/*----------------------------------------------------------------------------*/
int n,L,h,t;
ll pre[maxn],f[maxn],q[maxn];
inline ll calcK(int k,int j)
{
    return (f[j]+(pre[j]+L)*(pre[j]+L)-f[k]-(pre[k]+L)*(pre[k]+L))/(2*(pre[j]-pre[k]));
}
int main()
{
    n=read();L=read()+1;
    fer(i,1,n)pre[i]=read()+pre[i-1]+1;
    h=t=0;
    fer(i,1,n)
    {
        while(h<t&&calcK(q[h],q[h+1])<=pre[i])h++;
        f[i]=f[q[h]]+(pre[i]-pre[q[h]]-L)*(pre[i]-pre[q[h]]-L);
        while(h<t&&calcK(q[t-1],q[t])>calcK(q[t],i))t--;
        q[++t]=i;
    }
    cout<<f[n];
}