首页 > 其他分享 >2023 7月 dp做题记录

2023 7月 dp做题记录

时间:2024-04-20 11:44:21浏览次数:18  
标签:分割 sum 斜率 做题 复杂度 2023 我们 dp

目录

7月 dp做题记录

The Bakery

这道题的状态转移并不难列,经典的分段问题,设状态 \(dp_{i,j}\) 表示前 \(i\) 个数字分了 \(j\) 段的最大价值,转移可以写成

\(dp_{i,j}=\max(dp_{k,j-1}+val_{k+1,i})\)

这里的 \(val_{k+1,i}\) 指的是区间 \([k+1,i]\) 的价值,即不同数字的个数。如果不加优化,这个方程的转移复杂度是 \(O(n^2k)\) 的。

现在的瓶颈是,如何快速的求出先前 \(dp_{k,j-1}+val_{k+1,i}\) 的最大值。我们先分析 \(val\) 贡献的计算,对于相同数字,出现两次是不会增加贡献的,所以我们可以把相同数字的贡献分段,对于数字 \(a_i\),它在区间 \([pre_i,i]\) 中贡献,保证每个区间只存在一个 \(a_i\),这样贡献就不会重复计算。

求区间最大值,我们很容易想到用数据结构维护。需要区间修改和区间查询,求区间最大值,所以我们使用线段树。每枚举一个段数 \(j\),我们就给线段树赋上初值 \(dp_{k,j-1}\),再枚举 \(i\) 的同时更新贡献,并记录 \(dp_{i,j}\) 的答案,即区间最大值。

复杂度 \(O(kn\log n)\)

P5785 [SDOI2012] 任务安排

先考虑朴素的做法,当前费用与启动次数有关,并且与区间内的和有关,所以处理前缀和,设状态 \(dp_{i,j}\) 为到第 \(i\) 个任务,分了 \(j\) 批的最小费用,转移可以列成

\(dp_{i,j}=\min(dp_{k,j-1}+(s+sumt_i)\times (sumc_i-sumc_j))\)

这样的状态复杂度是 \(O(n^3)\) 的。

考虑转换状态表示,我们发现当前启动仅对之后的任务费用有影响,所以我们把费用提前,把启动多出的费用提前计算,这样就不用考虑先前的启动次数了,状态 \(dp_i\) 的意义也变成到第 \(i\) 个任务,分了若干批的费用再加上由于启动多出的费用的最小费用。

\(dp_i=\min(dp_j+s\times(sumc_n-sumc_j)+sumt_i\times(sumc_i-sumc_j))\)

复杂度变为 \(O(n^2)\),可以通过弱化版的 P2365 任务安排

这时候我们希望能更快找到先前的最小值。化简时发现当前枚举状态 \(i\) 会对先前的 \(j\) 贡献造成影响,所以不能用普通单调队列维护,同时也无法用数据结构维护。

化简式子,这里将 \(sumt\) 和 \(sumc\) 简写成 \(t\) 和 \(c\),

\(dp_{j}=(s+t_i)c_j+dp_i-sc_n-t_ic_i\)

这里为什么可以把 \(\min\) 去掉呢,是因为我们把 \(dp_j\) 和 \(c_j\) 看作一次函数中的 \(f(x)\) 和 \(x\),在平面直角坐标系中看成一个个点 \((c_j,dp_j)\) 组成的点集。

这样子,我们可以发现这是一次函数的形式,并且要使 \(dp_i\) 最小,即让这个一次函数的截距 \(dp_i-sc_n-t_ic_i\) 最小。所以我们把问题抽象化,取 \(j\) 的过程可以看作这条已知斜率的一次函数的直线不断向上移动的过程。可以发现,第一次碰到点集中的 \((c_j,dp_j)\) 时,截距最小。

所以问题就变成,维护点集中有价值的点,并快速求出第一次碰到的点。由于点是无序的,我们可以两两连线,发现对答案有贡献的点即为这些直线组成的下凸包。根据下凸包的性质,组成下凸包的直线的斜率不断增大,我们可以用单调队列维护不断增大的斜率,若出现一个点,它与队尾连线的斜率小于队尾与它的前一位连线的斜率时,我们弹出队尾(因为这时队尾会成为上凸包的顶点)。

在此题中,\(s+sumt_i\) 并不是单调递增,所以队头不一定是最优的,我们需要用二分查找找到第一个碰到的点。具体的,若队列中点 \(q_{mid}\) 与 \(q_{mid+1}\) 的连线大于斜率 \(s+sumt_i\) 时,右端点左移。

这样复杂度降到 \(O(n\log n)\),这个技巧叫做斜率优化。

P3195 [HNOI2008] 玩具装箱

经典的斜率优化动态规划问题。

首先考虑朴素做法,此题为典型的分段问题,并且不关心分了几段,可设状态 \(dp_{i}\) 表示到第 \(i\) 个玩具,分了若干段的最小费用。转移为

\(dp_i=\min(dp_j+(i-j-1-\sum\limits_{k=j+1}^jC_k-L)^2)\)

式子的变量太多,化简式子。可以发现求和部分可以用前缀和维护。将 \(sum_i=\sum\limits_{k=1}^iC_k\),有 \(\sum\limits_{k=j+1}^jC_k=sum_i-sum_j\)。同时我们可以把余下的 \(i\) 和 \(j\) 并入前缀和,设 \(f_i=sum_i+i\)。

式子最后化简为 \(dp_i=\min\limits_{0\le j< i}(dp_j+(f_i-f_j-1-L)^2)\)

复杂度是 \(O(n^2)\) 的,但很容易发现,我们只需要知道先前的最小值即可。这种式子并不可以用单纯的单调队列维护,因为化简出来后会带有类似 \(f_ig_j\) 的量,这说明当前枚举的 \(i\) 是会影响 \(j\) 的贡献的。对于这样的式子,我们考虑斜率优化。

这次用更为普遍的方法。考虑对于两个量 \(j_1\) 和 \(j_2\),有 \(j_1<j_2\),并且我们假设 \(j_2\) 的贡献更优,可以得出

\(dp_{j_2}+(f_i-f_{j_2}-1-L)^2<dp_{j_1}+(f_i-f_{j_1}-1-L)^2\)

\(dp_{j_2}+f_i^2-2f_i(f_{j_2}+1+L)+(f_{j_2}+1+L)^2<dp_{j_1}+f_i^2-2f_i(f_{j_1}+1+L)+(f_{j_1}+1+L)^2\)

\(dp_{j_2}+(f_{j_2}+1+L)^2-dp_{j_1}-(f_{j_2}+1+L)^2<2f_i(f_{j_2}-f_{j_1})\)

\(\frac{(dp_{j_2}+(f_{j_2}+1+L)^2)-(dp_{j_1}+(f_{j_2}+1+L)^2)}{f_{j_2}-f_{j_1}}<2f_i\)

将 \(g_i=(f_{i}+1+L)^2\)

\[\frac{(dp_{j_2}+g_{j_2}^2)-(dp_{j_1}+g_{j_1}^2)}{f_{j_2}-f_{j_1}}<2f_i \]

这就是经典的斜率优化的形式,即满足此不等式的 \(j_1\) 与 \(j_2\),有 \(j_2\) 优于 \(j_1\)。有了这个斜率式,我们就可以用单调队列维护一个斜率单调增加的点集,因为这个斜率式反映了我们需要维护下凸包

针对这一题,由于 \(2f_i\) 是单调增加的,所以答案只可能出现在队首,我们直接利用上面的斜率式维护队首,满足就弹出,最后取出队首更新答案即可。如果 \(2f_i\) 不是单调增加的,我们就只能用二分查找找到最优决策点了,复杂度会多一个 \(\log\)。

补充:

这里需要补充一下为什么对于点集,我们只需要维护上凸包或下凸包即可。

为了方便表述,我们定义 \(K(i,j)\) 表示 \(i,j\) 两点的斜率

接下来我们要证明: 如果 \(i,j,k\) 满足 \(k<j<i\) 且 \(K(i,j)<K(j,k)\),则 \(j\) 永远不可能成为最优解。

  • 当 \(K(i,j)\le 2\cdot sum_i\) 时,选择 \(i\) 优于选择 \(j\),排除 \(j\)。
  • 当 \(K(i,j)>2\cdot sum_i\) 时,选择 \(j\) 优于选择 \(i\),又因为 \(K(j,k)>K(i,j)>2\cdot sum_i\),所以选择 \(k\) 优于选择 \(j\),同样排除 \(j\)。

所以假设成立。

此处引用 「算法笔记」斜率优化 的证明

P3648 [APIO2014] 序列分割

把长度为 \(n\) 的序列分成 \(k + 1\) 块,也就是分 \(k\) 次,使得分最高。

对于这一类题,我们都可以探索一个性质。因为题目中的分块是无序的,可能同一分割方法不同分割顺序得到的答案不一样,这不利于我们转移,但我们可以猜测分割顺序是不影响答案的。

设有 \(abc\) 序列,

先分割 \(ab\),再分割 \(bc\) 的答案是 \(a(b+c)+bc=ab+ac+bc\)

先分割 \(bc\),再分割 \(ab\) 的答案是 \((a+b)c+ab=ab+ac+bc\)

所以顺序不影响答案,我们只需要求出如何分割即可。

从左到右按顺序考虑,可设状态 \(dp_{i,j}\) 为到第 \(i\) 块切割了 \(j\) 次的最大价值,转移我们枚举第一次分割(因为贡献更好计算),同时是最右边的分割,

\(dp_{i,j}=\max\limits_{1\le k<i}(dp_{k,j-1}+sum_k\times(sum_i-sum_k))\)

求先前最大值,并且带有 \(sum_ksum_i\) 的量,我们考虑斜率优化。首先空间上可以滚动数组,记 \(g_k=dp_{k,j-1}\)

原式为 \(dp_{i}=\max\limits_{1\le k<i}(g_{k}+sum_k\times(sum_i-sum_k))\)

设 \(k\) 和 \(k^{\prime}\),\(k<k^{\prime}\),有 \(k^{\prime}\) 的贡献优于 \(k\),即

\(g_{k}+sum_k\times(sum_i-sum_k)<g_{k{\prime}}+sum_{k^{\prime}}\times(sum_i-sum_{k^{\prime}})\)

\(g_k+sum_ksum_i-sum_{k}^2<g_{k^{\prime}}+sum_{k^{\prime}}sum_i-sum_{k^{\prime}}^2\)

\(g_k-sum_k^2-g_{k^{\prime}}+sum_{k^{\prime}}^2<sum_i(sum_{k^{\prime}}-sum_k)\)

\(\frac{(g_k-sum_k^2)-(g_{k^{\prime}}-sum_{k^{\prime}}^2)}{-sum_k-(-sum_{k^{\prime}}^2)}<sum_i\)

这可以看作点集 \((-sum_k,g_k-sum_k^2)\) 在第三象限组成的下凸包,维护不断增加的斜率即可。

对于询问方案,每次转移时,记录转移点,从后往前不断跳前继转移点即可。

标签:分割,sum,斜率,做题,复杂度,2023,我们,dp
From: https://www.cnblogs.com/FireRaku/p/18147522

相关文章

  • arc166D 做题小计
    线段树做法,拿下你谷最劣解。题意翻译很形象,就不说了。思路最大化最小值,我们很容易想到二分答案。很容易发现,答案具有单调性。我们二分一个答案\(x\),强制每次使用的区间长度都不小于\(x\),然后判断可行性。现在问题转化为怎么判断一个答案\(x\)是否可行。我们发现,如果枚......
  • [dpdk] rte_flow
     以下内容直接来自官网文档的整理。更精准的描述请阅读文档:https://doc.dpdk.org/guides/prog_guide/rte_flow.html一rte_flow是干嘛的一组用来创建自定义规则的api,该规则可以改变网络流量的命运,以及查询计数。 二规则啥样1match+actionmatch包括:两类,A报文内容(按......
  • 区间dp思想
    删数链接:https://www.luogu.com.cn/problem/P2426题目描述有\(N\)个不同的正整数\(x_1\),\(x_2\),...,\(x_N\)排成一排,我们可以从左边或右边去掉连续的\(i\)\((1\lei\len)\)个数(只能从两边删除数),剩下\(N-i\)个数,再把剩下的数按以上操作处理,直到所有的数都被删......
  • ICPC2023沈阳站(B C D E I J K M)
    本场金牌线为六题前一半,这里提供难度前八题的题解。本场真是个细节场,银牌金牌题细节都相当地多。ICPC2023沈阳站C:作为签到题,对这种赛制熟悉不熟悉直接区分了一血时间。谔谔。#include<algorithm>#include<iostream>#include<cstring>#include<cstdio>#include<cmat......
  • 概率dp四题(青蛙跳、吸血鬼、rating、k小姐的点赞之谜)
    青蛙跳Description有\(n\)个荷叶按顺序依次排列开,编号为\(1\)到\(n\),现在有只青蛙在编号为\(n\)的荷叶上。它现在自由愉快的跳跃,如果他在编号为\(i\)的荷叶上,它会等概率的跳到编号为\([1,i]\)的荷叶上,求它跳到编号为\(1\)的荷叶上的期望步数。Samples53.083333......
  • [dp 小计] wqs 二分
    天才算法!国外叫Alienstrick(外星人trick),真的太强了。其实是因为IOI2016Aliens这道题考了这个算法才开始普及。解决问题wqs二分一般用来解决如下问题。给定\(n\)个数,求强制选\(m\)个的价值最大。如果不是强制选\(m\)个,这类问题很好做。现在问题就是怎么取消......
  • Pyecharts制作动态GDP柱状图
    学习使用pyecharts制作动态柱状图使用csv模块进行csv数据文件处理importcsvfrompyecharts.chartsimportBar,Timelinefrompyecharts.optionsimport*frompyecharts.globalsimportThemeTypedefdealCSVFile():"""读取处理csv数据文件:retu......
  • DP10RF001一款200MHz~960MHz 低功耗(G)FSK/OOK无线收发芯片应用无线遥控工控设备无线
    产品概述.DP10RF001是一款工作于200MHz~960MHz范围内的低功耗、高性能、单片集成的(G)FSK/OOK无线收发机芯片。内部集成完整的射频接收机、射频发射机、频率综合器、调制解调器,只需配备简单、低成本的外围器件就可以获得良好的收发性能。芯片支持灵活可设的数据包格式,支持自动应......
  • 超低功耗Sub-1G收发芯片DP32RF002 M0内核(G)FSK/OOK 无线收发机的32位SoC芯
    产品概述DP32RF002是深圳市动能世纪科技有限公司研制的基于ARMCortex-MO+内核的超低功耗高性能的、单片集成(G)FSK/OOK无线收发机的32位SoC芯片。工作于200~960MHz范围内,支持灵活可设的数据包格式,支持自动应答和自动重发功能,支持跳频操作,支持FEC功能,同时内部集成了完整的射频......
  • 获取AWS lightsail Windows server RDP密码
    场景创建lightsail的linuxserver时已经生成SSHkey,建立Windows的实例(Instance)时,并未提示输入管理员密码。登录时,找密码登录,提示“DecipheryourpasswordYouusedthe"keyname"keywhenyoucreatedthisinstance.Seetheinstructionstodecipherthepasswordfromthe......