BZOJ3437 小p的牧场

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
2 4 2 4
3 1 4 2

Sample Output

9

Hint

【样例解释】
　　选取牧场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；

一类斜率优化问题，x单增，斜率单增，注意数形结合。

#include<bits/stdc++.h>
using namespace std;
#define Inc(i,L,r) for(register int i=(L);i<=(r);++i)
const int N = 1e6+10;
#define int long long
int n,a[N],b[N];
int s[N],ks[N];
inline void init(){
	scanf("%lld",&n);
	Inc(i,1,n)scanf("%lld",&a[i]);
	Inc(i,1,n)scanf("%lld",&b[i]);
	
	Inc(i,1,n)s[i]=s[i-1]+b[i];
	Inc(i,1,n)ks[i]=ks[i-1]+b[i]*i;
}
int f[N],q[N];
inline int y(int j){
	return f[j]+ks[j];
}
inline double slope(int j1,int j2){
	return (y(j1)-y(j2)-0.0)/(s[j1]-s[j2]-0.0);
}
inline void dp(){
	int L=1,r=0;
	q[++r]=0;
	Inc(i,1,n){
		while(L<r&&slope(q[L],q[L+1])<=i)++L;
		int j=q[L];f[i]=a[i]+f[j]+i*(s[i-1]-s[j])-(ks[i-1]-ks[j]);
		while(L<r&&slope(q[r-1],q[r])>=slope(q[r],i))--r;
		q[++r]=i;
//		Inc(j,0,i-1){
//			f[i]=min(f[i],a[i]+f[j]+i*(s[i-1]-s[j])-(ks[i-1]-ks[j]));
//		}
	}
	cout<<f[n]<<"\n";
}
signed main(){
	init();
	dp();
	return 0;
}