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(n2) ,不能承受。可以试试斜率优化:
把转移方程中和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;
}