DP--斜率优化--BZOJ1096--仓库建设

Description

　　L公司有N个工厂，由高到底分布在一座山上。如图所示，工厂1在山顶，工厂N在山脚。由于这座山处于高原内陆地区（干燥少雨），L公司一般把产品直接堆放在露天，以节省费用。突然有一天，L公司的总裁L先生接到气象部门的电话，被告知三天之后将有一场暴雨，于是L先生决定紧急在某些工厂建立一些仓库以免产品被淋坏。由于地形的不同，在不同工厂建立仓库的费用可能是不同的。第i个工厂目前已有成品Pi件，在第i个工厂位置建立仓库的费用是Ci。对于没有建立仓库的工厂，其产品应被运往其他的仓库进行储藏，而由于L公司产品的对外销售处设置在山脚的工厂N，故产品只能往山下运（即只能运往编号更大的工厂的仓库），当然运送产品也是需要费用的，假设一件产品运送1个单位距离的费用是1。假设建立的仓库容量都都是足够大的，可以容下所有的产品。你将得到以下数据：1：工厂i距离工厂1的距离Xi（其中X1=0）;2：工厂i目前已有成品数量Pi;:3：在工厂i建立仓库的费用Ci;请你帮助L公司寻找一个仓库建设的方案，使得总的费用（建造费用+运输费用）最小。

题解

首先不难写出转移方程

$$f[i]=min\{f[j]+c[i]+a[i]*(p\_sum[i]-p\_sum[j])+sum[j]-sum[i]\}$$
其中，a[i]为距离，p_sum表示p[i]的前缀和，sum表示a[i]*p[i]的前缀和。
显然，这样写的复杂度为 $O(n^2)$，不能承受。可以试试斜率优化：
把转移方程中和j有关的单独整理： $$M=min\{f[j]-a[i]*p\_sum[j]+sum[j]\}$$ 移项 $$f[j]-sum[j]=a[i]*p\_sum[j]+M$$
把f[j]-sum[j]（只和j有关）的看成y,a[i]（只和i有关）看成k,p_sum[j]看成（与i有关的那一项中与j有关的）看成x,M（结果）看成b，就得到了y=kx+b的形式，目的是要让b最小，那么，将一个决策点看成(x,y)的点，最优决策就是以当前的a[i]为斜率的直线从下往上推，遇到的第一个决策点，这样，所有上凸的决策点都是没用的，所以最终有用的决策点一定组成一个下凸壳，由于横坐标（即sum[i]）是递增的，所以每次加入一个决策点时在末尾加入，并且维护下凸性质。在挑选最优决策点时，由于斜率a[i]是递增的，所以每次的最优决策点前的所有决策点都是没用的，可以利用这个性质维护一个队列，队头就是最优决策点。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 1000006
#define LL long long
using namespace std;
LL n,hed,tal,p_x[maxn],p_y[maxn],que[maxn],a[maxn],f[maxn],p_sum[maxn],sum[maxn],c[maxn];
inline char nc(){
    static char buf[100000],*i=buf,*j=buf;
    return i==j&&(j=(i=buf)+fread(buf,1,100000,stdin),i==j)?EOF:*i++;
}
inline LL _read(){
    char ch=nc();LL sum=0;
    while(!(ch>='0'&&ch<='9'))ch=nc();
    while(ch>='0'&&ch<='9')sum=sum*10+ch-48,ch=nc();
    return sum;
}
double cross(int i,int j,int k){
    return (p_x[j]-p_x[i])*(p_y[k]-p_y[j])-(p_x[k]-p_x[j])*(p_y[j]-p_y[i]);
}
int main(){
    freopen("build.in","r",stdin);
    freopen("build.out","w",stdout);
    n=_read();
    for(int i=1,p;i<=n;i++)a[i]=_read(),p=_read(),c[i]=_read(),p_sum[i]=p_sum[i-1]+p,sum[i]=sum[i-1]+a[i]*p;
    f[0]=0;
    for(int i=1;i<=n;i++){
        while(hed<tal&&(p_y[que[hed]]-a[i]*p_x[que[hed]]>p_y[que[hed+1]]-a[i]*p_x[que[hed+1]]))hed++;
        f[i]=p_y[que[hed]]-a[i]*p_x[que[hed]]+c[i]+a[i]*p_sum[i]-sum[i];
        p_x[i]=p_sum[i];p_y[i]=f[i]+sum[i];
        while(hed<tal&&(cross(que[tal-1],que[tal],i)<=0))tal--;
        que[++tal]=i;
    }
    printf("%lld",f[n]);
    return 0;
}