仓库建设 [Bzoj 1096]

题目地址请点击——

---

仓库建设

---

Description

L 公司有 $N$ 个工厂，由高到底分布在一座山上。
如图所示，工厂 $1$ 在山顶，工厂 $N$ 在山脚。

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

---

Input

请你帮助L 公司寻找一个仓库建设的方案，使得总的费用（建造费用+运输费用）最小。
输入文件storage.in 第一行包含一个整数 $N$，表示工厂的个数。
接下来 $N$ 行每行包含两个整数 $X_i$, $P_i$,$C_i$, 意义如题中所述。

---

Output

输出文件storage.out 仅包含一个整数，为可以找到最优方案的费用。

---

Sample Input

3
0 5 10
5 3 100
9 6 10

---

Sample Output

32

---

Data Size

对于 $100\%$ 的数据， $N ≤1000000$。
所有的 $X_i$, $P_i$ , $C_i$ 均在 $32$ 位带符号整数 以内，保证中间计算结果不超过 $64$ 位带符号整数。

---

Solution

设 $f_i$ 表示在 $i$ 建仓库，解决 $1$~$i$ 的最小费用。

$$f_i=Max_{k=1}^{i-1}\{f_k+\sum_{j=k+1}^{i-1}(x_i-x_j)·p_j\}+c_i$$
用斜率优化和前缀和优化即可将复杂度降到 $\mathcal O(n)$

---

Code

#include <iostream>
#include <cstdio>
#include <stdint.h>

#define Min(x,y) ((x)<(y)?(x):(y))
#define LL long long

using namespace std;

LL m[1000010],y[1000010];
LL p[1000010],s[1000010];
double ks[1000010];
LL f[1000010],c[1000010];

LL n,Max_Right=0,lastcut=0;
const LL oo=9223372036854775807LL;

int main(){

    scanf("%lld",&n);
    for(int i=1;i<=n;i++){
        scanf("%lld%lld%lld",&s[i],&p[i],&c[i]);
        m[i]=m[i-1]+p[i];
        y[i]=y[i-1]+p[i]*s[i];
    }

    ks[Max_Right+1]=oo;ks[0]=-oo;

    for(int i=1;i<=n;i++){
        for(int j=lastcut;j<=Max_Right;j++)
            if(ks[j]<=s[i]&&ks[j+1]>=s[i]){
                lastcut=j;
                f[i]=f[j]+s[i]*m[i-1]-s[i]*m[j]-y[i-1]+y[j]+c[i];
                for(int k=Max_Right;k>=0;k--)
                    if((double)(f[i]+y[i]-(f[k]+y[k]))/(m[i]-m[k])>ks[k]){
                        Max_Right=k+1;
                        lastcut=Min(lastcut,Max_Right);
                        ks[Max_Right]=(double)(f[i]+y[i]-(f[k]+y[k]))/(m[i]-m[k]);
                        f[Max_Right]=f[i];
                        s[Max_Right]=s[i];
                        m[Max_Right]=m[i];
                        y[Max_Right]=y[i];
                        ks[Max_Right+1]=oo;
                        break;
                    }
                break;
            }
    }

    printf("%lld\n",f[n]);

    return 0;
}