小P的牧场
【题目描述】
背景:小P 是个特么喜欢玩MC 的孩纸。。。
小P 在MC 里有n 个牧场,自西向东呈一字形排列(自西向东用1…n 编号),于是他就烦恼了:为了控制这n 个牧场,他需要在某些牧场上面建立控制站,每个牧场上只能建立一个控制站,每个控制站控制的牧场是它所在的牧场一直到它西边第一个控制站的所有牧场(它西边第一个控制站所在的牧场不被控制)(如果它西边不存在控制站,那么它控制西边所有的牧场),每个牧场被控制都需要一定的花费(毕竟在控制站到牧场间修建道路是需要资源的嘛~),而且该花费等于它到控制它的控制站之间的牧场数目(不包括自身,但包括控制站所在牧场)乘上该牧场的放养量,在第i 个牧场建立控制站的花费是ai,每个牧场i 的放养量是bi,理所当然,小P 需要总花费最小,但是小P的智商有点不够用了,所以这个最小总花费就由你来算出啦。
【输入格式】
第一行一个整数n 表示牧场数目
第二行包括n 个整数,第i 个整数表示ai
第三行包括n 个整数,第i 个整数表示bi
【输出格式】
只有一行,包括一个整数,表示最小花费
【样例输入】
4
2 4 2 4
3 1 4 2
【样例输出】
9
【样例解释】
选取牧场1,3,4 建立控制站,最小费用为2+( 2 + 1 * 1 ) + 4 = 9。
数据范围与约定
对于10%的数据,1 <= n <= 10
对于30%的数据,1 <= n <= 1000
对于100%的数据,1 <= n <= 1000000 , 0 < ai,bi <= 10000

分析:

经过思考可以得到这样一个式子

f[i,1]表示在第i个牧场建站且右边的牧场已全部控制的最小费用,f[i,0]表示第i个牧场未建站而从i到n的牧场都被控制的最小代价,倒着来进行处理。

f[i,1]:=min(f[i+1,1],f[i+1,0])+a[i];

f[i,0]:=min(f[j,1]+t); (i:=n..1,i<j<n)

其中t表示i..j-1的牧场被j控制的总代价,在动规循环同时求出,这样要用三层循环,显然这个效率并不好。

对DP进行改进,将牧场排列倒过来,同时改进f,f[i]表示在第i个牧场建站并之前的牧场都被控制的最小总代价,我们增加一个牧场n+1,令其建站和被控制代价为0。则有:

f[i]:=min(f[j]+t)+a[i];(i:=1..n+1,1<=j<i) 其中t表示j+1..i-1的牧场都被控制的总代价,求t可以采用一个很好地方法:

令s[i]表示1..i的a数组元素和,u[i]表示1..i的s数组元素和:

t=s[i-1]-s[j]+s[i-1]-s[j+1]+...+s[i-1]-s[i-2]

=s[i-1]*(i-j-1)-(s[j]+s[j+1]+s[j+2]+..+s[i-2])

=s[i-1]*(i-j-1)-(u[i-2]-u[j-1])

于是得到:

f[i]:=min(f[j]+s[i-1]*(i-j-1)-(u[i-2]-u[j-1]))+a[i];

时间效率为O(n^2),这样还不够,于是进行斜率优化,要注意向队列添加元素时队内最少元素数目不低于2个,我就在这里处理错了导致之前一直没AC。

代码:

program pasture;
var
f:array[..]of int64;
a,b,u,s,q:array[..]of int64;
n,m,t,k,h,j:int64; i:longint;
function min(x,y:int64):int64;
begin
if x<y then min:=x else min:=y;
end;
function cal1(x,y:longint):int64;
begin
exit(f[x]-f[y]+u[x-]-u[y-]);
end;
function cal2(x,y:longint):int64;
begin
exit(x-y);
end;
begin
readln(n);
for i:=n downto do
read(a[i]);
readln;
for i:=n downto do
read(b[i]);
for i:= to n do
s[i]:=s[i-]+b[i];
for i:= to n do
u[i]:=u[i-]+s[i];
f[]:=a[]; q[]:=; h:=; t:=;
for i:= to n+ do
begin
while (h<t)and(cal1(q[h+],q[h])<s[i-]*cal2(q[h+],q[h])) do inc(h);
j:=q[h];
f[i]:=f[j]+s[i-]*(i-j-)-(u[i-]-u[j-])+a[i];
while (t-h>=)and(cal1(i,q[t])*cal2(q[t],q[t-])<=cal1(q[t],q[t-])*cal2(i,q[t])) do dec(t);
t:=t+; q[t]:=i;
end;
writeln(f[n+]);
end.

BZOJ 3437:小P的牧场(DP+斜率优化)的更多相关文章

  1. BZOJ 3437 小P的牧场(斜率优化DP)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3437 [题目大意] n个牧场排成一行,需要在某些牧场上面建立控制站, 每个牧场上只能建 ...

  2. bzoj 3437: 小P的牧场【斜率优化】

    emmm妹想到要倒着推 先假设只在n建一个控制站,这样的费用是\( \sum_{i=1}^{n} b[i]*(n-i) \)的 然后设f[i]为在i到n键控制站,并且i一定建一个,能最多节省下的费用, ...

  3. 【BZOJ-3437】小P的牧场 DP + 斜率优化

    3437: 小P的牧场 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 705  Solved: 404[Submit][Status][Discuss ...

  4. BZOJ 3437: 小P的牧场 斜率优化DP

    3437: 小P的牧场 Description 背景 小P是个特么喜欢玩MC的孩纸... 描述 小P在MC里有n个牧场,自西向东呈一字形排列(自西向东用1…n编号),于是他就烦恼了:为了控制这n个牧场 ...

  5. bzoj 3437: 小P的牧场 -- 斜率优化

    3437: 小P的牧场 Time Limit: 10 Sec  Memory Limit: 128 MB Description 小P在MC里有n个牧场,自西向东呈一字形排列(自西向东用1…n编号), ...

  6. bzoj3437 小P的牧场(斜率优化dp)

    3437: 小P的牧场 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2025  Solved: 1110[Submit][Status][Discu ...

  7. BZOJ3437 小P的牧场 【斜率优化dp】

    3437: 小P的牧场 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 1502  Solved: 836 [Submit][Status][Disc ...

  8. BZOJ 3437: 小P的牧场

    传送门 显然考虑 $dp$,设 $f[i]$ 表示前 $i$ 个牧场都被控制的最小代价 那么枚举所有 $j<i$ ,$f[i]=f[j]+val[i][j]+A[i]$ $val[i][j]$ ...

  9. BZOJ 1010: [HNOI2008]玩具装箱toy [DP 斜率优化]

    1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 9812  Solved: 3978[Submit][St ...

  10. BZOJ 1911: [Apio2010]特别行动队( dp + 斜率优化 )

    sum为战斗力的前缀和 dp(x) = max( dp(p)+A*(sumx-sump)2+B*(sumx-sump)+C )(0≤p<x) 然后斜率优化...懒得写下去了... ------- ...

随机推荐

  1. Handling Exceptions

    https://developer.apple.com/library/content/documentation/Cocoa/Conceptual/Exceptions/Tasks/Handling ...

  2. 【51nod1677】treecnt(树上数学题)

    点此看题面 大致题意: 给你一个节点从1~n编号的树,让你从中选择k个节点并通过选择的边联通,且要使选择的边数最少,让你计算对于所有选择k个节点的情况最小选择边数的总和. 题解 这道题乍一看很麻烦:最 ...

  3. python_60_装饰器3

    #嵌套函数 def foo(): print('in the foo') def bar(): print('in the bar') bar() #bar()#出错,无法在外边调用,bar函数的作用 ...

  4. 通过cmd查看环境变量名对应的环境变量值

    在VS环境中通常要添加路径,不过基本都是按照往上提供的方法添加变量名形如:$(VC_IncludePath),但是如何通过cmd命令找到真正的路径呢 未完待续……

  5. 多线程中使用HttpContext.Current为null的解决办法

    HttpContext.Current.Server.MapPath(logFile)   这个是得到具体路径的方法  正常情况下是可以的 多线程情况下就为null 下边的代码原本的作用是把网站的异常 ...

  6. 学习sqlserver的函数方法

    http://www.w3school.com.cn/sql/func_datediff.asp SQL Server DATEDIFF() 函数 SELECT DATEDIFF(day,'2008- ...

  7. Vscdoe技巧1

    vscdoe常用快捷键 主命令框 F1 或 Ctrl+Shift+P: 打开命令面板.在打开的输入框内,可以输入任何命令,例如: 按一下 Backspace 会进入到 Ctrl+P 模式 在 Ctrl ...

  8. 循环引用问题 -- dealloc方法不执行

    dealloc不执行 如果一个类在释放过后,dealloc方法没有执行,那么就代表着这个类还被其他对象所引用,引用计数不为0,这样就造成了内存泄露 昨天其他业务线开发告知他所依赖的我这边的父类VC的- ...

  9. Mysql--数据定义语言(DDL)

    DDL(Data Definition languages)数据定义语言,这些语句主要定义了不同的数据段,数据表.列.索引等操作,主要关键字有create.drop.alter. 一. 数据库的操作 ...

  10. [Wolfgang Mauerer] 深入linux 内核架构 第二章 进程管理与调度【未完】

     作为Linux开发爱好者,从事linux 开发有三年多时间.做过bsp移植,熟悉u-boot代码执行流程:看过几遍<linux 设备驱动程序开发>,分析过kernel启动流程,写过驱动, ...