当前位置:文档之家› K M P 算 法 详 解

K M P 算 法 详 解

KMP算法详解(转)

此前一天,一位MS的朋友邀我一起去与他讨论快速排序,红黑树,

字典树,B树、后缀树,包括KMP算法,唯独在讲解KMP算法的时候,言

语磕磕碰碰,我想,原因有二:1、博客内的东西不常回顾,忘了不少;2、

便是我对KMP算法的理解还不够彻底,自不用说讲解自如,运用自如了。

所以,特再写本篇文章。由于此前,个人已经写过关于KMP算法的两篇文

章,所以,本文名为:KMP算法之总结篇。

本文分为如下六个部分:

第一部分、再次回顾普通的BF算法与KMP算法各自的时间复杂度,并

两相对照各自的匹配原理;第二部分、通过我此前第二篇文章的引用,用

图从头到尾详细阐述KMP算法中的next数组求法,并运用求得的next数

组写出KMP算法的源码;第三部分、KMP算法的两种实现,代码实现一是

根据本人关于KMP算法的第二篇文章所写,代码实现二是根据本人的关于

KMP算法的第一篇文章所写;第四部分、测试,分别对第三部分的两种实

现中next数组的求法进行测试,挖掘其区别之所在;第五部分、KMP完整

准确源码,给出KMP算法的准确的完整源码;第六步份、一眼看出字符串

的next数组各值,通过几个例子,让读者能根据字符串本身一眼判断出

其next数组各值。 力求让此文彻底让读者洞穿此KMP算法,所有原理,

来龙去脉,让读者搞个通通透透(注意,本文中第二部分及第三部分的代

码实现一的字符串下标i

从0开始计算,其它部分如第三部分的代码实现二,第五部分,和第六部分的字符串下标i 皆是从1开始的)。

第一部分、KMP算法初解

1、普通字符串匹配BF算法与KMP算法的时间复杂度比较

KMP算法是一种线性时间复杂的字符串匹配算法,它是对BF算法

(Brute-Force,最基本的字符串匹配算法的)改进。对于给的原始串S

和模式串P,需要从字符串S中找到字符串P出现的位置的索引。

BF算法的时间复杂度O(strlen(S) * strlen(T)),空间复杂度O(1)。

KMP算法的时间复杂度O(strlen(S) + strlen(T)),空间复杂度

O(strlen(T))。

2、BF算法与KMP算法的区别

假设现在S串匹配到i位置,T串匹配到j位置。那么总的来说,两

种算法的主要区别在于失配的情况下,对的值做的处理:

BF算法中,如果当前字符匹配成功,即s[i+j] == T[j],令j++,

继续匹配下一个字符;如果失配,即S[i

+ j] != T[j],需要让i++,并且j= 0,即每次匹配失败的情况下,

模式串T相对于原始串S向右移动了一位。

而KMP算法中,如果当前字符匹配成功,即S[i]==T[j],令i++,

j++,继续匹配下一个字符;如果匹配失败,即S[i]

!= T[j],需要保持i不变,并且让j = next[j],这里next[j]

=j -1,即模式串T相对于原始串S向右移动了至少1位(移动的实际

位数j

- next[j]? =1), 同时移动之后,i之前的部分(即S[i-j+1 ~ i-1]),和j=next[j]

之前的部分(即T[0

~ j-2])仍然相等。显然,相对于BF算法来说,KMP移动更多的位数,

起到了一个加速的作用!

(失配的特殊情形,令j=next[j]导致j==0的时候,需要将i

++,否则此时没有移动模式串)。

3、BF算法为什么要回溯

首先说一下为什么BF算法要回溯。如下两字符串匹配(恰如上面所述:

BF算法中,如果当前字符匹配成功,即s[i+j]

== T[j],令j++,继续匹配下一个字符):

i+j(j随T中的j++变,而动)

S:aaaacefghij

T:aaac?

如果不回溯的话就是从下一位开始比起:

aaaacefghij

看到上面红颜色的没,如果不回溯的话,那么从a 的下一位c

比起。然而下述这种情况就漏了(正确的做法当然是要回溯:如果失

配,即S[i +

j] != T[j],需要让i++,并且j=

aaaacefghij

所以,BF算法要回溯,其代码如下:

view plain int?Index(SString?S,?SString?T,?int?pos)?{?--返回T在S中第

pos个字符之后的位

置?i=pos;?j=1;k=0;?while?(?i?=?S[0]?j?=?T[0]?)?{?if?(S[i+k]?=?=

T[j]){++k;++j;}--继续比较后续字

符?else?{i=i+1;?j=1;?k=0;}?--指针回溯到?下一首位,重新开

始?}?if(jT[0])?return?i;?--子串结束,说明匹配成

功?else?return?0;?}--Index?

不过,也有特殊情况可以不回溯,如下:abcdefghij(主

串)abcdefg(模式串)

即(模式串)没有相同的才不需要回溯。

算法思想

普通的字符串匹配算法必须要回溯。但回溯就影响了效率,回溯是

由T串本身的性质决定的,是因为T串本身有前后'部分匹配'的性质。像

上面所说如果主串为abcdef这样的,大没有回溯的必要。

改进的地方也就是这里,我们从T串本身出发,事先就找准了T自

身前后部分匹配的位置,那就可以改进算法。

如果不用回溯,那模式串下一个位置从哪里开始呢?

还是上面那个例子,T(模式串)为ababc,如果c失配,那就可以往

前移到aba最后一个a的位置,像这样:

这样i不用回溯,j跳到前2个位置,继续匹配的过程,这就是KMP

算法所在。这个当T[j]失配后,j

应该往前跳的值就是j的next值,它是由T串本身固有决定的,与S串(主串)无关。

5、next数组的含义

重点来了。下面解释一下next数组的含义,这个也是KMP算法中比较

不好理解的一点。

令原始串为: S[i],其中0=i=n;模式串为: T[j],其中0=j=m。

假设目前匹配到如下位置

S0,S1,S2.,Si-j,Si-j+1.,Si-1,

Si, Si+1.,Sn

T0,T1.,Tj-1,

S和T的绿色部分匹配成功,恰好到Si和Tj的时候失配,如果要保

持i不变,同时达到让模式串T相对于原始串S右移的话,可以更新j的

值,让Si和新的Tj进行匹配,假设新的j用next[j]表示,即让Si和

next[j]匹配,显然新的j值要小于之前的j值,模式串才会是右移的效

果,也就是说应该有next[j]

= j -1。那新的j值也就是next[j]应该是多少呢?我们观察如下的

匹配:

1)如果模式串右移1位(从简单的思考起,移动一位会怎么样),即

next[j] = j - 1, 即让蓝色的Si和Tj-1匹配

(注:省略号为未匹配部分)

S0,S1,S2.,Si-j,Si-j+1.,Si-1,

Si, Si+1.,Sn

T0,T1.,Tj-1, Tj, . (T的划线部分和S划线部分相等【1】)

T0,T1.Tj-2,Tj-1,

. (移动后的T的划线部分和S的划线部分相等【2】)

根据【1】【2】可以知道当next[j] =j -1,即模式串右移一位的时

候,有T[0

~ j-2] == T[1 ~ j-1],而这两部分恰好是字符串T[0 ~j-1]的前缀

和后缀,也就是说next[j]的值取决于模式串T中j前面部分的前缀和后

缀相等部分的长度(好好揣摩这两个关键字概念:前缀、后缀,或者再想

想,我的上一篇文章,从Trie树谈到后缀树中,后缀树的概念)。

2)如果模式串右移2位,即next[j] = j - 2, 即让蓝色的Si和

Tj-2匹配?

S0,S1.,Si-j,Si-j+1,Si-j+2.,Si-1,

Si, Si+1.,Sn

T0,T1,T2.,Tj-1,

Tj, .(T的划线部分和S划线部分相等【3】)

T0,T1.,Tj-3,Tj-2.(移动后的T的划线部分和S的划线部分相等【4】)

同样根据【3】【4】可以知道当next[j] =j -2,即模式串右移两位

的时候,有T[0

~ j-3] == T[2 ~ j-1]。而这两部分也恰好是字符串T[0 ~j-1]的前

缀和后缀,也就是说next[j]的值取决于模式串T中j前面部分的前缀和

后缀相等部分的长度。

3)依次类推,可以得到如下结论:当发生失配的情况下,j的新值next[j]取决于模式串中T[0

~ j-1]中前缀和后缀相等部分的长度, 并且next[j]恰好等于这个最

大长度。

!= T[j],需要保持i不变,并且让j

= next[j],这里next[j] =j -1,即模式串T相对于原始串S向右移

动了至少1位(移动的实际位数j

- next[j]? =1),

同时移动之后,i之前的部分(即S[i-j+1 ~ i-1]),和j=next[j]

之前的部分(即T[0 ~ j-2])仍然相等。显然,相对于BF算法来说,KMP

移动更多的位数,起到了一个加速的作用!

(失配的特殊情形,令j=next[j]导致j==0的时候,需要将i ++,否

则此时没有移动模式串)。”

于此,也就不难理解了我的关于KMP算法的第二篇文章之中:“当

匹配到S[i] != P[j]的时候有 S[i-j„i-1] = P[0„j-1]. 如果下面用

j_next去匹配,则有P[0„j_next-1] = S[i-j_next„i-1]

= P[j-j_next„j-1]。此过程如下图3-1所示。

当匹配到S[i] != P[j]时,S[i-j„i-1] = P[0„j-1]:

S: 0 „ i-j „ i-1 i „

P:? 0 „? j-1 j „

如果下面用j_next去匹配,则有P[0„j_next-1] =

S[i-j_next„i-1] = P[j-j_next„j-1]。

所以在P中有如下匹配关系(获得这个匹配关系的意义是用来求next

相关主题