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

3437: 小P的牧场

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 1493 Solved: 830
[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

2 4 2 4

3 1 4 2

Sample Output

很简单的斜率DP

首先最后一个牧场一定要放控制站，先假设其它地方都不放，求出价格

dp[i]表示在i处建站能节省的最多费用，有

sum[i] = ∑b[i]

dp[i] = max(dp[j]+sum[i]*(j-i), j∈(i, n] )-a[i]

ans = ∑b[i]*(n-i)-max(dp[i], i∈[1, n] )

#include<stdio.h>
#include<algorithm>
using namespace std;
#define LL long long
LL a[1000005], b[1000005], sum[1000005], dp[1000005], st[1000005];
double Jud(int j, int k)
{return 1.0*(dp[j]-dp[k])/(k-j);
}
int main(void)
{LL i, n, ans, bet, l, r;scanf("%lld", &n);for(i=1;i<=n;i++)scanf("%lld", &a[i]);for(i=1;i<=n;i++){scanf("%lld", &b[i]);sum[i] = sum[i-1]+b[i];}ans = a[n];for(i=1;i<=n;i++)ans += b[i]*(n-i);l = 1, r = 1;st[1] = n;bet = 0;for(i=n-1;i>=1;i--){while(r>l && Jud(st[l], st[l+1])>sum[i])l++;dp[i] = dp[st[l]]+sum[i]*(st[l]-i)-a[i];bet = max(bet, dp[i]);while(r>l && Jud(st[r], i)>Jud(st[r-1], st[r]))r--;st[++r] = i;}printf("%lld\n", ans-bet);return 0;
}