KMP算法的代码实现
上周算法班的BEN老师花了1个小时讲自动机和KMP的关系,结果failed...明天又要上课了,花了半天时间看了下KMP,暂且停留在利用next求模式中的跳跃长度,自动机那个还不能理解。。。
具体的可以百度阮一峰的KMP算法。
看着什么前缀后缀,突然想到上下文无关文法乔姆斯基范式了。。。。又想到了NFA和正则表达式的转换,是时候复习复习了。。
太晚了,直接上代码,明天继续看ML和统计学!加油!
1 #include <iostream> 2 #include <string> 3 #include <iterator>//输出 4 5 using namespace std; 6 void Next(string str, int next[]){ //自己跟自己匹配 7 int length=str.size(); 8 next[0]=-1; 9 int j=-1; 10 for(int i=1;i<length;i++){ 11 while(j>-1 && str[j+1]!=str[i]) j=next[j]; 12 if(str[j+1]==str[i]) j++; 13 next[i]=j; 14 } 15 } 16 void Match(string str1,string str2,int next[]){ 17 Next(str2,next); 18 int length1=str1.size(),length2=str2.size(); 19 int j=-1; 20 for(int i=0;i<length1;i++){ 21 while(j>-1 && str2[j+1]!=str1[i]) j=next[j]; 22 if(str2[j+1]==str1[i]) j++; 23 if(j==length2-1){ 24 cout<<"Pattern occurs with shift "<<i-length2+1<<endl; 25 j=next[j]; 26 } 27 } 28 } 29 int main(){ 30 string str1="bbcabcdababcdabcdabde"; 31 string str2="abcdabd"; 32 int next[20]; 33 Match(str1,str2,next); 34 copy(next,next+10,ostream_iterator<int>(cout," ")); 35 cout<<endl; 36 return 0; 37 }
输出结果
MacBook-Pro:Algorithm root# g++ kmpdemo2.cpp -o kmpdemo2
MacBook-Pro:Algorithm root# ./kmpdemo2
Pattern occurs with shift 13
-1 -1 -1 -1 0 1 -1 0 0 0