Typesetting math: 0%

SAM学习笔记

-1.参考资料

oi wiki

嗯,改了一些可能是错误的地方(迫真

0.定义

定义 ΣΣ 表示字符集。

定义 ab 表示 ab 的后缀。

定义 ab 表示 ab 的真后缀。

定义 |s| 表示 s 的长度。

定义 uv 表示 uv 的路径上节点组成的集合。

1.关于SAM

  • Q:SAM 是什么?可以吃吗?
  • A:SAM 全称后缀自动机(Suffix Automation),是一个接受 s 的所有后缀的最小 DFA ,不可以吃。
  • Q:DFA 是什么?可以吃吗?
  • A:DFA 全称确定有限状态自动机(Deterministic Finite Automaton),也不可以吃。DFA 由五个部分组成:
  1. 字符集(Σ),自动机中的所有字符均在字符集中。
  2. 状态集合(Q),如果将 DFA 看成 DAG ,那么 Q=V
  3. 起始状态(s),sQ,是一个特殊的状态。
  4. 接受状态(F),FQ,是一堆特殊的状态。
  5. 转移函数(δ),如果将 DFA 看成 DAG ,那么 δ=w ,且 δ:V×VΣ,此处乘法表示笛卡尔积。

总结一下:

  • SAM 是一个特殊的 DAG,结点称为状态,边称为转移。
  • 存在源点(初始状态) t0,其他结点均由 t0 可达。
  • 存在多个汇点(终止状态),满足任意一条由源点到汇点的路径上转移函数拼接组成的字符串均为原字符串的后缀。(这也是其被称为后缀自动机的原因)
  • 在所有满足以上条件的自动机中,SAM 是结点个数最少的。具体地,|V|=2n1,|E|=3n4

由定义我们可以得到一个简单的推论:由源点引出的所有路径上转移函数拼接组成的字符串(包括连到自己的路径,这种情况下我们视为空串)构成的集合即为该字符串的所有子串构成的集合,证明显然。

2.SAM 举例

盗来了几张 oi-wiki 上的图:

  1. 空串:
  2. s=a
  3. s=aa
  4. s=ab
  5. s=abb
  6. s=abbb

3.线性构造 SAM

Note:下面出现的一切字符串都是 s 的非空子串。

3.1 endpos

定义 endpos(t) 表示在字符串 st 的所有结束位置。

例如,当 s=abcbc 的时候,endpos(bc)={2,4}

我们称 s,t 为等价类,当且仅当 endpos(s)=endpos(t),于是我们可以将 s 的所有非空子串分为若干等价类。

则显然 SAM 上除了初始状态的每一个状态都对应一个等价类。

因为每个状态都能对应多个 endpos 相同的子串,也就是等价类,所以等价类的个数加上 1,就是SAM状态的个数!

下面我们会给出一些关于 endpos 的性质。

Lemma 1|u||w|,那么 uw 是等价类的充要条件是每次 u 出现都是作为 w 的后缀。

证明显然。

Lemma 2|u||w|,有:

{endpos(w)endpos(u)uwendpos(u)endpos(w)=otherwise

Proof

endpos(u)endpos(w),则 uw 在相同位置结束,则 uw,所以 endpos(w)endpos(u),反之亦然,得证。

Lemma 3
(1)对于同一个等价类中的两个子串 s,t,且 |s||t|,有 st
(2)一个等价类中的所有子串的长度,可以恰好不重复地覆盖一段连续区间。

Proof

如果等价类中只有一个子串,那么结论显然成立。

考虑有至少两个子串的情况,由 Lemma 1,短的子串是长的子串的后缀,所以等价类中子串长度两两不同。

设最长串为 u,最短串为 w,由 Lemma 1,wu

考虑长度在 [|u|,|w|] 之间的 w 的后缀,由 Lemma 1,这些字符串都在等价类中,所以(2)得证。

又因为等价类中子串长度两两不同,所以等价类中只有这些字符串,所以(1)得证。

取一个状态 vt0,我们目前已经得到了如下结论:

  1. v 对应一个等价类
  2. 如果设 w 是这个等价类中最长的一个,那么其他的字符串都是 w 的后缀。
  3. 如果将 w 的后缀按照长度降序排序,那么 w 的前几个后缀都在这个等价类中,且剩下的都不在。

我们设 t 是所有 w 的后缀中,最长的与 w 不在一个等价类中的后缀,那么定义 link(v)=t,即把 v 的后缀链接连到 t 上。

也就是说,一个后缀链接连接到对应于 w 的最长后缀的另一个等价类的状态。

下面给出一些性质。

P.S.为了方便起见,定义 endpos(t0)={1,0,,|S|1}

Lemma 4:所有后缀链接构成一棵根节点为 t0 的树,一般称之为 Parent Tree

Proof

证:vt0,由 Lemma 3,后缀链接连接的状态对应的字符串长度严格小于该状态字符串长度,因此必然会通向 t0,又因为恰有 |V|1 条边,所以构成树。

Lemma 5endpos(s)endpos(link(s))。亦即,parent tree的子节点的 endpos 是父节点的 endpos 的子集。

Proof

证:由定义 link(s)s 由 Lemma 2 得证。

这是一张 parent tree 对应自动机的图:

3.3 一些定义&小结

在学习算法之前,先来一些定义。

  • 定义 longest(v) 为状态 v 对应等价类中最长的字符串,len(v)=|longest(v)|
  • 定义 shortest(v) 为状态 v 对应等价类中最短的字符串,minlen(v)=|shortest(v)|

所以由 Lemma 3 以及 link 的定义立得:

minlen(v)=len(link(v))+1

下面是以上定理的总结:

  • s 的子串可以根据 endpos 划分为若干等价类。
  • SAM 由初始状态和等价类对应的状态构成。
  • 对于状态 vt0,有以下结论:
  1. 此状态对应的字符串均为 longest(v) 后缀,且长度恰好覆盖 [minlen(v),len(v)]
  2. endpos(v)endpos(link(v))
  3. parent tree 上 t0v 存在,且路径上每个节点对应的区间互不相交,并集恰为 [0,link(v)]

3.4 算法

这个算法是在线算法,也就是说可以每次加一个字符,同时更新SAM。

下面我们以字符串 abcbc 对应 SAM 为例,考虑加上一个 b(下面标记为字符 c)。

下面的图分别表示 SAM 和 parent tree。

注:节点上的字符串表示 longest(v),边上的字符串表示转移函数的值,蓝色点是初始状态,绿色点是终止状态。实际算法时,我们只会对每个点记录转移到的点,lenlink

在算法流程中,我们取 last 表示整个字符串对应的状态,每次更新完自动机后再更新 last,比如在上面的图中,last 就是 ABCBC 对应的点,后面会标记为橙色。

  1. 新建状态 cur,并令 len(cur)=len(last)+1longest(v)=s,后面会用红色标记。
  2. last 开始在 parent tree 上暴力跳祖先,只要跳到有到 c 的转移就停下,并将跳过的路径上的所有点
    (包括 last,不包括该点)全部向 cur 连转移为 c 的边,将此点设为 p。比如在这个例子中,找到的 p 就是 BC 对应的状态,用黄色标记,下面是操作好的自动机。
  3. 如果发现 p 不存在,那么令 link(cur)=t0,结束;否则,我们设它通过 c 转移到的状态为 q,标记为青色。比如在上图中,ABCB 对应的状态就是 q
  4. 如果发现 len(p)+1=len(q),那么令 link(cur)=q,结束。
  5. 否则我们新建一个状态 clone,复制 qlink 以及转移的信息,并且令 len(clone)=len(p)+1 ,将 clone 标记为紫色,下面是操作好的自动机和 parent tree。
  6. 然后我们令 link(cur)=link(q)=clone,下面是操作好的 parent tree。
  7. 最后我们用 parent tree 从 p 往上跳,跳到的点只要对 q 有转移,就重定向到 clone,下面是操作好的自动机。

注:易证这些有转移的点是连续的,实现时只要发现不能转移就终止。
8. 完成以上操作之后,我们将 last 的值更新为 cur,并转化为初始自动机的形式,如图所示。

这里补充一下怎么找到终止状态的:parent tree 上从 last 节点往上跳,跳的时候经过的所有状态都是终止状态,因为这样可以遍历所有 s 的后缀。

下面是代码实现:

Code
struct node{int len,link,to[30];};
struct Suffix_Automation{
	node state[N];
	int rt=1,last=1,tot=1;
	void extend(int ch){
		int cur=++tot,p,q;
		state[cur].len=state[last].len+1;
		for(p=last;p&&!state[p].to[ch];p=state[p].link)state[p].to[ch]=cur;
		if(!p)state[cur].link=rt;
		else{
			q=state[p].to[ch];
			if(state[p].len+1==state[q].len)state[cur].link=q;
			else{
				int clone=++tot;
				state[clone]=state[q];
				state[clone].len=state[p].len+1;
				state[cur].link=state[q].link=clone;
				for(;p&&state[p].to[ch]==q;p=state[p].link)state[p].to[ch]=clone;
			}
		}
		last=cur;
	}
}SAM;

4.应用

咕咕咕~

posted @   happydef  阅读(575)  评论(7编辑  收藏  举报
编辑推荐:
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
· .NET Core内存结构体系(Windows环境)底层原理浅谈
· C# 深度学习:对抗生成网络(GAN)训练头像生成模型
· .NET 适配 HarmonyOS 进展
阅读排行:
· 用 DeepSeek 给对象做个网站,她一定感动坏了
· DeepSeek+PageAssist实现本地大模型联网
· 手把手教你更优雅的享受 DeepSeek
· 腾讯元宝接入 DeepSeek R1 模型,支持深度思考 + 联网搜索,好用不卡机!
· 从 14 秒到 1 秒:MySQL DDL 性能优化实战
点击右上角即可分享
微信分享提示