2024.6
1. qoj8003 Anti-Plagiarism
不难想到 dp,即 dpx,a,bdpx,a,b 表示第一棵树以 xx 为根的子树能否和第二棵树中以 aa 为根时 bb 的子树匹配,其中 a,ba,b 相邻。
想想怎么转移,就是把子树的 dp 值算出来,然后对于 xx ,枚举一对 a,ba,b ,发现这里其实就是二分图匹配:右边的点是 xx 的儿子,左边的点是 bb 除了 aa 以外与之相邻的 cc 。然后如果 dpv,b,c=1dpv,b,c=1 就从 vv 往 cc 连边。看左边是否满流即可。最后再枚举一个 tt ,进行类似的过程来 check 以 xx 为根的子树是否和以 tt 为根时的第二棵树匹配。
虽然直接这样复杂度是对的,但是常数巨大,难以通过。考虑我们计算出以 tt 为根时做的匹配后,如果满了就直接输出 Yes 了;否则我们通过这个来计算所有 dpx,a,tdpx,a,t 的值:如果 a,xa,x 代表的左部点不是最大匹配的必经点,就说明 dpx,a,t=1dpx,a,t=1 了。这是可以在 dinic 结束后顺便求得的。这样,我们做匹配的次数变成原来的 1313 ,容易通过。
2. qoj8005 Crossing the Border
一个 trick 吧。首先 O(3n)O(3n) 的 dp 是容易的,即 sT≤WsT≤W 且 TT 最高位大于 SS 时用 dpSdpS 更新 dpS∪TdpS∪T 。考虑利用折半来优化这个过程,我们把状态表示成 LRLR 的形式,表示前 n2n2 个状态拼起来和后 n2n2 个物品的状态拼起来。核心思路就是:考虑在 LR′LR′ 处同时处理 LRLR 到 L′R′L′R′ 的转移,其中 LL 是 L′L′ 的子集,RR 是 R′R′ 的子集。
枚举 XX 和 YY ,如果 Y=0Y=0 就特殊处理一下;否则先枚举 YY 的子集 Y′Y′ ,其中 Y′Y′ 不能包含 YY 的最高位,把 (sY−Y′,dpXY′)(sY−Y′,dpXY′) 这个 pair 从小到大排序,并处理出前缀的信息和。
再枚举 XX 的超集 X′X′ ,找到第一个值 ≤W−sX′−X≤W−sX′−X 的前缀,用相应的信息更新 dpX′YdpX′Y 即可。注意到这里我们可以预处理出 sY−Y′/sX′−XsY−Y′/sX′−X 排序的结果,所以复杂度即
O(6n2+3n2logn)O(6n2+3n2logn) 。
3. qoj970 Best Subsequence
考虑二分答案 mm,把 ≤m2≤m2 的看成 11 ,>m2>m2 的看成 00 ,那可以发现 11 是会全部取的;对于 00 ,考虑相邻两个 11 之间只能取至多一个 00 ,我们尝试取其中最小的即可。
然后思考怎么快速的计算,首先 11 是容易的,可以主席树处理,对于 00 ,设第一个 11 在 xx ,最后一个 11 在 yy ,我们先看 x,yx,y 间的 00 有多少个被选的。仍然可以主席树处理:考虑 01 序列随着 mm 的变化只会变化 nn 次,考虑这个过程中所有存在过的相邻 11 ,设这两端更大的为 XX ,它们中间最小的数是 YY ,那么 m∈[X+Y,2Y)m∈[X+Y,2Y) 时这一对数就能产生贡献。我们把这一对的贡献放到左边,仍然用主席树查就好了,即查 mm 版本时 [x,y)[x,y) 的和,再加上 11 。最后是 [x,y][x,y] 外的 00 ,这里就单次的算即可。x,yx,y 的位置线段树二分查即可。于是复杂度 O(qlognlogV)O(qlognlogV) 。
4. qoj971 Binary Search Tree
首先离线下来,利用扫描线转化成:对于一个序列,我们会单点修改,以及查询:对于某个位置 xx ,找到序列建出的小根笛卡尔树中 xx 最深的权值 ≤v≤v 的祖先 ff ,并计算 ff 到根的点的权值和。
很 ez,我们找到 i≤xi≤x 中最大的 ii 使 ai≤vai≤v ,i≥xi≥x 同理,那么 ff 就是两个位置中 aa 更大的那个。这可以通过线段树二分计算;然后是计算 ff 到根的点权和,你发现某个点 ee 是 ff 的祖先等同于:ee 到 ff 的最小值就是 aeae 。用单侧递归线段树就能计算这个东西了。
5. qoj8686 Zoo Management
其实就是分类讨论一下。对于每一个边双,考虑其中的点双,如果有个点双不是简单环就可以乱排了;否则,如果只是一个单独的简单环,判一下循环同构即可;否则,如果存在偶环就能乱排了;否则看一下排列奇偶性是否符合即可。
6. qoj1812 Interesting Coloring
求出 dfs 树,注意到是边双,如果每个点到根路径的不同颜色 <8<8 ,那对于树边 e1e1,找到一个覆盖它的非树边 e2e2 ,把 e2e2 和 e2e2 两端在树上的路径除掉 e1e1 的部分拼起来即可;对于非树边 e1e1 ,取 e1e1 两端在树上的路径即可。
你发现在范围内这是可以构造出的:记当前 11 到 xx 有 cc 个不同颜色,把儿子按 siz 排序,先尽量用前 cc 种颜色,不能用了再用其他颜色,容易说明在 n≤5555n≤5555 的条件下这是一定合法的。
7. qoj1813 Joy with Permutations
对前 44 个数问 44 次中位数,若其排序后是 a1,a2,a3,a4a1,a2,a3,a4 ,我们就能得到 a2,a3a2,a3 的值(记为 l,rl,r) 以及 a1,a2a1,a2 下标的集合(记为 i1,i2i1,i2),a3,a4a3,a4 下标的集合(记为 i3,i4i3,i4)。如果 n=4n=4 的话你发现再用两次比较操作就得到答案了,我们尝试添加剩下的数并保持这个形式。现在加入 xx ,查询 i1,i3,xi1,i3,x 的中位数 mm,如果 l<m<rl<m<r 说明 x=mx=m 了,直接退出,否则我们不妨令 x≤lx≤l 。
如果 m=lm=l 说明 ai1=lai1=l ,让 i1:=xi1:=x ,再计算新的 ll ,算 i2,i3,xi2,i3,x 的中位数即可。
否则,说明 ai2=lai2=l ,让 i2:=xi2:=x ,l:=ml:=m 即可。这样就做完了。
8. P8478 「GLR-R3」清明
把 ∏a′i∏a′i 拆开。即,a′i=xi−m,i+xi−m+1,i+⋯+xi,ia′i=xi−m,i+xi−m+1,i+⋯+xi,i ,我们枚举其中一个 xpi,ixpi,i 即可。考虑所有情况下 ∏xpi,i∏xpi,i 的和怎么算,记 sisi 是有多少个 pxpx 的值等于 ii ,令 g=min(m,n−i)+1g=min(m,n−i)+1 ,把 (ai+g−1si+g−1)(ai+g−1si+g−1) 乘起来即可。
不难得到一个 O(2mpoly(n,m))O(2mpoly(n,m)) 的做法,即从 nn 到 11 枚举 ii ,记录 [i+1,i+m][i+1,i+m] 中哪些的 pp 没有确定,然后每次选其一个子集把 ii 填上,类似于子集卷积的优化即可做到 O(2mm2n)O(2mm2n) 。
尝试寻找 O(2n−mpoly(n,m))O(2n−mpoly(n,m)) 的做法平衡复杂度,你发现我们不需要记录 [1,m+1][1,m+1] 的状态,它们是本质相同的:转而记录其个数即可。复杂度就是 O(2n−mm3n)O(2n−mm3n) 。
平衡一下即可,复杂度 O(2n2n3.5)O(2n2n3.5),实际上 m≤16m≤16 时用算法一,否则用算法二即可。
9. qoj5478 Quiz Contest
经典问题。相当于有多项式 Fi,GiFi,Gi 满足 GiGi 项数 ≤Fi≤Fi 且所有 FiFi 的项数和 =m=m ,对于 i=1,2,…,ni=1,2,…,n 求 ∑hk[xk](Gi∏j≠iFj)∑hk[xk](Gi∏j≠iFj) 。
考虑分治 NTT,求出每个节点的 FF 的乘积。然后考虑 solve([l,r],H)solve([l,r],H) 表示 [l,r][l,r] 的答案都可以表示成 ∑hk[xk](Gi∏j≠i,l≤j≤rFj)∑hk[xk](Gi∏j≠i,l≤j≤rFj) 的形式。递归下去即可,即 solve([l,mid],H/F(mid+1,r)),sol([mid+1,r],H/F(l,mid))solve([l,mid],H/F(mid+1,r)),sol([mid+1,r],H/F(l,mid)) ,这里的 '/' 是差卷积。保留有效的项数即可。复杂度 O(mlog2m)O(mlog2m) 。
10. qoj2070 Heavy Stones
先想想暴力怎么写,这其实就是形如两条链的树做 exchange argument。
想想怎么优化,你考虑对两边分别把能往上合并的部分合并了,再合并两边的结果,即如果一个段的平均数小于其父亲的段的平均数,就给这两段合并了。考虑 ii 增大时前缀的段的变化,不难发现变化段数是均摊 O(1)O(1) 的,每次是删除栈尾的一些段再加入一段。那其实后缀的段的变化也是 O(1)O(1) 的。于是把所有出现过的段离散化下来,用线段树维护对两边合并的结果就好了。复杂度 O(nlogn)O(nlogn) 。
11. qoj8683 Bridging the Gap
先分析一下我们走的形式,有以下三种运法:
运满 mm 个人,其中有一些不动了,有一些以后为了拿手电筒会返回;
运 <m<m 个人,这些人之后全都会返回;
把所有人都运过去且不返回。
其中不动了的人我们会取权值尽量大的,回去的人会取尽量小的。现在假设我们已经确定了每次运的人数以及要回去的人数形成的二元组 (a1,b1),(a2,b2),…,(ak,bk)(a1,b1),(a2,b2),…,(ak,bk) ,那首先 ∑(bi−1)=−1∑(bi−1)=−1 ;其次我们直接令其费用为:k∑i=1sbi+tn−i−1∑j=1ajk∑i=1sbi+tn−i−1∑j=1aj ,其中 ss 是 tt 的前缀和。考虑这样为什么是对的,先说明其合法性,虽然确定了每次运的人,但我们完全可以编排这些二元组的顺序,令一个二元组的权值是 bi−1bi−1 ,我们在每个权值 x>0x>0 的二元组后接 xx 个 −1−1 即可,再把剩下的 −1−1 接到最后。每个 −1−1 就使用它前面那个二元组中产生的要返回的人。然后不难说明这是不劣的。
对着这个 dp 就好了,设 dpi,jdpi,j 是只考虑第一种,第三种运法使得 ∑a=n,∑(b−1)=j∑a=n,∑(b−1)=j 的最小代价,在令 fifi 是考虑第二种运法的最小代价,那么答案即为 minj≥0fj+dpn,−1−jminj≥0fj+dpn,−1−j 。直接算 dp 的值是 O(n3)O(n3) 的,注意到 |j|≤nm|j|≤nm 且每次转移 O(m)O(m) ,复杂度即降至 O(n2)O(n2) 。
12. CF1912C Cactus Transformation
题意是给一个仙人掌,每次可以删一条边加一条边,要求得到的还是仙人掌,要构造方案把它变成另一个给定的仙人掌。n≤1000n≤1000 ,要求步数 15000 以内。
先判无解,如果该仙人掌满足 nn 为奇数且 m=32(n−1)m=32(n−1) 就不能操作了,此时直接判两个图是否相同即可。
否则,由于操作可以可逆,尝试把这个仙人掌变成一个特定的结构较简单的图。考虑把其变成一个类似于菊花图的样子:点 11 上面挂了三元环 (1,2,3),(1,4,5),…,(1,2k,2k+1)(1,2,3),(1,4,5),…,(1,2k,2k+1) 以及一些边 (1,2k+2),(1,2k+3),…,(1,n)(1,2k+2),(1,2k+3),…,(1,n) 。
先把图简化,长度 >3>3 的环都容易改成三元环;对于不与 11 相连的割边 (u,v)(u,v) ,删 (u,v)(u,v) 连 (1,v)(1,v) 。然后考虑每次取一个与 11 相连的点 xx 使其度数 >2>2 ,那么有两种情况:
-
(1,x)(1,x) 为割边,取 xx 所在三元环 (x,y,z)(x,y,z) ,割 (y,z)(y,z) 连 (1,y)(1,y) ,割 (x,z)(x,z) 连 (1,z)(1,z) 。
-
不为割边,那么有三元环 (1,x,a)(1,x,a) ,考虑取一度数为 11 且与 11 相连的点 bb ,割 (x,a)(x,a) 连 (a,b)(a,b) 。注意这里就出现了割边 (1,x)(1,x) 。
我们优先找第一种点。那你发现两步以内就能使 11 度数 +2 。11 度数 =n-1 时肯定结束了,否则一定能找到这样的 xx 的,那么就做完了。
13. uoj700 可爱多的字符串
我们要算的大概是 r∑i=l+1wmin(r,i+lcp(i,l)−1)r∑i=l+1wmin(r,i+lcp(i,l)−1) 。
对反串建 sam 并在 fail 树上树剖。字符集是 {0,1}{0,1} 所以这个树是二叉树,在后面会帮到我们。
在每条重链上统计 i,li,l 代表的点的 lca 在这条链上时的贡献,令 didi 是 ii 代表的点往上遇到的第一个重链上的点的 len 值, di=dldi=dl 时,由于是二叉树,只有 ii 或 ll 为重链上的点才合法,这个直接算就好了;现在只考虑 di≠dldi≠dl 的情况。发现 lcp(i,l)=min(di,dl)lcp(i,l)=min(di,dl) ,那就可以把上面的式子写成 min(r,i+di−1,i+dl−1)min(r,i+di−1,i+dl−1) 。
先考虑 i+dl−1i+dl−1 为最小值的情况,这个时候 di>dldi>dl ,我们枚举重链上的每个点 uu,再枚举其轻子树上的挂的询问,此时 ii 一定在 uu 的重儿子 tt 的子树内,所以所有可能的 i+dl−1i+dl−1 其实是节点 tt 代表的最短串的所有出现位置的右端点 -1 的结果。
考虑建出正串的 sam 并找到这个最短串对应的节点,于是问题变成这个节点的所有在 (l+dl,min(r+1,r+dl)](l+dl,min(r+1,r+dl)] 的 endpos 的权值和,这容易排 dfn 转二维数点解决。
再考虑 i+dl−1≤ri+dl−1≤r 且 di<dldi<dl 的部分,直接二维数点做完。
再考虑 i+dl−1>ri+dl−1>r 的部分,相当于 wmin(r,i+di−1)wmin(r,i+di−1) 的和。仍然二维数点做完,不过注意要减去 di=dldi=dl 的贡献。
复杂度 O(nlog2n)O(nlog2n) 。
14. 6.4 T3
题意是给你 n,x,yn,x,y ,问所有可以表示成 axbyaxby 且 ab≤nab≤n 的数的和。x,y≤10,n≤1012x,y≤10,n≤1012 。
判掉 x=yx=y 后,不妨令 x<yx<y 。如果一对 a,ba,b 满足,在某个质数 pp 上 aa 的次数 ≥lcm(x,y)x≥lcm(x,y)x ,那它就是不优的:我们可以让 aa 除掉 plcm(x,y)/xplcm(x,y)/x 并让 bb 乘上 plcm(x,y)/yplcm(x,y)/y 。不难发现合法的 a,ba,b 对应了唯一的数,于是问题转化成,所有满足 ab≤nab≤n 且 aa 质因数分解后次数都小于 k=lcm(x,y)xk=lcm(x,y)x 的 axbyaxby 之和。
容斥处理这个东西,得到答案为 ∑ddkμ(d)F(ndk)∑ddkμ(d)F(ndk) ,其中 F(n)F(n) 是所有 ab≤nab≤n 的对的 axbyaxby 之和,这在预处理所有 [ni][ni] 处 ix,iyix,iy 的前缀和后即可 O(√n)O(√n) 得到。于是复杂度是 O(∑√ndk)O(∑√ndk) ,k=2k=2 时等同于 O(√nlogn)O(√nlogn) ,否则一定更快。预处理通过拉插可以做到 O((x+y)√n)O((x+y)√n) 。
15. qoj6299 Binary String
你发现 01 交替的部分都是比较稳定的,于是不妨分析所有大于 11 的连续段,手玩一下可以发现每次 0 连续段会往左移,1 连续段会往右移,撞上的时候长度会各自减 11 。剩下的过程就比较容易了。
16. qoj6304 Rectangle
只考虑一 xx 两 yy 的情况,其余都是对称或容易的。对 xx 轴扫描线,那么问题变成:每次加入/删除一个线段,同时查询取两个点覆盖所有线段的方案数。
分析一下,不妨计算第一个点位置 ≤A=minr≤A=minr ,第二个点 ≥B=maxl≥B=maxl 的方案数,若 A≥BA≥B 再减掉 (B−A+22)(B−A+22) 即可。
考虑第一个点的位置 xx 对第二个点 yy 的限制,对于一个线段 [l,r][l,r] 我们已经让 x≤r,y≥lx≤r,y≥l 了,所以只是当 x<lx<l 时会有 y≤ry≤r 的限制。
考虑维护一个序列 bxbx ,表示所有满足 x=l−1x=l−1 的线段的 rr 最小值,那么 xx 处 yy 就需要 ≥B≥B ,≤≤ bb 的后缀 min 。不难通过线段树二分算出合法的 xx 的区间。现在相当于是,设 cc 是 bb 的后缀 min 数组,要算 cc 的区间和,这可以用单侧递归的技巧维护。当然这里只是方便陈述,实际上得先离散化下来,计算的时候就稍微有点区别。
复杂度 O(nlog2n)O(nlog2n) 。
17.qoj7610 Bus Lines
记 gxgx 是覆盖 (x,fax)(x,fax) 的路径中最浅的 lca。考虑怎么计算 f(i,j)f(i,j) ,记它们的 lcalca 为 pp ,我们让两端各自从 uu 跳到 gugu ,直到深度 ≤p≤p ,但这样可能算大了一步:如果 i≠p,j≠pi≠p,j≠p ,记 i′i′ 是 ii 往上跳最后一个深度 >p>p 的点,j′j′ 同理,那如果有路径同时经过 (i′,j′)(i′,j′) 答案就减 11 了。
现在考虑统计答案,第一部分虽然比较繁琐但还是简单的,关键在第二部分,考虑从下往上枚举 pp 以及 pp 的儿子 vv ,计算 ii 在 vv 内,jj 在 pp 其他子树内的贡献。考虑此时所有 lca 为 pp 且一端在 vv 内的路径 (X,Y)(X,Y) ,建出 XX 的虚树,容易发现对虚树上的每一段我们只需计算:这一段子树内对应的所有 YY 到根的并的区域的权值和,其中此时一个点 xx 的权值是:有多少 uu 满足不跳到 pp 以上最高能跳到 xx 。写成一个 ∑F(ui)−∑F(lca(ui,ui+1))∑F(ui)−∑F(lca(ui,ui+1)) 的形式,其中 F(u)F(u) 是 uu 到根的权值和,这可以 BIT 维护;合并子树的信息可以线段树合并。这就做完了。复杂度 O(nlogn)O(nlogn) 。
18.学习 miller-rabin
先把 n−1n−1 写成 a∗2ka∗2k (aa 是奇数)。然后找一些质数 pp(OI 范围内取前 1212 个质数即可),如果 pa≡1(modn)pa≡1(modn) 就不管;否则我们看是否存在 j<kj<k 使得 pj≡n−1(modn)pj≡n−1(modn) ,不存在就寄了。
19.学习 pollard-rho
先判 nn 是不是质数,是就结束,否则尝试找到 1<d<n1<d<n 满足 d|nd|n ,递归的分解 dd 和 ndnd 。
考虑一个序列满足 fn=(f2n−1+c)modnfn=(f2n−1+c)modn ,其中 f0f0 和 cc 是随的。这个序列是会出现循环节的,不过在随机条件下大概是 O(√p)O(√p) 级别的,其中 pp 是 nn 的最小质因子。这里找循环节可以 floyd 判圈,其实就是找到第一个 fi=f2ifi=f2i 的地方。然后我们在过程不断看 gcd(n,|fi−f2i|)gcd(n,|fi−f2i|) 是否大于 11 ,大于 11 就可以结束了。否则我们重新随机 f0f0 和 cc 。
但这样因为要多次求 gcd 是不够快的,一个优化是以 B=1,2,4,…,128,128,128,…B=1,2,4,…,128,128,128,… 为块长,考虑对每一块把 ∏|fi−f2i|modn∏|fi−f2i|modn 算出来,再看它和 nn 的 gcd。这样就足够快了。要特判 n=4n=4 ,然后让 c,f0c,f0 在 [3,n)[3,n) 中随。
20.学习支配树
先跑出 dfs 树,不妨按 dfn 序重新编号。然后考虑求半支配点,点 xx 的半支配点是最小的 kk 满足存在 kk 到 xx 的路径使除了 kk 以外所有点 >x>x 。 求半支配点方法:找所有能一步到达 xx 的点 vv ,如果 v<xv<x 就直接更新,否则找 vv 在 dfs 树上所有 >x>x 的祖先(包括自己),用它们的半支配点更新 xx 的支配点。这个我们 xx 从大往小算,维护一个带权并查集即可。
找出半支配点后就容易找支配点了:保留 dfs 树的边,再从 xx 的半支配点向 xx 连边,按 DAG 求支配树的方法求即可,这里因为每个点就两个入度,所以是很好写的。复杂度 O(nlogn)O(nlogn) 。
21. loj3001「JOISC 2015 Day 3」AAQQZ
先讨论对称中心不在被排序的区间内的情况,枚举 [l,r][l,r] 使 [l,r][l,r] 回文且被排序的区间右端点在 l−1l−1 (在另一端可以 reverse 了做)。那么此时要做的就是,有两个序列 A,BA,B,会让 AA 的一个前缀从大往小排序,最后让 lcs(A,B)lcs(A,B) 最大。先考虑 lcs≤lcs≤ 排序前缀长度的情况,我们枚举这个 lcslcs 为 tt ,如果 BB 不单调就结束了,否则维护一个 zz 为最小的 zz 使得 A[1,z]A[1,z] 的可重集包含 B[1,t]B[1,t] ,若存在 >Bt>Bt 的数 pp 使得 A[1,z]A[1,z] 中 pp 的个数 >> B[1,t]B[1,t] 就也结束了。否则用 tt 更新答案。lcs>lcs> 排序前缀长度时,相当于枚举前缀长度,其实就是让上面的 zz 恰好等于 tt ,此时用 t+lcs(A[t+1:],B[t+1:])t+lcs(A[t+1:],B[t+1:]) 更新答案即可。看看上面的过程,容易发现 al−1=ar+1al−1=ar+1 时 [l,r][l,r] 贡献的答案不优,所以可能的 [l,r][l,r] 只有 O(n)O(n) 对。
再讨论对称中心在被排序的区间的情况,可以发现我们的对称中心一定是最大值/最小值构成的区间的中心,否则可以调整,不妨令中心为最大值的中心,枚举右端点 rr 。
从 rr 到 11 找到第一个 >ar+1>ar+1 的值,设其为 cc ,可以发现我们取的区间中 >ar+1>ar+1 的值只能有 cc ,否则就寄了,设 ll 必须大于 LL 。
类似上面的过程枚举 i=r+1,…,ni=r+1,…,n ,如果 ai−1<aiai−1<ai 就结束,否则仍然是维护一个 zz 为最大的 zz 使得 a[z,r]a[z,r] 可重集包含 a[r+1,i]a[r+1,i] 。仍然先看最后组成的回文串不超出 ll 的情况,那我们需要得知 ll 最小是多少,以使包含的 cc 尽量多,每次改变 ii 时对 j∈[ai+1,ai−1]j∈[ai+1,ai−1] ,用 <z<z 的最大位置更新 LL 即可,这容易 vector 实现。再看超出 ll 的情况,先判一下 aiai 在两边出现次数是否相等,此时的 ll 就需要 a[l,z)a[l,z) 都是 cc 。直接枚举这些 ll 就行,总量是 O(n)O(n) 的,预处理 lci,jlci,j 为 R(a[1,i−1])R(a[1,i−1]) 和 a[j+1,n]a[j+1,n] 的 lcp ,用 2(lcl−1,i+1+(i−r))+([l,r]内的c个数)2(lcl−1,i+1+(i−r))+([l,r]内的c个数) 更新答案即可。
22. loj3274. 「JOISC 2020 Day2」变色龙之恋
先想一想 O(n2)O(n2) 做法。对于两个点 u,vu,v ,我们称它们间有边当且仅当 u,vu,v 一起开会时会出现一种颜色。
思考得到这个图后如何计算答案,显然此时每个点度数为 11 或 33 ,为 11 时与之相连的点一定与之同色。否则设相连的点是 a,b,ca,b,c ,它自身是 xx ,三个点与它的关系形如:一个点与它同色,一个点爱 xx ,xx 爱一个点,且不存在双向恋爱关系。 分别询问 (x,a,b)(x,a,b) ,(x,b,c)(x,b,c) ,(x,a,c)(x,a,c) ,不妨令 xx 喜欢 cc ,可以发现 (x,b,c),(x,a,c)(x,b,c),(x,a,c) 都会返回 22 ,(x,a,b)(x,a,b) 会返回 11 。这样我们就找出了 xx 爱的点。找出所有不是双向恋的点所恋的点后,我们再来看 xx ,此时我们已得知哪个点爱 xx 了,我们便找出了与之同色的点。
于是问题变成如何在限定步数内得到这个图。可以发现很多点集,询问它得到的信息是很难用上的,原因是导致颜色个数变化的情况很多,难以分析。但你发现对于一个独立集 SS ,Query(S)Query(S) 一定等于 SS ;对于 x∉Sx∉S ,若 S∪{x}S∪{x} 不是独立集,则 Query(S∪{x})Query(S∪{x}) 一定不等于 |S|+1|S|+1 ,事实上设 xx 与 SS 的点间边数为 pp ,则 QueryQuery 的值基本是 |S|+1−p|S|+1−p ,只是喜欢 xx ,xx 喜欢的点同时在 SS 时多减了 11,但此时 p≥2p≥2 ,于是性质仍然成立。
由于每个点度数都 ≤3≤3 ,再加上上面对简单情况下 query 值的分析,我们尝试拉出一个极大独立集 TT ,可以发现 |T|≥14N|T|≥14N ,原因是 TT 和 U−TU−T 间至多有 3|U−T|3|U−T| 条边,而至少有 |T||T| 条边。找出 TT 和 U−TU−T 间的边,这个就枚举 U−TU−T 中的每个点 uu,把 TT 中的点排成一行,二分出最大的 xx 使得 T[1,x]T[1,x] 和 uu 构成独立集,那么 uu 和 Tx+1Tx+1 就存在一条边。再在 [x+2,|T|][x+2,|T|] 上二分即可。二分前先 check 一下 T[l,r]T[l,r] 和 uu 是否构成独立集。
再令 U:=U−TU:=U−T 做下去即可。
接下来计算一下总询问次数,每条边都来了一次二分,即为 32nlogn32nlogn ;最后还会多询问 3n3n 次,还是很够的了(这里 nn 是 2N2N) 。可以发现“二分前先 check 一下” 这一步很关键,省了很多次二分。
23.CF1616G Just Add an Edge
先想想怎么 check y→xy→x 是否合法,不妨在两侧加两个点 00 和 n+1n+1 ,那首先看 00 能否只通过 i→i+1i→i+1 走到 x−1x−1 ,y+1y+1 能否走到 n+1n+1 。
再看能否取出两条不交路径,分别是 x−1x−1 到 yy 及 xx 到 y+1y+1 ,且覆盖了中间所有点。那这个 dp 一下就好了,记 dpi,jdpi,j 表示 ii 在第 jj 条路径,i+1i+1 在另一条是否可行,转移是容易的。
考虑优化,令 Lu,RuLu,Ru 是只看 i→i+1i→i+1 时能到 uu 的最左点/uu 能到的最右点,于是我们有 x≤R0x≤R0 ,Ln−1≤yLn−1≤y 。可以发现 R0R0 和 R0+1R0+1 所在的路径一定是不相同的,于是考虑以 (R0,R0+1)(R0,R0+1) 为中心对两头分别做上面的 dp ,然后合并一下就好了。
24. qoj875 Arrange The Piranhas
观察 aa 的前缀和 ss,可以发现每次操作要么使某个位置 +1 ,要么使全局 -1。那其实在每一个空位的操作次数就可以固定下来了:先计算 didi 为 bb 的前 ii 个数减去 aa 的前 ii 个数,令其中最小的为 xx ,如果 x<0x<0 则我们就使 d0d0 到 dndn 都加上 −x−x 即可。
接下来只需要判方案合不合法。很简单,每次找一个能操作且 di>0di>0 的空位操作即可,首先这让两边的空位扩大了,其次虽然自己的空位变小了,但迟早也得变小的,严谨一点就是,对于一个合法操作序列,如果它第一步不是操作 ii ,那就把从前往后第一个 ii 移到最前面,显然依旧合法。现在只需要快速模拟这个过程,每次从前往后扫看有没有能操作的位置并尽全力操作。不太会分析这里的复杂度,但至少是很快的。
25. qoj2834 Nonsense
观察组合意义:对 n+1n+1 个球染色,使得中间有一个绿球,左边有 aa 个红球,右边有 bb 个白球,且左边每个白球都会带来 xx 的系数,右边每个白球带来 yy 的系数。
设 F(n,a,b)F(n,a,b) 是此时的答案,那如果 a>0a>0 ,枚举最左边的球颜色是红/白即可得到 F(n,a,b)=F(n−1,a−1,b)+xF(n−1,a,b)F(n,a,b)=F(n−1,a−1,b)+xF(n−1,a,b) ;否则枚举是绿/白即可得到
F(n,a,b)=(nb)yn−b+xF(n−1,a,b)F(n,a,b)=(nb)yn−b+xF(n−1,a,b) 。
同理尝试枚举最右边的球颜色,我们也有 b>0b>0 时 F(n,a,b)=F(n−1,a,b−1)+yF(n−1,a,b)F(n,a,b)=F(n−1,a,b−1)+yF(n−1,a,b) ,否则为 (na)xn−a+yF(n−1,a,b)(na)xn−a+yF(n−1,a,b) 。
对比两式,发现直接相减,我们即可由 F(n−1,a,b−1),F(n−1,a−1,b)F(n−1,a,b−1),F(n−1,a−1,b) 推到 F(n−1,a,b)F(n−1,a,b) 了,这个题就做完了。
另一种做法是注意到 F(n,a,b)=[zn−a−b](11−xz)a+1(11−yz)b+1F(n,a,b)=[zn−a−b](11−xz)a+1(11−yz)b+1 。将其乘上 (1−xz)+xz(1−xz)+xz 或 (1−yz)+yz(1−yz)+yz ,即可得到和上面类似的递推式。
26. qoj6842 Popcount Words
建出询问串的 AC 自动机,设 Sn=Sn−1+S′n−1Sn=Sn−1+S′n−1 ,其中 S′S′ 是指把 SS flip 得到的结果。
那么文本串可以容易拆成 nlogVnlogV 个 SiSi 或 S′iS′i 拼起来的形式。考虑设 f(u,n,0/1)f(u,n,0/1) 是节点 uu 通过 Sn/S′nSn/S′n 后转移到的节点,这个很容易算出来。再设 g(u,n,0/1)g(u,n,0/1) 为我们要对 uu 通过 Sn/S′nSn/S′n 的过程中的所有节点的权值加上这个 g(u,n,0/1)g(u,n,0/1) 。那这个也是容易算的。总之就是利用好文本串的递归性质。
值得一提的是这个文本串还有别的性质,一个很牛的性质是,对于 SnSn 来讲,其所有 00 位置的 kk 次方之和,与所有 11 位置的 kk 次方之和是相等的。(P7468)
27. 6.26模拟 T3
题意是给定 01? 组成的串 SS ,问有多少组 (T0,T1,T2,…,Tn)(T0,T1,T2,…,Tn) 使得 Ti 由 Ti−1 插入 0/1 得到,且 Tn 与 S 匹配。
很天才的题。思考如何计算一个固定串的答案,倒着看。如果我们看成:每次只能删一个连续段的末尾,这是不好处理的。一个天才的想法是,对 0 连续段只能删末尾,1 连续段只能删开头。那这样的话我们每次只能在 01 交界处进行操作,在这二者中选一个删掉。特别的在开头搞个 0,结尾搞个 1,且要求不能删除。
这还是不好处理。我们干脆不要只看固定串,就是说允许问号。那其实就是对于相邻两个字符,如果第一个是 0/? ,第二个是 1/? ,就能把它们合成为一个 ? 。因为此时我们确定了问号的值,所以这样弄是不重不漏的。
接下来就好处理了:设 pi 是第 i 个字符和第 i+1 个字符合并到一起的时间,那你发现 Si 就是对 pi−1 和 pi 的大小关系的限制!问题完全变成了一个经典老题:LOJ575。复杂度 O(nlog2n) 。
28. gym102201B Bohemian Rhaksody
不难转化成,取一个最大的在边框内的矩形,不严格包含给定点。把给定点对 x 排序,考虑 x 轴严格包含了点 [l,r] ,那贡献就是,(xr+1−xl−1) 乘上,把 0,yl,yl+1,…,yr,M 排序后相邻之差的最大值。这里就有个 O(n2) 做法:枚举 r ,l 从 1 到 r 推并利用链表更新这个值即可。
接下来考虑优化,我们分治,那么利用上面的方法就能从左右分别得到一个三元组集合 (l,r,h) ,我们要做的就是从左右集合中各取一个三元组,求 (hi+hj)(min(ri,rj)−max(li,lj)) 的最大值。但要注意这里的三元组集合有很好的性质:没有两组 (l,r) 相交不包含,且对于 A 包含 B 的情况,A 的 h 小于 B 的 h 。
现在我们讨论 (li,ri) 和 (lj,rj) 是否有包含关系。有包含关系的情况很简单,就是二维偏序算一算。考虑不包含的情况,不妨令 li<lj<ri<rj 。
进行第二次分治,计算 li<lj≤mid<ri<rj 的最优结果,然后你发现不跨过 mid 的三元组就直接分治下去,否则的话,考虑左集合的三元组,若某个右集合的与它交在 mid−1 ,那肯定同时交在 mid 了。于是只需分治到右边,右集合的分治到左边。
现在观察跨过 mid 的这些三元组,发现长得非常好看:对于随着 h 增大,l 增大,r 减小。
把左边的三元组按 h 从小到大排,右边的按 h 从大往小排,你发现两件事情:
如果不考虑包虑 li<lj<ri<rj 的限制,这是具有决策单调性的;
对于右边的某个三元组,其在左边能取的三元组形成了往左移的滑动窗口。
你考虑如果这些窗口都交于一点 x ,那其实问题就不难了:对于 x 两侧,分别单调栈二分即可。虽然它们现在并不交于一点,我们可以贪心的划分这些窗口,使每一段的窗口都会交于一点,复杂度仍然正确。为什么呢,考虑每一段的开头 i1,i2,…,ik 满足 Rij+1<Lij 。然后每一段的代价是 Rij−Lij+1−1<Lij−1−Lij+1−1 。求个和就能发现这个代价和是线性的了。
总复杂度 O(nlog3n) 。教训:传参传 vector 全都写成引用,速度显著提高。
29. qoj5312 Levenshtein Distance
朴素的 dp 可以算 T 所有前缀与 S 的编辑距离,那不妨对 T 的每个后缀分别进行 dp 。思考编辑距离 ≤k 带来了什么效果,首先 |i−j|≤k ,于是单次 dp 降为 O(nk) ,但这仍然是不优的,注意到这个 dp 的值有效范围也是 ≤k 的,经典的,尝试交换值域和定义域,你观察到对于固定的 i−j=x 而言 fi,j 是单调不降的,那我们就设 gx,y 是满足 fi,j≤y 且 i−j=x 的 (i,j) 中 i 最大者,转移就是令 a=gx,y,b=a+x+L−1 (L 是我们一开始枚举的后缀的左端点),枚举要不要跳过 a+1/b+1 ,通过算 lcp 转移给 gx+1,y+z (−1≤z≤1) 即可。这样状态就是 O(k2) ,单次转移在用 SA 预处理 lcp 后是 O(1) 的。总复杂度 O(nk2+nlogn) 。
30. qoj967 Rectangle Painting
类似于 loj3655 。复习一下这个东西:通过 ODT 来维护这个区间染色,观察每行的 0 连续段,每次会删掉一些加上一些,但均摊是 O(1) 的。然后你考虑对于每个段,把它拆到线段树上,就是在线段树每个节点维护一个 set 表示有哪些行的连续段会拆到这个节点上。我们可以设这个节点的权值是 set 中元素的 min ,然后再维护一个值是其子树内所有叶子到这个节点的权值 min 中最大的一个,这样查询就容易做了。
31. P9732 [CEOI2023] Trade
先想想怎么算答案,我们设一个区间的权值是前 k 大 a 的和减掉 b 的和,我们要找最大的 f(l,r) ,不难发现 f 是满足四边形不等式的,那我们分治一下就算出来答案了。
接下来看怎么判一个位置在不在最优方案内,只需要找出所有 f(l,r)=ans 的 [l,r] ,算出区间内第 k 大值 z ,那 [l,r] 内那些 ai≥z 的位置就都合法了,问题变成找 [l,r] 。但你发现这样的 [l,r] 数量可以卡到 O(n2) ,怎么办呢。
神秘观察:考虑 f(a,d)=f(b,c)=ans ,其中 a<b<c<d 。那由于 f(a,c)+f(b,d)≥f(a,d)+f(b,c) ,f(a,c),f(b,d)≤ans ,可以发现 f(a,c)=f(b,d)=ans 。而 [a,c] 的第 k 大值肯定 ≥ [a,d] 的第 k 大值,[b,d] 同理,于是可以发现我们不需要用 [a,d] 更新答案。
考虑分治的时候在值最大的前提下取 l 最大的,设其为 r。可以发现对于 r ,设上一个能得到的最大值的是 r′ ,那你只需要扫一遍 l∈[pr′,pr] 然后更新答案即可。复杂度 O(nlog2n) 。
32.qoj2562 Fake Plastic Trees 2
暴力 dp 就是设 dpx,i,j 是以 x 为根,子树内割了 i 条边,当前联通块权值和为 j 是否可能。
考虑优化,我们真的需要记下所有可能的 j 吗?想一想,我们只需要每个联通块在 [L,R] 内就行了,那有没有方法能把一个 j 支配掉呢?
考虑三个数 x<y<z 。你发现如果 z−x≤R−L+1 ,那 y 就能支配掉了:设选择了 y 之后最后形成的大小是 a ,那如果 a−y+z≤R 就替换成 z ,否则的话把 a−y+z≥R+1 和 x−z≥L−R−1 两个式子加起来就能得到 a+x−y≥L 了,所以就能把 y 替换成 x 。
再分析一下 j 的值域,由于已经割出的联通块大小总和在 [iL,iR],可以得到可能的 j 的极差是 ≤i(R−L) 的。
于是可以得到可能的 j 只有 O(K) 种,总复杂度 O(nK3)。
33.qoj5016 Range Minimum Element
不妨假设每个位置都被至少一个区间包含,设 hi 是确定了 b 后包含 i 的区间对应的 b 的 max 。考虑知道 h 后能不能确定 b 的权值,发现是能的:找到 h 最小的所有位置,设其为 i1,i2,…,ik ,那包含其中任意一个位置的区间的 b 就等于这个最小值了。剩下的区间权值一定大于最小值,递归的填 b 即可。当然这组 h 还可能不合法, 我们要求递归到的每个 [l,r] 都满足:使用被 [l,r] 包含的区间能覆盖完 [l,r] 。
现在考虑 dp ,设 fl,r,k 是有多少组合法的 h ,其最大值 ≤k 。从小到大枚举 k ,再枚举 l ,转移是容易的,复杂度 O(n3m) ,插一下就是 O(n4) 。
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列1:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
· 【杂谈】分布式事务——高大上的无用知识?