AhoCorasickAutomata AC自动机

问题

在文本texttext中查找kk个模式patternpattern出现的所有位置。其中texttext长度为nnpatternipattern_{i}的长度为mim_{i},其中最长的模式长度为mmaxm_{max}n>mmaxn \gt m_{max},所有模式长度之和为msum=i=1kmim_{sum} = \sum_{i=1}^{k} m_{i},且所有模式两两互不重复。

简单字符串匹配算法

将简单字符串匹配SimpleMatch应用在本问题上,搜索所有模式需要重复kk次,每次的时间复杂度为O(nmi)O(n \cdot m_i)。直接应用SimpleMatch算法的时间复杂度为O(nmsum)O(n \cdot m_sum)

前缀树

能否在一次匹配texttext的过程中就同时找出所有模式呢?即并行算法(算法的并行,而非多线程/多进程的并行)。

首先用所有patternpattern构造一个前缀树ptpt,如图所示:

(1)(1)texttext的首个字符text[i=0]text[i = 0]开始,将其与前缀树中的节点依次向下匹配,可知text[01]=pt[25]text[0 \dots 1] = pt[2 \dots 5],但text[2]pt[8]text[2] \ne pt[8]因此匹配失败;

(2)(2) 匹配位置向右移动一位,从text[i=1]text[i = 1]开始,可知text[1]=pt[1]text[1] = pt[1]text[2]pt[3],text[2]pt[4]text[2] \ne pt[3], text[2] \ne pt[4]

(3)(3) 匹配位置向右移动一位,从text[i=2]text[i = 2]开始匹配,可知text[23]=pt[13],text[25]=pt[19]text[2 \dots 3] = pt[1 \dots 3], text[2 \dots 5] = pt[1 \dots 9]匹配成功;

\cdots

构建前缀树的时间复杂度为O(msum)O(m_{sum})。利用前缀树,文本上的每个字符匹配前缀树即可,该算法的时间复杂度为O(msum)+O(nmmax)O(m_{sum}) + O(n \cdot m_{max})。当nn远大于mmaxm_{max}时显然构造第二种算法更优。

失败指针

下图中,当匹配到text[i=5]text[i = 5]时有text[56]=pt[25]text[5 \dots 6] = pt[2 \dots 5],但text[7]pt[8]text[7] \ne pt[8]匹配失败。我们不希望从text[i=6]text[i = 6]处从前缀树的根节点重新开始匹配,显然text[i=6]text[i = 6]在前缀树中已经存在。因为pt[11]pt[1 \dots 1]pt[25]pt[2 \dots 5]的后缀字符串,这时将前缀树的匹配位置调整到pt[1]pt[1],那么pt[1]pt[1]text[i=6]text[i = 6]可以继续匹配,尝试找到一个成功的匹配。图中红色的连线称为失败链接/失败指针failurelinkfailure link

设字符串α\alpha的末尾字符为pt[i]pt[i],尝试在前缀树中寻找α\alpha的最长后缀字符串β\beta(设pt[j]pt[j]β\beta的末尾字符)。若找到这样一个合适的β\beta,建立从pt[i]pt[i]pt[j]pt[j]的指针,否则建立从pt[i]pt[i]pt[root]pt[root]的指针。显然前缀树中每个节点只有一个失败指针。失败指针的出发节点是前缀树中最后一个成功匹配的字符,其实质是后缀字符串,也称后缀链接/后缀指针suffixlinksuffix link

失败指针的核心思路在于匹配文本失败时,希望避免从前缀树的根部重新开始匹配。失败指针要么指向一个与当前位置上字符串相同的最长的后缀字符串(这样的指针就是后缀指针),要么指向前缀树的根节点。比如下图中pt[25]="sh"pt[2 \dots 5] = "sh"的最长后缀字符串是pt[11]="h"pt[1 \dots 1] = "h"pt[28]="she"pt[2 \dots 8] = "she"的最长后缀字符串是pt[13]="he"pt[1 \dots 3] = "he"pt[14]="hi"pt[1 \dots 4] = "hi"找不到最长后缀字符串(也可以认为最长后缀字符串为空),因此有失败指针pt[4]pt[root]pt[4] \rightarrow pt[root]

前缀树中的失败指针联系的两个节点可能在同一个字符串上,比如下图中有失败指针pt[11]pt[5]pt[11] \rightarrow pt[5]pt[10]pt[2]pt[10] \rightarrow pt[2],这两个失败指针在前缀树中构成了环形图。

输出指针

下图中,当匹配到text[10]=pt[8]text[10] = pt[8]时(即使在该位置没有text[810]=pt[28]text[8 \dots 10] = pt[2 \dots 8]匹配成功也同样适用),我们发现不管前缀树当前位置pt[8]pt[8]匹配成功与否,一定存在成功的匹配pt[13]pt[1 \dots 3]。利用这一特性避免了从text[i=9]text[i = 9]和前缀树的根节点重新开始匹配。显然这也是失败指针,但并非在匹配失败时才跳转,这类跳转称为输出指针/输出链接outputlinkoutput link,用红色虚线表示。

再给一个特别情况,当匹配到text[02]=pt[27]text[0 \dots 2] = pt[2 \dots 7]时,有失败指针pt[7]pt[6]pt[7] \rightarrow pt[6],输出指针pt[6]pt[1]pt[6] \rightarrow pt[1]。因此对于前缀树中的节点pt[7]pt[7],需要递归的沿着所有失败指针,找出一次成功匹配text[22]=pt[11]text[2 \dots 2] = pt[1 \dots 1]。当匹配到text[03]=pt[29]text[0 \dots 3] = pt[2 \dots 9]时,有输出指针pt[9]pt[8],pt[8]pt[4]pt[9] \rightarrow pt[8], pt[8] \rightarrow pt[4]。如图所示:

仔细观察可以发现,输出指针pt[i]pt[j]pt[i] \rightarrow pt[j]有几个特性:

(1)(1) 两个节点不在同一个字符串分支上,pt[i]pt[i]是前缀树中的任意节点;

(2)(2) 输出指针是一种特殊的失败指针,pt[j]pt[root]pt[j] \ne pt[root]。显然每个节点上只有至多一个输出指针;

(3)(3) pt[j]pt[j]是前缀树中的叶子节点;

在匹配过程中,尝试递归的沿着前缀树上当前节点的失败指针,找出所有输出指针,这些输出指针都是(在其他分支的字符串上的)成功匹配。

最终得到AC自动机算法:对于文本texttext上的任意字符text[i]text[i],从前缀树ptpt的根部开始匹配:

(1)(1) 沿着前缀树完成一次成功匹配,text[i]text[i]上的位置ii向右移动一位,从前缀树ptpt的根节点重新开始匹配;

(2)(2) 匹配失败时,若前缀树上的当前节点上有非pt[root]pt[root](非前缀树根节点)的失败指针,则跳到失败指针处继续匹配;若没有这样的失败指针,则文本texttext的匹配位置向右移动一位,从前缀树ptpt的根节点重新开始匹配;

(3)(3) 匹配途中若遇到输出指针,立刻找到一次输出指针所处的成功匹配,但不影响当前字符串分支上的匹配,当前的匹配仍然继续;

AC自动机的匹配时间复杂度为O(n+msum+z)O(n + m_{sum} + z)。其中zz是所有模式patternpattern在文本texttext上出现的次数。

构建AC自动机

构建AC自动机需要三步:构建前缀树;构建失败指针;构建输出指针。

构造前缀树的过程详见本书的DataStructure-PrefixTree。

构造失败指针的过程是一种类似BFS/层序遍历树的算法。初始时令根节点的失败指针指向自己,首先将前缀树的第一层节点加入空队列QueueQueue中,所有的失败指针指向根节点。然后依次从QueueQueue中取出头节点pt[i]pt[i],对于头节点的某个孩子节点pt[child]pt[child],寻找它的失败指针,并将pt[child]pt[child]推入QueueQueue中,直到QueueQueue为空:

(1)(1) 对于前缀树根节点pt[root]pt[root],其失败指针指向自己;

(2)(2) 对于前缀树第一层节点,其失败指针指向pt[root]pt[root]

(3)(3) 对于前缀树中其他的节点pt[i]pt[i],设该节点的字符为xx,其父节点为pt[father]pt[father],且pt[fail]pt[fail]pt[father]pt[father]的失败指针。若pt[fail]pt[fail]有字符为xx的孩子节点pt[child]pt[child],则显然pt[child]pt[child]所在的字符串为pt[i]pt[i]所在字符串的最长后缀字符串。因此有失败指针pt[i]pt[child]pt[i] \rightarrow pt[child]。若不存在这样的孩子节点,则递归的再考虑pt[j]pt[j]的失败指针,直到失败指针本身是pt[root]pt[root],则有失败指针pt[i]pt[root]pt[i] \rightarrow pt[root],递归结束;如图所示:

在构造失败指针的同时构造输出指针,若pt[i]pt[i]的失败指针pt[j]pt[j]不是前缀树根节点pt[root]pt[root],又是前缀树的叶子节点,则有输出指针pt[i]pt[j]pt[i] \rightarrow pt[j]。显然pt[root]pt[root]不存在输出指针,前缀树第一层节点也都不存在输出指针。

AC自动机的构造时间复杂度为O(msum)O(m_{sum}),加上匹配的时间,AC自动机算法的时间复杂度为O(n+msum+z)O(n + m_{sum} + z)

Aho Corasick Automata

源码

AhoCorasickAutomata.h

AhoCorasickAutomata.cpp

测试

AhoCorasickAutomataTest.cpp

Last updated