本文共 3800 字,大约阅读时间需要 12 分钟。
字符串匹配是极为常见的一种模式匹配。简单地说,就是判断主串TT中是否出现该模式串PP,即PP为TT的子串。特别地,定义主串为T[0…n−1]T[0…n−1],模式串为P[0…p−1]P[0…p−1],则主串与模式串的长度各为nn与pp。
暴力匹配方法的思想非常朴素:
下图给出了暴力匹配的例子,主串T="ababcabcacbab"
,模式串P="abcac"
,第一次匹配:
第二次匹配:
第三次匹配:
C代码实现:
int brute_force_match(char *t, char *p) { int i, j, tem; int tlen = strlen(t), plen = strlen(p); for(i = 0, j = 0; i <= tlen - plen; i++, j = 0) { tem = i; while(t[tem] == p[j] & j < plen) { tem++; j++; } // matched if(j == plen) { return i; } } // [p] is not a substring of [t] return -1;}
时间复杂度:i
在主串移动次数(外层的for循环)有n−pn−p次,在失配时j
移动次数最多有p−1p−1次(最坏情况下);因此,复杂度为O(n∗p)O(n∗p)。
我们仔细观察暴力匹配方法,发现:失配后下一次匹配,
P[0]
。如此未能利用已经匹配上的字符的信息,造成了重复匹配。举个例子,比如:第一次匹配失败时,主串、模式串失配位置的字符分别为 a
与 c
,下一次匹配时主串、模式串的起始位置分别为T[1]
与P[0]
;而在模式串中c
之前是ab
,未有重复字符结构,因此T[1]
与P[0]
肯定不能匹配上,这样造成了重复匹配。直观上,下一次的匹配应从T[2]
与P[0]
开始。
根据暴力方法的缺点,而引出KMP算法的思想。首先,一般化匹配失败,如下图所示:
在暴力匹配方法中,下一次匹配开始时,主串指针会回溯到i+1
,模式串指针会回退到0
。那么,如果不让主串指针发生回溯,模式串的指针应回退到哪个位置才能保证正确匹配呢?首先,我们从上图中可以得到已匹配上的字符:
T[i…i+j−1]=P[0…j−1]T[i…i+j−1]=P[0…j−1]
KMP算法思想便是利用已经匹配上的字符信息,使得模式串的指针回退的字符位置能将主串与模式串已经匹配上的字符结构重新对齐。当有重复字符结构时,下一次匹配如下图所示:
从图中可以看出,下一次匹配开始时,主串指针在失配位置i+j
,模式串指针回退到m+1
;模式串的重复字符结构:
T[i+j−m−1…i+j−1]=P[j−m−1…j−1]=P[0…m](1)(1)T[i+j−m−1…i+j−1]=P[j−m−1…j−1]=P[0…m]
且有
T[i+j]≠P[j]≠P[m+1]T[i+j]≠P[j]≠P[m+1]
那么应如何选取mm值呢?假定有满足式子(1)的两个值m1>m2m1>m2,如下图所示:
如果选取m=m2m=m2,则会丢失m=m1m=m1的这一种字符匹配情况。由数学归纳法容易知道,应取所有满足式子(1)中最大的mm值。
KMP算法中每一次的匹配,
模式串P="abcac"
匹配主串T="ababcabcacbab"
的KMP过程如下图:
根据上面的讨论,我们定义部分匹配函数(Partial Match,在数据结构书[2]称之为失配函数):
f(j)={max{m}−1P[0…m]=P[j−m…j],0≤m<jelsef(j)={max{m}P[0…m]=P[j−m…j],0≤m<j−1else
其表示字符串P[0…j]P[0…j]的前缀与后缀完全匹配的最大长度,也表示了模式串中重复字符结构信息。KMP中大名鼎鼎的next[j]
函数表示对于模式串失配位置j+1
,下一轮匹配时模式串的起始位置(即对齐于主串的失配位置);则
next[j]=f(j)+1next[j]=f(j)+1
如何计算部分匹配函数呢?首先来看一个例子,模式串P="ababababca"
的部分匹配函数与next函数如下:
j | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | |
---|---|---|---|---|---|---|---|---|---|---|---|
P[j] | a | b | a | b | a | b | a | b | c | a | |
f(j) | -1 | -1 | 0 | 1 | 2 | 3 | 4 | 5 | -1 | 0 | |
next[j] | 0 | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 0 | 1 |
模式串的f(j)
满足P[0…f(j)]=P[j−f(j)…j]P[0…f(j)]=P[j−f(j)…j],在计算f(j+1)
分为两类情况:
f(j+1)=f(j)+1
。P[f(j+1)]=P[j+1]
的f(j+1)
,从而得到P[0…f(j+1)]=P[j+1−f(j+1)…j+1]P[0…f(j+1)]=P[j+1−f(j+1)…j+1]其中,根据f(j)
的定义有:
P[j]=P[f(j)]=P[f(f(j))]=⋯=P[fk(j)]P[j]=P[f(j)]=P[f(f(j))]=⋯=P[fk(j)]
其中,fk(j)=f(fk−1(j))fk(j)=f(fk−1(j))。通过上面的例子可知,函数fk(j)fk(j)是随着kk递减的,并最后收敛于-1
。此外,P[j]
与p[j+1]
相邻;因此若存在P[f(j+1)]=P[j+1]
,则必有
f(j+1)=fk(j)+1f(j+1)=fk(j)+1
为了求满足条件的最大的f(j+1)
,因fk(j)fk(j)是随着kk递减的,故应为满足上式的最小kk值。
综上,部分匹配函数的计算公式如下:
f(j)={fk(j−1)+1−1minkP[fk(j−1)+1]=P[j]elsef(j)={fk(j−1)+1minkP[fk(j−1)+1]=P[j]−1else
部分匹配函数(失配函数)的C实现代码:
int *fail(char *p) { int len = strlen(p); int *f = (int *) malloc(len * sizeof(int)); f[0] = -1; int i, j; for(j = 1; j < len; j++) { for(i = f[j-1]; ; i = f[i]) { if(p[j] == p[i+1]) { f[j] = i + 1; break; } else if(i == -1) { f[j] = -1; break; } } } return f;}
KMP的C实现代码:
int kmp(char *t, char *p) { int *f = fail(p); int i, j; for(i = 0, j = 0; i < strlen(t) && j < strlen(p); ) { if(t[i] == p[j]) { i++; j++; } else if(j == 0) i++; else j = f[j-1] + 1; } return j == strlen(p) ? i - strlen(p) : -1;}
时间复杂度:fail
函数的复杂度为O(p)O(p),kmp
函数的复杂度为O(n)O(n),所以整个KMP算法的复杂度为O(n+p)O(n+p)。
[1] dekai, .
[2] E. Horowitz, S. Sahni, S. A. Freed, 《Fundamentals of Data Structures in C》. [3] Jake Boxer, .转载地址:http://motgn.baihongyu.com/