专栏文章

DP solution set

算法·理论参与者 7已保存评论 6

文章操作

快速查看文章及其快照的属性,并进行相关操作。

当前评论
6 条
当前快照
1 份
快照标识符
@mibyrpbp
此快照首次捕获于
2025/11/24 01:00
3 个月前
此快照最后确认于
2025/12/01 20:18
3 个月前
查看原文

传统计数

传统的分析性质+计数 DP,DP 本身不算难。需要明确判据并根据判据 DP。

[AGC002F] Leftmost Ball 紫(计数 DP)

白球对其他球的限制是:任意后缀(或前缀)的白球个数大于等于颜色种类数。
只限制后缀也能保证正确性,因此从后往前 DP,并且需要把白球个数和颜色种类数计入状态。
具体地,我们钦定 ii 第二次出现的位置小于 i+1i+1 第二次出现的位置。
dpi,jdp_{i,j} 表示填入了 ini\sim n 的数字,填入了 jj 个白球的方案数,每次钦定填入的第一个球在最前面以避免重复。
转移是简单的。

树上 DP

树上 DP 的本质在于合并子树,明确 DP 的状态定义后通常转移并不困难。

联考 NOIP 模拟 T3(计数,随机选数转排列,树上 DP)

形式化题意

考虑这样的一个算法:
给定树 TT,初始时图上所有节点均为白色,计数器 c=0c = 0。不断进行下列操作直到 TT 上不存在白色点:
  • 从树上所有白色点中随机选择一个点 uu,将 uu 以及与其有边直接相连的所有白色点 vv 全部染黑,并令 cc+1c \leftarrow c + 1
最后返回 cc 作为答案。
林尼想知道,对于给定的树 TT,算法的正确率是多少,即有多大概率使得算法返回的答案确实为原图的最大独立集大小。
你只需要输出概率对 109+710 ^ 9 + 7 取模的结果。
nn 为树的大小,n5000n\leq 5000
部分分:最大独立集唯一。

题解

首先按照经典套路,将白点选数改成随机一个排列,按排列选数,然后跳过不合法的选择。
这样最终方案和排列构成双射,至此我们完成了概率转计数。
考虑部分分,分析排列性质,首先考虑限制所有独立集的点都被选,但是发现这样不好限制。
考虑限制所有不在独立集的点都不被选,由于是独立集,所以这样限制也是充要的。
记排列的逆为 rkrk,独立集为 SS。一个不是 SS 的点 uu 不被选,当且仅当存在一个邻域 vvrkv<rkurk_v<rk_u,且 vSv\in S,这样我们就得到了充要条件。
fu,i,0/1f_{u,i,0/1} 表示子树 uu 内,离散化后的 rku=irk_u=i,且如果 uSu\notin S,是否有 uu 的儿子 vv,满足 vSrkv<rkuv\in S\land rk_v<rk_u
转移时如果 uSu\in S,则直接合并 fv,x,1f_{v,x,1},对于 fv,x,0f_{v,x,0},我们需要做的问题如下:
  • 合并两个数列,第一个序列的第 xx 个位置在第二个序列的第 yy 个位置前。
转移时进行树上背包:
fu,i+k,1=fu,i,1×fv,j,0×(i+k1k)×(szui+szvkszui)(j>k)f_{u,i+k,1}'=f_{u,i,1}\times f_{v,j,0}\times \binom{i+k-1}{k}\times \binom{sz_u-i+sz_v-k}{sz_u-i}(j> k)
如果 uSu\notin S,大致同理,当第三维由 0011 时做上面的转移,只是把 >> 改为 \leq
可以前缀和优化做到 O(n2)O(n^2)
考虑正解。
类似 DP of DP,先使用最大独立集 DP(gu,0/1g_{u,0/1} 表示 uu 子树 uu 是否选择)建立出最优转移 DAG,然后在 ff 后加一维表示 uu 是否属于 SS,只进行 DAG 上边的转移即可。
我们无需显式建立 DAG,只需要在往上做独立集 DP 时只转移最优情况即可。

状压 DP

既然都是状压 DP 了,状态设计和转移都十分暴力,主要难点在于把题目转化为可以 DP 的形式,即有限的状态。

QOJ4314(多项式奇偶分离,(多项式)状压 DP)

形式化题意

给定一个所有系数都 {0,1}\in \{0,1\}n1n-1 次多项式 PP,求 PkP^{k} 中有多少项的系数为奇数。
n20n\leq 20k1016k\leq 10^{16}

题解

考虑快速幂,有两个方向:
  • 计算 f2kf^{2^k},卷积。
  • 通过 fif^{i} 递推到 f2if^{2i}f2i+1f^{2i+1}
第一个方向的卷积不保证项数,倒闭。
考虑第二个方向,根据模二性质,有:
f2i(x)=fi(x2)f2i+1(x)=fi(x2)f(x)f^{2i} (x)=f^{i}(x^2) \\ f^{2i+1} (x)=f^{i}(x^2)f(x)
考虑 DP,我们无法处理第二个式子后面的 f(x)f(x)
正常 DP 思路,把不能处理的东西记到状态里,记 dpi,gdp_{i,g} 表示 kk 的高 ii 位的 f?(x)f^{?}(x) 乘上 g(x)g(x)popcnt\text{popcnt}gg 是一个多项式(额好像记多项式为状态不是很正常,但确实能做)。
推式子。
f2i(x)g(x)=fi(x2)g(x)f^{2i}(x)g(x)=f^{i}(x^2)g(x)
由于是 popcntpopcntf(x2)f(x^2)f(x)f(x) 等价,考虑将 gg 奇偶分离,将式子写成 f(x2)f(x^2) 的形式。
f2i(x)g(x)=fi(x2)g(x)=fi(x2)[F(x2)+xG(x2)]fi(x)F(x)+xfi(x)G(x)f^{2i}(x)g(x)=f^{i}(x^2)g(x)=f^{i}(x^2)[F(x^2)+xG(x^2)]\equiv f^{i}(x)F(x)+xf^{i}(x)G(x)
奇偶项相加互不干扰,即 dpi,gdpi1,F+dpi1,Gdp_{i,g}\leftarrow dp_{i-1,F}+dp_{i-1,G},其中 F,GF,Ggg 奇偶分离后的两个多项式。
f2i+1f^{2i+1} 同理,每次奇偶拆分会将项数除以二,因此 gg 的项数不超过 nn,状压存 gg 即可。
需要预处理所有 gg 所拆成的 FFGG

AT_agc020_f [AGC020F] Arcs on a Circle 黑(状压 DP,实数处理技巧,DP 转移处理)

小数处理技巧:暴力枚举 NN 个数小数部分的排名。
枚举完把小数部分离散化一下变为 NCNC 个点,状压 DP 即可。
注意要先让最长的线段放置,要不然可能存在线段绕环大半圈后回来的情况,此时 DP 有后效性。
DP 时按照线段的左端点为顺序转移,避免计重。

子问题划分

有些 DP 题的难点在于子问题划分,可以通过手玩样例或者考虑去掉一个元素后解的形态来确定子问题。

CF1792F1 Graph Coloring (easy version) *2700(补图性质,计数 DP)

红蓝互为补图,因此红蓝都不连通的情况不存在。
第一反应是对红蓝都连通的情况容斥,但是发现这样根本做不了。
不妨转化为有一个颜色不连通,对恰有一种颜色不连通的情况计数,此时另一种颜色必定连通。
fif_i 表示 ii 个点,每个子集满足蓝色不连通的方案数,枚举与 11 相连极大的蓝色连通块 SS,保证跨过 SS[1,i]/S[1,i]/S 的集合中不存在红蓝都连通的集合,这个连通块需要满足红色不连通,用“极大”保证不重。
由于是"极大连通块",所以两个集合之间的连边方案是固定的,只能是红色边。
恰好,这样的连边方式满足【跨过 SS[1,i]/S[1,i]/S 的集合中不存在红蓝都连通的集合】,因此逻辑以一种很奇怪的方式自洽了。
由于对称性,红色不连通方案数 == 蓝色不连通方案数。
fi=j=1n1(i1j1)×fj×(2fnj[nj=1])f_i=\sum\limits_{j=1}^{n-1}\binom{i-1}{j-1}\times f_{j}\times (2f_{n-j}-[n-j=1])
乘二是因为当集合有边时可以任意交换红蓝。

CF251E Tree and Table *3000(子问题划分,分类讨论)

注意到树的度数不超过 44
当树是一条链时,分类讨论,一定是从某点开始走到头,拐回来走一段之后开始走折线,简单讨论后发现答案是 2(n2n+2)2(n^2-n+2),注意链可以翻转,翻转后算不同的方案。
对于一般情况,先放入一个三度点,设其为树根。
当根存在一个儿子是叶子时,两侧被划分为了两个子问题,否则枚举儿子放在哪边,也是两个子问题。
容易发现我们有两种不同的子问题:
  • 矩形左上(或左下)有一个点,要用这个点的子树铺成一个矩形。
  • 矩形左上和左下各有一个点,要用这两个点的子树铺成一个矩形。
设第一种为 F(u)F(u),第二种为 G(u,v)G(u,v)
先讨论 F(x)F(x)(假设 xx 在左上)。
  • 如果 xx 没有儿子,则无解。
  • 如果 szx=2sz_x=2,答案为 11
  • 如果 xx 有两个儿子,则枚举两个儿子是怎么放置的,设放在下面的为 uu,右边的为 vv
    • uu 的儿子个数为 00,答案加上 F(v)F(v)
    • 否则,uu 的儿子个数必为 11,答案加上 G(v,sonu)G(v,son_u)
  • 如果 xx 有一个儿子:
    • 若子树 xx 为一条链,答案为 szx2\frac{sz_{x}}{2}
    • sonxson_x 放在下面,答案加上 F(sonsonx)F(son_{son_x})
    • 否则找到最浅的有两个儿子的点,设其为 yy,一直直走直到遇到这个点停止,枚举哪个点向下拐。
      • 如果向下拐的点(设其为 z=sony,0z=son_{y,0})有一个儿子,判断一下它能不能把第二行的前缀填满,如果可以,答案加上 F(sony,1)F(son_{y,1})
      • 否则,枚举一下哪个儿子(设其为 sonz,0son_{z,0})能把第二行的前缀填满,如果可以,答案加上的 G(sony,1,sonz,1)G(son_{y,1},son_{z,1})
再讨论 G(x,y)G(x,y)
  • 如果 x,yx,y 没有儿子,答案为 11
  • 如果 x,yx,y 各有一个儿子,答案为 G(sonx,sony)G(son_x,son_y)
  • 否则,设 xx 有一个儿子,答案为 F(sonx)F(son_x)
复杂度看似 O(n2)O(n^2),但是提交后直接过了。
因为 G(u,v)G(u,v) 被调用时,两个点的深度差不超过 11,每次 GG 往后转移只有一个后继,而这颗树每层节点个数是 O(1)O(1) 的(似乎最多 66),因此 GG 的总状态是 O(n)O(n) 的。

对解的生成方式 DP

这类 DP 的难点在于设计一个解的生成方式,并根据生成方式 DP。
通常常用的技巧有分布转移和操作后置。

CF2068D Morse Code *3100(DP 贡献计算技巧,分步转移)

感觉还是对 DP 过程中算贡献和分布转移的理解不是很深刻。
对二叉树形态 DP,由于节点放置不连续,显然要按深度 DP。
fi,j,k,lf_{i,j,k,l} 表示 ii 个叶子,jj 个最大深度的叶子,kk 个次大深度的叶子,最大深度是 ll,可以做到 O(n5)O(n^5)
我尝试了在转移时算贡献,但是没有成功,因为最终排名不确定。
事实上,DP 过程中算贡献和分布转移本质上是让每个点实时有贡献,并且固定一些点停止贡献。
对于扩展深度为 22 的点,我们并不用管,因为我们每次只会把点的最浅深度集体下移 11,且点只有在作为最浅深度时才会被固定。
因此这个点的实际深度就是这个点(这个点不存在时我们认为这个点的祖先就是这个点)在状态里的时间 1-1
fi,j,kf_{i,j,k} 表示固定了 ii 个叶子不再扩展,在未固定的叶子中有 jj 个最大深度、kk 个次大深度的叶子。
则要么固定一个次大深度的叶子,要么把次大深度的叶子全部扩展。
fi,j,kfi+1,j,k1fi,j,k+wfi,k,j+kf_{i,j,k}\to f_{i+1,j,k-1} \\ f_{i,j,k}+w\to f_{i,k,j+k}
ww 表示 pp 的第 i+1ni+1\sim n 大的和,答案为 fn,0,0f_{n,0,0},构造方案是简单的。
O(n3)O(n^3)
double 要加上 0.50.5 再转为 int 啊啊啊。

P11802 【MX-X9-T6】『GROI-R3』Graph 紫(分析性质,计数 DP,分步转移)

模拟赛赛时认为 gcd\gcd 必须等于 11,爆零了/ll。
看来考试时要把发现结论的依据写一写(周期结论)。
首先由题意,GG 上任意一个点能走到他自己,因此 GG 是若干环的并。
得到 gcd=1\gcd=1 的错误结论是由于周期结论的直觉,让我们仔细想想这个结论能带给我们什么信息。
对于 GG 上一个长为 lenlen 的环,我们从一个点出发,每次走 kk 步,能得到 gcd(k,len)\gcd(k,len) 个长为 lengcd(k,len)\frac{len}{\gcd(k,len)} 的环。
因此,我们还可以选择把若干个相同长度的环在 GG 中拼在一起
具体地,用若干个长度为 ii 的环拼成长为 i×ti\times t 的大环,则有:
itgcd(k,len)=it=gcd(k,it)gcd(kt,i)=1\frac{it}{\gcd(k,len)}=i \\ t=\gcd(k,it) \\ \gcd(\frac{k}{t},i)=1
因此,设 cntloopicntloop_i 表示环长为 ii 的环数,找到最小的 tt 使得 gcd(kt,i)=1\gcd(\frac{k}{t},i)=1,则 tcntloopit|cntloop_i。这个条件是充要的。
考虑计数 DP。
按链 id DP要记录每个环长的 cntcnt,不可接受。
因此考虑按需要记录的值 DP,这样的好处是我们可以在进行 ii+1i\to i+1 的转移时顺便判断 ii 是否合法。
fi,j,kf_{i,j,k} 表示当前处理长度 ii,长度为 ii 的链有 jj 个,有 kk 个剩余孤点(记录长度 j\leq j 的链个数不容易保证不重不漏,因此选择 == 的定义)。
有一个技巧:我们可以把没有闭合的链各加上一个孤点,以保证 ii 增加时链长 =i=i 的定义。
然而 ww 个孤点在环上的贡献为 w!w!,这要求了我们不能实时计算这个贡献,只能在最后计算。
具体地,我们在转移时不区分孤点,最后乘上孤点个数的阶乘,根据组合意义,这个方法在不存在全是孤点组成的环时正确。
在孤点成环时,我们发现是阶乘和圆排列的区别。
具体地,对于一个 xx,设长为 xx 的全孤点环为 yy 个,则我们多算了 xyy!x^y y! 倍,因为我们要:
  • 除掉圆排列的循环。
  • 当孤点都相同时,这些环没有顺序之分,除掉【最后乘上孤点个数的阶乘】对这些环顺序的错误贡献。
我们尝试用分步的思想来推导 DP 转移方程。
  • 加入原图上的长为 ii 的链,这一步是模拟题意。fi,j,kfi,j+cntlinki,kf_{i,j,k}\to f'_{i,j+cntlink_i,k}cntlinkicntlink_i 指原图上长为 ii 的链个数)。
  • 求出最小的 tt 使得 gcd(kt,i)=1\gcd(\frac{k}{t},i)=1,加入一些环,满足长为 ii 的总环数是 tt 的倍数,形成环有如下两种方式:
    • 选出若干孤点构成环。
    • 将已有的链闭合形成环。
    即枚举 c,lc,lfi,j,k(jl)icc!fi,jl,kic\frac{f'_{i,j,k}\binom{j}{l}}{i^c c!}\to f''_{i,j-l,k-ic},转移条件是 (cntloopi+l+c)modt=0(cntloop_i+l+c)\bmod t=0
  • 把没有闭合的链加上一个孤点,fi,j,kfi+1,j,kjf''_{i,j,k}\to f_{i+1,j,k-j}
记孤点个数为 cnt1cnt1,答案即为 fn+1,0,0×cnt1!f_{n+1,0,0}\times cnt1!
分析复杂度,有 ijn,ickn,ljij\leq n,ic\leq k\leq n,l\leq j
ij,ic,ilij,ic,il 各贡献 ni\sum \frac{n}{i}kk 贡献 nn
总复杂度 O(nn3i3)=O(n41i3)=O(n4)O(n\sum\frac{n^3}{i^3})=O(n^4\sum\frac{1}{i^3})=O(n^4)
在推导时,我们发现没有转移条件时完全可以把【选出若干孤点构成环】和【将已有的链闭合形成环】分步转移做到 O(n3)O(n^3)
把条件改成 cntloopil=c(modt)-cntloop_i-l=c\pmod t 仍然可以分步转移,枚举 cntloopi+lmodtcntloop_i+l\bmod t 的值做若干遍即可。

P9523 [JOIST 2022] 复制粘贴 3 / Copy and Paste 3 紫(区间 DP,简化转移,操作后置,调和级数)

这个操作看起来就很符合区间 DP 的形式。
dpl,rdp_{l,r} 表示 [l,r][l,r] 的生成代价,一个区间形如 ???[Y]???[Y]??[Y]?????? [Y] ??? [Y] ?? [Y] ???,即若干个字符加上若干个剪切板字符串。
转移时枚举第一个 [Y][Y] 的位置可以做到 O(n4)O(n^4),考虑简化转移,要求转移到的区间末尾是 [Y][Y],再用区间扩展转移实现往后加 ??,可以做到 O(n3)O(n^3)
我们的贡献和字符串 YY[l,r][l,r] 的出现次数有关,考虑反过来刷表,枚举出现次数以做到调和级数的复杂度。
dpl,r+B+C(occur(s[x:r],s[l:r]))+A(rx+1occur(s[x:r],s[l:r])×(rl+1))dpx,r  (xl)dpl,rdpl,r+1dp_{l,r}+B+C(\text{occur}(s[x:r],s[l:r]))+A(r-x+1-\text{occur}(s[x:r],s[l:r])\times (r-l+1))\to dp_{x,r}\ \ (x\leq l) \\ dp_{l,r}\to dp_{l,r+1}
我们发现我们还要处理 [x,r][x,r] 开头 ?? 段,可以上数据结构做到双 log\log,但是无法通过。
**Hint:**将操作 A 后置,即允许在字符串开头添加字符。
dpl,rdpl1,rdp_{l,r}\to dp_{l-1,r},这相当于把操作 A 在操作序列中插入到了头部。
这样我们只需要枚举 O(nlnn)O(n\ln n)xx,复杂度 O(n2lnn)O(n^2\ln n)
后置操作 A 本质上就是我们新加入了一个操作 D 是在串前加入字符,然后发现和原问题等价。
实际写的时候写的是 dpl,rdpl,x  (x>r)dp_{l,r}\to dp_{l,x}\ \ (x>r),总之是等价的。

匹配型 DP 与贡献后置

匹配型 DP 即对一个匹配方案 DP,较为常见的对排列 DP 本质上也是匹配型 DP。
通常转移是无法保证已经扫到的元素全部被匹配的,需要记录失配数,即把失配点的贡献后置处理。
当然贡献后置的技巧不只局限于记录失配数,所有无法当前处理的转移 case 都可以考虑使用贡献后置。

P4859 已经没有什么好害怕的了 紫(二项式反演,匹配型 DP)

t=nk2t=\frac{n-k}{2}
二项式反演,钦定 iia>ba>b,其余不管的方案数为 gig_i,则 ans=i=tn(it)(1)itgians=\sum\limits_{i=t}^{n} \binom{i}{t}(-1)^{i-t}g_i
匹配型 DP 求出 gig_i,具体地令 fi,jf_{i,j} 表示 b[1i]b[1\sim i] 中有 jj 个匹配了比他小的 aa 的方案数,限制由严到松,因此转移是容易的。
最后 gi=fn,i×(ni)!g_i=f_{n,i}\times \color{red}(n-i)!,不要忘了这个系数,表示其余的匹配任意。

P14424 [JOISC 2014] 邮戳收集 / Collecting Stamps 紫(贡献后置,状态简化)

题意配图

题解

经过一个邮戳台有四种方式:
  1. u+eu+e
  2. d+vd+v
  3. u+vu+v
  4. e+de+d
考虑在东向电车行走的路线,可以跳过若干邮戳台不走,其余走 u+vu+v,或者走 u+eu+e 上去,把没走的邮戳台挑若干个走 e+de+d,挑一个走 d+vd+v 下去。
注意 u+eu+ed+vd+v 可能被走多次。
e+de+dd+vd+v 的贡献后置,在转移时加上额外贡献,令 fi,j,kf_{i,j,k} 表示到了站台 ii,有 jj 个邮戳台被指定走 d+ed+ekk 个邮戳台被指定走 d+vd+v,转移时加上 2T(j+k)2T(j+k) 的贡献,枚举当前有 k1k_1u+eu+e 和前面匹配,k2k_2d+vd+v 和后面匹配。可以做到 O(poly(n))O(\text{poly}(n))
简化状态,由于在上面只能往左走,因此只有存在没有匹配的 d+vd+v 时才能匹配存在未匹配的 d+ed+e,另一方面,我们做 u+eu+ed+vd+v 的匹配时实际上顺便把 d+ed+e 走了,可以把向上和向下理解称括号匹配,当存在未匹配的左括号时可以选择 d+ed+e
dpi,jdp_{i,j} 表示到了站台 ii,有 jj 个未匹配的 d+vd+v 的最小代价。
转移时,可以选择 u+vu+v,当 j>0j>0 时可以选择 d+ed+e。枚举这个站台有 k1k1u+eu+e 和前面匹配,k2k2d+vd+v 和后面匹配,转移到 dpi+1,jk1dp_{i+1,j-k1}dpi+1,j+k2dp_{i+1,j+k2}
前缀和优化可以做到 O(n2)O(n^2)

[CSP-S 2025] employ 紫(贡献后置,匹配型 DP)

会了贡献后置但是不会推式子,要练组合数学了。
首先这是一个匹配问题,具体匹配方案如下:
  • 当这个点值为 11,且有人被录取时,匹配一个 cc 大于【之前失败人数】的人。
  • 当这个点值为 11,且无人被录取时,匹配一个 cc 小于等于【之前失败人数】的人。
  • 当这个点值为 00 时,匹配任意一个人。
匹配问题考虑 DP + 贡献后置,按 cc 的值 DP,令 fi,j,kf_{i,j,k} 表示 [1,i][1,i] 拒绝 jj 个人,匹配了 kkc>jc>j 的人的方案数。
jj 既代表当前拒绝的人数,也代表 cc 匹配到的值,在本题中可以合并成一维状态。
cnti=x[cx=i],prei=x=1icntxcnt_i=\sum\limits_x [c_x=i],pre_i=\sum\limits_{x=1}^{i} cnt_x
[si=1]18\sum [s_i=1]\leq 18 的部分分启示我们暴力枚举每个点的状态:
  • si=1s_i=1,实际有人录取,则有:fi,j,kfi+1,j,k+1f_{i,j,k}\to f_{i+1,j,k+1}
以上是我自己会的部分 qwq\tiny{\text{以上是我自己会的部分 qwq}}
  • si=1s_i=1,实际无人录取,有一个 j\leq jcc 与当前位匹配,有 jj+1j\to j+1,那么我们要枚举 c=j+1c=j+1 的有多少和前面匹配,则有:w=1min(k,cntj+1)fi,j,k(kw)(cntj+1w)w!(prej(ik))fi+1,j,kw\sum\limits_{w=1}^{\min(k,cnt_{j+1})} f_{i,j,k}\binom{k}{w}\binom{cnt_{j+1}}{w}w!(pre_j-(i-k))\to f_{i+1,j,k-w}
  • si=1s_i=1,有 jj+1j\to j+1,那么我们要枚举 c=j+1c=j+1 的有多少和前面匹配,并且我们要讨论与这个位置匹配的人的 cc 是否 >j+1>j+1
    • c>j+1c>j+1 时:w=1min(k,cntj+1)fi,j,k(kw)(cntj+1w)w!fi+1,j+1,kw+1\sum\limits_{w=1}^{\min(k,cnt_{j+1})} f_{i,j,k}\binom{k}{w}\binom{cnt_{j+1}}{w}w!\to f_{i+1,j+1,k-w+1}
    • cj+1c\leq j+1 时:w=1min(k,cntj+1)fi,j,k(kw)(cntj+1w)w!(prej+1(i(kw)))fi+1,j+1,kw\sum\limits_{w=1}^{\min(k,cnt_{j+1})} f_{i,j,k}\binom{k}{w}\binom{cnt_{j+1}}{w}w!(pre_{j+1}-(i-(k-w)))\to f_{i+1,j+1,k-w}
讨论到这里,我们发现已经没有必要枚举每个点状态了,直接转移所有可能的情况就做完了。
注意统计答案时必须要求 k=nprejk=n-pre_j

P7213 [JOISC 2020] 最古の遺跡 3 黑(分析性质,贡献后置,子问题划分)

当然贡献后置的技巧不只局限于记录失配数,所有无法当前处理的转移 case 都可以考虑使用贡献后置。
bib_i 表示第 ii 个石柱最后的高度,根据题意有且仅有 bAi0b_{A_i}\neq 0
考虑如何从 hhbb,首先我们发现最后 bb 会变成一个排列,因为每次都会固定 1N1\sim N 恰一个数不动。排列继续进行操作不会受任何影响,所以进行 NN 次操作等价于不断进行操作。
模拟一下发现每个时刻相同高度的石柱最多有两个,且下一步前面的石柱高度会减一。
一个石柱从不变到减少的状态转变,当且仅当:
  • 后面某个石柱上一步减少到了这个高度,之后由于相同高度的石柱最多有两个,因此后面的石柱被固定,前面的石柱开始减少。
  • 石柱减少到某个位置,后面不存在高度相同的石柱,石柱被固定。
容易发现,高度为 ii 的石柱位置会不断前移,严谨的说:
  • 一个后缀中如果出现了高度为 ii 的石柱,则这个后缀在接下来的时间里必定含有高度为 ii 的石柱。
  • 反之,一个后缀中没有出现高度为 ii 的石柱,则这个后缀在之前也没有出现过高度为 ii 的石柱。
从后往前考虑,如果 ii 后面的 bb 出现了 1hi1\sim h_i 的所有整数,那么 bi=0b_i=0,否则,记 xx 表示满足以下条件的最大整数:
  • x[1,hi]x\in[1,h_i]
  • 不存在 bj=xj>ib_j=x\land j>i
那么 bi=xb_i=x,因为上述条件表示 ii 后面高度为 xx 的石柱从未出现过,那么 ii 必定变为高度为 xx 的石柱。
现在我们知道了 NNb0b\neq 0 的位置,我们需要对 hh 计数。
但是我们不知道如何根据 bbhh 的方案数,所以直接对 bb 的排列进行计数的方法是没有前途的。
考虑直接对 hh DP,模拟 bb 的生成方式,具体地令 fi,jf_{i,j} 表示 [i,n][i,n] 中,bb 出现了 1n1\sim n 的所有整数,所对应 hh 的方案数。
我们要干的事情是对 hh 进行匹配计数,1N1\sim N 的每个数可以被匹配两次,我们记 cnt0icnt0_i 表示 1i1\sim ib=0b=0 的位置个数。
核心思想:贡献前置、不同高度之间的独立性。
为了计算可以填数的数量,我们把两个相同的 ii 视作本质不同的数字,最后除以 2n2^n
讨论 bib_i
  • bi=0b_i=0 时,hijh_i\leq jfi,jfi+1,j×(2j(j+cnt0ncnt0i))f_{i,j}\leftarrow f_{i+1,j}\times (2j-(j+cnt0_{n}-cnt0_{i})),表示匹配一个 j\leq j 的数,j\leq j 的中 jj 个数已经被用掉并且 b0b\neq 0cnt0ncnt0icnt0_{n}-cnt0_{i} 个数被用掉且 b=0b=0。由于 jj 在从后往前的模拟中递增,限制由紧到松,因此转移逻辑自洽。
  • bi0b_i\neq 0 时:
    • bi>j+1b_i>j+1 时,考虑延迟贡献这一部分,先不管,fi,jfi+1,jf_{i,j}\leftarrow f_{i+1,j}
    • bi=j+1b_i=j+1 时,如果 [j+2,k][j+2,k] 在后面都出现过且 k+1k+1 没有出现过,则 fi,kfi+1,jf_{i,k}\leftarrow f_{i+1,j}。在这里,我们需要处理之前 hi>j+1h_i>j+1 所带来的贡献。有 j+1bikj+1\leq b_i\leq k
      • 枚举 kk,计算 [j+2,k][j+2,k] 都出现的方案数,这部分与之前的状态独立,因此可以直接乘起来。
      • suf1=(ni)(cnt0ncnt0i)suf_1=(n-i)-(cnt0_n-cnt0_i),我们需要在 jsuf1j-suf_1 个数中选 kj1k-j-1 个数,要求这 kj1k-j-1 个数排成连续段,则有 fi,kfi+1,j(suf1jkj1)(2(kj+1)(kj1))gkj1f_{i,k}\leftarrow f_{i+1,j}\binom{suf_1-j}{k-j-1}(2(k-j+1)-(k-j-1))g_{k-j-1},其中 gig_i 表示 ii 个数最后形成的 bb 排成连续段的方案数。容易发现在这次转移后所有未确定的石柱高度都 >k>k,所以这些石柱不会对当前 kj1k-j-1 个数排成连续段造成任何影响,即 gg 完全独立。
考虑如何求 gig_i。考虑类似区间 DP 的子问题划分结构,枚举最后一个被固定的石柱的高度,gi=j=1igj1gij(i1j1)(2(ij+1)(ij))g_{i}=\sum\limits_{j=1}^{i} g_{j-1}g_{i-j}\binom{i-1}{j-1}(2(i-j+1)-(i-j))
时间复杂度 O(N3)O(N^3),常数极小可以通过。
本题的启示:
  1. 当 DP 遇到无法处理的 case 时(通常是对另一种转移有后效性),不管这个贡献,放到受影响的转移 case 中后置处理。
  2. 不知道如何拆分子问题时,可以考查个别数,例如如果把第一个/最后一个满足某条件的元素拿走会怎么样。

数据结构/贪心思想辅助 DP

这类问题难点不在 DP,主要在于使用贪心思想简化问题或使用数据结构大力优化。
传统的 DS 优化 DP 太常见了,这里放一个不是很常见的套路。

CF1610G AmShZ Wins a Bet *3300(贪心,trie+hash,DP)

形式化题意

给定一个括号串 SS,每次可以删去一个 ()() 子序列,求能得到的最小字典序字符串。
1S3×1051\leq |S|\leq 3\times 10^5

题解

对于本题,只有后缀的局部最优等于整体最优。
首先删除一个左括号时,贪心地,肯定删除最前的右括号。
并且我们可以把删去一个左括号等价为删除这个左括号连续段中最后一个左括号。
那么我们可以把题目中的子序列改为子串。
那么我们删去一个左括号时,一定会删除后面极短的合法括号串。
删一个合法括号串后变为子问题,考虑 DP,令 fif_i 表示 [1,i][1,i] 的答案,记 lstilst_i 表示第 ii 个字符向前极短的括号串为 [lsti,i][lst_i,i](若 Si=(S_i=(lsti=1lst_i=-1)。
Si=(S_i=( 时转移是唯一的。
Si=)S_i=) 时有:
fi=min(flsti1,fi1+Si)f_{i}=\min(f_{lst_i-1},f_{i-1}+S_i)
把所有 ff 放到 trie 上,相当于动态加叶子的 trie,可以求 lca\text{lca} 快速比较字典序。
但是,这是错的。
因为我们对前缀进行 DP 时,局部最优不等于整体最优,比如 (((()))))(((())))),则应该不删任何东西,但是我们会令 f8=f_8=\emptysetf9=)f_9=)
本质上是我们删完前缀的后缀时不知道后面拼接的是啥,后面拼的字符可能影响最优性,局部最优不等于整体最优。
考虑对后缀 DP,则:
fi=min(fnxti+1,Si+fi+1)f_{i}=\min(f_{nxt_i+1},S_i+f_{i+1})
需要在反 trie 上加叶子,那么在反 trie 上倍增倍增维护哈希即可。
后缀 DP 的好处是删完后缀的前缀时知道后面拼什么,而前面拼什么不影响后缀串的最优性,因此是正确的。

评论

6 条评论,欢迎与作者交流。

正在加载评论...