题目链接

http://www.lydsy.com/JudgeOnline/problem.php?id=1010

思路


斜率优化DP

我们知道,有些DP方程可以转化成DP[i]=f[j]+x[i]的形式,其中f[j]中保存了只与j相关的量。这样的DP方程我们可以用单调队列进行优化,从而使得O(n^2)的复杂度降到O(n)。

可是并不是所有的方程都可以转化成上面的形式,比如dp[i]=dp[j]+(x[i]-x[j])*(x[i]-x[j])。如果把右边的乘法化开的话,会得到x[i]*x[j]的项。这就没办法使得f[j]里只存在于j相关的量了。于是上面的单调队列优化方法就不好使了,所以我们需要一种新的优化方法,叫做斜率优化。

拿本题来说,我们设dp[i]表示装到第i个玩具的时候最少的花费,sum[i]表示前i个玩具的长度和。于是方程就是:

dp[i]  = min{ dp[j] + [i+sum[i]-(j+1+sum[j])-L]2 };

题目的N<=50000,二维铁定超时了。我们就来看看斜率优化如何做到从O(n^2)复杂度降到O(n)。

分析:

我们假设k<j<i。如果在j的时候决策要比在k的时候决策好,那么也是就是dp[j] + [i+sum[i]-(j+1+sum[j])-L]< dp[k] + [i+sum[i]-(j+1+sum[k])-L]2。(因为是最小花费,所以优就是小于)

两边移项一下得到:[dp[j]+(j+1+sum[j])^2-(dp[k]+(k+1+sum[k])^2)]/(2*( (j+1+sum[j])-(k+1+sum[k]) )) < (i + sum[i]-L)。我们把dp[j]-num[j]^2看做是yj,把2*num[j]看成是xj。

那么不就是(yj-yk)/(xj-xk) < A[i]么?   左边是不是斜率的表示?

那么(yj-yk)/(xj-xk) < A[i]说明了什么呢?  我们前面是假设j的决策比k的决策要好才得到这个表示的,那么g[j,k]=(yj-yk)/(xj-xk) < A[i]代表这j的决策比k的决策要更优

斜率优化在于:①设k<j<i,如果g[i,j] < g[j,k],那么j点便永远不可能成为最优解,可以直接将它踢出我们的最优解集。‘

为什么呢?我们假设g[i,j] < A[i],那么就是说i点要比j点优,排除j点。如果g[i,j] >= A[i],那么j点此时是比i点要更优,但是同时g[j,k]>g[i,j]>A[i]。这说明还有k点会比j点更优,同样排除j点。排除多余的点,这便是一种优化!

由于我们排除了g[i,j] < g[j,k]的情况,所以整个有效点集呈现一种下凸性质,即k j的斜率要小于j i的斜率。这样,从左到右的斜率之间就是单调递增的了,所以我们就可以对g维护一个单调队列。

同时,函数A[i]也要具有单调递增的特性,所以第二个优化就在于:②如果单调队列的头两个点为i, j,斜率g[j, i] < A[i],则说明j优于i,并且由于A单调递增,所以g[j,i]恒小于A[],所以可以直接把i排除掉。

于是对于这类题斜率优化做法可以总结如下:

1,用一个单调队列来维护解集。

2,假设队列中从头到尾已经有元素a b c。那么当d要入队的时候,我们维护队列的下凸性质,即如果g[d,c]<g[c language=",b"][/c],那么就将c点删除。直到找到g[d,x]>=g[x,y]为止,并将d点加入在该位置中。

3,求解时候,从队头开始,如果已有元素a b c,当i点要求解时,如果g[b,a]<A[i],那么说明b点比a点更优,a点可以排除,于是a出队。最后dp[i]=getDp(q[head])。


代码

[cpp]
#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <string>
#include <cstring>
#include <vector>
#include <queue>
#define MID(x,y) ((x+y)/2)
#define MEM(a,b) memset(a,b,sizeof(a))
#define REP(i, begin, end) for (int i = begin; i <= end; i ++)
using namespace std;

long long dp[50005], sum[50005];
vector <int> Q;
int n, L;
inline long long getDP(int i, int j){
return dp[j] + (i+sum[i]-(j+1+sum[j])-L) * (i+sum[i]-(j+1+sum[j])-L);
}
inline double getUP(int j, int k){
return (double)dp[j] + (double)(j+1+sum[j]) * (double)(j+1+sum[j]) - ((double)dp[k] + (double)(k+1+sum[k]) * (double)(k+1+sum[k]));
}
inline double getDOWN(int j, int k){
return 2.0 * ((j+1+sum[j]) - (k+1+sum[k]));
}
inline double getRIGHT(int i){
return (double)(sum[i] + i - L);
}
inline double slope(int j, int k){
return getUP(j, k) / getDOWN(j, k);
}
int main(){
scanf("%d %d", &n, &L);
sum[0] = 0, dp[0] = 0;
for (int i = 1; i <= n; i ++){
scanf("%lld", &sum[i]);
sum[i] += sum[i-1];
}
Q.clear();
Q.push_back(0);
for (int i = 1; i <= n; i ++){
while((int)Q.size() > 1 && slope(Q[1], Q[0]) <= getRIGHT(i))
Q.erase(Q.begin());
int tmp = Q[0];
dp[i] = getDP(i, Q[0]);
while( (int)Q.size() > 1 && slope(Q[(int)Q.size()-1], Q[(int)Q.size()-2]) >= slope(i, Q[(int)Q.size()-1]) )
Q.pop_back();
Q.push_back(i);
}
printf("%lld\n", dp[n]);
return 0;
}
[/cpp]

BZOJ 1010 [HNOI2008]玩具装箱 (斜率优化DP)的更多相关文章

  1. BZOJ 1010 HNOI2008 玩具装箱 斜率优化

    题目链接: http://www.lydsy.com/JudgeOnline/problem.php?id=1010 Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的 ...

  2. bzoj 1010 [HNOI2008]玩具装箱toy(DP的斜率优化)

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

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

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

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

    1010: [HNOI2008]玩具装箱toy Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再 ...

  5. Bzoj 1010: [HNOI2008]玩具装箱toy(斜率优化)

    1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MB Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定 ...

  6. BZOJ 1010 [HNOI2008]玩具装箱toy

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

  7. BZOJ 1010: [HNOI2008]玩具装箱toy(斜率优化dp)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1010 题意: 思路: 容易得到朴素的递归方程:$dp(i)=min(dp(i),dp(k)+(i-k ...

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

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1010 题意: 有n条线段,长度分别为C[i]. 你需要将所有的线段分成若干组,每组中线段的 ...

  9. BZOJ.1010.[HNOI2008]玩具装箱toy(DP 斜率优化/单调队列 决策单调性)

    题目链接 斜率优化 不说了 网上很多 这的比较详细->Click Here or Here //1700kb 60ms #include<cstdio> #include<cc ...

随机推荐

  1. cpu与寄存器,内核态与用户态及如何切换

    cpu:相当于计算机的大脑负责运算和发送命令: 寄存器:寄存器是cpu当中的一个有限存储部件,cpu从内存调用数据时,寄存器会将从内存调用的数据进行更新在寄存器中以一个字或变量进行存储. 寄存器总共分 ...

  2. delphi webbrowser 跨域访问

    procedure IterateFrames(const AWB: IWebBrowser2);var Doc: IHTMLDocument2; Container: IOleContainer; ...

  3. 实参相依查找[条款25]----《C++必知必会》

    名字空间对现代C++编程和设计有着深远的影响.其中有些影响直接而明显,例如using声明和using 指令以及采用名字空间作用域加以限定的名字.然而,名字空间还有一些在语法上不那么明显但任然很有基础. ...

  4. SQL Server DDL触发器

    DDL 触发器作用: DDL 触发器主要用于防止对数据库架构.视图.表.存储过程等进行的某些修改. DDL 触发器事件: DDL 触发器在创建用来监视并响应该数据库或服务器实例中的活动的事件通知时,可 ...

  5. ECMAScript与JavaScript

    ECMAScript发展史: 1997 1998.6 1999.12 2008.7 2009.12 2015.61996年11月 javaScript创造者Netscape公司将javaScript提 ...

  6. hdu6215 Brute Force Sorting

    地址:http://acm.split.hdu.edu.cn/showproblem.php?pid=6215 题目: Brute Force Sorting Time Limit: 1000/100 ...

  7. $思维导图——numpy基本知识

  8. [Linux 001]——计算机和操作系统的基础知识

    在正式开始学习 Linux 操作系统之前,有必要先回顾/学习一下计算机和操作系统的基本知识,为我们在后续的学习中铺路搭桥,在了解计算机一些基础原理的条件下再去进行学习,理解应该会更透彻一些.我会从一个 ...

  9. 20145312 《Java程序设计》第四周学习总结

    20145312 <Java程序设计>第四周学习总结 学习笔记 Chapter 6 6.1何为继承 1.定义:面向对象中子类继承父类,避免重复的行为定义. 6.1.1 继承共同行为 1.以 ...

  10. nmap与Nessus扫描特定靶机分析

    打开装载Metasploitable2虚拟机的靶机,并获取靶机ip: 使用nmap+ip初步扫描靶机 PORT为端口,STATE为端口开放状态,SERVICE为端口的提供的服务.靶机的MAC地址为: ...