很详尽KMP算法(厉害)ZzUuOo666

本KMP原文最初写于2年多前的2011年12月,因当时初次接触KMP,思路混乱导致写也写得混乱。所以一直想找机会重新写下KMP,但苦于一直以来对KMP的理解始终不够,故才迟迟没有修改本文。

假设现在我们面临这样一个问题:有一个文本串S,和一个模式串P,现在要查找P在S中的位置,怎么查找呢?

如果用暴力匹配的思路,并假设现在文本串S匹配到i位置,模式串P匹配到j位置,则有:

举个例子,如果给定文本串S“BBCABCDABABCDABCDABDE”,和模式串P“ABCDABD”,现在要拿模式串P去跟文本串S匹配,整个过程如下所示:

1.S[0]为B,P[0]为A,不匹配,执行第②条指令:“如果失配(即S[i]!=P[j]),令i=i-(j-1),j=0”,S[1]跟P[0]匹配,相当于模式串要往右移动一位(i=1,j=0)

2.S[1]跟P[0]还是不匹配,继续执行第②条指令:“如果失配(即S[i]!=P[j]),令i=i-(j-1),j=0”,S[2]跟P[0]匹配(i=2,j=0),从而模式串不断的向右移动一位(不断的执行“令i=i-(j-1),j=0”,i从2变到4,j一直为0)

3.直到S[4]跟P[0]匹配成功(i=4,j=0),此时按照上面的暴力匹配算法的思路,转而执行第①条指令:“如果当前字符匹配成功(即S[i]==P[j]),则i++,j++”,可得S[i]为S[5],P[j]为P[1],即接下来S[5]跟P[1]匹配(i=5,j=1)

4.S[5]跟P[1]匹配成功,继续执行第①条指令:“如果当前字符匹配成功(即S[i]==P[j]),则i++,j++”,得到S[6]跟P[2]匹配(i=6,j=2),如此进行下去

5.直到S[10]为空格字符,P[6]为字符D(i=10,j=6),因为不匹配,重新执行第②条指令:“如果失配(即S[i]!=P[j]),令i=i-(j-1),j=0”,相当于S[5]跟P[0]匹配(i=5,j=0)

6.至此,我们可以看到,如果按照暴力匹配算法的思路,尽管之前文本串和模式串已经分别匹配到了S[9]、P[5],但因为S[10]跟P[6]不匹配,所以文本串回溯到S[5],模式串回溯到P[0],从而让S[5]跟P[0]匹配。

而S[5]肯定跟P[0]失配。为什么呢?因为在之前第4步匹配中,我们已经得知S[5]=P[1]=B,而P[0]=A,即P[1]!=P[0],故S[5]必定不等于P[0],所以回溯过去必然会导致失配。那有没有一种算法,让i不往回退,只需要移动j即可呢?

答案是肯定的。这种算法就是本文的主旨KMP算法,它利用之前已经部分匹配这个有效信息,保持i不回溯,通过修改j的位置,让模式串尽量地移动到有效的位置。

比如对于字符串aba来说,它有长度为1的相同前缀后缀a;而对于字符串abab来说,它有长度为2的相同前缀后缀ab(相同前缀后缀的长度为k+1,k+1=2)。

比如对于aba来说,第3个字符a之前的字符串ab中有长度为0的相同前缀后缀,所以第3个字符a对应的next值为0;而对于abab来说,第4个字符b之前的字符串aba中有长度为1的相同前缀后缀a,所以第4个字符b对应的next值为1(相同前缀后缀的长度为k,k=1)。

下面,咱们就结合之前的《最大长度表》和上述结论,进行字符串的匹配。如果给定文本串“BBCABCDABABCDABCDABDE”,和模式串“ABCDABD”,现在要拿模式串去跟文本串匹配,如下图所示:

通过上述匹配过程可以看出,问题的关键就是寻找模式串中最大长度的相同前缀和后缀,找到了模式串中每个字符之前的前缀和后缀公共部分的最大长度后,便可基于此匹配。而这个最大长度便正是next数组要表达的含义。

由上文,我们已经知道,字符串“ABCDABD”各个前缀后缀的最大公共元素长度分别为:

而且,根据这个表可以得出下述结论

把next数组跟之前求得的最大长度表对比后,不难发现,next数组相当于“最大长度值”整体向右移动一位,然后初始值赋为-1。意识到了这一点,你会惊呼原来next数组的求解竟然如此简单:就是找最大对称长度的前缀后缀,然后整体右移一位,初值赋为-1(当然,你也可以直接计算某个字符对应的next值,就是看这个字符之前的字符串中有多大长度的相同前缀后缀)。

换言之,对于给定的模式串:ABCDABD,它的最大长度表及next数组分别如下:

根据最大长度表求出了next数组后,从而有

而后,你会发现,无论是基于《最大长度表》的匹配,还是基于next数组的匹配,两者得出来的向右移动的位数是一样的。为什么呢?因为:

所以,你可以把《最大长度表》看做是next数组的雏形,甚至就把它当做next数组也是可以的,区别不过是怎么用的问题。

接下来,咱们来写代码求下next数组。

基于之前的理解,可知计算next数组的方法可以采用递推:

举个例子,如下图,根据模式串“ABCDABD”的next数组可知失配位置的字符D对应的next值为2,代表字符D前有长度为2的相同前缀和后缀(这个相同的前缀后缀即为“AB”),失配后,模式串需要向右移动j-next[j]=6-2=4位。

向右移动4位后,模式串中的字符C继续跟文本串匹配。

对于P的前j+1个序列字符:

用代码重新计算下“ABCDABD”的next数组,以验证之前通过“最长相同前缀后缀长度值右移一位,然后初值赋为-1”得到的next数组是否正确,计算结果如下表格所示:

从上述表格可以看出,无论是之前通过“最长相同前缀后缀长度值右移一位,然后初值赋为-1”得到的next数组,还是之后通过代码递推计算求得的next数组,结果是完全一致的。

下面,我们来基于next数组进行匹配。

还是给定文本串“BBCABCDABABCDABCDABDE”,和模式串“ABCDABD”,现在要拿模式串去跟文本串匹配,如下图所示:

在正式匹配之前,让我们来再次回顾下上文2.1节所述的KMP算法的匹配流程:

匹配过程一模一样。也从侧面佐证了,next数组确实是只要将各个最大前缀后缀的公共元素的长度值右移一位,且把初值赋为-1即可。

我们已经知道,利用next数组进行匹配失配时,模式串向右移动j-next[j]位,等价于已匹配字符数-失配字符的上一位字符所对应的最大长度值。原因是:

但为何本文不直接利用next数组进行匹配呢?因为next数组不好求,而一个字符串的前缀后缀的公共元素的最大长度值很容易求。例如若给定模式串“ababa”,要你快速口算出其next数组,乍一看,每次求对应字符的next值时,还得把该字符排除之外,然后看该字符之前的字符串中有最大长度为多大的相同前缀后缀,此过程不够直接。而如果让你求其前缀后缀公共元素的最大长度,则很容易直接得出结果:00123,如下表格所示:

然后这5个数字全部整体右移一位,且初值赋为-1,即得到其next数组:-10012。

next负责把模式串向前移动,且当第j位不匹配的时候,用第next[j]位和主串匹配,就像打了张“表”。此外,next也可以看作有限状态自动机的状态,在已经读了多少字符的情况下,失配后,前面读的若干个字符是有用的。

行文至此,咱们全面了解了暴力匹配的思路、KMP算法的原理、流程、流程之间的内在逻辑联系,以及next数组的简单求解(《最大长度表》整体右移一位,然后初值赋为-1)和代码求解,最后基于《next数组》的匹配,看似洋洋洒洒,清晰透彻,但以上忽略了一个小问题。

比如,如果用之前的next数组方法求模式串“abab”的next数组,可得其next数组为-1001(0012整体右移一位,初值赋为-1),当它跟下图中的文本串去匹配的时候,发现b跟c失配,于是模式串右移j-next[j]=3-1=2位。

右移2位后,b又跟c失配。事实上,因为在上一步的匹配中,已经得知p[3]=b,与s[3]=c失配,而右移两位之后,让p[next[3]]=p[1]=b再跟s[3]匹配时,必然失配。问题出在哪呢?

问题出在不该出现p[j]=p[next[j]]。为什么呢?理由是:当p[j]!=s[i]时,下次匹配必然是p[next[j]]跟s[i]匹配,如果p[j]=p[next[j]],必然导致后一步匹配失败(因为p[j]已经跟s[i]失配,然后你还用跟p[j]等同的值p[next[j]]去跟s[i]匹配,很显然,必然失配),所以不能允许p[j]=p[next[j]]。如果出现了p[j]=p[next[j]]咋办呢?如果出现了,则需要再次递归,即令next[j]=next[next[j]]。

所以,咱们得修改下求next数组的代码。

利用优化过后的next数组求法,可知模式串“abab”的新next数组为:-10-10。可能有些读者会问:原始next数组是前缀后缀最长公共元素长度值右移一位,然后初值赋为-1而得,那么优化后的next数组如何快速心算出呢?实际上,只要求出了原始next数组,便可以根据原始next数组快速求出优化后的next数组。还是以abab为例,如下表格所示:

只要出现了p[next[j]]=p[j]的情况,则把next[j]的值再次递归。例如在求模式串“abab”的第2个a的next值时,如果是未优化的next值的话,第2个a对应的next值为0,相当于第2个a失配时,下一步匹配模式串会用p[0]处的a再次跟文本串匹配,必然失配。所以求第2个a的next值时,需要再次递归:next[2]=next[next[2]]=next[0]=-1(此后,根据优化后的新next值可知,第2个a失配时,执行“如果j=-1,或者当前字符匹配成功(即S[i]==P[j]),都令i++,j++,继续匹配下一个字符”),同理,第2个b对应的next值为0。

对于优化后的next数组可以发现一点:如果模式串的后缀跟前缀相同,那么它们的next值也是相同的,例如模式串abcabc,它的前缀后缀都是abc,其优化后的next数组为:-100-100,前缀后缀abc的next值都为-100。

然后引用下之前3.1节的KMP代码:

接下来,咱们继续拿之前的例子说明,整个匹配过程如下:

1.S[3]与P[3]匹配失败。

2.S[3]保持不变,P的下一个匹配位置是P[next[3]],而next[3]=0,所以P[next[3]]=P[0]与S[3]匹配。

3.由于上一步骤中P[0]与S[3]还是不匹配。此时i=3,j=next[0]=-1,由于满足条件j==-1,所以执行“++i,++j”,即主串指针下移一个位置,P[0]与S[4]开始匹配。最后j==pLen,跳出循环,输出结果i-j=4(即模式串第一次在文本串中出现的位置),匹配成功,算法结束。

“KMP的算法流程:

BM算法定义了两个规则:

下面举例说明BM算法。例如,给定文本串“HEREISASIMPLEEXAMPLE”,和模式串“EXAMPLE”,现要查找模式串是否在文本串中,如果存在,返回模式串在文本串中的位置。

1.首先,"文本串"与"模式串"头部对齐,从尾部开始比较。"S"与"E"不匹配。这时,"S"就被称为"坏字符"(badcharacter),即不匹配的字符,它对应着模式串的第6位。且"S"不包含在模式串"EXAMPLE"之中(相当于最右出现位置是-1),这意味着可以把模式串后移6-(-1)=7位,从而直接移到"S"的后一位。

2.依然从尾部开始比较,发现"P"与"E"不匹配,所以"P"是"坏字符"。但是,"P"包含在模式串"EXAMPLE"之中。因为“P”这个“坏字符”对应着模式串的第6位(从0开始编号),且在模式串中的最右出现位置为4,所以,将模式串后移6-4=2位,两个"P"对齐。

3.依次比较,得到“MPLE”匹配,称为"好后缀"(goodsuffix),即所有尾部匹配的字符串。注意,"MPLE"、"PLE"、"LE"、"E"都是好后缀。

4.发现“I”与“A”不匹配:“I”是坏字符。如果是根据坏字符规则,此时模式串应该后移2-(-1)=3位。问题是,有没有更优的移法?

5.更优的移法是利用好后缀规则:当字符失配时,后移位数=好后缀在模式串中的位置-好后缀在模式串中上一次出现的位置,且如果好后缀在模式串中没有再次出现,则为-1。所有的“好后缀”(MPLE、PLE、LE、E)之中,只有“E”在“EXAMPLE”的头部出现,所以后移6-0=6位。可以看出,“坏字符规则”只能移3位,“好后缀规则”可以移6位。每次后移这两个规则之中的较大值。这两个规则的移动位数,只与模式串有关,与原文本串无关。

6.继续从尾部开始比较,“P”与“E”不匹配,因此“P”是“坏字符”,根据“坏字符规则”,后移6-4=2位。因为是最后一位就失配,尚未获得好后缀。

由上可知,BM算法不仅效率高,而且构思巧妙,容易理解。

Sunday算法由DanielM.Sunday在1990年提出,它的思想跟BM算法很相似:

下面举个例子说明下Sunday算法。假定现在要在文本串"substringsearchingalgorithm"中查找模式串"search"。

substringsearchingalgorithmsearch^3.结果第一个字符就不匹配,再看文本串中参加匹配的最末位字符的下一位字符,是'r',它出现在模式串中的倒数第3位,于是把模式串向右移动3位(r到模式串末尾的距离+1=2+1=3),使两个'r'对齐,如下:substringsearchingalgorithmsearch^

4.匹配成功。

回顾整个过程,我们只移动了两次模式串就找到了匹配位置,缘于Sunday算法每一步的移动量都比较大,效率很高。完。

THE END
1.最浅显易懂的KMP算法讲解易易易苼学习资料KMP算法-超细超全讲解(上)原理篇 __阿岳__ 1.1万 65 03:51 KMP算法Next和NextVal 数组求解过程 诗岸梦行舟 1937 4 06:37 KMP 学一遍忘一遍?ACM 金牌选手用可视化直击本质,理解了内核后想忘记都难! NotOnlySuccess 3.0万 361 35:23 KMP算法-超细超全讲解(下)代码篇 __阿岳__ 9217 47 03https://www.bilibili.com/list/ml3236443846
2.KMP算法c语言源代码gubojun的技术博客KMP算法--c语言源代码 KMP算法 首先kmp算法的核心问题就是求next数组,next数组是为了得到匹配字符串中重复的位置。 假如 文本内容为abdaaeabdaaeaeaeffd,匹配字符串为aeabdaaeaea 1.根据匹配字符串p求出next 首先next[0]赋值为-1,next[1]赋值为0,循环遍历p中每个字符,如果2个字符相同则next[i+1]=next[i]https://blog.51cto.com/u_10101161/7177098
3.数据结构—kmp算法和strstr函数腾讯云开发者社区KMP原理 三、复杂度分析 四、KMP算法代码 引言 现实生活中,字符串匹配在很多的应用场景里都有着极其重要的作用,包括生物信息学、信息检索、拼写检查、语言翻译、数据压缩、网络入侵检测等等,至此诞生了很多的算法,那么我们今天就来探索这两种经典的算法。 https://cloud.tencent.com/developer/article/2103095
4.一文带你入木三分地理解字符串KMP算法以及C++实现C语言归纳:至此我们已经可以对于任意索引index,求出匹配串[1,index]区间的最大公共前后缀长度,结合上部分指针移动公式即可完成KMP算法。 4. 用c++代码实现 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 #include<bits/stdc++.h> using nahttps://www.jb51.net/article/269926.htm
5.动态规划之KMP算法详解(配代码版)PS:本文的代码参考《算法4》,原代码使用的数组名称是dfa(确定有限状态机),因为我们的公众号之前有一系列动态规划的文章,就不说这么高大上的名词了,本文还是沿用dp数组的名称。 一、KMP 算法概述 首先还是简单介绍一下 KMP 算法和暴力匹配算法的不同在哪里,难点在哪里,和动态规划有啥关系。 http://www.360doc.com/content/19/0925/14/9482_863133021.shtml
6.刷题29天贪心算法那么又要贪心了,局部最优:取candyVec[i + 1] + 1 和 candyVec[i] 最大的糖果数量,保证第i个小孩的糖果数量既大于左边的也大于右边的。全局最优:相邻的孩子中,评分高的孩子获得更多的糖果。 c++ classSolution{public:intcandy(vector<int>&ratings){vector<int>candyVec(ratings.size(),1);//1从前往后https://zhuanlan.zhihu.com/p/11941030486
7.KMP算法详解及next数组代码解释kmpnext数组代码KMP算法理解起来并不算太过于困难,从图像实例可以很直观得明晰算法原理,难点在于理解KMP算法中生成next数组的代码: voidnext(char*s,int*next){intk=-1;intj=0;next[0]=-1;while(j<strlen(s)){if(k==-1||s[k]==s[j]){k++;j++;next[j]=k;}elsej=next[j];}} https://blog.csdn.net/Amahisa/article/details/105137652