标签:后缀自动机 == 定义 搜索 bsp 使用 结构 方便 else
后缀自动机 (suffix automaton, SAM) 是一个能解决许多字符串相关问题的有力的数据结构。(否则我们也不会用它)
举几个例子,以下的字符串问题都可以在线性时间内通过 SAM 解决
1.在另一个字符串中搜索一个字符串的所有出现位置。(诶?KMP好像能做)
2.计算给定的字符串中有多少个不同的子串。(诶?好像也能做)
3.统计一个子串所有不同子串的总长度 (额,线性复杂度?以前的知识好像就不行了)
4.字典序第k大/小子串(字典树好像可以做,但线性?不会!)
5.给定一个字符串 。找出字典序最小的循环移位。(这个我会!最小表示法!)
6.对于一个给定的文本串T ,有多组询问,每组询问给一个模式串S ,回答模式串S在字符串T中作为子串出现了多少次。(AC自动机!也可以)
7.给定一个字符串S 和一个特定的字符集,我们要找一个长度最短的没有在S中出现过的字符串。(以前的知识做不到我太菜了)
8.两个字符串的最长公共子串(线性??!以前的知识做不到我太菜了)
9.多个字符串的最长公共子串(我不会!!)
这么多东西都可以线性时间内解决,厉害吧;
厉害的前提是你先要学会后缀自动机
接下来进入正题:
1.SAM 是一张有向无环图。结点被称作 状态 ,边被称作状态间的 转移 。
2.图存在一个源点root称作 初始状态 ,其它各结点均可从root出发到达。每个 转移 都标有一些字母。从一个结点出发的所有转移均 不同 。
3.存在一个或多个 终止状态 。如果我们从初始状态出发,最终转移到了一个终止状态,则路径上的所有转移连接起来一定是字符串S的一个后缀。
4.S的每个后缀均可用一条从到某个终止状态的路径构成。
5.在所有满足上述条件的自动机中,SAM 的结点数是最少的。
具体可以证明,SAM的点数最多是2*n-1,边数最多时3*n-4; (暂时为了方便理解并不证明,你只要露出"哇!这是上限吗","我靠,这SAM的复杂度也太***的强了"的表情就好了);
SAM 最简单、也最重要的性质是,它包含关于字符串的所有子串的信息。任意从初始状态开始的路径,如果我们将转移路径上的标号写下来,都会形成S的一个子串 。反之每个S的子串对应从初始状态开始的某条路径。
为了更好的理解SAM(背模板,水经验),我们定义一些奇怪的东西;(不要一看到不认识的概念就跳过去,否则建图会看不懂的)
1.endpos等价类
考虑字符串S的任意非空子串t,我们记endpos(t)为在字符串S中t的所有结束位置.
让我们举个栗子:(假设对字符串中字符的编号从零开始)
对于S=“cbacbaba” endpos(ba)="2,5,7" endpos(a)="2,5,7" endpos(cba)="2,5";
我们发现,对于"ba"和"a"的endpos等价类的数值是完全相同的,那么这对与SAM来说是个好东西,我们可以将所有endpos相同的子串放到SAM的一个节点上,大大节省了空间与时间复杂度;
那么一个重要的结论出现了:SAM上的每个节点都代表一个endpos等价类;(比如刚才的栗子中,子串"ba","a"在SAM上用一个节点表示)
我们发现一个不是引理的性质:对于一个子串S,不断的删除它开头的元素,那么它的endpos一定是在已有的旧元素基础上加入一些新元素(或者不加),大致理解一下对于下面的引理很好理解
引理1:字符串S的两个非空子串x和y的endpos相同,当且仅当字符串x在S中的每次出现,都是以y后缀的形式存在的。(x<=y) (很显然的对吧?)
引理2:两个非空子串x,y, 他们的endpos的不可能存在交集 (很显然的对吧)
引理3:考虑一个endpos等价类,将类中的所有子串按长度递减的顺序排序。每个子串是它前一个子串的后缀。
换句话说,对于同一等价类的任一两子串,较短者为较长者的后缀,且该等价类中的子串长度恰好覆盖整个区间;
粗略证明感性理解:对于一个endpos等价类中的最长子串x,和最短子串y,若y的开头加一个元素,使得新形成的字符串z是S的子串,那么如果z的长度<=x,则endpos(z)
一定=endpos(y)=endpos(x); 使用数学归纳法可以证明所有的y<=z<=x;z的endpos都与x,y相同;
2.后缀链接 link:
考虑 SAM 中某个不是root的一个状态 。我们已经知道,状态v对应于具有相同endpos的等价类。我们如果定义x为这些字符串中最长的一个,则所有其它的字符串都是x的后缀。
一个后缀链接link连接到对应于状态v的最短子串的所有后缀中长度最大的那个endpos等价类。
引理1:所有后缀链接构成一棵根节点为root的树 (和造一颗随机树的原理相同:对于一个endpos等价类v,只会连到比v的最短子串小的endpos等价类)
下面放两张图:对于字符串"abcbc",前者是表示边的转移状态,后者表示link的连接方式;
好了,定义暂时告一段落了;接下来就是构造环节了(代码):(别人家的构造过程 (逃~)
!!极度建议结合之前的定义来阅读!!
代码实现就是:
class node{ public: int ch[28]; int link; int len; }SAM[3000010]; int ans[3000010]; int last=0,size=1; void SAM_build() { SAM[0].len=0; n=strlen(s+1); SAM[0].link=-1; inc(u,1,n){ int c=s[u]-‘a‘; int cur=size++; ans[cur]=1; SAM[cur].len=SAM[last].len+1; int p=last; while(p!=-1&&!SAM[p].ch[c]){ SAM[p].ch[c]=cur; p=SAM[p].link; } if(p==-1){ SAM[cur].link=0; } else{ int q=SAM[p].ch[c]; if(SAM[q].len==SAM[p].len+1){ SAM[cur].link=q; } else{ int newq=size++; SAM[newq]=SAM[q]; SAM[newq].len=SAM[p].len+1; while(p!=-1&&SAM[p].ch[c]==q){ SAM[p].ch[c]=newq; p=SAM[p].link; } SAM[q].link=newq; SAM[cur].link=newq; } } last=cur; } }
标签:后缀自动机 == 定义 搜索 bsp 使用 结构 方便 else
原文地址:https://www.cnblogs.com/kamimxr/p/12056278.html