BZOJ3437 小P的牧场【斜率优化dp】

3437: 小P的牧场

Time Limit: 10 Sec Memory Limit: 128 MB

Submit: 1502 Solved: 836

[Submit][Status][Discuss]

Description

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

Input

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

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

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

Output

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

Sample Input

4

2424

3142

Sample Output

9

样例解释

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

1<=n<=1000000, 0 < a i ,bi < = 10000

与P1096仓库建设很像，不过路径权值变成了1

我们用一个比较牛的前缀和:s[i]为b[i]前缀和，c[i]为b[i]前缀和

那么从j + 1全部搬到i的代价就是(s[i] - s[j]) * i - (c[i] - c[j])【想象一下】

那么我们设f[i]表示在i建厂的最小代价

f[i] = min{f[j] + (s[i] - s[j]) * i - (c[i] - c[j])} + A[i]

去掉常量化简得到(f[j] + c[j]) = i * s[i] + f[i]

我们就的到了y = i * x + f[j]这样的直线求截距最大，维护凸包就好了

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#define eps 1e-9

#define LL long long int

#define REP(i,n) for (int i = 1; i <= (n); i++)

#define fo(i,x,y) for (int i = (x); i <= (y); i++)

#define Redge(u) for (int k = head[u]; k != -1; k = edge[k].next)

using namespace std;

const int maxn = 1000005,maxm = 100005,INF = 1000000000;

inline LL read(){

	LL out = 0,flag = 1;char c = getchar();

	while (c < 48 || c > 57) {if (c == '-') flag = -1; c = getchar();}

	while (c >= 48 && c <= 57) {out = out * 10 + c - 48; c = getchar();}

	return out * flag;

}

LL n,A[maxn],B[maxn],s[maxn],c[maxn],f[maxn],q[maxn],head,tail;

inline double slope(int u,int v){

	return (double)(f[u] + c[u] - f[v] - c[v]) / (s[u] - s[v]);

}

inline LL getf(int i,int j){

	return f[j] + (s[i] - s[j]) * i - (c[i] - c[j]) + A[i];

}

int main()

{

	n = read();

	REP(i,n) A[i] = read();

	REP(i,n) B[i] = read(),s[i] = s[i - 1] + B[i],c[i] = c[i - 1] + B[i] * i;

	head = tail = 0;

	for (int i = 1; i <= n; i++){

		while (head < tail && slope(q[head],q[head + 1]) < i + eps) head++;

		f[i] = getf(i,q[head]);

		while (head < tail && slope(q[tail],q[tail - 1]) + eps > slope(i,q[tail])) tail--;

		q[++tail] = i;

	}

	cout<<f[n]<<endl;

	return 0;

}