【bzoj3437】小P的牧场

3437: 小P的牧场

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

【题解】

f[i]表示前i个牧场被控制的最小代价，sumb表示b数组的前缀和，sumkb表示 b[i]*i 的前缀和

则状态转移方程：f[i]=min{f[j]+a[i]+i*(sumb[i]-sumb[j])-(sumkb[i]-sumkb[j])} (1<=j<i)

设j<k，且从k转移到i比j更优。

得斜率表达式：[(f[j]-sumkb[j])-(f[k]-sumkb[k])]/(sumb[j]-sumb[k])>i

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstdlib>
 4 #include<cstring>
 5 #include<ctime>
 6 #include<cmath>
 7 #include<algorithm>
 8 using namespace std;
 9 typedef long long ll;
10 #define FILE "read"
11 #define MAXN 1000100
12 #define up(i,j,n) for(ll i=j;i<=n;i++)
13 namespace INIT{
14     char buf[1<<15],*fs,*ft;
15     inline char getc(){return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;}
16     inline ll read(){
17         ll x=0,f=1;  char ch=getc();
18         while(!isdigit(ch))  {if(ch=='-')  f=-1;  ch=getc();}
19         while(isdigit(ch))  {x=x*10+ch-'0';  ch=getc();}
20         return x*f;
21     }
22 }using namespace INIT;
23 ll n,head,tail,a[MAXN],b[MAXN],sumb[MAXN],sumkb[MAXN],q[MAXN],f[MAXN];
24 void init(){
25     n=read();
26     up(i,1,n) a[i]=read();
27     up(i,1,n) b[i]=read(),sumb[i]=sumb[i-1]+b[i],sumkb[i]=sumkb[i-1]+b[i]*i;
28 }
29 inline double slop(ll j,ll k) {return (double)((f[j]+sumkb[j])-(f[k]+sumkb[k]))/(double)(sumb[j]-sumb[k]);}
30 void solve(){
31     up(i,1,n){
32         while(head<tail&&slop(q[head],q[head+1])<i)  head++;
33         ll t=q[head];
34         f[i]=f[t]+a[i]+i*(sumb[i]-sumb[t])-(sumkb[i]-sumkb[t]);
35         while(head<tail&&slop(q[tail-1],q[tail])>slop(q[tail],i))  tail--;
36         q[++tail]=i;
37     }
38     printf("%lld
",f[n]);
39 }
40 int main(){
41     freopen(FILE".in","r",stdin);
42     freopen(FILE".out","w",stdout);
43     init();
44     solve();
45     return 0;
46 }