BZOJ3437:小P的牧场(斜率优化DP)
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
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
样例解释
选取牧场1,3,4建立控制站,最小费用为2+(2+1*1)+4=9。
1<=n<=1000000, 0 < a i ,bi < = 1000
Source
Solution
1A就很开心了
非常裸的一眼斜率优化DP
设sumi为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<iostream> 2 #include<cstring> 3 #include<cstdio> 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<tail && Y(i,q[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<tail && cover(i,q[tail],q[tail-1])) tail--; 37 q[++tail]=i; 38 } 39 printf("%lld",f[n]); 40 }
标签:
动态规划
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 聊一聊 操作系统蓝屏 c0000102 的故障分析
· SQL Server 内存占用高分析
· .NET Core GC计划阶段(plan_phase)底层原理浅谈
· .NET开发智能桌面机器人:用.NET IoT库编写驱动控制两个屏幕
· 用纯.NET开发并制作一个智能桌面机器人:从.NET IoT入门开始
· 我干了两个月的大项目,开源了!
· 推荐一款非常好用的在线 SSH 管理工具
· 聊一聊 操作系统蓝屏 c0000102 的故障分析
· 千万级的大表,如何做性能调优?
· .NET周刊【1月第1期 2025-01-05】