小P 的牧场(pasture)题解

题目描述

小P是个特么喜欢玩MC 的孩纸。。。

小P在MC 里有n 个牧场，自西向东呈一字形排列(自西向东用1…n 编号)，于是他就烦恼了：为了控制这n 个牧场，他需要在某些牧场上面建立控制站，每个牧场上只能建立一个控制站，每个控制站控制的牧场是它所在的牧场一直到它西边第一个控制站的所有牧场(它西边第一个控制站所在的牧场不被控制)(如果它西边不存在控制站，那么它控制西边所有的牧场)，每个牧场被控制都需要一定的花费(毕竟在控制站到牧场间修建道路是需要资源的嘛~)，而且该花费等于它到控制它的控制站之间的牧场数目(不包括自身，但包括控制站所在牧场)乘上该牧场的放养量，在第i 个牧场建立控制站的花费是ai，每个牧场i 的放养量是bi，理所当然，小P需要总花费最小，但是小P的智商有点不够用了，所以这个最小总花费就由你来算出啦。

输入

第一行一个整数n 表示牧场数目

第二行包括n个整数，第i个整数表示ai

第三行包括n个整数，第i个整数表示bi

输出

只有一行，包括一个整数，表示最小花费

样例输入

4

2 4 2 4

3 1 4 2

样例输出

9

提示

【样例解释】

选取牧场1，3，4 建立控制站，最小费用为2+( 2 + 1 * 1 ) + 4 = 9。

【数据范围与约定】

对于10%的数据，1 <= n <= 10

对于30%的数据，1 <= n <= 1000

对于100%的数据，1 <= n <= 1000000 , 0 < ai,bi <= 10000

想法

同仓库建设（DP+斜率优化）

算法

见代码

代码

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <string>
#include <cmath>
#define MAXN 1000005
using namespace std;
typedef long long LL;
int n,Q[MAXN],l,r;
LL F[MAXN],X[MAXN],P[MAXN],C[MAXN],S[MAXN],B[MAXN];
inline double cal(int x,int y)
{
return double(F[x]-F[y]+B[x]-B[y])/double(S[x]-S[y]);
}
int main()
{
//freopen("pasture.in","r",stdin);
//freopen("pasture.out","w",stdout);
scanf("%lld",&n);
for (int i=1;i<=n;i++)scanf("%d",&C[i]);
for (int i=1;i<=n;i++)
{
scanf("%lld",&P[i]);
X[i]=i-1;
S[i]=S[i-1]+P[i];
B[i]=B[i-1]+P[i]*X[i];
}
for (int i=1;i<=n;i++)
{
while(l<r&&cal(Q[l+1],Q[l])<=X[i])++l;
int num=Q[l];
F[i]=F[num]+(S[i-1]-S[num])*X[i]-B[i-1]+B[num]+C[i];
while(l<r&&cal(Q[r],Q[r-1])>=cal(i,Q[r]))--r;
Q[++r]=i;
}
printf("%lld",F
);
return 0;
}