CF917E Upside Down

我永远喜欢数据结构。

大家好,我喜欢暴力数据结构,所以使用分块通过了此题。

3400 是你的谎言。我觉得这题远大于 CF1608G 啊。


洛谷题目传送门

CodeForces 题目传送门

  • 给出 n 个点的树,第 i 条边上有字母 ci。有 m 个字符串 s1sm 以及 q 组询问。每次询问给出 x,y,k

  • str(x,y)x,y 简单有向路径边上的字母按顺序拼接得到的字符串,形式化地,若 x,y 简单有向路径上一共有 E 条边,记 eix,y 有向路径上的第 i 条边,则 str(x,y)=ce1ce2ceE

  • skstr(x,y) 中出现了多少次。形式化地,求有多少个正整数 i[1,|str(x,y)||sk|+1] 使得 j[0,|sk|1]Z,(sk)i+j=[str(x,y)]i+j

  • M=i=1m|si|n,m,M,q105

  • 3 s / 512 MB

约定:

  • 本文中所有下标均从 1 开始。钦定 1 为根。用打印机字体(\texttt)表示具体的字符 / 字符串。

  • 默认 O(n)=O(m)=O(M)=O(q)O(n)>O(log2n)

  • 记一个点 u 的父亲为 fau,深度(到根的边数)为 depu

  • xy 表示 xy 的简单有向路径,ufau 这条边上的字符为 valu。特别地,val1=1

  • lca(x,y) 表示 x,y 两点的最近公共祖先。lcp(x,y) 表示两个字符串 x,y 的最长公共前缀。

  • 记一个串 s 的反串为 sR。形式化地,|sR|=|s|siR=s|s|i+1

  • anc(k,x) 表示 x 向上走 k 条边到达的点,即 x 的树上 k 级祖先。

  • 若字符参与运算,则其值等于其 ASCII 值。


考虑弱化版 CF1045J 的做法,sk 出现的位置要么完全包含在 xlca(x,y)lca(x,y)y 两条直链内,要么跨过了 lca(x,y)

Case 1:完全包含在直链内的情况

在这部分中考虑使用哈希实现字符串匹配。我们的哈希方式为多项式哈希,即对于字符串 s,其哈希值为:

H(s)=i=1|s|(sibase|s|i)modMOD

其中 base 为乘数,MOD 为模数。本题卡乘数和模数,我使用了【】生日的日期做乘数,109+9 做模数。不能自然溢出,因为后面需要用到逆元。

在弱化版中,我们运用 |S|100 的条件,对于每一种长度的字符串单独处理。在此题中我们也可以如法炮制,需要运用到一个引理:

Lemma 1

在字符串总长度为 n 的长为 m 的字符串序列 s1sm 中,本质不同的字符串长度种数为 O(n) 级别。

Proof 1

考虑 s1sm 种出现了 k 种本质不同的长度,从小到大依次是 len1lenk,记 cnti 表示第 i 种长度的出现次数,形式化地,cnti=j=1m[|sj|=leni]

那么有:i=1k(lenicnti)=n

可以发现 lenii,cnti1。后面那个很显然,因为这种长度出现时一定存在一个字符串满足其长度为 leni

至于前面那个使用归纳法证明:

  • i=1len11 显然成立。

  • 假设对于 i[1,p]Z 时成立,则对于 i[1,p+1]Z 时,由于 len1lenk 中的每一个数都是一种本质不同的长度,且从小到大排列,所以 lenp+1>lenpp。由于都是整数,所以 lenp+1p+1

由于这里涉及到的量都是正的,所以 lenicntii,因此 i=1kii=1k(lenicnti)=n,因此有 k2+k2n。可以得到 k2n=O(n)。注意这里不是在解不等式,由 k2+k2n 推导出一个成立的条件。

证毕。

那么我们可以对于这 O(n) 种长度分别求解。在求解一种长度 len 的询问时,我们对于每个点 u 预处理 str(u,anc(len,u))str(anc(len,u),u) 的哈希值 ukudku(若不存在则不处理)。需要先预处理 upudwnu 表示 str(u,1)str(1,u) 的哈希值。容易得到:

  • upuupfau+valubasedepu1(modMOD)

  • dwnudwnfaubase+valu(modMOD)

由于这个做法比较垃圾,我们不能在求解每种长度时重新遍历树计算哈希值,否则会超时。可以考虑牺牲空间,在第一次遍历树时就对于每个点存下这些哈希值。这样可以省去 O(n) 次遍历树的时间。

ukudku 都可以通过与 anc(len,u) 的哈希值差分得到,注意差分时的移位操作。

具体地:

  • ukuupuupanc(len,u)basedepulen(modMOD)

  • dkudwnudwnanc(len,u)baselen(modMOD)

可以预处理 base 的若干次方以及对应的逆元。考虑怎么快速求 anc(len,u)。由于查询次数为 O(nn),所以单次查询必须为 O(1)。可以考虑维护 1u 构成的序列 stk,使得 stki 为路径上的第 i 个点。则所求即为 stkdepulen。每次遍历到一个点 u 时,将 u 加入序列末尾。结束 u 的遍历时,将 u 从序列末尾删除。

在求解每种长度时再考虑对于每种询问串的哈希值 hsh 单独求解。考虑记 num0,u 表示 1u 上有多少个点 v 满足 ukv=hshnum1,u 表示 1u 上有多少个点 v 满足 dkv=hsh。这个可以考虑从 1n 扫描 i,依次维护 v[1,i] 的情况,则每次新扫到一个位置,需要修改(+1)的值拍平成 dfn 序后形如一段区间。

至于询问,对于 xlca(x,y) 那条链上的贡献,考虑匹配的起点在哪个位置,容易发现链上的一个点 u 能够匹配当前询问串当且仅当 depudeplca(x,y)len,且 uku=hsh。因为这样才能包含在直链内。进一步发现满足这个条件的点位于 xanc(depxlendeplca(x,y),x) 上,那么拿 num0,xnum0,anc(depxlendeplca(x,y),x) 差分即可。lca(x,y)y 的查询方式类似,注意此时匹配的方向是自上而下,用 num1 值差分计算。

此时,一共有 O(nn) 次修改,O(n) 次查询,维护以 dfn 序为下标的差分数组,那么只需要分块支持 O(1) 单点修改,O(n) 前缀查询即可。

值得注意的是,为了将同种哈希值的询问一起做,考虑使用排序将它们排在一个连续的区间内时,需要使用基数排序确保排序复杂度线性,才能保证 O(n) 次排序的总复杂度为 O(nn)


Case 2:跨过直链的情况

考虑分别处理每种串 sk 的询问。

假设跨过直链的匹配发生在 uv 上,其中 u,vxy 上的节点且 uv 前。此时一定满足,str(u,lca(x,y))sk 的前缀,str(lca(x,y),v)sk 的后缀。

同时,str(u,lca(x,y))str(x,lca(x,y)) 的后缀,str(lca(x,y),v)str(lca(x,y),y) 的前缀。

考虑找到最长的长度 P,Q,使得 str(x,lca(x,y)) 存在长度为 P 的后缀为 sk 的前缀;str(lca(x,y),y) 存在长度为 Q 的前缀为 sk 的后缀。

考虑一个基础问题:

给出字符串 a,b,找到 b 的最长前缀使得它是 a 的后缀。求出这个最长长度。

解决方法是:找到 a 的一个后缀 a[i,|a|] 使得 |lcp(a[i,|a|],b)| 最大。记 L=|lcp(a[i,|a|],b)|,则 a[i,|a|] 最长的长度不超过 Lborder 的长度即为所求。

接下来证明正确性。

Proof 2

首先,这个 border 一定是同时是 b 的前缀和 a 的后缀。因为它是 a[i,|a|] 的前缀又是它的 border,说明 a[i,|a|] 存在这个 border 作为后缀。自然 a 也存在这个 border 作为后缀。记这里求出来的长度为 len

考虑是否存在更长的答案。假设存在更长的答案长度为 tmp,其一定不超过 L,不然 a[i,|a|] 就不是使得 lcp 长度更大的后缀了。此时,a[i,|a|]b 存在长度为 tmplcp。这时候 a[i,|a|] 开头的 tmp 个字符形成的字符串与结尾的 tmp 个字符形成的字符串相等。此时 tmpa[i,|a|] 的一个更长的、长度不超过 Lborder,矛盾。

因此不存在更长的答案,len 即为所求。

证毕。

将原问题转化成上述形式,那么 P 就是 str(lca(x,y),x) 最长的前缀长度满足它是 skR 的后缀。这个和原问题显然是等价的,因为两个询问串都是原问题询问串的反串。不妨令新问题答案为 w,则反过来后原串对应的位置也相等,原问题的答案至少w;若原问题存在更长的答案 z,则反串中这些对应的位置也相等,z 就是新问题的一个更长的答案,与 w 的定义矛盾。

因此,P 就是在这个基础问题中 b=str(lca(x,y),x),a=skR 的情况;Q 就是在这个基础问题中 b=str(lca(x,y),y),a=sk 的情况。

先对 sk 以及 skR 进行后缀排序。

最长的长度不超过 Lborder 很好求,由于要求的是某个后缀的 border,在其反串上就变成了前缀的 border,这两个问题也是等价的,和上面的证明类似。

因此,对于 skskR 建立失配树 TkTkR,并进行轻重链剖分。找到满足条件的后缀后,先看一下它的反串对应的是哪一棵失配树,然后在失配树上一条一条重链向上跳。失配树的根链是单调递减的(从节点到根)。若链顶大于 L,就整条链跳过,否则在链上二分,单次询问时间复杂度为 O(logn)

接着考虑如何找到使得 L 最大的后缀。这个后缀一定满足,它要么是 a 中最大的字典序大小不超过 b 的后缀,要么是 a 中最小的字典序大小个超过 b 的后缀。换句话说,设这两个后缀与 blcp 长度分别为 A,B,则 L=max{A,B}

接下来给出证明:

Proof 3

设它们的排名分别为 i,j。则一定有 j=i+1。因为根据定义,排名为 i+1 的后缀字典序大小已经超过了 b,但是排名在 [1,i]Z 内的后缀字典序大小都不超过 b

考虑反证,假设排名为 rnk(rnkirnki+1) 的后缀会得到更大的 lcp 长度。

记这个更大的 lcp 长度为 len。分两种情况讨论:

  • rnk[1,i)Z,则排名为 rnk 的后缀的前 len 位均与 b 的前 A 位相同。根据 len 的定义可知其第 A+1 位也与 b 的这一位相同。根据定义,排名为 i 的后缀的第 A+1 位小于 b 的这一位,或者说这一位不存在(空字符)。此时,排名为 i,rk 的两个后缀前 A 位相同都等于 b 的前 A 位。且后者的第 A+1 位大于前者的这一位。说明后者比前者字典序大,这与 rnk[1,i)Z 矛盾。

  • rnk(i+1,|a|]Z,与上一种情况类似推导得到字典序大小上的矛盾即可证明。

证毕。

于是考虑求得排名 i。由于经过后缀排序,即这些后缀的字典序递增,所以答案有单调性,直接二分这个排名即可。

考虑如何求一条链上的字符串和序列上的字符串的最长公共前缀长度。对原树进行轻重链剖分,将边权转化为深度较深的端点的点权,则这条链会被表示成 O(logn) 条连续的重链区间。对于 dfn 序列形成的字符串维护哈希,对 skskR 也维护哈希。

一条一条重链匹配,若能全部匹配上,就算上这些长度,否则二分第一个不同的位置。只有第一条不匹配的重链需要二分,因此时间复杂度为 O(logn)

这部分细节比较多,尤其是一方匹配完的边界情况,具体看代码中的 qlcp 部分。

此时,这个过程已经求出了 lcp 长度,顺带比较一下大小配合套在外面的二分。

那么 P,Q 均被我们求出来了。

求出来之后,我们只需要考虑 xlca(x,y) 的后 P 个位置为开头处形成的匹配,因为不能匹配 sk 更长的前缀了。

记这 P 个位置依次为 u1uP,满足它们按照 lca(x,y)x 路径上的顺序排列;记后 Q 个位置依次为 v1vQ,满足它们按照 lca(x,y)y 路径上的顺序排列。

ui 为开头处能形成合法的匹配,当且仅当一下三点同时满足:

  • sk[1,P] 存在长度为 iborder
  • skR[1,Q] 存在长度为 |sk|iborder
  • i|sk|

证明:

Proof 4

  • 充分性:

    因为 i|sk|,所以跨过了 lca(x,y)。因为 sk[1,P] 存在长度为 iborder,根据 P 的定义可以得到 str(ui,lca(x,y))=sk[Pi+1,P]=sk[1,i];因为 skR[1,Q] 存在长度为 |sk|iborder,类似地,str(v|sk|i,lca(x,y))=skR[1,|sk|i],根据反串的定义得到 str(lca(x,y),v|sk|1)=sk[i+1,|sk|]。两者拼接恰好是 sk

  • 必要性:

    ui 开头处可以形成合法的匹配,首先一定有 i|sk|,因为要跨过 lca(x,y)。其次 str(ui,lca(x,y))=sk[1,i],根据 P 的定义,str(ui,lca(x,y))=sk[Pi+1,P],因此 sk[1,i]=sk[Pi+1,P],即 sk[1,P] 存在长度为 iborder;类似地,str(v|sk|i,lca(x,y))=skR[1,|sk|i]=skR[Q|sk|+i+1,Q],因此 skR[1,Q] 存在长度为 |sk|iborder

证毕。

所以,我们要统计有多少 i 合法,就是要统计有多少 i 满足这三个条件。

转化成失配树上的限制,就是要求有多少 i 满足 iTkP 的根链上,且 |sk|iTkRQ 的根链上。

考虑离线 + 扫描线。对于所有 sk 的询问,将它挂在 TkP 节点上。考虑深度优先搜索 Tk,在过程中一并维护数组 a0a|sk|。其中 aj 表示有多少 i,满足:

  • iTk 的当前搜到的点 u 的根链上。
  • |sk|iTkR 中点 j 的根链上。
  • i|sk|

则只要在 P 处单点查 aQ 即可。

每次新扫到一个点 u,则和上一层深搜相比根链上增加了一个点 u 的贡献。考虑 |sk|u 的贡献,此时首先满足 u|sk|,发现只有它子树内的点的根链经过它,即只要这些点的 aj 值要增加 u 的贡献。拍平成 dfn 序后将 a 映射过去,再进行差分,则需要支持的操作形如区间加、单点查,由于此处修改、询问同阶,树状数组维护即可。每次结束一个点的深搜时,删去它的贡献。

这部分就做完了,时间复杂度为 O(nlog2n)


综上,这个做法时空复杂度均为 O(nn),可以接受。前面也说过空间可以做到线性,只是需要一些精湛的卡常技艺。

AC Link & Code

码量 12.5 KB。

总结一下,本题考察了一下几方面的知识:

  • 思想:离线、扫描线、平衡复杂度等。
  • 树论:轻重链剖分、最近公共祖先等。
  • 数据结构:树状数组、分块等。
  • 字符串:后缀数组、失配树、字符串哈希等。
  • 数论:快速幂、乘法逆元等。

题目思路巧妙、性质极多、码量极大,对于选手的思维能力和代码能力都是一种不小的挑战,算是一道不可多得的字符串神仙 & 毒瘤题目。

完结撒花!

posted @   lzyqwq  阅读(26)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· 没有Manus邀请码?试试免邀请码的MGX或者开源的OpenManus吧
· 【自荐】一款简洁、开源的在线白板工具 Drawnix
· 园子的第一款AI主题卫衣上架——"HELLO! HOW CAN I ASSIST YOU TODAY
· Docker 太简单,K8s 太复杂?w7panel 让容器管理更轻松!
点击右上角即可分享
微信分享提示