By zghtyarecrenj
本文包括:Lyndon Words & Significant Suffixes & Lyndon Array & Runs & Lyndon Tree & Three Squares Lemma。
禁止转载全文,转载部分需要注明出处。
如果你发现笔者有写错的地方,请联系笔者。(尽量不要用博客园评论,我不会经常看,建议用洛谷私信我或者加 QQ 3360155688)。
截至 2021-06-14 的 To-do List:
目录
关于字符串部分的基础详见 oi-wiki 字符串,这里略提。
我们定义两个字符串 \(a\) 和 \(b\),如果 \(a\) 的字典序 \(<b\),则我们称 \(a < b\)。
如果 \(a\) 是 \(b\) 的前缀且 \(a \ne b\),则我们称 \(a \sqsubset b\)。
如果 \(a\) 是 \(b\) 的前缀,则我们称 \(a \sqsubseteq b\)。
如果 \(a < b\) 且 \(a\) 不是 \(b\) 的前缀,则我们称 \(a \triangleleft b\)。即 \(a \triangleleft b \Longleftrightarrow (a < b) \wedge (a \not\sqsubseteq b)\)。Fact:如果 \(a \triangleleft b\),则 \({au} < {bv}\)。
\({abc}\) 表示拼接 \(a, b, c\) 三个字符串。
\(a^n\) 表示 \(n\) 个 \(a\) 拼接在一起。e.g. \({a^2b} = {aab}\)
\(\epsilon\) 表示空串。
我们定义字符集为 \(\Sigma\),组成的字符串为 \(\Sigma^*\),\(\Sigma^+ = \Sigma^* \setminus \{\epsilon\}\)
\(\operatorname{pref}(a)\) 表示所有 \(a\) 的前缀的集合,\(\operatorname{suf}(a)\) 表示所有 \(a\) 的后缀的集合(包含 \(a\) 和 \(\epsilon\))
\(\operatorname{pref}^+(a) = \operatorname{pref}(a) \setminus \{a,\epsilon\},\ \operatorname{suf}^+(a) = \operatorname{suf}(a) \setminus \{a, \epsilon\}\)
若无特殊定义,字符串 \(s\) 是从 \(0\) 开始。
\(|s|\) 表示 \(s\) 的长度,\(s[i..j]\) 表示 \(s\) 的一个子串,第一个字符的标号为 \(i\),最后一个字符的标号为 \(j\)。
$\hat{w}=w$ $,其中 $$ $ 是一个比字符集里面任何数小的字符。
这部分这位老哥写的可能比我的详细= = link
Lyndon Word:一个字符串是一个 Lyndon Word 当且仅当 \(\forall a\) 的后缀 \(b\),有 \(a < b\)。
还有一个定义:对于一个 \(n\) 的串有 \(n\) 个循环同构,则其中严格最小的那个是一个 Lyndon Word。
比如:\(\text{ab}\) 是一个 Lyndon Word,但是 \(\text{ba}\) 不是。
\(\mathcal L\) 表示所有 Lyndon Word 组成的集合。
又称 Lyndon Decomposition。
我们定义 \(\operatorname{CFL}(s)\) 是一个对于 \(s\) 串的划分,即划分成了 \({w_1w_2\cdots w_k} = s\),使得所有 \(w_i\) 是 Lyndon Word,并且 \(w_1 \ge w_2 \ge \cdots \ge w_n\)。
比如:串 \(\text{bbababaabaaabaaaab}\) 的 Lyndon 分解是 \(\text{b|b|ab|ab|aab|aaab|aaaab}\)。
这是一个很显然的结论。
如果 \(a, b \in \mathcal L\),且 \(a < b\),则 $ {ab} \in \mathcal L$。
由于 \(a < b\),我们有 \({ab} < b\)。接下来我们分两种情况讨论。
当 \(a \not \sqsubseteq b\) 时:根据 \(a < b\),我们有 \(a \triangleleft b\)。所以 \({ab} \triangleleft b \implies {ab} < b\)。
当 \(a \sqsubseteq b\) 时:令 \(b={ac}\),则 \({ab} = {a^2c}\)。因为 \(b \in \mathcal L\),所以 \({ab} < b \implies {a^2c} < {ac} \implies {ac} < c\),所以 \(b < c\)。
所以,\(\forall d \in \operatorname{suf}^+(b), \ {ab} < b < d \implies \forall c \in \operatorname{suf}^+(a),\ a \triangleleft e \implies {ab} \triangleleft {eb}\)。\(\blacksquare\)
这个结论和 [Theorem 1.2.3] Uniqueness of CFL 是两个很有趣的结论。
对于任意的串 \(s\),\(\operatorname{CFL}(s)\) 一定存在。
构造法。我们考虑,单个的字母一定是 Lyndon Word。
根据 [Theorem 1.2.1 Lyndon Concatanation],我们可以把字典序小的两个 Lyndon Word 并起来,所以我们把所有的字典序单增的序列都并起来,剩下的就是一个合法的 CFL。\(\blacksquare\)
对于任意的串 \(s\),\(\operatorname{CFL}(s)\) 一定唯一。
反证法,假设有两种方案。我们考虑第一个不同的位置的情况,可以很容易地得到矛盾,和 CFL 的定义矛盾。\(\blacksquare\)
然后我们就得到了 CFL 存在且唯一。由此有两个推论:
\(w_1\) 是最长的 Lyndon 前缀且 \(w_k\) 是最长的 Lyndon 后缀。
反证法。因为如果 \(w_1\) 不是最长,那么还能再拼,产生了两个合法的 CFL,和 [Theorem 1.2.3 Uniqueness of CFL] 矛盾。所以 \(w_1\) 是最长的 Lyndon 前缀。
\(w_k\) 同理。\(\blacksquare\)
一个字符串 \(s\) 的最小后缀是 \(w_k\)。
首先,我们有这样的一个 CFL:
首先,我们记 \(w_n\) 的起始位置为 \(pos\),则显然
如图,最小后缀的其实位置不可能 \(>pos\),因为根据 Lyndon Word 的定义,\(w_n\) 的每个后缀都大于他自身。
接下来我们考虑最小后缀在另一个位置的情况,即他在另一个 \(w_i\) 之中
根据 \(w_i \ge w_{i+1}\ge \cdots \ge w_n\),而 \(w_i\) 的一个后缀 \(> w_i\),所以这个后缀大于 \(w_n\)。
所以唯一可能的最小后缀就是 \(w_n\)。
简单来说,假设最小后缀是 \({xw_{i+1}w_{i+2}\cdots w_{k}}\) 而不是 \(w_k\) 且 \(|x| < |w_i|\)。我们有 \({x w_{i + 1} \dots w_k} \geq x > w_i \ge w_k\),矛盾。\(\blacksquare\)
就是求出 CFL 的算法啦~
我们有一个非常优美的算法
简要思想:
$ {uav} \in \mathcal{L}\(,\)u,v,h \in \Sigma^*\(,\)a<b\in\Sigma\(,\)k\ge1$
换成代码实现就是:
我们需要维护两个部分:\(ubv\) 和 \(u\)。
如果还是不太懂可以移步 oi-wiki,那里写的比较详细。
模拟即可,显然空间复杂度 \(\mathcal O(1)\)。接下来证明一下时间复杂度。
接下来证明一下复杂度为什么是对的。
最优情况为一个分解走到底,\(\mathcal O(n)\)。
最坏情况为不停地在重新找,由于至多回退 \(n\) 次,每次回退的距离不超过前进的距离,所以是 \(\mathcal O(n)\)。
CFL 板子,欢迎来蒯。char *s
从 \(0\) 开始编号,int *ans
返回 CFL 的右端点(注意从 \(0\) 开始),函数返回的 int
是 Lyndon 串的个数。
inline int CFL(char *s, int *ans) {
int res = 0, n = strlen(s), i = 0, j, k;
while (i < n) {
j = i, k = i + 1;
while (k < n && s[j] <= s[k]) j = s[j] == s[k++] ? j + 1 : i;
while (i <= j) ans[res++] = i + k - j - 1, i += k - j;
}
return res;
}
给定串 \(S\),求前缀的最小后缀。要求 \(\mathcal O(n)\)。
这是很裸的 Lyndon 的题吧……
不用 Lyndon 的话传统做法是后缀树,维护一下最右边的路径即可。然而是带 \(\log\) 的。
根据 [Theorem 1.2.5 Theorem of Minsuf],\(w_n\) 就是最小后缀,所以我们可以记忆化一下 CFL 分解:
我们知道,Lyndon 分解的时候是一个 Lyndon 串不停地在重复,如果碰到冲突会重新跳,我们考虑把这个东西记忆化,记为 \(\operatorname{minsuf}(i)\)。
然后这个问题就迎刃而解了。
求前缀的最大后缀。
你可能看到这个问题会觉得把字典序反一下就好了,那也太 naive 了,试试自己举一个反例。
揭晓答案:\(\color{white}\text{我们在最小后缀之中需要的是字典序最短的,但是最大后缀要字典序最长的}\)。(试试复制一下,你就能看了)
但是好消息是这个问题可以用另一种非常 naive 的方法解决,把 Lyndon 串在分解的时候直接记一下,输出当前的串就好了。(感性理解一下
最小循环表示
对着 \({x}\) 求一遍 CFL,找到前半串开始的最长 Lyndon Word 就是了。时间复杂度 \(\mathcal O(n)\)。
FAQ:诶,这用最小表示不是更好做吗?A:但是你用 Lyndon 不就能少学一个算法吗= =
代码可以看 xht 的。
给你一个串 \(S\ (|S|\le 5\times 10^5)\)。给你 \(q\) 个询问,每次询问 \(\operatorname{maxsuf}(s[l..r])\) 的下标。
这个显然可以SA,但是这个问题其实是可以使用 Lyndon 的知识解决的。
首先我们以 \(r\) 为关键词离线所有的询问,所以现在这个东西变成这样了:
我们定义一个位置 \(i\) 是有用的当且仅当 \(f(i\cdots r)\) 是 Lyndon 的。
我们有一堆 \(i\),他们有这些特点:他们满足一种后缀关系(废话
我们可以维护一个类似栈的东西,把相邻的两个点 LCP 算出来,比如:
C 和 B 在某个位置不同,B 和 A 在某个位置不同。我们把小的删掉,一直删光。
如果我们要回答询问,我们可以 std::lower_bound
求出 \(l\) 右边的第一个。
如果你不想离线,我们把一个 \(r\) 对应的集合用主席树记下来,也可以做的。
时间复杂度 \(\mathcal O(n \log n)\)。
题意:给你一个字符串 \(S\),你可以选择任意一个区间进行翻转,使操作后的字符串字典序最小。\(|S| \le 10^7\)
来观察性质。
Lemma:如果我们在开头有一个不大于其他所有的字符的字符,则我们一定不动这个字符,否则一定会反转某一个前缀。
假设翻转一个前缀 \(s[1..t]\) 后的第一个字符比当前的小,则显然 \(s[t] < s[1]\),矛盾。
如果第一个字符不是最小,那么我们应该找到最小的字符 \(s[t]\) 并且翻转前缀 \(s[1..t]\),使得第一个字符成为最小。\(\blacksquare\)
既然我们要考虑一个翻转的区间,则我们不妨把原创翻转的结果记为 \(S^R\)。
我们考虑翻转过后的一个前缀是 \(S^R\) 中的一个后缀,而我们要求翻转之后的字典序最小,所以我们需要挑选 \(S^R\) 中的最小后缀(这里的最小后缀和一般说的有所不同,只要比其他不存在前缀关系都要小即可,所以可能不止一个)。
可以发现,所以符合条件的后缀都是最长的一个的前缀。
我们注意到,一个串 \(S'\) 是可以表示成若干的 \(w\) 的幂和一个 \(w\) 的前缀。
接下来我们考虑什么时候需要翻转。
如果 \(|A|=|B|\),则 \(A + A + C < A + C + B < C + B + B\) 或 \(A+A+C>A+C+B>C+B+B\) 或 \(A+A+C=A+C+B=C+B+B\)。
我们考虑 \(A=w, B=w^R,C=w'\) 的情况,发现满足条件的后缀之中只有最长的和最短的后缀可能对答案产生贡献。
接下来我们应该使用 Duval 算法来找到 \(\operatorname{CFL}(S^R)\) 来找到这些后缀。
时间复杂度 \(\mathcal O(n)\)。
题意:给你一个字符串 \(S\),你可以把他分成 \(k\) 段,并且翻转其中的若干段,使操作后的字符串字典序最小。\(1\le k\le |S| \le 5\times 10^6\)
留作习题。
我们令 \(\operatorname{minsuf}(u)\) 为 \(u\) 的最小后缀,且 \(\operatorname{minsuf}(u, v) = \min _{w \in \operatorname{suf}(u)} wv\)。
Significant Suffixes:\(\Lambda(u) = \arg\min_{w\in \operatorname{suf}(u)} wv\)。
\(\operatorname{minsuf}(s)\) 表示 S 字典序最小的后缀,且
由 \(\operatorname{minsuf}\) 的性质可知,\(\operatorname{minsuf}(u, \epsilon) = \operatorname{minsuf}(u)\)。\(\implies \operatorname{minsuf}(u) \in \Lambda(u)\)。
所以,显然 \(\forall u \in \Lambda(u),\ \operatorname{minsuf}(u) \sqsubseteq u\)。
我们注意到一个 CFL 分解中的 Lyndon Words 是存在一定的循环的。因此,我们可以记一个 CFL 为次方的形式。
\[\operatorname{CFL}(u) = {{w_1}^{k_1}{w_2}^{k_2}\cdots {w_n}^{k_n}}
\]
我们记 \(s_i\) 为一个后缀,即 \(s_i = {{w_i}^{k_i}{w_{i+1}}^{k_{i+1}}\cdots{w_n}^{k_n}}\)。边界:\(s_{n+1} = \epsilon\)。
首先,我们需要一个引理。
在 [Theorem 2.2.2 Significant Suffixes Theorem] 里面会用到,建议先食用下一个 Theorem。
如果 \(u^\infty < v\),则 \(v > {uv} > {u^2v} > \cdots\)。
\(u^\infty < v \implies u^\infty < {uv}\)。
令 \(u = {xay}\),\(v = {(xay)^k xbh}\),其中 \(x,y,h\in\Sigma^*\),\(a,b\in\Sigma\),\(a<b\)。
我们有 \(v \succ uv \Longleftrightarrow (xay)^{k - 1} xbh \succ (xay)^k xbh \Longleftrightarrow xbh \succ (xay) xbh\)。
\(v>{uv} \implies {u^iv} > u^{i+1} \implies \blacksquare\)
同理如果 \(u^\infty > v\),则 \(v < {uv} < {u^2v} < \cdots\)。
\[\Lambda(u)\subseteq \{s_i | i \in [1,n]\}
\]
反证法:如果这个命题不成立,则我们分类讨论
i) 假设有一个串 \(v = {b{w_i}^ks_{i+1}} \in \Lambda(u)\),\(|b| < |w_i|,\ 0 \le k < k_i\)。
\(w_i \in \mathcal L \implies w_i \triangleleft b \implies s_i = {w_is_{i+1}} < {bs_{i+1}}\),矛盾。
ii) 假设有一个串 \(v = {{w_i}^ks_{i+1}} \in \Lambda(u)\),\(1 < k < k_i\)。
根据 [Theorem 2.2.1 Infinite Theorem],如果 \({w_i}^\infty < s_{i+1}\),则 \({{w_i}^{k_i}s_{i+1}} < {{w_i}^{k_i - 1}s_{i+1}}<\cdots
我们令 \(\lambda = \min \{i : s_{i+1} \sqsubset s_i\}\)。\(\forall i \ge \lambda, \ w_i = {s_{i+1}y_i},\ x_i = {y_is_{i+1}}\)。\(\implies s_i = {{w_i}^{k_i}s_{i+1}}= {(s_{i+1}y_i)^{k_i}s_{i+1}} = {s_{i+1}{x_i}^{k_i}}\)。
根据 CFL 的性质,\(s_{\lambda} \triangleleft w_{\lambda - 1}\)。所以 \(\Lambda(u)\subseteq \{s_i | i \in [1,n]\}\)。\(\blacksquare\)
如果有 \(2\) 个串 \(u\) 和 \(v\),满足 \(|u| \le |v|\),则我们有
\[\begin{aligned}\Lambda(uv) &\subseteq \Lambda(v) \cup \{\operatorname{maxsuf}^R(u, v)\} \\&= \Lambda(u) \cup{\max _{s \in \Lambda(u)}}^R \{sv\}\end{aligned}
\]
理由很简单,因为 \(\{\operatorname{maxsuf}^R(u, v)\}\) 也是一个 Significant Suffix,随意我们就可以把它展成第二行的式子的形式。\(\blacksquare\)
一个字符串 \(S\) 的 Significant Suffixes 至多有 \(\log n\) 个。
原命题可以很容易地转化为:
如果两个 Significant Suffixes \(u\),\(v\) 满足 \(|u| < |v|\),那么 \(2|u| < |v|\)。
反证法。设存在 \(|u| < |v| < 2|u|\)。因为 \(u, v \in \operatorname{suf}^+(u)\),所以 \(u \in \operatorname{suf}^+(v)\)。
所以我们可以非常容易地知道,\(u \triangleleft v\)。\(\implies v\) 有一个长度为 \(|v| - |u| < \frac {|v|} 2\) 的周期,记为 \(T\)。
所以,\(u = {Tw}, v = {T^2w}\)。
由于 \(u\) 是一个 Significant Suffix,因此存在串 \(t\),满足 \(vt>ut\),即 \({T^2wt} > {Twt} \implies {Twt} > {wt}\)。
而 \(w \in \operatorname{suf}^+(s)\),所以与 \(u\) 是 Significant Suffix 矛盾。\(\blacksquare\)
我们知道 \(\Lambda(S)\) 中有很多串,其中最短的是 \(\operatorname{minsuf}(S)\),而最长的是 \(\operatorname{maxsuf}^R(S)\)。这里的 \(^R\) 代表 reverse。
这里的 Proof 先咕着吧。贴个 Reference:
Tomohiro, I., Nakashima, Y., Inenaga, S., Bannai, H., & Takeda, M. (2016). Faster Lyndon factorization algorithms for SLP and LZ78 compressed text. Theor. Comput. Sci., 656, 215-224.
Kociumaka, T. (2016). Minimal Suffix and Rotation of a Substring in Optimal Time. ArXiv, abs/1601.08051
留作习题。
简要题意:维护动态字符串,值域 \(10^9\),区间加,求区间 \(\operatorname{minsuf}\) 的位置。\(|S| \le 2 \times 10 ^ 5,m \le 3 \times 10 ^ 4\)
前置芝士:Lyndon 分解,Significant Suffixes,线段树,字符串哈希,分块。
我们先考虑不带修的情况。
由 [Theorem 2.3.1 Lambda Subset Theorem],我们可以很容易地想到考虑建一棵线段树来维护 Significant Suffixes。
细节:如果线段树的 mid = l + r >> 1
,则左边的区间比右边长一些。但是上面的这个结论对于 \(|u| \le |v|\) 有效,所以我们需要调整一下,使得左儿子比右儿子要长一些(即:mid = l + r + 1 >> 1
,使得左儿子总不比右儿子短)
可以存一下当前代表的串的所有 Significant Suffixes,然后直接考虑合并(把右边的所有的直接加进来,左边的都循环一遍,字典序最长的加进去)得到父节点的 Significant Suffixes 即可。(看不懂的看代码)
由 [Theorem 2.3.2 Lambda Log Theorem] 可知,每一个集合都是 \(O(\log n)\) 大小的。这样的话,我们求出了每一个线段树上的区间的 Significant Suffixes。然后查询就在这 \(O(\log n)\) 个区间内求 Significant Suffixes 的并,暴力比较即可。所以我们需要一个 \(O(1)\) 比较两个串的方法(否则复杂度就挂了)。所以如果不带修的话我们可以考虑 SA。
接下来考虑带修的情况。
我们需要快速地求两个串的 LCP,又有一个线段树,所以可以很自然地想到一个线段树+字符串哈希+二分LCP的算法。复杂度 \(O(q \log^4 n)\),慢了点,我这种人傻常数大的就不用想了。
我们考虑分块维护一些哈希,分 \(\sqrt n\) 的块。我们维护一下每个点到块的末端的哈希值,然后维护一下每个块到串的末尾的哈希值。然后我们可以记一个块的全局的偏移量,就可以算了。每次查询的时候,我们只需要查 \(2\) 次即可,\(O(1)\) 查找。最终是 \(O(q \log ^3 n + q\sqrt n)\) 的复杂度。
Q: 道理我都懂,但是为什么我挂了? A: 你是用了自然溢出哈希吧,换个双哈希试试。
我们有一个字符串 \(s\),则
Lyndon Array:\(\mathcal L[i] = \max \{j : s_i \cdots s_{j-1} \in L\}\),其中 \(L\) 表示 Lyndon 串的集合。在 \(\prec_l\) 意义下的 \(\mathcal L\) 记为 \(\mathcal L_l\)。
这有啥子用?别急,先看性质。
最长的 Lyndon 子串是无交集的,即 \(i < j < \mathcal L[i]\),我们有 \(\mathcal L[j] \le \mathcal L[i]\)。
我们假设存在 \(i,j\) 使得 \(\mathcal L[i] < \mathcal L[j]\)。
我们假设 \(u = s_i \cdots s_{j-1}\),\(v = s_j \cdots s_{\mathcal L[i] - 1}\),\(w = s_{\mathcal L[i]} \cdots s_{\mathcal L[j] - 1}\),且 \(u,v,w\) 满足 \({uv}, {vw} \in L\)。
对于所有的 \(s \in \operatorname{suf}^+({uvw})\),且满足 \(|s| \le |v| + |w|\),有 \(s \triangleleft {vw} \sqsupseteq v \triangleleft {uv}\)。\(\implies {uvw} \triangleleft s\)。
对于所有的 \(s \in \operatorname{suf}^+({uvw})\),且满足 \(|s| > |v| + |w|\),有 \({svw} \sqsupseteq {sv} \triangleleft {uv}\)。\(\implies {uvw} \triangleleft {svw}\)。
所以 \({uvw} \in L\),矛盾。\(\blacksquare\)
我们设 \(\operatorname{suf}(i) = s_i \cdots s_{n-1}\),即一个后缀。
而我们有 \(s_i \cdots s_{\mathcal L[i] - 1} \triangleleft s_j \cdots s_{\mathcal L[i] - 1}\)。
\(\implies \operatorname{suf}(i) \triangleleft \operatorname{suf}(j)\quad(i<j<\mathcal L[i])\)
于是我们设
\[\operatorname{NSV}(i) = \min\{\{j > i : \neg(\operatorname{suf}(i) \triangleleft \operatorname{suf}(j))\} \cup \{n\}\}
\]
显然 \(\mathcal L[i] \le \operatorname {NSV}(i)\)。
我们还有 \(\neg(\operatorname{suf}(i) \triangleleft \operatorname{suf}(j)) \Longleftrightarrow \operatorname{suf}(j) \sqsubseteq \operatorname{suf}(i) \vee \operatorname{suf}(j) \triangleleft \operatorname{suf}(i) \Longleftrightarrow \operatorname{rank}(i) > \operatorname{rank}(j)\)。
有了上述定义,证这个是不是非常简单呢 XD
\(\mathcal L[i] = \operatorname{NSV}(i)\)
原命题可以很方便地转化为 \(s_i \cdots s_{\operatorname{NSV}(i) - 1} \in L\)。
分类讨论:
如果 \(\operatorname{NSV}(i) = n\),\(s_i \cdots s_{\operatorname{NSV}(i) - 1} = \operatorname{suf}(i)\)
否则的话我们肯定有一些 \(j\) 使得 \(\operatorname{suf}(i) \triangleleft \operatorname{suf}(j)\)。
然后我们继续来讨论:
i) 如果 \(s_1 \cdot s_{\operatorname{NSV}(i) - 1} \triangleleft s_{j} \cdots s_{\operatorname{NSV}(i) - 1}\),易证。
ii) 反之,结合 \(\operatorname{suf}(i + (\operatorname{NSV}(i) - j) - 1) \triangleright \operatorname{suf}(i)>\operatorname{suf}(\operatorname{NSV}(i))\),易证矛盾。(你看不出来?明显与 \(\operatorname{suf}(i) \triangleright \operatorname{suf}(j)\) 矛盾)。\(\blacksquare\)
留作习题。
这个东西的英文名是 runs,他的中文名是顶天立地串……(好中二啊)
我们有一个串,runs 是他的一些子串,满足:
\(p = \operatorname{per}(s_i\cdots s_{j-1})\le \dfrac {j-i}2\),\(s_{i-1} \ne s_{i-1}+p\),\(s_{j-p}=s_{j}\)
更好理解的定义:
定义一个字符串 \(|S|\) 里的一个 run,指其内部一段两侧都不能扩展的周期子串,且周期至少完整出现两次。
严格地说,一个 run 是一个 三元组 \((i,j,p)\),满足 \(p\) 是 \(S[i..j]\) 的最小周期,\(j-i+1 \ge 2p\),且满足如下两个条件:
例如:\(S = \text{aababaababb}\) 之中有 7 个 runs:\(S[1..2] = \text a^2\),\(S[1..10] = (\text{aabab})^2\),\(S[2..6] = (\text{ab})^{2.5}\),\(S[4..9] = (\text{aba})^2\),\(S[6..7] = \text a^2\),\(S[7..10] = (\text{ab})^2\),\(S[10..11] = \text b^2\)。
(实际上是 LOJ #173 的题面,题目是我造的,这一段是 EtoainWu 的文字)
定义 \(Runs(w)\) 表示字符串 \(w\) 的所有 runs 的集合。
\(\rho(n)\) 表示了在一个长为 \(n\) 的字符串之中至多有多少组 runs,而 \(\sigma(n)\) 表示了在一个长为 \(n\) 的字符串之中所有 runs 的幂之和的最大值。
Lyndon Root:令 \(r=(i,j,p)\) 是一个run,则他的 Lyndon Root 是一个 \(s[i..j]\) 的长度为 \(p\) 的 Lyndon 子串。
每一个 run 都有一个 Lyndon root。
我们假设 \(\prec^0\) 表示 \(<\),而 \(\prec^1\) 表示 \(<^R\)。(此处的 \(^R\) 表示 reverse,给 \(\prec\) 标号是为了方便)
\(\prec^0\) 和 \(\prec^1\) 的对应的 Lyndon Array 是 \(\mathcal L^0\) 和 \(\mathcal L^1\).
\(\rho(n) \le 2n\)
原命题可以转化为
对于每个 runs,我们有存在 \(i\) 和 \(t\) 使得 \(s[i..\mathcal L^t[i] - 1]\) 是 Lyndon root。
我们令 \(w\) 是 Lyndon root,\(w=s[k..s-1]\)。
分类讨论:
如果 \(j=|S|\),
我们可以把 \(s[k .. |S| - 1]\) 表示成 \({w^pw'} \ (p \in N, w' \in \operatorname{pref}(w))\)。
因为 \(\operatorname{CFL}({w^pw'}) = {w^p\operatorname{CFL}(w')}\),所以 \(w\) 是从 \(k\) 开始的最长 Lyndon 前缀。
如果 \(j<|S|\),
我们可以把 \(w\) 表示成 \({uab}\),其中 \(a \ne b\)。
所以我们可以把 \(s_k \cdots s_{|S| - 1}\) 表示成为 \({(uav)^pub}\)。
我们不妨假设 \(b \prec^t a\)、
因为我们有 \(\operatorname{CFL}^t({(uav)^pubh}) = (uav)^p\operatorname{CFL}^t({ubh})\),所以 \({uav}\) 是 \(\prec^t\) 下的最长 Lyndon 前缀。\(\blacksquare\)
\(\rho(n)<n,\sigma(n)\leq3n-3\)
几乎从 WC2019 课件搬运的证明
定义 \(Beg(I)\) 表示 \(I\) 中所有区间的起始端点的集合。
对于一个串的 Lyndon Array \(\mathcal L^0[i]\) 和 \(\mathcal L^1[i]\),总有 \(\mathcal L^{l}[i] = [i..i], \mathcal L^{1-l}[i] = [i..j] (j\ne i)\),其中 \(l\in \{0,1\}\)。
令 \(k=\max\{k'\ |\hat{w}_{k'}\ne \hat{w}_i,k'>i\}\)。
由 [Theorem 1.2.1 Lyndon Concatanation] 可得:
若 \(r=(i,j,p)\) 为一个run,则对于 \(\hat{w}[j+1]\prec_l \hat{w}[j+1-p]\) 的 \(l\),\(\forall r\) 的 \(\prec_l\) 意义下的 Lyndon Root \(\hat w[i_{\lambda}..j_{\lambda}]\) 都与 \(\mathcal L^l(i_{\lambda})\)相等。
\(\because \hat{w}[j+1]\ne\hat{w}[j+1-p]\),令 \(l\in\{0,1\}\) 满足 \(\hat{w}[j+1]\prec_l\hat{w}[j+1-p]\)。
令 \(\lambda=[i_{\lambda}…j_{\lambda}]\) 为 \(r\) 的 \(\prec_l\) 意义下的一个 Lyndon Root,由 [Theorem 1.2.1 Lyndon Concatanation],\([i_{\lambda}…j_{\lambda}]=\mathcal L^l(i_{\lambda})\)。\(\blacksquare\)
对于一个run \(r=(i,j,p)\),令 \(B_r=\{\lambda=[i_{\lambda}…j_{\lambda}]|\lambda\) 为 \(r\) 的 \(\prec_l\) 意义下的一个 Lyndon Root 且 \(i_{\lambda}\ne i\}\)。即 \(B_r\) 表示所有 \(r\) 的关于 \(\prec_l\) 的 Lyndon Root 构成的集合,但要除去开头位置 \(i\) 处开始的 Lyndon Root。有 \(|Beg(B_r)|=|B_r|\geq \lfloor e_r-1\rfloor\geq 1\),其中 \(e_r\) 为 \(r\) 的指数。
两个不同的 run \(r,r'\),\(Beg(B_r)\cap Beg(B_{r'})\) 为空。
反证,假设存在 \(i\in Beg(B_r)\cap Beg(B_{r'})\),并且 \(\lambda=[i…j_{\lambda}]\in B_r\),\(\lambda'=[i…j_{\lambda'}]\in B_{r'}\)。
令 \(l\in\{0,1\}\) 满足 \(\lambda=\mathcal L^l[i]\),由于 \(\lambda\ne \lambda'\),有 \(\lambda'=\mathcal L^{1-l}[i]\)。
由 Lemma A,\(\lambda\) 和 \(\lambda'\) 中有且只有一个为 \([i..i]\)。
不妨设 \(\lambda=[i..i]\),那么 \(j_{\lambda'}>i\)。
由于 \(w[i…j_{\lambda'}]\) 为一个 Lyndon Word,有 \(w[i]\ne w[j_{\lambda'}]\)。
由 \(B_r\) 和 \(B_{r'}\) 的定义,\(r\) 和 \(r'\) 的开始位置均小于 \(i\),这意味着 \(w[i-1]=w[i]\)(由 \(r\) 的周期性),并且 \(w[i-1]=w[j_{\lambda'}]\)(由 \(r'\) 的周期性)。矛盾 \(\blacksquare\)
任意的一个 run \(r\) 可以被赋予一个两两不交的非空位置集合 \(Beg(B_r)\)。并且,由于 \(1\notin Beg(B_r)\) 对于任意的一个 \(r\) 均成立,有 \(\sum_{r\in Runs(w)}|B_r|=\sum_{r\in Runs(w)}|Beg(B_r)|\leq |w|-1\)。
考虑字符串 \(w\),由于对于任意 \(r\in Runs(w)\),有 \(|B_r|\geq1\),由 Lemma C,有 \(|Runs(w)|\leq\sum_{r\in Runs(w)}|B_r|\leq |w|-1\)。
考虑字符串 \(w\),令 \(e_r\) 表示 \(r\) 的指数。由于对于任意 \(r\in Runs(w)\),有 \(|B_r|\geq \lfloor e_r-1\rfloor>e_r-2\),由 Lemma C,有 \(\sum_{r\in Runs(w)}(e_r-2)<\sum_{r\in Runs(w)}\lfloor e_r-1\rfloor\leq\sum_{r\in Runs(w)}|B_r|\leq |w|-1\)。因为 \(|Runs(w)|\leq |w|-1\),可得 \(\sum_{r\in Runs(w)}e_r\leq3n-3\)。\(\blacksquare\)
在实现 Runs 之前,你需要会字符串哈希或者后缀数组或者其他后缀数据结构。
根据以上证明中的 Lemma B,每一个 runs 都会对应一个 Lyndon root,所以如果我们把 Lyndon Array 算出来了,就可以把每个 runs 对应的 Lyndon root 求出来。
那么我们怎么求 Lyndon Array \(\mathcal L^0\) 和 \(\mathcal L^1\)?
\[\mathcal{L}^0[i] = \mathrm{NSV}(i) = \min\{\{j > i : \neg(\mathrm{suf}(i) \triangleleft \mathrm{suf}(j))\} \cup \{n\}\}
\]
而
\[\begin{aligned}\mathcal{L}^1[i] &= \mathrm{NSV}^R(i) \\&= \min\{\{j > i : \neg(\mathrm{suf}(i) \triangleleft^R \mathrm{suf}(j))\} \cup \{n\}\} \\&= \min\{\{j > i : \mathrm{suf}(i) \triangleleft \mathrm{suf}(j) \vee \mathrm{suf}(i) \sqsupseteq \mathrm{suf}(j))\} \cup \{n\}\}\end{aligned}
\]
所以我们考虑对于字符串的每个后缀都维护他的 CFL,方法是在头上插入一个新字符,然后判断是否合法。根据 [Theorem 1.2.1 Lyndon Concatanation],如果遇到一个 Lyndon word 大于下一个的情况,合并即可。可以保证正确性。
这里有一个实现上的细节,可能会好写一点。根据 [Theorem 1.2.4 Lyndon Suffixes and Lyndon Prefixes],\(w_i\) 是 \(w_iw_{i+1}\cdots w_k\) 的最小前缀,所以比较两个 Lyndon word 的字典序相当于比较两个后缀的大小,而这个是比做一个 lcp 要简单多的。
所以至此 Lyndon Array 已经求完了,具体实现细节可以看代码。
接下来我们只要使用 Lyndon Array 扩展出 runs 就可以了,具体的做法是求出 lcp,即如果当前的 Lyndon Array \(\mathcal L[i] = (l..r)\),则我们 lcp 求出最长的 \(s[l..l+l_1-1]=s[r+1..r+l_1], s[l-l_2..l-1]=s[r-l_2..r-1]\)。
根据 runs 的定义,如果 \(l_1 + l_2 \ge 2(l - r + 1)\),那么我们就找到了一个 run \((l-l_2, r+l_1-1, r-l+1)\)。
如果我们使用 \(\mathcal O(n \log n)\) 的 SA 和 \(\mathcal O(n \log n)-\mathcal O(1)\) 的 rmq 算法,我们就可以在 \(\mathcal O(n \log n)\) 的时间复杂度之内求出所有的 runs。
如果我们使用 SAIS 和 \(\mathcal O(n) - \mathcal O(1)\) 的 rmq 算法,我们就可以 \(\mathcal O(n)\) 求出所有的 runs。
比较优秀的 \(\mathcal O(n) - \mathcal O(1)\) 的 rmq 方法:叉姐的 和 hqztrue的。
是我出的啦owo
标准划分:把一个 Lyndon Word 划分为他的字典序最小的后缀和其余的部分。两部分一定是 Lyndon Word。
Lyndon Tree:每一个节点都是一个 Lyndon Word,左儿子和右儿子对应一个标准划分。
\(\operatorname {Ltree}(S)\) 表示一个根节点为 \(S\) 的 Lyndon Tree。
如图是一个串 \(\text{aababbaababa}\) 的 Lyndon Tree:(论文图真香)
我们先来观察性质,他是一棵笛卡尔树,所以如果你不嫌麻烦的话可以用 Lyndon Tree 实现线性 runs。
性质:如果我们的 Lyndon Tree 有一个 Lyndon word,记为 \(s[i..j]\),那么从 \(i\) 到 \(j\) 所有叶节点的 LCA 是 \([i..k]\),其中 \(i \le j \le k\)。如果这个 Lyndon Word 是从 \(i\) 开始的最长 Lyndon 串,则 LCA 一定是一个右子节点。如果 \(i=1\),则 LCA 为根节点。
我们有一个串 \(S\) 有 \(p,q\) 两个周期,满足 \(p<q\) 且 \(p + q \le |S|\),则 \(\gcd(p,q)\) 也是 \(S\) 的周期。
这个应该很显然把……
显然,\(S_i = S_{i+q} = S_{i-p+q} = S_{i-p} = S_{i-p+q}\)。所以 \(q-p\) 也是 \(S\) 的周期。
然后辗转相除法可以证得 \(\gcd(p,q)\) 也是 \(S\) 的周期。\(\blacksquare\)
事实上有一个更强的结论。
我们有一个串 \(S\) 有 \(p,q\) 两个周期,满足 \(p<q\) 且 \(p + q - \gcd(p,q) \le |S|\),则 \(\gcd(p,q)\) 也是 \(S\) 的周期。
经典问题,可能是 Lyndon tree 的唯一应用= =
\(S\) 是一个串,询问 \(S\) 的某个子串有没有不超过长度一半的周期,如果有就求出最小周期。
我们定义一个 run 是一个 \(\operatorname{exrun}(i,j)\) 当且仅当他是 \((i', j', p)\) 满足 \(i' \le i, j' \ge j, p \le (j-i+1)/2\)。
根据 [Theorem 5.2.1 Weak Periodicity Lemma],exrun 唯一。如果我们找出了 \(\operatorname{exrun}(i,j)\),则我们就可以回答询问。
我们先构造出两棵 Lyndon Tree \(\operatorname {Ltree}^0(S)\) 和 \(\operatorname {Ltree}^1(S)\)。
然后我们在两棵 Lyndon Tree 上面分别寻找从 \(i\) 到 \((i+j)/2\) 的 LCA,判断其右子节点的 runs 是否符合 exrun。
证明留作习题。
听说可以搬到树上什么的……还是不要出这种东西吧……
Squares: 能表示成 \(x^2\) 的串。
Primitive Squares: 不能再拆的 Squares,即最小周期为自身长度的一半的字符串。如 \(x^2+x\) 一定是一个 Primitive Square。
我们有 3 个 Primitive Squares,为 \(u^2\),\(v^2\) 和 \(w^2\),满足 \(|u|<|v|<|w|\)。
则我们有 \(|u| + |v| \le |w|\)。
反证法,假设 \(|u|+|v| > |w|\) 成立,则 \(|w|-|v|\) 为 \(u\) 和 \(v\) 的周期。
如果 \(2|u| \le |v|\),则 \(|w|-|v| < |u|\) 是 \(u^2\) 的周期,故 \(u^2\) 不是一个 Primitive Square,矛盾。
否则,\(|u|\) 和 \(|w|-|v|\) 都是 \(v\) 的周期,而 \(u^2\) 是一个 Primitive Square,所以 \(u\) 不是一个周期串,则 \(|u|+|w|>2|v|\)。
我们令 \(w=vs_1\),\(u=s_1s_3\),\(v=ws_2=s_1s_3s_2\)。显然,\(|s_2| < |s_1|\)。
考虑串 \(s_3 s_2\),它有周期 \(|s_2|\)。由于 \(|s_1|\) 是 \(u\) 的周期,可得 \(s_3s_2\) 是 \(u\) 的前缀,所以 \(|s_3|\) 也是它的周期。
根据 [Theorem 5.2.1 Weak Periodicity Lemma],\(r=\gcd(|s_2|,|s_3|)\) 是它的周期。而 \(|s_2|\) 本身同时是 \(u\) 的周期,因此可得 \(r\) 是 \(u\) 的周期。
接着考虑串 \(u=s_1s_3\)。它的周期有 \(|s_1|\) 和 \(r\),而 \(r\le|s_3|\),根据 [Theorem 5.2.1 Weak Periodicity Lemma],\(r'=\gcd(r,|s_1|)\) 也是 \(u\) 的周期。然而 \(|s_1|\) 和 \(|s_3|\) 都是 \(r'\) 的倍数,这表示 \(u^2\) 也有周期 \(r'\),矛盾。\(\blacksquare\)
由此我们有 2 个显然的结论:
Primitive Squares 的数量级为 \(\mathcal O(n\log n)\)。
显然。
本质不同的 Primitive Squares 不超过 \(2n\) 个。
每个位置最多出现 2 个以这个位置为结尾的 Primitive Squares。\(\blacksquare\)
首先我们可以想到一种十分显然的方法,因为数量级为 \(\mathcal O(n \log n)\),而 primitive squares 一定是 runs 的一部分,所以所以我们可以在每一个 runs 上面暴力求是否有 runs。时间复杂度 \(\mathcal O(n \log n)\)。
方法二,枚举平方串的长度,lcp。
我在出 LOJ #173 的时候,发现有一个非常强的暴力,就是暴力地去做 lcp 和 lcs。
后来我经过不懈努力使用 aaaaa...ab
这个串卡掉了他,但是其实考场上如果时间不够的话写这个暴力其实是非常优秀的,毕竟我见过的两个 runs 题都没有卡这个暴力。
给你个串 \(S\ (|S| \le 10^5)\),切成若干段,其中每一段是 primitive 的,且相邻两个串是不相等的。求方案数模 \(998244353\)。
首先考虑 dp。\(f_i\) 表示 \(\operatorname{suf}(i)\) 的划分方案数。可以得到这么一个柿子:(当然,如果考虑划分前缀也是可以的
\[f_i = \sum _{j<i}f_j \cdot (-1)^{\sigma(j\cdots i) - 1}
\]
其中,\(\sigma(s)\) 是这样定义的:对于字符串 \(s\),他的最小循环节的循环的长度为 \(\sigma(s)\)。
比如,对于串 \(\color{blue}\text{abcabcabc}\color{red}\text{defg}\),他的 \(\sigma = 3\),因为 \(\color{blue}\text{abc}\) 为最小循环节,循环了 \(3\) 次。
这个复杂度是 \(\mathcal O(n^2)\),无法接受。
考虑优化:把循环串分成一个初始的循环串和一些循环,分开计算贡献。
我们认为一个循环串初始的时候的最小周期一定是大于等于长度的一半的(否则会重复的)。
\[f_i = \left[\sum_{j<i}f_j\right] + \sum _{j<i} f_j\left[(-1)^{\sigma(j\cdots i)-1}-1\right]
\]
然于是我们把 Primitive squares 和 runs 求出来转移就好了。时间复杂度为 \(\mathcal O(n \log n)\)。
留作习题。
手机扫一扫
移动阅读更方便
你可能感兴趣的文章