https://www.lydsy.com/JudgeOnline/problem.php?id=1233

数据结构优化dp的代码总是那么抽象

题意:奶牛们讨厌黑暗。 为了调整牛棚顶的电灯的亮度,Bessie必须建一座干草堆使得她能够爬上去够到灯泡 。一共有N大包的干草(1<=N<=100000)(从1到N编号)依靠传送带连续的传输进牛棚来。第i包干草有一个 宽度W_i(1<=w_i<=10000)。所有的干草包的厚度和高度都为1. Bessie必须利用所有N包干草来建立起干草堆,并且按照他们进牛棚的顺序摆放。她可以相放多少包就放 多少包来建立起tower的地基(当然是紧紧的放在一行中)。接下来他可以放置下一个草包放在之前一级 的上方来建立新的一级。注意:每一级不能比下面的一级宽。她持续的这么放置,直到所有的草包都被安 置完成。她必须按顺序堆放,按照草包进入牛棚的顺序。说得更清楚一些:一旦她将一个草包放在第二级 ,她不能将接下来的草包放在地基上。 Bessie的目标是建立起最高的草包堆。

Input

第1行:一个单一的整数N。 第2~N+1行:一个单一的整数:W_i。

奥妙重重dp题。

第一眼觉得是个反向贪心,从上往下取,每次取最小的可以形成下一层的干草集合,后来发现是个假算法。

反例


6 2 3 7 10 11

既然贪心不行,就自然而然想到dp,考虑dp[i]表示到这一层的最大层数,发现不可行,因为这样的dp方法不能固定一个状态。

但是这一题有一个很强的证明:对于一座干草堆而言,基底最小的状态一定为最优状态。

也就是说,每一个状态中,地基最小的状态就是高度最高的状态,我们可以简单的推断一下,每一个高度最高的状态必然可以从下往上贪心的堆干草最终达到基底最小的状态。

这样我们原来的dp就变成了dp[i]表示i -- n种所有干草组成的集合形成的最小的干草堆的基底。

我们反向从N - > 1遍历,对于i,需要找到一个最小的大于i的j满足 sum[j - 1] - sum[i - 1] >= dp[j], 这样就变成了一个O(N * N)的dp,但是这还不够,到这一步需要考虑优化。

我们将转移方程变形 变成sum[i - 1] <= sum[j - 1] - dp[j],注意到我们每次需要的考虑的集合在每一次询问i的时候都多加入1个,以及当存在k < j && sum[k - 1] - dp[k]  >= sum[j - 1] - dp[j]的时候,j就是一个多余的状态,sum[i - 1]由于其前缀和的性质,事实上是递减的,对于每一次集合内的查找,我们在查找到最小的j的时候,可以将其余的比sum[i - 1]大的状态全部舍弃。

以上三个性质中,显然可以发现用来维护这个集合的数据结构就是单调队列,最终时间复杂度O(n)

#include <map>
#include <set>
#include <ctime>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
#define For(i, x, y) for(int i=x;i<=y;i++)
#define _For(i, x, y) for(int i=x;i>=y;i--)
#define Mem(f, x) memset(f,x,sizeof(f))
#define Sca(x) scanf("%d", &x)
#define Sca2(x,y) scanf("%d%d",&x,&y)
#define Scl(x) scanf("%lld",&x);
#define Pri(x) printf("%d\n", x)
#define Prl(x) printf("%lld\n",x);
#define CLR(u) for(int i=0;i<=N;i++)u[i].clear();
#define LL long long
#define ULL unsigned long long
#define mp make_pair
#define PII pair<int,int>
#define PIL pair<int,long long>
#define PLL pair<long long,long long>
#define pb push_back
#define fi first
#define se second
typedef vector<int> VI;
const double eps = 1e-9;
const int maxn = 1e5 + 10;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
int N,M,tmp,K;
int a[maxn];
int sum[maxn];
int dp[maxn];
int f[maxn];
int Queue[maxn];
int main()
{int N; Sca(N);  //sum[i - 1] <= sum[j - 1] - dp[j]For(i,1,N){Sca(a[i]);sum[i] = sum[i - 1] + a[i];} int head = 1,tail = 1; Queue[1] = N + 1;_For(i,N,1){while(head < tail && sum[i - 1] <= sum[Queue[head + 1] - 1] - dp[Queue[head + 1]]) head++;dp[i] = sum[Queue[head] - 1] - sum[i - 1];f[i] = f[Queue[head]] + 1;while(head <= tail && sum[i - 1] - dp[i] >= sum[Queue[tail] - 1] - dp[Queue[tail]]) tail--;Queue[++tail] = i;}printf("%d",f[1]);#ifdef VSCodesystem("pause");#endifreturn 0;
}

转载于:https://www.cnblogs.com/Hugh-Locke/p/9671948.html

bzoj1233 单调队列优化dp相关推荐

  1. 算法笔记--单调队列优化dp

    单调队列:队列中元素单调递增或递减,可以用双端队列实现(deque),队列的前面和后面都可以入队出队. 单调队列优化dp: 问题引入: dp[i] = min( a[j] ) ,i-m < j ...

  2. poj 2373(单调队列优化dp)

    在长为L(<=1000000)的草地(可看成线段)上装喷水头,喷射是以这个喷水头为中心,喷水头的喷洒半径是可调节的调节范围为[a,b].要求草地的每个点被且只被一个喷水头覆盖,并且有些连续区间必 ...

  3. poj 1821(单调队列优化dp)

    题意:有一道线性篱笆由N个连续的木板组成.有K个工人,你要叫他们给木板涂色.每个工人有3个参数:L 表示 这个工人可以涂的最大木板数目,S表示这个工人站在哪一块木板,P表示这个工人每涂一个木板可以得到 ...

  4. 洛谷P3195 [HNOI2008]玩具装箱TOY(单调队列优化DP)

    题目描述 P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维容器中.P教授有编号为1...N的N件玩具, ...

  5. 【计蒜客 - 蓝桥训练】蒜厂年会(单调队列优化dp,循环数列的最大子段和)

    题干: 在蒜厂年会上有一个抽奖,在一个环形的桌子上,有 nn 个纸团,每个纸团上写一个数字,表示你可以获得多少蒜币.但是这个游戏比较坑,里面竟然有负数,表示你要支付多少蒜币.因为这些数字都是可见的,所 ...

  6. 单调队列以及单调队列优化DP

    单调队列定义: 其实单调队列就是一种队列内的元素有单调性的队列,因为其单调性所以经常会被用来维护区间最值或者降低DP的维数已达到降维来减少空间及时间的目的. 单调队列的一般应用: 1.维护区间最值 2 ...

  7. POJ 1821 Fence(单调队列优化DP)

    题解 以前做过很多单调队列优化DP的题. 这个题有一点不同是对于有的状态可以转移,有的状态不能转移. 然后一堆边界和注意点.导致写起来就很难受. 然后状态也比较难定义. dp[i][j]代表前i个人涂 ...

  8. 【单调队列优化DP】烽火传递 LibreOJ - 10180

    题目来源 点我进入提交题目 反思 因为目前在学习单调队列优化DP,所以会往单调队列上面想.然后犯了一个错误就是,认为这个题目只要用单调队列就可以完成,单调队列只是用来减少时间复杂度的,遇到了求最优解的 ...

  9. AcWing 1089 烽火传递 题解(动态规划—DP—单调队列优化DP)

    AcWing 1089 烽火传递 单调队列优化DP,思路比较简单,维护一个保持元素单调递增的单调队列,队首就是第i座烽火台能接收到的,代价最小的方案,加上第i座烽火台的代价就是这座烽火台的最小值 #i ...

最新文章

  1. 大学计算机专业副修课,计算机学院举行本科课程教学大纲修订工作研讨会
  2. python源文件保存在哪里_python数据爬下来保存在哪里
  3. [Beta]第五次 Scrum Meeting
  4. 光刷题不参加这些算法竞赛?太亏了!
  5. 'webpack-dev-server' 不是内部或外部命令,也不是可运行的程序 或批处理文件。
  6. C# 8.0 的默认接口方法
  7. linux 文件inode,linux文件系统-inode学习整理
  8. 【华为云技术分享】数据管理服务DAS 之 数据库自动化运维功能展播4:慢SQL
  9. Selenium UI自动化测试(六)WebDriver常用操作方法
  10. oracle帮助文档_Spring Boot Config文档,使用IntelliJ IDEA的两种方法
  11. JMeter的编码与HTTP请求
  12. Linux-install-mysql5.6
  13. 1007 计算(a+b)*c的值
  14. 2分钟部署人生模拟器,解锁人生新剧情
  15. linux中常用打开pdf文档指令软件
  16. 如何在指定网站搜索内容
  17. pk 与fk mysql_什么是MySQL FK的正确命名约定?
  18. 365赚钱宝养猫小程序程序源码下载
  19. 逆向世界杯直播App 央视影音-iOS客户端
  20. java实现2进制转16进制,16进制转2进制(通过查表法和强转实现)

热门文章

  1. 计算机二级考试Access教程
  2. 为何不精通C? 03 深入剖析声明
  3. 各种分页存储过程 (转)
  4. ASIHTTPRequest 常用方法讲解
  5. 老陈学 C++ 序列之二: 友元函数
  6. C++之STL理论基础
  7. 用友云开发者中心助你上云系列之在线调试
  8. Win32-Application的窗口和对话框
  9. 3D打印技术在医疗领域能做些什么?帮助精确完成手术
  10. iOS逆向工程(简单利用dumpdecrypted给ipa砸壳)