【HNOI2008】玩具装箱
P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京。他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维容器中。P教授有编号为1...N的N件玩具,第i件玩具经过压缩后变成一维长度为Ci.为了方便整理,P教授要求在一个一维容器中的玩具编号是连续的。同时如果一个一维容器中有多个玩具,那么两件玩具之间要加入一个单位长度的填充物,形式地说如果将第i件玩具到第j个玩具放到一个容器中,那么容器的长度将为 x=j-i+Σ(Ck) i<=K<=j 制作容器的费用与容器的长度有关,根据教授研究,如果容器长度为x,其制作费用为(X-L)^2.其中L是一个常量。P教授不关心容器的数目,他可以制作出任意长度的容器,甚至超过L。但他希望费用最小
第一行输入两个整数N,L.接下来N行输入Ci.1<=N<=50000,1<=L,Ci<=10^7
输出最小费用
提取信息:
1.n件玩具
2.可以将某段玩具压缩,需要付出代价
3.压缩次数不限
4.要求计算最小费用
我们首先设计状态:
直观的:f[i]表示到了第i件玩具
这时我们发现剩余的2个信息无需也无法在向状态空间中添加,因为压缩次数是不限的,或者说,这样就足够描述状态空间了。
所以我们就用f[i]表示处理到第i件玩具时的最小费用
接着思考转移:
我们能够用来连接转移的信息只有将某段玩具压缩了,所以我们很容易地想到:
设j表示i之前的某个点,我们要压缩j+1到i之间的玩具。
我们再记一个前缀和数组S表示数组C的前缀和
这样我们就能写出状态转移方程
f[i] = min{f[j] + (i - j + S[i] - S[j] - L - 1) ^ 2}
复杂度O(n^2)
肯定无法通过这道题,思考如何优化,首先我们将循环变量i看做一个常量
将min函数去掉,得到
f[i] = f[j] + (i - j + S[i] - S[j] - L - 1)^2
为了方便,我们将多项式分离成: i + S[i]和 j + S[j] + L + 1
用a[i]存储i+S[i],即a[i] = i + S[i],同理:b[i] = i + S[i] + L + 1
这样我们将方程变为
f[i] = f[j] + (a[i] - b[j])^2
==> f[i] = f[j] + a[i]^2 - 2 * a[i] * b[j] + b[j]^2
==> f[j] + b[j]^2 = 2 * a[i] * b[j] + f[i] - a[i]^2
我们将b[j]看做横坐标,2*a[i]看做斜率,f[j] + b[j]^2看做y
这个式子就可以看做以2*a[i]为斜率的直线,纵截距为f[i] - a[i]^2,且截距未知
于是我们很容易能想到,我们要在某个位置取到最优决策,就要求截距最小化
我们很容易看出2*a[i]和f[i] - a[i]^2是常量,而对于自变量b[j]本身具有单调性,所以显然该方程具有决策单调性,所以我们使用斜率优化
对于任意三个决策点(f[j1] + b[j1]^2, b[j]), (f[j2] + b[j2]^2, b[j2]),(f[j3] + b[j3]^2, b[j3]), 设b[j1] < b[j2] < b[j3],若j2为最优决策,则有
f[j1] + b[j1]^2 - b[j1] > f[j2] + b[j2]^2 - b[j2]
同理:f[j2] + b[j2]^2 - b[j2] < f[j3] + b[j3]^2 - b[j3]
可以变形为:f[1] + b[j1]^2 - f[j2] - b[j2]^2 / (b[j1] - b[j2]) > 1
f[j2] + b[j2] ^ 2 - f[j3] - b[j3]^2 / (b[j2] - b[j3]) < 1
也就是说,加入某个中间点j2是最优决策,则必须满足:
f[j1] + b[j1]^2 - f[j2] - b[j2]^2 / (b[j1] - b[j2]) > f[j2] + b[j2] ^ 2 - f[j3] - b[j3]^2 / (b[j2] - b[j3])
//2*a[i]被约分掉了
然后我们可以发现等式两边实际分别表示的是直线j1j2和直线j2j3的斜率
也就是说我们可以通过维护斜率的单调递增来保证每次我们取到的是最优解,也就是维护一个下凸壳
对于下凸壳的维护,我们需要建立一个支持如下操作的单调队列:
1.每次取出队头q[l], q[l + 1],若斜率 (f[q[l]] + b[q[l]]^2 - f[q[l+1]] - b[q[l+1]]^2) / (b[q[l]] - b[q[l+1]]) <= 2 * a[i],则将队头出列,直至不满足
2.此时对于队头,就是当前状态的最优决策,取出队头进行转移
3.检查队尾,对于队尾q[r - 1], q[r],以及新决策i,若不满足斜率单调递增,则出队,直到满足
4.将决策i插入队尾
//可能有部分出入和错误,因为题解太长而我又是记事本写的....
//建议读代码
- #include<bits/stdc++.h>
- #define ll long long
- #define ld long double
- #define db double
- using namespace std;
- const int maxn = ;
- ld f[maxn], s[maxn];//, a[maxn], b[maxn];
- int n, L;
- int l, r, q[maxn];
- inline ld a(int i) {
- return s[i] + i;}
- inline ld b(int i) {
- return a(i) + L + ;}
- inline ld X(int i) {
- return b(i);}
- inline ld Y(int i) {
- return f[i] + b(i) * b(i);}
- inline ld ask(int i,int j) {
- return (Y(i) - Y(j)) / (X(i) - X(j));}
- int main() {
- cin >> n >> L;
- for(int i = ; i <= n; ++i) {
- cin >> s[i];
- s[i] += s[i - ];
- }
- l = r = ;
- for(int i = ; i <= n; ++i) {
- while(l < r && ask(q[l], q[l + ]) < * a(i)) l++;
- f[i] = f[q[l]] + (a(i) - b(q[l]))* (a(i) - b(q[l]));
- while(l < r && ask(i, q[r]) < ask(q[r], q[r - ])) r--;
- q[++r] = i;
- }
- cout << (ll)f[n] << '\n';
- return ;
- }
【HNOI2008】玩具装箱的更多相关文章
- 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教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再 ...
- 【bzoj1010】[HNOI2008]玩具装箱toy
1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 9281 Solved: 3719[Submit][St ...
- 【BZOJ 1010】 [HNOI2008]玩具装箱toy (斜率优化)
1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 9330 Solved: 3739 Descriptio ...
- 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 ...
- bzoj1010[HNOI2008]玩具装箱toy 斜率优化dp
1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 11893 Solved: 5061[Submit][S ...
- P3195 [HNOI2008]玩具装箱TOY(斜率优化dp)
P3195 [HNOI2008]玩具装箱TOY 设前缀和为$s[i]$ 那么显然可以得出方程 $f[i]=f[j]+(s[i]-s[j]+i-j-L-1)^{2}$ 换下顺序 $f[i]=f[j]+( ...
- [luogu P3195] [HNOI2008]玩具装箱TOY
[luogu P3195] [HNOI2008]玩具装箱TOY 题目描述 P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆, ...
- cogs 1330 [HNOI2008]玩具装箱toy
cogs 1330 [HNOI2008]玩具装箱toy 瞎扯,急忙AC的请跳过 感觉数据结构写的太多了有点晕=+ 发现还没学斜率优化+- 于是来学一学QwQ 上次这题打了个决策优化直接水过了..理论O ...
随机推荐
- [洛谷P2482][SDOI2010]猪国杀
题目大意:猪国杀,又一道大模拟题 题解:模拟,对于一个没有玩过三国杀的人来说,一堆细节不知道,写的十分吃力 卡点:无数,不想说什么了,这告诉我要多玩游戏 C++ Code: #include < ...
- [TJOI2018]数学计算 线段树
---题面--- 题解: ,,,考场上看到这题,没想到竟然是省选原题QAQ,考场上把它当数学题想了好久,因为不知道怎么处理有些数没有逆元的问题....知道这是线段树后恍然大悟. 首先可以一开始就建出一 ...
- 遇到问题---java---git下载的maven项目web用tomcat发布时不带子项目
遇到的情况是用git下载maven项目,然后用mvn eclipse:eclipse命令标记为eclipse项目之后,使用maven插件导入之后用tomcat发布运行,发现maven关联的几个子项目没 ...
- MySQL里执行SHOW INDEX结果中Cardinality的含义
今天在写一个Perl脚本,想自动化查找出MySQL数据库中可能无效的索引,于是根据朝阳的书上提到的一些规则,我来设计了一些判断方法,其中发现某个我想要的值就是SHOW INDEX FROM table ...
- mysql删除id最小的条目
DELETE FROM 表1 WHERE Mid in (select Mid from (SELECT Min(Mid) Mid FROM 表1 c1) t1);
- APP兼容性测试
一.APP兼容性范围以及问题 1.硬件 各个硬件结构 2.软硬件之间 硬件dll库(C++) 软硬件之间的通信,各个厂商提供的ROM 3.软件 浏览器.操作系统.数据库.手机.功能兼容性(功能修改,二 ...
- Spring + Mybatis - 原始dao开发整合 与 Mapper代理整合
Spring + Mybatis - 原始dao开发整合 与 Mapper代理整合 标签: mybatisSpringbeanApplicationContextMapper 2015-12-31 1 ...
- MDK stm32 仿真
直接选择simulator,仿真时报错 *** error 65: access violation at 0x40021000 : no 'read' permission 修改系统配置,原配置如下 ...
- canvas动态绘制饼状图,
当我们使用Echrts很Highcharts的时候,总是觉得各种统计图表是多么神奇,今天我就用现代浏览器支持的canvas来绘制饼状统计图,当然仅仅是画出图并没什么难度,但是统计图一般都有输入,根据不 ...
- 网络基础(osi、协议)
*互联网协议 人和人沟通需要一套共同的标准,英语就是普遍的一种,计算机如果需要进行联网互通,也需要一种统一的标准,如果所有的计算机都遵守这种标准,就会实现网络的互联. 1.一系列统一的标准,这些标准称 ...