博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ3437:小P的牧场(斜率优化DP)
阅读量:6193 次
发布时间:2019-06-21

本文共 1781 字,大约阅读时间需要 5 分钟。

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
样例解释
选取牧场1,3,4建立控制站,最小费用为2+(2+1*1)+4=9。
1<=n<=1000000, 0 < a i ,bi < = 1000

Source

Solution

1A就很开心了

非常裸的一眼斜率优化DP

设$sum_i$为$b_i$的后缀和,$s_i$为$(n-i)*b_i$的后缀和

那么DP式子很明显可以写为

$f_i=f_j+s_{j+1}-s_i-(n-i)*(sum_{j+1}-sum_{i+1})+a_i$

感性理解一下其实式很子好懂

把式子拆开后发现满足斜率单调性,直接上斜率优化就好了

Code

1 #include
2 #include
3 #include
4 #define N (1000000+1000) 5 #define LL long long 6 using namespace std; 7 8 LL n,a[N],b[N],sum[N],s[N],q[N],f[N]; 9 10 LL K(LL j){
return sum[j+1];}11 LL B(LL j){
return f[j]+s[j+1]-n*sum[j+1];}12 LL Y(LL i,LL j){
return K(j)*i+B(j);}13 14 bool cover(LL x1,LL x2,LL x3)15 {16 LL w1=(B(x2)-B(x1))*(K(x1)-K(x3));17 LL w2=(B(x3)-B(x1))*(K(x1)-K(x2));18 return w1<=w2;19 }20 21 int main()22 {23 scanf("%lld",&n);24 for (int i=1; i<=n; ++i) scanf("%lld",&a[i]);25 for (int i=1; i<=n; ++i) scanf("%lld",&b[i]);26 for (int i=n; i>=1; --i)27 {28 sum[i]=sum[i+1]+b[i];29 s[i]=s[i+1]+(n-i)*b[i];30 }31 int head=0,tail=0;32 for (int i=1; i<=n; ++i)33 {34 while (head
=Y(i,q[head+1])) head++;35 f[i]=Y(i,q[head])-s[i+1]+n*sum[i+1]-i*sum[i+1]+a[i];36 while (head

转载于:https://www.cnblogs.com/refun/p/9279361.html

你可能感兴趣的文章
Linux 下 项目管理
查看>>
Steve Souders 高性能网站建设指南14条提升网站页面速度的最佳实践(1)
查看>>
ORA-00932: 数据类型不一致: 应为 NUMBER, 但却获得 BINARY
查看>>
YMP开发框架快速上手(一)
查看>>
Oracle11g RAC下ASM 的管理与维护
查看>>
centos 安装ftp
查看>>
android 简单的 左右上下手势 判断
查看>>
CSS重置, 批量设置指定所有类型控件的CSS风格
查看>>
Go 1.4 src/pkg → src
查看>>
成为JavaGC专家Part I — 深入浅出Java垃圾回收机制
查看>>
linux 配置软连接的需要注意的一个问题
查看>>
如何对system.img和userdata.img解包,再重新打包
查看>>
都是 HBase 上的 SQL 引擎,Kylin 和 Phoenix 有什么不同?
查看>>
使用jMeter构造大量并发HTTP请求进行微服务性能测试
查看>>
实现RabbitMQ 延时消息
查看>>
系统目录结构、ls命令、文件类型、alias命令
查看>>
京东八年架构师: Redis 如何分布式,金融的设计原理
查看>>
完美Python入门基础知识点总结,看完你的Python就完全入门了!
查看>>
hadoop副本放置策略
查看>>
FastDFS-图片服务器
查看>>