仓库建设
Description
L 公司有
N
个工厂,由高到底分布在一座山上。
如图所示,工厂
由于这座山处于高原内陆地区(干燥少雨),L 公司一般把产品直接堆放在露天,以节省费用。突然有
一天,L 公司的总裁L 先生接到气象部门的电话,被告知三天之后将有一场暴雨,于是L 先生决定紧急在
某些工厂建立一些仓库以免产品被淋坏。
由于地形的不同,在不同工厂建立仓库的费用可能是不同的。第
对于没有建立仓库的工厂,其产品应被运往其他的仓库进行储藏,而由于L 公司产品的对外销售处设置在山脚的工厂
N
,故产品只能往山下运(即只能运往编号更大的工厂的仓库),当然运送产品也是需要费用的,假设一件产品运送
假设建立的仓库容量都都是足够大的,可以容下所有的产品。
你将得到以下数据:
工厂
工厂
i
目前已有成品数量
在工厂
i
建立仓库的费用
Input
请你帮助L 公司寻找一个仓库建设的方案,使得总的费用(建造费用+运输费用)最小。
输入文件storage.in 第一行包含一个整数
N
,表示工厂的个数。
接下来
Output
输出文件storage.out 仅包含一个整数,为可以找到最优方案的费用。
Sample Input
3
0 5 10
5 3 100
9 6 10
Sample Output
32
Data Size
对于
100%
的数据,
N≤1000000
。
所有的
Xi
,
Pi
,
Ci
均在
32
位带符号整数 以内,保证中间计算结果不超过
64
位带符号整数。
Solution
设
fi
表示在
i
建仓库,解决
用斜率优化和前缀和优化即可将复杂度降到
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;
}