【暑假】[实用数据结构]KMP

KMP算法

 

KMP算法是字符串匹配算法，可以在O(n)的时间完成，算法包含两部分，分别是：构造适配函数与两串匹配。

失配边的使用大大提高了算法效率，可以理解为已经成功匹配的字符不在重新匹配，因为我们已经知道它是什么，对应到算法中 匹配失败后应该在最大前缀之后继续匹配，因为某后缀已与最大前缀匹配成功而不用重新比较。

 

以下为代码实现：

 

 1 const int maxn = 1000 + 5; 2  3 void getFail(char* P,int* f){  //构造失配边  4     int n=strlen(P);     5     f[0]=f[1]=0; 6     for(int i=0;i<n;i++){  //自己和自己匹配  7         int j=f[i]; 8         while(j && P[i]!=P[j]) j=f[j]; 9         f[i+1]= P[i]==P[j]? j+1 : 0;10     }11 }12 13 //可以如是理解f数组，f[i]表示到i与后缀成功匹配的最大前缀的下一指针14 //f数组节约了算法时间，对于已知(已比较)的字符不在重新比较15  16 void KMP(char* T,char* P){17 int f[maxn];18 int n=strlen(T), m=strlen(P);19   getFail(P,f);20   int j=0;21   for(int i=0;i<n;i++){22       while(j && T[i] != P[j]) j=f[j];  //沿失配边走 23       if(T[i]==P[j]) j++;       24       if(j==m) { printf("%d\n",i-m+1); j=0; }   //成功匹配P，输出匹配点 25       //将j重调为0代表重新检查，否则P[j]调用出错 26   }27 }

 

 

有趣的是，在这里发现了作者的一个小错误，如果使字符串成功匹配所有位置且KMP算法可以完成的话，需要加入语句j=0;而在课本中这条语句是没有的。