bzoj 3437: 小P的牧场（斜率优化DP）

3437: 小P的牧场

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 994  Solved: 566
[Submit][Status][Discuss]

Description

 背景

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

 描述

    小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。

数据范围与约定


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

HINT

Source

KpmCup#0 By Greens

[Submit][Status][Discuss]

#include<iostream>#include<cstdio>#include<algorithm>#include<cstring>#include<cmath>#define N 1000003#define LL long long using namespace std;LL f[N],a[N],b[N],sum[N],sum1[N],c[N];int head,tail,q[N];int n;/*LL calc(int x,int y){if (x>y) return 0;LL ans=0;for (int i=x;i<=y;i++) ans+=b[i]*(LL)(y-i+1);return ans;}*/LL calc(int j,int i){return -sum1[j]*(LL)i+(LL)n*sum1[j]-sum[j]+f[j];}LL K(int i){return -sum1[i];}LL B(int j){return (LL)n*sum1[j]-sum[j]+f[j];}int pd(int x1,int x2,int x3){LL w1=(K(x3)-K(x2))*(B(x3)-B(x1));LL w2=(K(x1)-K(x3))*(B(x2)-B(x3));return w1<=w2;}int main(){freopen("a.in","r",stdin);scanf("%d",&n);for (int i=1;i<=n;i++) scanf("%I64d",&a[i]);for (int i=1;i<=n;i++) scanf("%I64d",&b[i]);for (int i=1;i<=n;i++) sum1[i]=sum1[i-1]+b[i];for (int i=1;i<=n;i++) c[i]=b[i]*(LL)(n-i);for (int i=1;i<=n;i++) sum[i]=sum[i-1]+c[i];//for (int i=1;i<=n;i++) cout<<sum[i]<<" "<<sum1[i]<<endl;memset(f,127/3,sizeof(f));f[0]=0;/*for (int i=1;i<=n;i++)    for (int j=0;j<i;j++)   f[i]=min(f[i],f[j]+calc(j+1,i-1)+a[i]);*/head=tail=1; q[head]=0;for (int i=1;i<=n;i++){while (head<tail&&calc(q[head],i)>calc(q[head+1],i)) head++;//cout<<calc(q[head],i)<<" "<<a[i]+sum[i-1]-(LL)n*sum1[i-1]+(LL)i*sum1[i-1]<<endl;f[i]=calc(q[head],i)+a[i]+sum[i-1]-(LL)n*sum1[i-1]+(LL)i*sum1[i-1];//cout<<f[i]<<" "<<q[head]<<endl;while (head<tail&&pd(q[tail-1],q[tail],i)) tail--;q[++tail]=i;}printf("%I64d\n",f[n]);}