【bzoj3437】【小p的牧场】【斜率优化dp】
2016-03-23 19:37
387 查看
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
42424
3142
Sample Output
9样例解释
选取牧场1,3,4建立控制站,最小费用为2+(2+1*1)+4=9。
数据范围与约定
对于100%的数据,1<=n<=1000000,0<ai,bi<=10000
题解:同bzoj1096.
代码:
#include<iostream> #include<cstdio> #include<cstring> #define N 1000010 using namespace std; int q ,n,l,r; long long s ,a ,p ,sp ,c ,f ; long long G(int a,int b){return f[b]-f[a]+sp[b]-sp[a];} long long S(int a,int b){return p[b]-p[a];} double work(int a,int b){return G(a,b)*1.0/S(a,b)*1.0;} int main(){ scanf("%d",&n); for (int i=1;i<=n;i++) scanf("%lld",&c[i]); for (int i=1;i<=n;i++) scanf("%lld",&a[i]); for (int i=1;i<=n;i++) p[i]=p[i-1]+a[i]; s[1]=0;for (int i=2;i<=n;i++) s[i]=s[i-1]+1; for (int i=1;i<=n;i++) sp[i]=sp[i-1]+s[i]*a[i]; for (int i=1;i<=n;i++){ while (l<r&&work(q[l],q[l+1])<s[i]) l++; f[i]=f[q[l]]+s[i]*(p[i]-p[q[l]])-(sp[i]-sp[q[l]])+c[i]; while (l<r&&work(q[r-1],q[r])>work(q[r],i)) r--; q[++r]=i; } cout<<f <<endl; }
相关文章推荐
- consul系列二之集群
- const int * 与 int *const
- 如何省略显示过多的child views
- Struts2中 Result类型配置详解
- 为什么std::cin>>int input_value是不合法的
- 共享数据的包含const
- Common Lisp 笔记(1)
- jQuery学习笔记(一)
- linux中libevent-2.0.so.5 not found 解决 ldd 方法
- 第4周项目5(4) 用递归方法求解
- Delphi 中 FindWindow 和 FindWindowEx 找到外部进程,然后发送消息(比如最大化)
- 记录Hibernate的缓存知识
- Linux内核分析第五周学习总结:扒开系统调用的三层皮(下)
- Linux内核分析作业 NO.5
- 比特承诺的概念
- 第4周项目3:随机数函数应用于游戏-猜数字游戏
- fast-rcnn配置运行demo.py(Ubuntu14.04),不成功的朋友请与我(lee)联系,后面附带邮箱
- 张量的展开与matlab下的工具包操作
- Linux内核分析第五章读书笔记
- 关闭磁感应,光感,,