BZOJ 1010: [HNOI2008]玩具包装toy
职务地址:http :// www . lydsy . com / JudgeOnline / problem . php ? id = 1010
题目大意:见原题。
算法分析:
设s[i]为c[i]的前缀和,f[i]表示第1个物品到第i个物品的最小代价。
易得DP方程为f[i]=min(f[j]+(s[i]-s[j]+i-j+1-l)^2)。
设t[i]=s[i]+i。则f[i]=min(f[j]+(t[i]-t[j]+1-l)^2。
设m=t[i]-l-1。则f[i]=min(f[j]+(m-t[j])^2)。
对i做决策时,设j,k为2个一般的决策点。设j<k。
若k比j优,则有
f[k]+(m-t[k])^2<f[j]+(m-t[j])^2
展开,得
f[k]+m^2-2*m*t[k]+t[k]^2<f[j]+m^2-2*m*t[j]+t[j]^2
化简,得
((f[k]+t[k]^2)-(f[j]+t[j]^2))/(t[k]-t[j])<2*m=2*(s[i]+i-l-1)
于是能够用单调队列来优化DP。
单调队列优化DP的伪代码:
for (int i=1;i<=n;++i){
while (h+1<t && slope(q[h],q[h+1])<...) h++;
f[i]=...;
while (h+1<t && slope(q[t-2],q[t-1])>slope(q[t-1],i)) t--;
q[t++]=i;
}
斜率优化的2个重要结论(对于此题):
1)j<k,若slope(j,k)<2*(s[i]+i-l-1),则k比j优。
(PS:不等式左边为仅仅与j,k有关的式子。不等式右边为仅仅与i有关的式子)
2)j<k<l,若slope(j,k)>slope(k,l),则k能够被舍去。
Code:
#include <cstdio> #define N 50000 using namespace std; int n,l,h,t=1,c[N+10],q[N+10];
long long s[N+10],f[N+10]; template <typename T>
inline T sqr(T x){return x*x;}
inline double slope(int i,int j){return (f[i]+sqr((double)s[i]+i)-f[j]-sqr((double)s[j]+j))/(s[i]+i-s[j]-j);} int main(){
#ifndef ONLINE_JUDGE
freopen("1010.in","r",stdin);
freopen("1010.out","w",stdout);
#endif
scanf("%d%d",&n,&l);
for (int i=1;i<=n;++i) scanf("%d",&c[i]);
for (int i=1;i<=n;++i) s[i]=s[i-1]+c[i];
for (int i=1;i<=n;++i){
while (h+1<t && slope(q[h],q[h+1])<2*(s[i]+i-l-1)) h++;
f[i]=f[q[h]]+sqr(s[i]-s[q[h]]+i-q[h]-1-l);
while (h+1<t && slope(q[t-2],q[t-1])>slope(q[t-1],i)) t--;
q[t++]=i;
}
printf("%lld\n",f[n]);
return 0;
}
By Charlie Pan
Jul 18,2014
版权声明:本文博客原创文章。博客,未经同意,不得转载。
BZOJ 1010: [HNOI2008]玩具包装toy的更多相关文章
- BZOJ 1010: [HNOI2008]玩具装箱toy [DP 斜率优化]
1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 9812 Solved: 3978[Submit][St ...
- BZOJ 1010: [HNOI2008]玩具装箱toy 斜率优化DP
1010: [HNOI2008]玩具装箱toy Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再 ...
- bzoj 1010 [HNOI2008]玩具装箱toy(DP的斜率优化)
1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 7874 Solved: 3047[Submit][St ...
- BZOJ 1010 [HNOI2008]玩具装箱toy
1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 7184 Solved: 2724[Submit][St ...
- Bzoj 1010: [HNOI2008]玩具装箱toy(斜率优化)
1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MB Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定 ...
- BZOJ 1010: [HNOI2008]玩具装箱toy(斜率优化dp)
http://www.lydsy.com/JudgeOnline/problem.php?id=1010 题意: 思路: 容易得到朴素的递归方程:$dp(i)=min(dp(i),dp(k)+(i-k ...
- BZOJ 1010: [HNOI2008]玩具装箱toy | 单调队列优化DP
原题: http://www.lydsy.com/JudgeOnline/problem.php?id=1010 题解: #include<cstdio> #include<algo ...
- BZOJ 1010 [HNOI2008]玩具装箱toy:斜率优化dp
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1010 题意: 有n条线段,长度分别为C[i]. 你需要将所有的线段分成若干组,每组中线段的 ...
- BZOJ.1010.[HNOI2008]玩具装箱toy(DP 斜率优化/单调队列 决策单调性)
题目链接 斜率优化 不说了 网上很多 这的比较详细->Click Here or Here //1700kb 60ms #include<cstdio> #include<cc ...
随机推荐
- Linux 下开启ssh服务(转)
二.SSH SSH 为 Secure Shell 的缩写,由 IETF 的网络工作小组(Network Working Group)所制定:SSH 为建立在应用层和传输层基础上的安全协议.SSH 是目 ...
- poj 2478 Farey Sequence(欧拉函数是基于寻求筛法素数)
http://poj.org/problem?id=2478 求欧拉函数的模板. 初涉欧拉函数,先学一学它主要的性质. 1.欧拉函数是求小于n且和n互质(包含1)的正整数的个数. 记为φ(n). 2. ...
- A WPF File ListView and ComboBox
源码下载: Download FileListView_Version_2.zip Download FileListView_Version_2_Binaries.zip Download File ...
- Gradle cookbook(转)
build.gradle apply plugin:"java" [compileJava,compileTestJava,javadoc]*.options*.encoding ...
- 深入理解spring中的各种注解(转)
Spring中的注解大概可以分为两大类: 1)spring的bean容器相关的注解,或者说bean工厂相关的注解: 2)springmvc相关的注解. spring的bean容器相关的注解,先后有:@ ...
- linux文件打开模式
文件打开 int open(const char *pathname, int flags, mode_t mode); 普通方式(Canonical mode) flags中没有设置O_SYN ...
- 最简单的视频编码器:基于libx265(编码YUV为H.265)
===================================================== 最简单的视频编码器系列文章列表: 最简单的视频编码器:编译 最简单的视频编码器:基于libx ...
- UVA434 - Matty's Blocks
option=com_onlinejudge&Itemid=8&page=show_problem&category=457&problem=375&mosms ...
- 返璞归真 asp.net mvc (4) - View/ViewEngine
原文:返璞归真 asp.net mvc (4) - View/ViewEngine [索引页] [源码下载] 返璞归真 asp.net mvc (4) - View/ViewEngine 作者:web ...
- java ClassLoader static
package init; class Person { private static Person person = new Person(); public static int count2 = ...