后缀数组简介

字符串 / sa

本地源文件:docs/string__sa.md

后缀数组简介

一些约定

字符串相关的定义请参考 字符串基础

字符串下标从 11 开始.

字符串 𝑠s 的长度为 𝑛n

" 后缀 𝑖i" 代指以第 𝑖i 个字符开头的后缀,存储时用 𝑖i 代表字符串 𝑠s 的后缀 𝑠[𝑖…𝑛]s[i…n]

后缀数组是什么?

后缀数组(Suffix Array)主要关系到两个数组:𝑠𝑎sa 和 𝑟𝑘rk

其中,𝑠𝑎[𝑖]sa[i] 表示将所有后缀排序后第 𝑖i 小的后缀的编号,也是所说的后缀数组,后文也称编号数组 𝑠𝑎sa

𝑟𝑘[𝑖]rk[i] 表示后缀 𝑖i 的排名,是重要的辅助数组,后文也称排名数组 𝑟𝑘rk

这两个数组满足性质:𝑠𝑎[𝑟𝑘[𝑖]] =𝑟𝑘[𝑠𝑎[𝑖]] =𝑖sa[rk[i]]=rk[sa[i]]=i

解释

后缀数组示例:

后缀数组怎么求?

O(n^2logn) 做法

相信这个做法大家还是能自己想到的:将盛有全部后缀字符串的数组进行 sort 排序,由于排序进行 𝑂(𝑛log⁡𝑛)O(nlog⁡n) 次字符串比较,每次字符串比较要 𝑂(𝑛)O(n) 次字符比较,所以这个排序是 𝑂(𝑛2log⁡𝑛)O(n2log⁡n) 的时间复杂度.

O(nlog^2n) 做法

这个做法要用到倍增的思想.

首先对字符串 𝑠s 的所有长度为 11 的子串,即每个字符进行排序,得到排序后的编号数组 𝑠𝑎1sa1 和排名数组 𝑟𝑘1rk1

倍增过程:

  1. 用两个长度为 11 的子串的排名,即 𝑟𝑘1[𝑖]rk1[i] 和 𝑟𝑘1[𝑖 +1]rk1[i+1],作为排序的第一第二关键字,就可以对字符串 𝑠s 的每个长度为 22 的子串:{𝑠[𝑖…min(𝑖 +1,𝑛)] | 𝑖 ∈[1, 𝑛]}{s[i…min(i+1,n)] | i∈[1, n]} 进行排序,得到 𝑠𝑎2sa2 和 𝑟𝑘2rk2
  1. 之后用两个长度为 22 的子串的排名,即 𝑟𝑘2[𝑖]rk2[i] 和 𝑟𝑘2[𝑖 +2]rk2[i+2],作为排序的第一第二关键字,就可以对字符串 𝑠s 的每个长度为 44 的子串:{𝑠[𝑖…min(𝑖 +3,𝑛)] | 𝑖 ∈[1, 𝑛]}{s[i…min(i+3,n)] | i∈[1, n]} 进行排序,得到 𝑠𝑎4sa4 和 𝑟𝑘4rk4
  1. 以此倍增,用长度为 𝑤/2w/2 的子串的排名,即 𝑟𝑘𝑤/2[𝑖]rkw/2[i] 和 𝑟𝑘𝑤/2[𝑖 +𝑤/2]rkw/2[i+w/2],作为排序的第一第二关键字,就可以对字符串 𝑠s 的每个长度为 𝑤w 的子串 𝑠[𝑖…min(𝑖 +𝑤 −1, 𝑛)]s[i…min(i+w−1, n)] 进行排序,得到 𝑠𝑎𝑤saw 和 𝑟𝑘𝑤rkw.其中,类似字母序排序规则,当 𝑖 +𝑤 >𝑛i+w>n 时,𝑟𝑘𝑤[𝑖 +𝑤]rkw[i+w] 视为无穷小;
  1. 𝑟𝑘𝑤[𝑖]rkw[i] 即是子串 𝑠[𝑖…𝑖 +𝑤 −1]s[i…i+w−1] 的排名,这样当 𝑤 ⩾𝑛w⩾n 时,得到的编号数组 𝑠𝑎𝑤saw,也就是我们需要的后缀数组.

过程

倍增排序示意图:

显然倍增的过程是 𝑂(log⁡𝑛)O(log⁡n),而每次倍增用 sort 对子串进行排序是 𝑂(𝑛log⁡𝑛)O(nlog⁡n),而每次子串的比较花费 22 次字符比较;

除此之外,每次倍增在 sort 排序完后,还有额外的 𝑂(𝑛)O(n) 时间复杂度的,更新 𝑟𝑘rk 的操作,但是相对于 𝑂(𝑛log⁡𝑛)O(nlog⁡n) 被忽略不计;

所以这个算法的时间复杂度就是 𝑂(𝑛log2⁡𝑛)O(nlog2⁡n)

实现

---|---

### O(nlogn) 做法

在刚刚的 𝑂(𝑛log2⁡𝑛)O(nlog2⁡n)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 做法中,单次排序是 𝑂(𝑛log⁡𝑛)O(nlog⁡n)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 的,如果能 𝑂(𝑛)O(n)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 排序,就能 𝑂(𝑛log⁡𝑛)O(nlog⁡n)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 计算后缀数组了.

前置知识:[计数排序](../../basic/counting-sort/),[基数排序](../../basic/radix-sort/).

由于计算后缀数组的过程中排序的关键字是排名,值域为 𝑂(𝑛)O(n)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),并且是一个双关键字的排序,可以使用基数排序优化至 𝑂(𝑛)O(n)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).

实现

---|---

一些常数优化

如果你把上面那份代码交到 LOJ #111: 后缀排序 上:

这是因为,上面那份代码的常数的确很大.

第二关键字无需计数排序

思考一下第二关键字排序的实质,其实就是把超出字符串范围(即 𝑠𝑎[𝑖] +𝑤 >𝑛sa[i]+w>n)的 𝑠𝑎[𝑖]sa[i] 放到 𝑠𝑎sa 数组头部,然后把剩下的依原顺序放入:

---|---

#### 优化计数排序的值域

每次对 𝑟𝑘rk![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 进行更新之后,我们都计算了一个 𝑝p![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),这个 𝑝p![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 即是 𝑟𝑘rk![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 的值域,将值域改成它即可.

#### 若排名都不相同可直接生成后缀数组

考虑新的 𝑟𝑘rk![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 数组,若其值域为 [1,𝑛][1,n]![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 那么每个排名都不同,此时无需再排序.

实现

---|---

O(n) 做法

在一般的题目中,常数较小的倍增求后缀数组是完全够用的,求后缀数组以外的部分也经常有 𝑂(𝑛log⁡𝑛)O(nlog⁡n) 的复杂度,倍增求解后缀数组不会成为瓶颈.

但如果遇到特殊题目、时限较紧的题目,或者是你想追求更短的用时,就需要学习 𝑂(𝑛)O(n) 求后缀数组的方法.

SA-IS

可以参考 诱导排序与 SA-IS 算法,另外它的 评论页面 也有参考价值.

DC3

可以参考[[2009] 后缀数组——处理字符串的有力工具 by. 罗穗骞](https://github.com/OI-wiki/libs/blob/master/%E9%9B%86%E8%AE%AD%E9%98%9F%E5%8E%86%E5%B9%B4%E8%AE%BA%E6%96%87/%E5%9B%BD%E5%AE%B6%E9%9B%86%E8%AE%AD%E9%98%9F2009%E8%AE%BA%E6%96%87%E9%9B%86/11.%E7%BD%97%E7%A9%97%E9%AA%9E%E3%80%8A%E5%90%8E%E7%BC%80%E6%95%B0%E7%BB%84%E2%80%94%E2%80%94%E5%A4%84%E7%90%86%E5%AD%97%E7%AC%A6%E4%B8%B2%E7%9A%84%E6%9C%89%E5%8A%9B%E5%B7%A5%E5%85%B7%E3%80%8B/%E5%90%8E%E7%BC%80%E6%95%B0%E7%BB%84%E2%80%94%E2%80%94%E5%A4%84%E7%90%86%E5%AD%97%E7%AC%A6%E4%B8%B2%E7%9A%84%E6%9C%89%E5%8A%9B%E5%B7%A5%E5%85%B7.pdf "\[2009\] 后缀数组——处理字符串的有力工具 by. 罗穗骞").

后缀数组的应用

寻找最小的循环移动位置

将字符串 𝑆S 复制一份变成 𝑆𝑆SS 就转化成了后缀排序问题.

例题:「JSOI2007」字符加密

在字符串中找子串

任务是在线地在主串 𝑇T 中寻找模式串 𝑆S.在线的意思是,我们已经预先知道知道主串 𝑇T,但是当且仅当询问时才知道模式串 𝑆S.我们可以先构造出 𝑇T 的后缀数组,然后查找子串 𝑆S.若子串 𝑆S 在 𝑇T 中出现,它必定是 𝑇T 的一些后缀的前缀.因为我们已经将所有后缀排序了,我们可以通过在 𝑝p 数组中二分 𝑆S 来实现.比较子串 𝑆S 和当前后缀的时间复杂度为 𝑂(|𝑆|)O(|S|),因此找子串的时间复杂度为 𝑂(|𝑆|log⁡|𝑇|)O(|S|log⁡|T|).注意,如果该子串在 𝑇T 中出现了多次,每次出现都是在 𝑝p 数组中相邻的.因此出现次数可以通过再次二分找到,输出每次出现的位置也很轻松.

从字符串首尾取字符最小化字典序

例题:「USACO07DEC」Best Cow Line

题意:给你一个字符串,每次从首或尾取一个字符组成字符串,问所有能够组成的字符串中字典序最小的一个.

题解

暴力做法就是每次最坏 𝑂(𝑛)O(n) 地判断当前应该取首还是尾(即比较取首得到的字符串与取尾得到的反串的大小),只需优化这一判断过程即可.

由于需要在原串后缀与反串后缀构成的集合内比较大小,可以将反串拼接在原串后,并在中间加上一个没出现过的字符(如 #,代码中可以直接使用空字符),求后缀数组,即可 𝑂(1)O(1) 完成这一判断.

参考代码

---|---

## height 数组

### LCP(最长公共前缀)

两个字符串 𝑆S![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 和 𝑇T![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 的 LCP 就是最大的 𝑥x![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)(𝑥 ≤min(|𝑆|,|𝑇|)x≤min(|S|,|T|)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)) 使得 𝑆𝑖 =𝑇𝑖 (∀ 1 ≤𝑖 ≤𝑥)Si=Ti (∀ 1≤i≤x)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).

下文中以 𝑙𝑐𝑝(𝑖,𝑗)lcp(i,j)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 表示后缀 𝑖i![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 和后缀 𝑗j![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 的最长公共前缀(的长度).

### height 数组的定义

ℎ𝑒𝑖𝑔ℎ𝑡[𝑖] =𝑙𝑐𝑝(𝑠𝑎[𝑖],𝑠𝑎[𝑖 −1])height[i]=lcp(sa[i],sa[i−1])![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),即第 𝑖i![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 名的后缀与它前一名的后缀的最长公共前缀.

ℎ𝑒𝑖𝑔ℎ𝑡[1]height[1]![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 可以视作 00![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).

### O(n) 求 height 数组需要的一个引理

ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖]] ≥ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖 −1]] −1height[rk[i]]≥height[rk[i−1]]−1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)

证明

当 ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖 −1]] ≤1height[rk[i−1]]≤1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 时,上式显然成立(右边小于等于 00![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)).

当 ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖 −1]] >1height[rk[i−1]]>1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 时:

根据 ℎ𝑒𝑖𝑔ℎ𝑡height![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 定义,有 𝑙𝑐𝑝(𝑠𝑎[𝑟𝑘[𝑖 −1]],𝑠𝑎[𝑟𝑘[𝑖 −1] −1]) =ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖 −1]] >1lcp(sa[rk[i−1]],sa[rk[i−1]−1])=height[rk[i−1]]>1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).

既然后缀 𝑖 −1i−1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 和后缀 𝑠𝑎[𝑟𝑘[𝑖 −1] −1]sa[rk[i−1]−1]![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 有长度为 ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖 −1]]height[rk[i−1]]![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 的最长公共前缀,

那么不妨用 𝑎𝐴aA![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 来表示这个最长公共前缀.(其中 𝑎a![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 是一个字符,𝐴A![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 是长度为 ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖 −1]] −1height[rk[i−1]]−1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 的字符串,非空)

那么后缀 𝑖 −1i−1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 可以表示为 𝑎𝐴𝐷aAD![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),后缀 𝑠𝑎[𝑟𝑘[𝑖 −1] −1]sa[rk[i−1]−1]![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 可以表示为 𝑎𝐴𝐵aAB![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).(𝐵 <𝐷B<D![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),𝐵B![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 可能为空串,𝐷D![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 非空)

进一步地,后缀 𝑖i![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 可以表示为 𝐴𝐷AD![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),存在后缀(𝑠𝑎[𝑟𝑘[𝑖 −1] −1] +1sa[rk[i−1]−1]+1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7))𝐴𝐵AB![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).

因为后缀 𝑠𝑎[𝑟𝑘[𝑖] −1]sa[rk[i]−1]![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 在大小关系的排名上仅比后缀 𝑠𝑎[𝑟𝑘[𝑖]]sa[rk[i]]![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 也就是后缀 𝑖i![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),小一位,而 𝐴𝐵 <𝐴𝐷AB<AD![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).

所以 𝐴𝐵 ⩽AB⩽![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 后缀 𝑠𝑎[𝑟𝑘[𝑖] −1] <𝐴𝐷sa[rk[i]−1]<AD![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),显然后缀 𝑖i![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 和后缀 𝑠𝑎[𝑟𝑘[𝑖] −1]sa[rk[i]−1]![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 有公共前缀 𝐴A![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).

于是就可以得出 𝑙𝑐𝑝(𝑖,𝑠𝑎[𝑟𝑘[𝑖] −1])lcp(i,sa[rk[i]−1])![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 至少是 ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖 −1]] −1height[rk[i−1]]−1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),也即 ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖]] ≥ℎ𝑒𝑖𝑔ℎ𝑡[𝑟𝑘[𝑖 −1]] −1height[rk[i]]≥height[rk[i−1]]−1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).

### O(n) 求 height 数组的代码实现

利用上面这个引理暴力求即可:

---|---

𝑘k 不会超过 𝑛n,最多减 𝑛n 次,所以最多加 2𝑛2n 次,总复杂度就是 𝑂(𝑛)O(n)

height 数组的应用

两子串最长公共前缀

𝑙𝑐𝑝(𝑠𝑎[𝑖],𝑠𝑎[𝑗]) =min{ℎ𝑒𝑖𝑔ℎ𝑡[𝑖 +1..𝑗]}lcp(sa[i],sa[j])=min{height[i+1..j]}

感性理解:如果 ℎ𝑒𝑖𝑔ℎ𝑡height 一直大于某个数,前这么多位就一直没变过;反之,由于后缀已经排好序了,不可能变了之后变回来.

严格证明可以参考[[2004] 后缀数组 by. 许智磊](https://github.com/OI-wiki/libs/blob/master/%E9%9B%86%E8%AE%AD%E9%98%9F%E5%8E%86%E5%B9%B4%E8%AE%BA%E6%96%87/%E5%9B%BD%E5%AE%B6%E9%9B%86%E8%AE%AD%E9%98%9F2004%E8%AE%BA%E6%96%87%E9%9B%86/%E8%AE%B8%E6%99%BA%E7%A3%8A--%E5%90%8E%E7%BC%80%E6%95%B0%E7%BB%84.pdf "\[2004\] 后缀数组 by. 许智磊").

有了这个定理,求两子串最长公共前缀就转化为了 RMQ 问题

比较一个字符串的两个子串的大小关系

假设需要比较的是 𝐴 =𝑆[𝑎..𝑏]A=S[a..b] 和 𝐵 =𝑆[𝑐..𝑑]B=S[c..d] 的大小关系.

若 𝑙𝑐𝑝(𝑎,𝑐) ≥min(|𝐴|,|𝐵|)lcp(a,c)≥min(|A|,|B|),𝐴 <𝐵 ⟺ |𝐴| <|𝐵|A<B⟺|A|<|B|

否则,𝐴 <𝐵 ⟺ 𝑟𝑘[𝑎] <𝑟𝑘[𝑐]A<B⟺rk[a]<rk[c]

不同子串的数目

子串就是后缀的前缀,所以可以枚举每个后缀,计算前缀总数,再减掉重复.

「前缀总数」其实就是子串个数,为 𝑛(𝑛 +1)/2n(n+1)/2

如果按后缀排序的顺序枚举后缀,每次新增的子串就是除了与上一个后缀的 LCP 剩下的前缀.这些前缀一定是新增的,否则会破坏 𝑙𝑐𝑝(𝑠𝑎[𝑖],𝑠𝑎[𝑗]) =min{ℎ𝑒𝑖𝑔ℎ𝑡[𝑖 +1..𝑗]}lcp(sa[i],sa[j])=min{height[i+1..j]} 的性质.只有这些前缀是新增的,因为 LCP 部分在枚举上一个前缀时计算过了.

所以答案为:

𝑛(𝑛+1)2 −𝑛∑𝑖=2ℎ𝑒𝑖𝑔ℎ𝑡[𝑖]n(n+1)2−∑i=2nheight[i]

出现至少 k 次的子串的最大长度

例题:「USACO06DEC」Milk Patterns

题解

出现至少 𝑘k 次意味着后缀排序后有至少连续 𝑘k 个后缀以这个子串作为公共前缀.

所以,求出每相邻 𝑘 −1k−1 个 ℎ𝑒𝑖𝑔ℎ𝑡height 的最小值,再求这些最小值的最大值就是答案.

可以使用单调队列 𝑂(𝑛)O(n) 解决,但使用其它方式也足以 AC.

参考代码

---|---

### 是否有某字符串在文本串中至少不重叠地出现了两次

可以二分目标串的长度 |𝑠||s|![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),将 ℎh![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 数组划分成若干个连续 LCP 大于等于 |𝑠||s|![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 的段,利用 RMQ 对每个段求其中出现的数中最大和最小的下标,若这两个下标的距离满足条件,则一定有长度为 |𝑠||s|![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 的字符串不重叠地出现了两次.

### 连续的若干个相同子串

我们可以枚举连续串的长度 |𝑠||s|![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7),按照 |𝑠||s|![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 对整个串进行分块,对相邻两块的块首进行 LCP 与 LCS 查询,具体可见[[2009] 后缀数组——处理字符串的有力工具](https://github.com/OI-wiki/libs/blob/master/%E9%9B%86%E8%AE%AD%E9%98%9F%E5%8E%86%E5%B9%B4%E8%AE%BA%E6%96%87/%E5%9B%BD%E5%AE%B6%E9%9B%86%E8%AE%AD%E9%98%9F2009%E8%AE%BA%E6%96%87%E9%9B%86/11.%E7%BD%97%E7%A9%97%E9%AA%9E%E3%80%8A%E5%90%8E%E7%BC%80%E6%95%B0%E7%BB%84%E2%80%94%E2%80%94%E5%A4%84%E7%90%86%E5%AD%97%E7%AC%A6%E4%B8%B2%E7%9A%84%E6%9C%89%E5%8A%9B%E5%B7%A5%E5%85%B7%E3%80%8B/%E5%90%8E%E7%BC%80%E6%95%B0%E7%BB%84%E2%80%94%E2%80%94%E5%A4%84%E7%90%86%E5%AD%97%E7%AC%A6%E4%B8%B2%E7%9A%84%E6%9C%89%E5%8A%9B%E5%B7%A5%E5%85%B7.pdf "\[2009\] 后缀数组——处理字符串的有力工具 by. 罗穗骞").

例题:[「NOI2016」优秀的拆分](https://loj.ac/p/2083).

### 结合并查集

某些题目求解时要求你将后缀数组划分成若干个连续 LCP 长度大于等于某一值的段,亦即将 ℎh![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 数组划分成若干个连续最小值大于等于某一值的段并统计每一段的答案.如果有多次询问,我们可以将询问离线.观察到当给定值单调递减的时候,满足条件的区间个数总是越来越少,而新区间都是两个或多个原区间相连所得,且新区间中不包含在原区间内的部分的 ℎh![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 值都为减少到的这个值.我们只需要维护一个并查集,每次合并相邻的两个区间,并维护统计信息即可.

经典题目:[「NOI2015」品酒大会](https://uoj.ac/problem/131)

### 结合线段树

某些题目让你求满足条件的前若干个数,而这些数又在后缀排序中的一个区间内.这时我们可以用归并排序的性质来合并两个结点的信息,利用线段树维护和查询区间答案.

### 结合单调栈

例题:[「AHOI2013」差异](https://loj.ac/problem/2377)

题解

被加数的前两项很好处理,为 𝑛(𝑛 −1)(𝑛 +1)/2n(n−1)(n+1)/2![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)(每个后缀都出现了 𝑛 −1n−1![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 次,后缀总长是 𝑛(𝑛 +1)/2n(n+1)/2![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)),关键是最后一项,即后缀的两两 LCP.

我们知道 𝑙𝑐𝑝(𝑖,𝑗) =𝑘lcp(i,j)=k![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 等价于 min{ℎ𝑒𝑖𝑔ℎ𝑡[𝑖 +1..𝑗]} =𝑘min{height[i+1..j]}=k![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7).所以,可以把 𝑙𝑐𝑝(𝑖,𝑗)lcp(i,j)![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 记作 min{𝑥|𝑖 +1 ≤𝑥 ≤𝑗,ℎ𝑒𝑖𝑔ℎ𝑡[𝑥] =𝑙𝑐𝑝(𝑖,𝑗)}min{x|i+1≤x≤j,height[x]=lcp(i,j)}![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 对答案的贡献.

考虑每个位置对答案的贡献是哪些后缀的 LCP,其实就是从它开始向左若干个连续的 ℎ𝑒𝑖𝑔ℎ𝑡height![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 大于它的后缀中选一个,再从向右若干个连续的 ℎ𝑒𝑖𝑔ℎ𝑡height![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7) 不小于它的后缀中选一个.这个东西可以用 [单调栈](../../ds/monotonous-stack/) 计算.

单调栈部分类似于 [Luogu P2659 美丽的序列](https://www.luogu.com.cn/problem/P2659) 以及 [悬线法](../../misc/hoverline/).

参考代码

---|---

类似的题目:「HAOI2016」找相同字符

习题

参考资料

本页面中(4070a9b 引入的部分)主要译自博文 Суффиксный массив 与其英文翻译版 Suffix Array.其中俄文版版权协议为 Public Domain + Leave a Link;英文版版权协议为 CC-BY-SA 4.0.

论文:

  1. [[2004] 后缀数组 by. 许智磊](https://github.com/OI-wiki/libs/blob/master/%E9%9B%86%E8%AE%AD%E9%98%9F%E5%8E%86%E5%B9%B4%E8%AE%BA%E6%96%87/%E5%9B%BD%E5%AE%B6%E9%9B%86%E8%AE%AD%E9%98%9F2004%E8%AE%BA%E6%96%87%E9%9B%86/%E8%AE%B8%E6%99%BA%E7%A3%8A--%E5%90%8E%E7%BC%80%E6%95%B0%E7%BB%84.pdf "\[2004\] 后缀数组 by. 许智磊")
  1. [[2009] 后缀数组——处理字符串的有力工具 by. 罗穗骞](https://github.com/OI-wiki/libs/blob/master/%E9%9B%86%E8%AE%AD%E9%98%9F%E5%8E%86%E5%B9%B4%E8%AE%BA%E6%96%87/%E5%9B%BD%E5%AE%B6%E9%9B%86%E8%AE%AD%E9%98%9F2009%E8%AE%BA%E6%96%87%E9%9B%86/11.%E7%BD%97%E7%A9%97%E9%AA%9E%E3%80%8A%E5%90%8E%E7%BC%80%E6%95%B0%E7%BB%84%E2%80%94%E2%80%94%E5%A4%84%E7%90%86%E5%AD%97%E7%AC%A6%E4%B8%B2%E7%9A%84%E6%9C%89%E5%8A%9B%E5%B7%A5%E5%85%B7%E3%80%8B/%E5%90%8E%E7%BC%80%E6%95%B0%E7%BB%84%E2%80%94%E2%80%94%E5%A4%84%E7%90%86%E5%AD%97%E7%AC%A6%E4%B8%B2%E7%9A%84%E6%9C%89%E5%8A%9B%E5%B7%A5%E5%85%B7.pdf "\[2009\] 后缀数组——处理字符串的有力工具 by. 罗穗骞")
本页面最近更新: 2026/1/7 08:56:54,更新历史 发现错误?想一起完善?在 GitHub 上编辑此页! 本页面贡献者:Ir1d, ksyx, Early0v0, ouuan, countercurrent-time, StudyingFather, Tiphereth-A, H-J-Granger, NachtgeistW, CCXXXI, Enter-tainer, hsfzLZH1, minghu6, AngelKitty, Chrogeek, cjsoft, diauweb, ezoixx130, GekkaSaori, Henry-ZHR, HeRaNO, iamtwz, Konano, LovelyBuggies, Makkiy, Marcythm, mgt, P-Y-Y, PotassiumWings, SamZhangQingChuan, sshwy, Suyun514, weiyong1024, alphagocc, billchenchina, ChungZH, Clouder0, Eletary, gavinliu266, GavinZhengOI, Gesrua, gi-b716, Great-designer, hly1204, Hukeqing, ImpleLee, interestingLSY, Joob1n, kenlig, kxccc, LingeZ3z, lychees, Menci, Nanarikom, Peanut-Tang, r-value, renbaoshuo, SukkaW, TOMWT-qwq, WAAutoMaton, werner-yin, yafngzh, yjl9903 本页面的全部内容在CC BY-SA 4.0SATA 协议之条款下提供,附加条款亦可能应用