2023.01.25 题目选讲 学习笔记
讲课人:罗恺。
这讲都是好题,只可惜有的题我感觉很难用文字表达清楚做法,把人讲懂,比如 \(10,15\) 两题。之后再接再厉!
1.CF1361E James and the Chase
首先一个点是好的当且仅当以这个点为根搜出一棵 dfs 树,dfs 树中只有树边和返祖边。由后者推前者是容易的,由前者推后者只需从反证法的角度考虑,发现不可能有前向边或横叉边。
于是我们可以先随 \(100\) 个点并判断它们是不是好点,如果都不是我们可以认为好点数没有 \(20\%\),否则我们找到了一个好点。
以我们找到的好点为根搜出一棵 dfs 树,现在尝试判断每个点是否是好点。对于点 \(u\)(\(u\) 不为根),首先 \(u\) 到子树内的点都恰有一条路径,其次 \(u\) 子树内指向子树外的边至少有一条(因为强连通),如果 \(u\) 子树内指向子树外的边有至少两条,那么 \(u\) 一定存在至少两条简单路径到 \(u\) 的父亲;否则,\(u\) 子树内指向子树外的边恰有一条,设指向了 \(v\)。
接下来我们说明 \(u\) 到每个点的简单路径数和 \(v\) 到每个点的简单路径数相同,从而 \(u\) 是好的当且仅当 \(v\) 是好的。首先 \(u,v\) 到 \(v\) 子树内每个点的路径都恰有一条;其次如果 \(v\) 不是根,则 \(v\) 到子树外的点一定要经过 \(v\) 子树内的某条返祖边,并且该路径一定不经过 \(u\) 子树,原因是 \(u\) 子树内没有能到 \(v\) 子树外的返祖边。所以按深度从小到大更新即可。
设随机一个点并判断其是否是好的的次数为 \(B\),则时间复杂度 \(O(Bn)\),空间复杂度 \(O(n)\)。
2.Gym104090E Oscar is All You Need
首先 \(n=3\) 时只有两种可能的结果。不妨设 \(n\ge 4\)。
首先 \(1\) 可以在两步操作内被换到位置 \(1\)。现在假设 \(1\sim k-1\) 已经分别被换到前 \(1\sim k-1\) 个位置了,我们尝试把 \(k\) 换到位置 \(k\)。
此时序列一定形如:\(1\sim k-1,A,k,B\)。其中 \(k\ge 2\),\(A,B\) 含有若干个数(可能为空)。
如果 \(|A|=0\),那就做完了。
如果 \(|A|\ge1,|B|\ge1\),那么可以 \([1\sim k-1][A,k][B]\rightarrow [B,A][k][1\sim k-1]\rightarrow 1\sim k-1,k,B,A\)。
如果 \(|A|\ge2,|B|=0\),将 \(A\) 分成 \(A_1,A_2\) 两部分且两部分都不为空,那么可以 \([1\sim k-1][A_1][A_2,k]\rightarrow [A_2][k,A_1][1\sim k-1]\rightarrow 1\sim k-1,k,A_1,A_2\)。
剩下 \(|A|=1,|B|=0\) 的情况,此时序列前 \(n-2\) 个数一定都归位了。注意到上述过程每将一个数移动到对应位置,花费次数不超过 \(2\),因此该部分不超过 \(2n-4\) 步。
此时,只剩下 \(1\sim n-2,n-1,n\) 和 \(1\sim n-2,n,n-1\) 两种情况,前者已经做完了,后者我们把 \(1\sim n-2\) 划分成 \(A,B\) 两部分且两部分都不为空,那么可以 \([A][B,n][n-1]\rightarrow [n-1][B][n,A]\rightarrow [n][A][B,n-1]\rightarrow [B][n-1,A][n]\rightarrow [n][n-1][A,B]\rightarrow A,B,n-1,n\),花费 \(5\) 步。总步数不超过 \(2n+1\)。
时间复杂度 \(O(n^2)\),可以用 splay 优化到 \(O(n\log n)\),空间复杂度 \(O(n)\)。
3.CF843D Dynamic Shortest Path
首先第一次求最短路就直接用堆优化的 dijkstra 求。求完之后设 \(1\) 到 \(i\) 的最短路长度为 \(dis_i\),那么对于一条边 \(i\rightarrow j\) 我们把它的边权 \(w(i,j)\) 做一个调整得到一张新图:\(w'(i,j)\leftarrow w(i,j)+dis_i-dis_j\)。这个调整的作用是,如果一条边在最短路图上那它的边权为 \(0\),否则边权一定大于 \(0\)。同时,新图上 \(i\) 到 \(j\) 的任何一条路径都比原图上对应路径的长度大了 \(dis_i-dis_j\),所以新图上求出的最短路即为原图求出的最短路。有了这个调整后我们每进行一次 \(2\) 操作就把新图建出来并求出 \(1\) 号点到每个点的最短路。注意到新图上 \(1\) 到每个点最短路长度都不超过 \(n\),因此我们就可以用 \(n\) 个队列来代替 priority_queue,从而在 \(O(n+m)\) 时间内求出新的最短路(即 \(O(n+m+V)\) 的 dijkstra,\(V\) 是源点到每个点最短路的最大值)。
时间复杂度 \(O(m\log m+q(n+m))\),空间复杂度 \(O(n+m)\)。
4.UOJ#703. 赵云八卦阵
首先发现 \(a_i\) 是可以异或上它前面的任意一个数 \(a_j\) 的,只要依次操作 \(j+1,\cdots,i,i-1,\cdots,j+1,j+2,\cdots,i,i-1,\cdots,j+2\)。所以第 \(i\) 个数的取值范围为 \(a_1\sim a_{i-1}\) 的线性基能表示出的数异或上 \(a_i\) 得到的集合,记该集合为 \(S_i\)。
我们从前往后求出每个前缀的线性基。如果 \(a_i\) 已经能被 \(i-1\) 前缀的线性基表示出来了,那么 \(S_i\) 就是 \(i-1\) 前缀的线性基能表示出的数的集合,此时对于任意 \(j<i\) 都有 \(S_j\subseteq S_i\)。那么如果有一种最优解中第 \(i\) 个位置不在 LIS 中但是它前面有数在 LIS 中,我们可以把第 \(i\) 个位置上的数变成它前面最后一个在 LIS 中的数,把前面最后一个在 LIS 中的数从 LIS 中删去并把第 \(i\) 个位置加入 LIS 中,得到的还是最优解。这样,对于每个 \(a_i\) 能被 \(i-1\) 前缀的线性基表示出来的位置 \(i\),要么位置 \(i\) 在 LIS 中,要么前 \(i-1\) 个位置没有数在 LIS 中。
所以,存在一种最优解,满足 \(a_i\) 能被 \(i-1\) 前缀的线性基表示出来的位置选的是一段后缀。那么我们就可以把 \(a_i\) 能被 \(i-1\) 前缀的线性基表示出来的位置依据它们的 \(S\) 集合分成若干个连续段,满足每个连续段的 \(S\) 完全相同,不同连续段的 \(S\) 不同,可以发现一共有 \(O(\log a)\) 个连续段。同时,\(a_i\) 不能被 \(i-1\) 前缀的线性基表示出来的位置也有 \(O(\log a)\) 个。
设 \(a_i\) 不能被 \(i-1\) 前缀的线性基表示出来的位置有 \(m\) 个,依次为 \(1=k_1<k_2<\cdots<k_m\le n\)。从后往前 dp,设 \(dp_{i,j}\) 表示 \([k_i,n]\) 中所有 \(a_i\) 能被 \(i-1\) 前缀的线性基表示出来的位置都选,\([k_i,n]\) 中 \(a_i\) 不能被 \(i-1\) 前缀的线性基表示出来的位置选了 \(j\) 个,此时 LIS 开头的最大值。从 \(dp_i\) 转移到 \(dp_{i-1}\) 时首先要在 \([k_{i-1}+1,k_i-1]\) 这些位置的 \(S\) 中定位出小于某个数的最大数,接下来贪心能选就选,有两种可能,一种是贪心地选的过程中直接把 \(S\) 中剩余能选的元素都选完了,此时直接更新答案;另一种可能是没有选完,此时再枚举 \(k_{i-1}\) 是否选,并在对应的线性基中定位出小于某个数的最大数即可。
\(dp\) 状态数是 \(O(\log^2 a)\) 的,每次在线性基中定位是 \(O(\log a)\) 的,从单个状态转移到之后的状态要定位 \(O(1)\) 次,所以该部分复杂度是 \(O(\log^3 a)\) 的。
总时间复杂度 \(O(n\log a+\log^3 a)\),空间复杂度 \(O(n+\log^2 a)\)。
5.Gym104090B Useful Algorithm
我们对每个 \(0\le k<m\) 求答案,实际上就是求 \(\max_{i,j}\{b_i+b_j|(a_i\bmod 2^k)+(a_j\bmod 2^k)\ge 2^k\}\)。令 \(w_i=\max_j\{b_j|(a_j\bmod 2^k)=i\}\)(如果集合为空则 \(w_i=-\inf\)),也就是求 \(\max_{0\le i,j\lt 2^k}\{w_i+w_j|i+j\ge 2^k\}\)。令 \(v_i=w_{2^k-i}\),变成了求 \(\max_{0\le i,j\lt 2^k}\{w_i+v_j|i\ge j\}\)。\(w,v\) 可以用 multiset 维护,所求式子可以用线段树维护。
时间复杂度 \(O(nm\log n+2^m+qm^2+qm\log n)\),空间复杂度 \(O(n+2^m)\)。
6.Gym104076I Shortest Path
钦定一条边,表示它是我们经过的路径上长度最小的边。不妨设这条边是 \((u,v)\),由于我们经过了它所以我们走的路线一定形如 \(1\) 到 \((u,v)\) 的某个端点,经过 \((u,v)\) 并从另一个端点离开。在 \(1\) 到 \((u,v)\) 的某个端点的路上,如果有 \(u,v\) 外的点 \(x\) 在路线中被经过了至少 \(3\) 次,那么第一次到 \(x\) 后走了一个环又回到了 \(x\),第二次到 \(x\) 后走了一个环又回到了 \(x\),这两个环中如果有偶环可以把偶环删掉,并替换成等量边数的最小边;如果都是奇环可以把两个环都删掉,并替换成等量边数的最小边,此时路径不会变长。所以在 \(1\) 到 \((u,v)\) 的某个端点的路上经过 \(u,v\) 外的每个点的次数都不超过 \(2\),同理从 \((u,v)\) 另一个端点离开到 \(n\) 经过 \(u,v\) 外的每个点的次数也都不超过 \(2\)。也就是说,在所有经过 \(c\) 条边的从 \(1\) 出发到 \(n\) 的路径外,除了 \(u,v\) 外的点都经过不超过 \(4\) 次,那么剩下的次数一定是在 \((u,v)\) 这条边上来回走。
因此在经过的路径边数足够大时,我们一定是从 \(1\) 出发走到一条边 \((u,v)\) 的一个端点,在 \((u,v)\) 上来回走,再从 \((u,v)\) 的某个端点出发走到 \(n\)。同时注意到按照上述流程,真实答案一定会被算上,且如果钦定的边不是路径上最小的边答案并不会更优(因为也是一种合法的走法)。
我们先求出从 \(1\) 和 \(n\) 出发经过 \(0\sim 2n\) 条边到每个点的最短路。接下来枚举一条边 \((u,v)\),枚举从 \(1\) 走到 \((u,v)\) 的哪个端点以及走过的路径边数 \(l_1\),设此时最短路径长度 \(v_1\),枚举从 \((u,v)\) 的哪个端点离开走到 \(n\) 以及走过的路径边数 \(l_2\),设此时最短路径长度 \(v_2\),设总共经过 \(c\) 条边,那么在 \((u,v)\) 上来回走的次数就是 \(c-l_1-l_2\),设 \((u,v)\) 长度为 \(w(u,v)\),则路径长度为 \(v_1+v_2+(c-l_1-l_2)w(u,v)=(v_1-l_1w(u,v))+(v_2-l_2w(u,v))+cw(u,v)\),
显然要选 \(v_1-l_1w(u,v)\) 最小和 \(v_2-l_2w(u,v)\) 最小的两条路径。同时注意到路径长度的形式可以表示成 \(A+Bc\) 的形式,不难想到可以用凸包来维护。
时间复杂度 \(O(nm)\),空间复杂度 \(O(n+m)\)。
7.CF833E Caramel Clouds
我们尝试对每个位置 \(x\) 求出 \(f_x\) 表示 \([0,x]\) 内最长可能的未被覆盖的长度,这样询问就可以通过二分快速回答。进一步地,我们发现只需要求出 \(x\in\{l_i\}\cup\{r_i\}\) 这些位置的 \(f\) 值即可。
离散化后,设第 \(i\) 个位置对应的原坐标为 \(v_i\)。假设我们已经求出了 \(f_{i-1}\),如果 \((v_{i-1},v_i)\) 没被覆盖,那么 \(f_i=f_{i-1}+v_i-v_{i-1}\)。如果 \((v_{i-1},v_i)\) 被覆盖了至少 \(3\) 次,那么 \(f_i=f_{i-1}\)。
如果 \((v_{i-1},v_i)\) 被覆盖了 \(2\) 次,一种情况是 \((v_{i-1},v_i)\) 最终被覆盖,\(f_i=f_{i-1}\);另一种情况是 \((v_{i-1},v_i)\) 最终未被覆盖,\(f_{i}=v_i-v_{i-1}+h_{i-1}\),其中 \(h_i\) 表示只考虑右端点(离散化前)不超过 \(v_i\) 的区间,未被覆盖的长度有多长。
如果 \((v_{i-1},v_i)\) 被覆盖了 \(1\) 次,一种情况是 \((v_{i-1},v_i)\) 最终被覆盖,\(f_i=f_{i-1}\);另一种情况是 \((v_{i-1},v_i)\) 最终未被覆盖,此时要求的东西形如,在右端点不超过 \(v_{i-1}\) 的区间中,去掉某个区间后,新增的未被覆盖的长度,也即只被这个区间覆盖的长度。记当前右端点不超过 \(v_{i-1}\) 的区间中,只被区间 \(j\) 覆盖的长度为 \(g_j\),那么 \(f_{i}\) 还可以为 \(v_i-v_{i-1}+h_{i-1}\)(前提是唯一覆盖 \((v_{i-1},v_i)\) 的区间花费小等于 \(C\))或 \(v_i-v_{i-1}+h_{i-1}+g_j\)(前提是唯一覆盖 \((v_{i-1},v_i)\) 的区间花费和区间 \(j\) 的花费之和小等于 \(C\),且区间 \(j\) 的右端点不超过 \(v_{i-1}\))。
所以我们需要动态维护出每个时刻 \(t\),只考虑右端点不超过 \(v_t\) 的区间时的 \(h\) 和 \(g_i\)。注意到 \(h\) 只在 \(t\) 增加或有位置被覆盖的次数变化为 \(1\) 时才会变化,\(g_i\) 只在有位置被覆盖的次数变化为 \(1\) 或 \(2\) 时才会变化,所以 \(h\) 和 \(g_i\) 的总变化次数都是 \(O(n)\) 的,可以用珂朵莉树维护。
动态维护出 \(g_i\) 后,我们还要在一个时刻查询 \(c_i\) 不超过某个值的 \(g_i\) 的最大值,这可以用单点修改,前缀最大值线段树维护。
时间复杂度 \(O((n+m)\log n)\),空间复杂度 \(O(n)\)。
8.Gym104076B Torch
称走在前面的人为 A,后面的人为 B。设不考虑 B 会被 A 挡住的情况下,\(t\) 时刻 A 的位置为 \(1+u_t\),B 的位置为 \(v_t\),则实际上 \(t\) 时刻 B 的位置为 \(v_t-\max_{0\le i\le t}\{v_i-u_i\}\)。\(u,v\) 是可以 \(O(1)\) 计算的,因此我们只需快速计算 \(\max_{0\le i\le t}\{v_i-u_i\}\)。
注意到可能成为 \(v_i-u_i\) 最大的时刻 \(i\) 一定满足下述条件之一:\(i=t\),或 \(i\) 时 A 的灯暗且 \(i+1\) 时 A 的灯亮,或 \(i\) 时 B 的灯亮且 \(i+1\) 时 B 的灯暗。但这还不够,我们发现还可以把 \(t\ge (a_1+b_1)(a_2+b_2)\) 的情况转成 \(t<(a_1+b_1)(a_2+b_2)\) 的情况,那么所有可能成为 \(v_i-u_i\) 最大的时刻 \(i\) 只有 \(O(a_1+b_1+a_2+b_2)\) 种了。所以把询问的时间(要先转成 \(t< (a_1+b_1)(a_2+b_2)\))和可能成为 \(v_i-u_i\) 最大的时刻 \(i\) 一起排序,就做完了。
时间复杂度 \(O((a_1+b_1+a_2+b_2+n)\log n)\) 或 \(O(a_1+b_1+a_2+b_2+n)\)(取决于排序),空间复杂度 \(O(a_1+b_1+a_2+b_2+n)\)。
9.Gym104076J Skills
注意到开始学一个技能后则这个技能的熟练度不会被减到 \(0\) 及以下,不然不如把之前学习该技能的时间拿去学其它技能,所以在开始学一个技能后我们可以忽略熟练度和 \(0\) 取 \(\max\) 的条件。
我们还发现,开始学一个技能后,中断的时间不会太长。具体来说,不会超过 \(2\sqrt a+O(1)\) 天。比如说我们找到一个技能,它是所有技能中最早中断了 \(2\sqrt a+O(1)\) 天的,那我们在其第一次中断了至少 \(2\sqrt a+O(1)\) 天的这段中断期间中的第 \(\sqrt a\) 天放弃学习别的技能,转而学习这个技能,此时损失是原本要学习的技能能获得的熟练度和那个技能中断的天数,收益至少是接下来 \(\sqrt a+O(1)\) 天每天少损失 \(\sqrt a\) 熟练度。损失中前者是 \(O(a)\) 的,对于后者由于我们找到的技能是最早中断了 \(2\sqrt a+O(1)\) 天的,所以此时原本要学习的技能中断的天数不超过 \(2\sqrt a+O(1)\),因此调整完会更优。
于是就可以 dp 了,记一下当前的天数,这天学习的技能,以及另外两个技能上次学习的时间或者还未开始学,状态数是 \(O(na)\) 的,转移是 \(O(1)\) 的。
总时间复杂度 \(O(na)\),空间复杂度 \(O(a)\)。
10.CF1616G Just Add an Edge
先特判掉 \(1\) 到 \(n\) 本来就有哈密顿路径的情况。
其余情况一定会经过新加的边 \((x,y)(x>y)\)。为了统一起点和终点,我们添加两个点 \(0,n+1\),\(0\) 向所有点都连边,所有点都向 \(n\) 连边。于是最终路径一定形如从 \(0\) 走编号升序且连续的点到 \(y-1\),从 \(y-1\) 走编号升序的点到 \(x\),从 \(x\) 走到 \(y\),从 \(y\) 走编号升序的点到 \(x+1\),从 \(x+1\) 走编号升序且连续的点到 \(n+1\)。
我们考虑把每个点染色:如果走到该点时还未走过新加的边,就染成 \(0\),否则染成 \(1\)。我们把 \(0\sim n+1\) 从左往右从小往大排成一个序列,于是整个序列中最靠右的 \(0\) 就是 \(x\),最靠左的 \(1\) 就是 \(y\)。一种染色方案是合法的,当且仅当每个点有恰好一种颜色,且同种颜色的点按编号升序排序后,编号相邻的点在图中都有边相连。
首先枚举最靠左的 \(1\)(即 \(y\))的位置。从小往大考虑每个数染成什么颜色,可以得到一个比较显然的 dp:设 \(dp_{i,j}\) 表示目前染成 \(0\) 的点中最大的是 \(i\),染成 \(1\) 的点中最大的是 \(j\),此时是否有合法的染色方案。转移枚举下一个点的颜色。总时间复杂度是 \(O(n^3)\) 的。
一个优化是,颜色相同的段我们可以合并转移,具体而言,判断 \(i\sim j\) 能否都染成某种颜色,只需判断 \(i\rightarrow i+1,\cdots,j-1\rightarrow j\) 这些边是否都存在即可。也就是说,在上述的 dp 中,我们只用考虑形如 \(dp_{i,i+1}\) 和 \(dp_{i+1,i}\) 这些状态。于是可以重新设 \(dp_{i,0/1}\) 表示考虑了前 \(i+1\) 个点的颜色,\(i\) 和 \(i+1\) 的颜色不同,\(i\) 的颜色为 \(0/1\),此时是否有合法的染色方案。时间复杂度降为 \(O(n(n+m))\)。
由于此时图中没有哈密顿路径,所以一定可以找到一个最小的 \(l\) 满足 \(l\rightarrow l+1\) 的边不存在,那么此时 \(l\) 和 \(l+1\) 的颜色就一定不相同!所以我们先得到一个结论,\(x\ge l,y\le l+1\)。
注意到我们对 \(x,y\) 有个限制 \(x>y\),实际上我们 dp 的过程中可能出现 \(x\le y\) 的情况只有 \(x=l,y=l+1\),而这种情况能使图产生哈密顿路径当且仅当 \(1\rightarrow2,\cdots,n-1\rightarrow n\) 中只有 \(l\rightarrow l+1\) 不存在。这是很好特判的,即如果有这种情况,我们先把答案 \(-1\)。于是 dp 时就没有 \(x>y\) 的限制了!
再来看看 dp,我们是先枚举了 \(y\),再 dp 求答案。但是实际上,对于不同的 \(y\),转移都是一样的,唯一不同的只有初值,我们必须减少无用的重复转移。
这就又要用到我们前面说的「\(l\) 和 \(l+1\) 的颜色一定不相同」了,可以发现,\(dp_{1}\sim dp_{l-1}\) 不可能一步转移到 \(dp_{l+1}\sim dp_{n}\)。这启发我们可以不枚举 \(y\),因为 \(dp_{1}\sim dp_{l-1}\) 都要先转移到 \(dp_{l}\),再通过 \(dp_{l}\) 转移到 \(dp_{l+1}\sim dp_{n}\)。所以 \(dp_{1}\sim dp_{l-1}\) 和 \(dp_{l+1}\sim dp_{n}\) 其实是相对独立的,我们只关心 \(dp_{1}\sim dp_{l-1}\) 能使 \(dp_l\) 变成什么样,和在只给定的 \(dp_l\) 的情况下 \(dp_{l+1}\sim dp_n\) 能变成什么样。
后者可以直接枚举所有可能的 \(dp_l\) 初值,用上述 dp \(O(n+m)\) 求出最后 \(dp_{l+1}\sim dp_n\) 的形态。前者我们可以从后往前求出 \(f_{i,o1,o2}\) 表示如果一开始只有 \(dp_{i,o1}=1\),能否让 \(dp_{l,o2}=1\),这也可以 \(O(n+m)\) 求。最后把两边 dp 的信息合并一下就好了。
时间复杂度 \(O(n+m)\),空间复杂度 \(O(n+m)\)。
11.CF643G Choosing Ads
如果 \(p\ge 51\) 就是要求区间严格众数,可以用线段树维护摩尔投票的方法做。具体的,线段树上每个节点维护对应区间摩尔投票后留下的数,以及它的次数,那么合并两个节点可以 \(O(1)\) 合并。
拓展到 \(p\ge 20\) 的情况,我们时刻维护不超过 \(\lfloor\frac{100}{p}\rfloor\) 个数,表示摩尔投票筛出来的数,每加入一个数 \(x\),如果 \(x\) 出现在了这些数中,则将 \(x\) 的次数 \(+1\);否则,如果当前维护的数个数小于 \(\lfloor\frac{100}{p}\rfloor\),则把 \(x\) 加进维护的数中并将 \(x\) 的次数设为 \(1\);否则,将这 \(\lfloor\frac{100}{p}\rfloor\) 个数的次数全部 \(-1\),并把次数变为 \(0\) 的数从我们维护的数中删掉。接下来我们说明为什么出现次数至少占数的总数的 \(p\%\) 的数一定会出现在我们最终维护的数中。
考虑什么时候我们维护的数的次数会减少,只有当我们维护的数恰有 \(\lfloor\frac{100}{p}\rfloor\) 个,且加入一个不同于它们的数时才会发生。那么我们从某种数 \(x\) 的角度看,如果 \(x\) 某次加入时 \(x\) 已经存在或维护的数个数小于 \(\lfloor\frac{100}{p}\rfloor\),则 \(x\) 的次数会 \(+1\);如果 \(x\) 某次加入时维护的数恰有 \(\lfloor\frac{100}{p}\rfloor\) 个,则 \(x\) 的次数不会 \(+1\),但是可以认为是 \(x\)「抵消」掉了 \(\lfloor\frac{100}{p}\rfloor\) 个不同于 \(x\) 的数;如果加入某个不同于 \(x\) 的数 \(y\) 时维护的数恰有 \(\lfloor\frac{100}{p}\rfloor\) 个,且 \(x\) 在其中,则 \(x\) 的次数会 \(-1\),但是也可以认为是 \(x\)「抵消」掉了 \(\lfloor\frac{100}{p}\rfloor\) 个不同于 \(x\) 的数。假设数的总数为 \(l\),\(x\) 出现了 \(t(\ge p\%\cdot l)\) 次,那么不同于 \(x\) 的数有 \(l-t\) 个,\(x\) 最多会「抵消」\(\lfloor\frac{l-t}{\lfloor\frac{100}{p}\rfloor}\rfloor\) 次,故 \(x\) 最终的次数至少为 \(t-\lfloor\frac{l-t}{\lfloor\frac{100}{p}\rfloor}\rfloor\),下面我们说明 \(t-\lfloor\frac{l-t}{\lfloor\frac{100}{p}\rfloor}\rfloor>0\)。
由于 \(t-\lfloor\frac{l-t}{\lfloor\frac{100}{p}\rfloor}\rfloor=\lceil\frac{t\lfloor\frac{100}{p}\rfloor-l+t}{\lfloor\frac{100}{p}\rfloor}\rceil\),只需说明 \(\lceil\frac{t\lfloor\frac{100}{p}\rfloor-l+t}{\lfloor\frac{100}{p}\rfloor}\rceil>0\),即 \(t\lfloor\frac{100}{p}\rfloor-l+t>0\)。
\(t\lfloor\frac{100}{p}\rfloor-l+t\ge \frac{pl}{100}\lfloor\frac{100}{p}\rfloor-l+\frac{pl}{100}=\frac{1}{100}(pl\lfloor\frac{100}{p}\rfloor-100l+pl)=\frac{1}{100}(l(100-100\bmod p)-100l+pl)>0\)。
所以至少占数的总数的 \(p\%\) 的数一定会出现在我们最终维护的数中。
上述过程还是可以用线段树维护,只不过合并的复杂度变成了 \(O(\lfloor\frac{100}{p}\rfloor^2)\)。
总时间复杂度 \(O(n\log n\lfloor\frac{100}{p}\rfloor^2)\),空间复杂度 \(O(n\lfloor\frac{100}{p}\rfloor)\)。
顺带一提,如果不允许输出询问的答案可以包含错的数或重复输出,那就变成了 CF840D Destiny 这题,求出可能答案后用主席树算一下真实出现次数即可。
12.UOJ#749. 【UNR #6】稳健型选手
称先手为 A,后手为 B。先考虑单组询问怎么做。对于每个长度为 \(l\) 的前缀,B 至少会选 \(\lfloor\frac{l}{2}\rfloor\) 个元素,并且对于满足该条件的所有选法,A 每次选他要选的编号最小的元素,最终一定能选到他想要的所有数。于是可以反悔贪心。
同时由上述条件知道如果区间长度是奇数,那么最后一个数一定是归 A 的,所以我们只用考虑区间长度为偶数的情况。
假设我们已经用堆维护出了只看长度为偶数的区间 \([l,r]\) 中的数,A 选的数的集合和 A 不选的数的集合。如果往区间末尾加两个数,区间内原先不选的一定不选,原先选的可能不选,那么根据贪心我们应该把这两个数加入 A 选的数的堆,并把 A 选的数的堆中最小的数移到 A 不选的数的堆中。如果往区间开头加两个数,区间内原先选的一定选,原先不选的可能选,那么根据贪心我们应该把这两个数加入 A 不选的数的堆,并把 A 不选的数的堆中最大的数移到 A 选的数的堆中。这个过程如果用不删莫队实现可以做到 \(O(n\sqrt q\log n)\)。
实际上,我们可以把上述做法拓展一下。考虑分治,假设当前分治区间为 \([l,r]\),令 \(mid=\lfloor\frac{l+r}{2}\rfloor\),根据上述做法可以发现 \([l,mid]\) 的任意两个长度差为 \(1\) 的后缀,选的数至多有 \(1\) 个不同,所以我们可以用主席树维护 \([l,mid]\) 的每个后缀选了哪些数,同理也可以维护 \([mid+1,r]\) 的每个前缀选了哪些数。然后考虑回答所有 \(l\le l1\le mid\lt r1\le r\) 的询问,我们已知 \([l1,mid],[mid+1,r1]\) 中选的数和不选的数的集合,那么合并在一起后,\([l1,mid]\) 中不选的数一定仍不选,选的数可能变不选,\([mid+1,r1]\) 中选的数一定仍选,不选的数可能变选。同时可以发现,\([l1,mid]\) 中新的不选的数一定是之前选的数中最小的几个,\([mid+1,r1]\) 中新的选的数一定是之前不选的数中最大的几个,于是可以主席树上二分。
时间复杂度 \(O(n\log^2n+q\log n)\),空间复杂度 \(O(n\log n)\)。
13.CF1764G3 Doremy's Perfect DS Class (Hard Version)
设询问 \((l,r,k)\) 得到的结果为 \(f(l,r,k)\)。
如果 \(n\) 是奇数,设当前已经确定 \(1\) 在 \([l,r]\) 中,令 \(mid=\lfloor\frac{l+r}{2}\rfloor\),询问 \((1,mid,2)\) 和 \((mid+1,n,2)\),考虑 \(f(1,mid,2)-f(mid+1,n,2)\),显然除以 \(2\) 下取整结果相等的只有 \(2\) 和 \(3\),\(4\) 和 \(5\),……,\(n-1\) 和 \(n\) 这些对。讨论可知,如果 \(1\) 在 \([1,mid]\) 中,则 \(f(1,mid,2)-f(mid+1,n,2)=1+\frac{mid-1-(n-mid)}{2}\),否则 \(f(1,mid,2)-f(mid+1,n,2)=\frac{mid-(n-mid-1)}{2}-1\),于是可以根据 \(f(1,mid,2)-f(mid+1,n,2)\) 的值在 \(20\) 次内二分出 \(1\) 的位置。
如果 \(n\) 是偶数,我们仿照 \(n\) 是奇数的情况做,发现除以 \(2\) 下取整结果相等的对有 \(2\) 和 \(3\),\(4\) 和 \(5\),……,\(n-2\) 和 \(n-1\),而 \(n\) 和 \(1\) 一样是独立的,没有数和它配对(除以 \(2\) 下取整相等)。同样询问 \((1,mid,2)\) 和 \((mid+1,n,2)\),讨论可知,我们可以根据 \(f(1,mid,2)-f(mid+1,n,2)\) 的值,判断出三种情况:1. \(1,n\) 都在 \([1,mid]\);2. \(1,n\) 都在 \([mid+1,n]\);3. \(1,n\) 一个在 \([1,mid]\),一个在 \([mid+1,n]\)。前两种情况直接得到了 \(1\) 的位置,但是第三种还没法区分哪个是 \(1\) 哪个是 \(n\)。于是对于第三种情况我们再问一次 \((1,mid,n)\)(或者 \((mid+1,n,n)\))便可知道 \(1\) 和 \(n\) 谁在前面了。注意到我们只需要求一次 \(1,n\) 的相对前后关系,所以这么做最坏要询问 \(21\) 次。
考虑把多出来的 \(1\) 次优化掉。当 \(n\) 是偶数且 \(l+1=r\) 时,我们原本要花 \(2\sim 3\) 次才能询问出 \(1\) 在哪,现在我们直接讨论。如果此时还没确定 \(1,n\) 的相对前后关系,那么 \(1,n\) 必定都在 \([l,r]\) 中,那么询问 \((1,l,n)\) 或者 \((l+1,n,n)\) 可以确定出哪个是 \(1\) 哪个是 \(n\),最坏花费 \(19\) 次。如果已经确定了 \(1,n\) 的相对前后关系,那么 \([l,r]\) 中的 \(2\) 个数一定是一个 \(1\) 和一个 \([2,n-1]\) 中的数 \(x\),注意到此时我们已经知道了 \(f(1,l-1,2),f(l,n,2),f(1,r,2),f(r+1,n,2)\),如果和 \(x\) 配对的数在 \([1,l-1]\) 中,那么我们询问 \((1,l,2)\) 看 \(f(1,l,2)\) 是否等于 \(f(1,l-1,2)\) 便可知道 \(1\) 的位置;同理,如果和 \(x\) 配对的数在 \([r+1,n]\) 中,那么询问 \((r,n,2)\) 看 \(f(r,n,2)\) 是否等于 \(f(r+1,n,2)\) 即可。最坏询问次数 \(20\) 次。
14.UOJ#708. 【UER #10】磁球与磁棍
一个连通块合法,当且仅当连通块内所有叶子到块内深度最浅的点的距离都是偶数,或者连通块内所有叶子到块内深度最浅的点的距离都是奇数且深度最浅的点不是叶子。
不难想到树形 dp。设 \(dp_{i,c,0/1,0/1}\) 表示 \(i\) 子树内,划分成了 \(c\) 个连通块,\(i\) 所在的连通块中所有叶子到 \(i\) 的距离模 \(2\) 都为 \(0/1\),\(i\) 到父亲的边是否有断掉,此时是否有合法划分方案。时间复杂度是 \(O(n^2)\) 的。
我们称一个点为自由点,当且仅当它到连通块内所有叶子的距离都是奇数。
可以发现,我们划分出的连通块要么大小为 \(1\),要么大小至少为 \(3\),后者中一定有自由点。
对于一个 \(2\) 度的自由点,我们可以将其两条出边都断掉,连通块数相比之前增加 \(2\)。
对于一个至少 \(4\) 度的自由点,我们可以将其任意两条出边断掉,连通块数相比之前增加 \(2\)。
对于一个 \(3\) 度的自由点,我们可以将其任意一条出边断掉,连通块数相比之前增加 \(1\)。
如果有至少两个 \(3\) 度的自由点,那么我们任选两个 \(3\) 度点各断一条出边,连通块数相比之前增加 \(2\)。
也就是说,只要有自由点,并且自由点不是一个 \(3\) 度点,那么就存在方案使得连通块数相比之前增加 \(2\)。
注意 \(k=n\) 时有解,\(k=n-1\) 时无解,\(k\le n-2\) 时合法划分方案中一定有自由点。自由点只有一个 \(3\) 度点的情况,当且仅当 \(k=n-3\) 时可能出现。
也就是说,\(k\le n-2\) 且 \(k\neq n-3\) 时,如果存在划分方案有 \(k\) 个连通块,则也存在划分方案有 \(k+2\) 个连通块。
现在我们只关心已知连通块数的奇偶,最少能划分出几个连通块。设 \(dp_{i,0/1,0/1,0/1}\) 表示 \(i\) 子树内,划分成了偶数/奇数个连通块,\(i\) 所在的连通块中所有叶子到 \(i\) 的距离模 \(2\) 都为 \(0/1\),\(i\) 到父亲的边是否有断掉,此时能划分出的最少连通块数。
时间复杂度 \(O(n)\),空间复杂度 \(O(n)\)。
15.CF1007E Mini Metro
我们先来理清一下题目的时间线:对于自然数 \(i\),在 \([i,i+1)\) 这些时间里发生的事有,\(i+0.5\) 时刻来车接人,\(i+1-\epsilon\) 时刻站台来人。
首先发现有车接第 \(i\) 站的人时,前 \(i-1\) 站的人都已被接完。其次可以在最后补一个站 \(n+1\),\(a_{n+1}=c_{n+1}=\inf,b_{n+1}=0\),这样每辆车都会满载。
设 \(sa_i=\sum_{j=1}^i a_j\),\(sb_i=\sum_{j=1}^i b_j\),\(sc_i=\sum_{j=1}^i c_j\)。
设 \(f_{i,j,o}\) 表示只考虑前 \(i\) 站,初始第 \(x(1\le x\le i)\) 站有 \(oa_x\) 个人,存活到 \(j\) 时刻,每辆车都必须是满载的(只能载前 \(i\) 站的人),需要的最少车数(如果无解则为 \(\inf\))。
设 \(g_{i,j,o}\) 表示只考虑前 \(i\) 站,初始第 \(x(1\le x\le i)\) 站有 \(oa_x\) 个人,存活到 \(j+0.5\) 时刻,每辆车都必须是满载的(只能载前 \(i\) 站的人),且 \(j+0.5\) 时刻来完车后前 \(i-1\) 站的人都已被接完,需要的最少车数(如果无解则为 \(\inf\))。
答案就是 \(f_{n+1,t,1}\)。下面考虑怎么转移。
情况 \(1\):\(j\) 时刻前没有车从站台 \(i\) 接走人。
此时需满足 \(oa_i+jb_i\le c_i\),即在 \(j\) 时刻第 \(i\) 站的人数不会超限。并且,前 \(i-1\) 站能独立存活到 \(j\) 时刻,即 \(f_{i-1,j,k}\neq \inf\)。如果上述条件都满足,便有转移:
令 \(x=\lceil\frac{osa_{i-1}+jb_{i-1}}{k}\rceil\),这是清空前 \(i-1\) 站至少要的车数。如果 \(xk\gt osa_i+jsb_i\),就有车没满载没法更新 \(g\) 了。否则,现在至少要用 \(x\) 辆车,事实上,用 \(x\) 辆车是能清空前 \(i-1\) 站的,原因是 \(j\) 时刻前没有车从站台 \(i\) 接走人,那么在 \(j+0.5\) 时刻把剩余的车来了即可。
情况 \(2\):\(j\) 时刻前有车从站台 \(i\) 接走人。枚举时刻 \(r(0\le r\lt j)\),表示上次从站台 \(i\) 接走人是在 \(r+0.5\) 时刻。
这种情况的流程如下:
-
前 \(i\) 站独立存活 \(r\) 时刻,并在 \(r+0.5\) 时刻来车让前 \(i-1\) 站全空。
-
如果在 \(r+0.5\) 时刻来车后第 \(i\) 站剩余人数加 \((j-r)b_i\) 比 \(c_i\) 大,即如果 \((r+0.5,j]\) 时间都不来车的话,\(i\) 站台的人数会超限,则在 \(r+0.5\) 时刻继续来车直到 \((r+0.5,j]\) 时间都不来车 \(i\) 站台的人数也不会超限为止。注意要判断车能否满载。
-
从 \(r+0.5\) 时刻来完车开始到 \(j\) 时刻为止来的车都不从站台 \(i\) 接走人。我们要保证在这段时间中前 \(i-1\) 站能独立存活。
-
存活到了 \(i\) 时刻,就可以直接更新 \(f\) 了。对于 \(g\),我们可能要再来一些车使前面的站都清空。此时算出 \(g\) 的下界,类似情况 \(1\) 分析可得下界依然能取到。
这便是大致的流程了。如果你完全理解了情况 \(1\),那么情况 \(2\) 根据上述流程也不难推出(至少我理解流程后已经能写出代码了)。最后情况 \(2\) 转移的式子长成这样(式子没判转移是否合法):
令 \(x=osa_i+rsb_i-g_{i,r,o}*k\),\(y=\lceil\frac{\max(x+(j-r)b_i-c_i,0)}{k}\rceil\),\(z=g_{i,r,o}+y+\lceil\frac{(j-r)sb_{i-1}}{k}\rceil\)。
时间复杂度 \(O(nt^2)\),空间复杂度 \(O(nt)\)。