本文主要是介绍看到一篇博客写的kmp算法,由于才疏学浅,难以理解,干脆自己实现个kmp算法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
KMP字符串模式匹配详解
来自CSDN A_B_C_ABC 网友
KMP字符串模式匹配通俗点说就是一种在一个字符串中定位另一个串的高效算法。简单匹配算法的时间复杂度为O(m*n);KMP匹配算法。可以证明它的时间复杂度为O(m+n).。
一. 简单匹配算法
先来看一个简单匹配算法的函数:
int Index_BF ( char S [ ], char T [ ], int pos )
{
/* 若串 S 中从第pos(S 的下标0≤pos<StrLength(S))个字符
起存在和串 T 相同的子串,则称匹配成功,返回第一个
这样的子串在串 S 中的下标,否则返回 -1 */
int i = pos, j = 0;
while ( S[i+j] != '\0'&& T[j] != '\0')
if ( S[i+j] == T[j] )
j ++; // 继续比较后一字符
else
{
i ++; j = 0;// 重新开始新的一轮匹配
}
if ( T[j] == '\0')
return i; // 匹配成功 返回下标
else
return -1; // 串S中(第pos个字符起)不存在和串T相同的子串
} // Index_BF
此算法的思想是直截了当的:将主串S中某个位置i起始的子串和模式串T相比较。即从 j=0 起比较 S[i+j] 与 T[j],若相等,则在主串 S 中存在以 i 为起始位置匹配成功的可能性,继续往后比较( j逐步增1 ),直至与T串中最后一个字符相等为止,否则改从S串的下一个字符起重新开始进行下一轮的"匹配",即将串T向后滑动一位,即 i 增1,而 j 退回至0,重新开始新一轮的匹配。
例如:在串S=”abcabcabdabba”中查找T=” abcabd”(我们可以假设从下标0开始):先是比较S[0]和T[0]是否相等,然后比较S[1] 和T[1]是否相等…我们发现一直比较到S[5] 和T[5]才不等。如图:
当这样一个失配发生时,T下标必须回溯到开始,S下标回溯的长度与T相同,然后S下标增1,然后再次比较。如图:
这次立刻发生了失配,T下标又回溯到开始,S下标增1,然后再次比较。如图:
这次立刻发生了失配,T下标又回溯到开始,S下标增1,然后再次比较。如图:
又一次发生了失配,所以T下标又回溯到开始,S下标增1,然后再次比较。这次T中的所有字符都和S中相应的字符匹配了。函数返回T在S中的起始下标3。如图:
二. KMP匹配算法
还是相同的例子,在S=”abcabcabdabba”中查找T=”abcabd”,如果使用KMP匹配算法,当第一次搜索到S[5]和T[5]不等后,S下标不是回溯到1,T下标也不是回溯到开始,而是根据T中T[5]==’d’的模式函数值(next[5]=2,为什么?后面讲),直接比较S[5]和T[2]是否相等,因为相等,S和T的下标同时增加;因为又相等,S和T的下标又同时增加。。。最终在S中找到了T。如图:
KMP匹配算法和简单匹配算法效率比较,一个极端的例子是:
在S=“AAAAAA…AAB“(100个A)中查找T=”AAAAAAAAAB”,简单匹配算法每次都是比较到T的结尾,发现字符不同,然后T的下标回溯到开始,S的下标也要回溯相同长度后增1,继续比较。如果使用KMP匹配算法,就不必回溯.
对于一般文稿中串的匹配,简单匹配算法的时间复杂度可降为O (m+n),因此在多数的实际应用场合下被应用。
KMP算法的核心思想是利用已经得到的部分匹配信息来进行后面的匹配过程。看前面的例子。为什么T[5]==’d’的模式函数值等于2(next[5]=2),其实这个2表示T[5]==’d’的前面有2个字符和开始的两个字符相同,且T[5]==’d’不等于开始的两个字符之后的第三个字符(T[2]=’c’).如图:
也就是说,如果开始的两个字符之后的第三个字符也为’d’,那么,尽管T[5]==’d’的前面有2个字符和开始的两个字符相同,T[5]==’d’的模式函数值也不为2,而是为0。
前面我说:在S=”abcabcabdabba”中查找T=”abcabd”,如果使用KMP匹配算法,当第一次搜索到S[5]和T[5]不等后,S下标不是回溯到1,T下标也不是回溯到开始,而是根据T中T[5]==’d’的模式函数值,直接比较S[5]和T[2]是否相等。。。为什么可以这样?
刚才我又说:“(next[5]=2),其实这个2表示T[5]==’d’的前面有2个字符和开始的两个字符相同”。请看图 :因为,S[4] ==T[4],S[3] ==T[3],根据next[5]=2,有T[3]==T[0],T[4] ==T[1],所以S[3]==T[0],S[4] ==T[1](两对相当于间接比较过了),因此,接下来比较S[5]和T[2]是否相等。。。
有人可能会问:S[3]和T[0],S[4]和T[1]是根据next[5]=2间接比较相等,那S[1]和T[0],S[2]和T[0]之间又是怎么跳过,可以不比较呢?因为S[0]=T[0],S[1]=T[1],S[2]=T[2],而T[0] != T[1], T[1] != T[2],==> S[0] != S[1],S[1] != S[2],所以S[1] != T[0],S[2] != T[0]. 还是从理论上间接比较了。
有人疑问又来了,你分析的是不是特殊轻况啊。
假设S不变,在S中搜索T=“abaabd”呢?答:这种情况,当比较到S[2]和T[2]时,发现不等,就去看next[2]的值,next[2]=-1,意思是S[2]已经和T[0]间接比较过了,不相等,接下来去比较S[3]和T[0]吧。
假设S不变,在S中搜索T=“abbabd”呢?答:这种情况当比较到S[2]和T[2]时,发现不等,就去看next[2]的值,next[2]=0,意思是S[2]已经和T[2]比较过了,不相等,接下来去比较S[2]和T[0]吧。
假设S=”abaabcabdabba”在S中搜索T=“abaabd”呢?答:这种情况当比较到S[5]和T[5]时,发现不等,就去看next[5]的值,next[5]=2,意思是前面的比较过了,其中,S[5]的前面有两个字符和T的开始两个相等,接下来去比较S[5]和T[2]吧。
总之,有了串的next值,一切搞定。那么,怎么求串的模式函数值next[n]呢?(本文中next值、模式函数值、模式值是一个意思。)
三.怎么求串的模式值next[n]
定义:
(1)next[0]= -1 意义:任何串的第一个字符的模式值规定为-1。
(2)next[j]= -1 意义:模式串T中下标为j的字符,如果与首字符
相同,且j的前面的1—k个字符与开头的1—k
个字符不等(或者相等但T[k]==T[j])(1≤k<j)。
如:T=”abCabCad”则 next[6]=-1,因T[3]=T[6]
(3)next[j]=k 意义:模式串T中下标为j的字符,如果j的前面k个
字符与开头的k个字符相等,且T[j] != T[k] (1≤k<j)。
即T[0]T[1]T[2]。。。T[k-1]==
T[j-k]T[j-k+1]T[j-k+2]…T[j-1]
且T[j] != T[k].(1≤k<j);
(4) next[j]=0 意义:除(1)(2)(3)的其他情况。
举例:
01)求T=“abcac”的模式函数的值。
next[0]= -1 根据(1)
next[1]=0 根据 (4) 因(3)有1<=k<j;不能说,j=1,T[j-1]==T[0]
next[2]=0 根据 (4) 因(3)有1<=k<j;(T[0]=a)!=(T[1]=b)
next[3]= -1 根据 (2)
next[4]=1 根据 (3) T[0]=T[3]且 T[1]=T[4]
即
下标 | 0 | 1 | 2 | 3 | 4 |
T | a | b | c | a | c |
next | -1 | 0 | 0 | -1 | 1 |
若T=“abcab”将是这样:
下标 | 0 | 1 | 2 | 3 | 4 |
T | a | b | c | a | b |
next | -1 | 0 | 0 | -1 | 0 |
为什么T[0]==T[3],还会有next[4]=0呢,因为T[1]==T[4],根据 (3)”且T[j] != T[k]”被划入(4)。
02)来个复杂点的,求T=”ababcaabc”的模式函数的值。
next[0]= -1 根据(1)
next[1]=0 根据(4)
next[2]=-1 根据 (2)
next[3]=0 根据 (3)虽T[0]=T[2]但T[1]=T[3]被划入(4)
next[4]=2 根据 (3) T[0]T[1]=T[2]T[3]且T[2] !=T[4]
next[5]=-1 根据 (2)
next[6]=1 根据 (3) T[0]=T[5]且T[1]!=T[6]
next[7]=0 根据 (3)虽T[0]=T[6]但T[1]=T[7]被划入(4)
next[8]=2 根据 (3) T[0]T[1]=T[6]T[7]且T[2] !=T[8]
即
下标 | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
T | a | b | a | b | c | a | a | b | c |
next | -1 | 0 | -1 | 0 | 2 | -1 | 1 | 0 | 2 |
只要理解了next[3]=0,而不是=1,next[6]=1,而不是= -1,next[8]=2,而不是= 0,其他的好象都容易理解。
03) 来个特殊的,求 T=”abCabCad”的模式函数的值。
下标 | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
T | a | b | C | a | b | C | a | d |
next | -1 | 0 | 0 | -1 | 0 | 0 | -1 | 4 |
next[5]= 0 根据 (3)虽T[0]T[1]=T[3]T[4],但T[2]==T[5]
next[6]= -1 根据 (2)虽前面有abC=abC,但T[3]==T[6]
next[7]=4 根据 (3)前面有abCa=abCa,且 T[4]!=T[7]
若T[4]==T[7],即T=” adCadCad”,那么将是这样:next[7]=0,而不是= 4,因为T[4]==T[7].
下标 | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
T | a | d | C | a | d | C | a | d |
next | -1 | 0 | 0 | -1 | 0 | 0 | -1 | 0 |
如果你觉得有点懂了,那么
练习:求T=”AAAAAAAAAAB”的模式函数值,并用后面的求模式函数值函数验证。
意义:
next函数值究竟是什么含义,前面说过一些,这里总结。
设在字符串S中查找模式串T,若S[m]!=T[n],那么,取T[n]的模式函数值next[n],
1. next[n]= -1表示S[m]和T[0]间接比较过了,不相等,下一次比较 S[m+1]和T[0]
2. next[n]=0表示比较过程中产生了不相等,下一次比较 S[m]和T[0]。
3. next[n]= k >0但k<n,表示,S[m]的前k个字符与T中的开始k个字符已经间接比较相等了,下一次比较S[m]和T[k]相等吗?
4. 其他值,不可能。
楼主根据以上讲解,实现了找出next数组方法,和字符串匹配方法
求next【j】数组伪算法:
next[j]=?;如何求出next数组下标为j的值呢?,且听我慢慢分析。。。
1. 首先求出s数组从0到J的前缀和后缀有多少个相同的。此函数为
the same=int GetSame(string&s,int j);
举例:s="abcacd"
a | b | c | a | c | |
j | 0 | 1 | 2 | 3 | 4 |
same | 0 | 0 | 0 | 0 | 1 |
if 前缀和后缀相同
if s[i]!=s[same]
next[j]=same;
else
if s[i]==s[0]
next[j]=-1;
else
next[j]
endif
endif
if 前缀和后缀不相同
if s[j]==s[0]
next[j]=-1;
else
next[j]=0;
#include "iostream"
using namespace std;
int GetSame(string&s,int j);
void GetNextArry(int *arry,string &s);
bool IsMatch(string&t,string&s,int*next);
int main()
{
string s,t;
t="dfafsfdabcacdfj";
s="abcacdf";
/* GetSame(s,7);*/
int next[20];
GetNextArry(next,s);
if(IsMatch(t,s,next))
cout<<"YES"<<endl;
else
{
cout<<"NO"<<endl;
}
return 0;
}
void GetNextArry(int *arry,string &s)
{
int lens=s.size();
int i=1,j=0;
arry[0]=-1;
while (i<lens)
{
int same;
same=GetSame(s,i);
if(same)
{
if(s[i]==s[same])
{
if(s[i]==s[0])
arry[i]=-1;
else
{
arry[i]=0;
}
}
else
{
arry[i]=same;
}
}
else
{
if(s[i]==s[0])
{
arry[i]=-1;
}
else
{
arry[i]=0;
}
}
i++;
}
}
int GetSame(string&s,int j)
{
if(j==1)
return 0; //一个字符串的 下标j==1 直接返回
int index; // j为偶数时匹配最大个数j/2,j为基数时匹配最大个数为(j+1)/2
if(j%2==0)
index=j/2;
else
{
index=(j+1)/2;
}
int max=0;
string L,R;
for(int i=0;i<index;i++)
{
L=s.substr(0,i+1); //前缀
R=s.substr(j-i-1,i+1); //后缀
/* cout<<L.c_str()<<" "<<R.c_str()<<endl;*/
if(L.compare(R)==0)
max=i+1;
}
/* cout<<max<<endl;*/
return max;
}
bool IsMatch(string&t,string&s,int*next)
{
int lent=t.size();
int lens=s.size();
if(lens>lent)
return 0;
int i=0;
int j=0;
while (i<lent&&j<lens)
{
if(t[i]==s[j])
{
i++;
j++;
}
else
{
j=next[j];
if(j==-1)
{
i++;
j=0;
}
}
}
if(j==lens)
return 1;
else
{
return 0;
}
}
这篇关于看到一篇博客写的kmp算法,由于才疏学浅,难以理解,干脆自己实现个kmp算法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!