字符串匹配之KMP---全力解析

转载

mb5ff980b461ced 2014-07-21 15:57:00

文章标签 next数组 kmp算法字符串匹配微信公众号搜索 文章分类 代码人生

字符串匹配方法大概有：BF(暴力破解法), 简化版的BM，KMP，BM，普通情况下，大家听说最多的应该就是KMP算法了。之前学习过，因为时间间隔比較大，记不太清楚了，今天上网查了下，发现写KMP的文章是不少，可是真正清楚简洁就没有了(july的文章太繁琐)，所以自己就研究了一晚上，弄清楚了kmp的计算过程，也就在此分享下。

KMP算法最难理解的是就是next数组的计算过程，在此分享下我所理解的kmp算法以及next数组的计算过程(假设看前面理论比較头大，能够先看后面样例的计算过程，在回过头来看理论就会释怀)：

1. next数组的计算过程:

申明:next数组下标从0算起, 定义next[0]=-1, next[1]=0; 模式串记为T[ ]

假如求 T中 j+1 位的next[j+1]：

将其前一位(模式字符)的内容与其前一位的next值(next[j])的内容(T[next[j]])进行比較：

假设它们相等(T[j]==T[next[j]])，则next[j+1] = next[j]+1;

假设他们不相等，则继续向前寻找，直到找到next值相应的内容与前一位相等为止，则在这个next值上加一；

假设直到第一位都没有与之相等，则next[j+1] = 0;

例: 有模式串 "abaababc"

j=0时，next[0] = -1 ; j=1时，next[1] = 0;

j=2时，t1!=t0, k=next[0]=-1, next[2]=0;

j=3时，t2==t0, next[3] = next[2]+1 = 1;

j=4时，k=next[3]=1, t3!=T[1], k=next[1]=0, T[3]==T[0], next[4]=next[1]+1 = 1;

j=5时，k=next[4]=1, T[4]==T[1], next[5]=next[4]+1=2;

j=6时，k=next[5]=2, T[5]==T[2], next[6]=next[5]+1=3;

j=7时，k=next[6]=3, T[6]!=T[3], k=next[3]=1, T[6]==T[1], next[7]=next[3]+1 = 2;

字符串匹配之KMP---全力解析_字符串匹配

2. 上述算法的实现：

//update 2014-04-19 10:08
       void calNext(const char *T, int *next){
          int n = strlen(T);
          if(n<=0) return ;
          next[0] = -1;
          next[1] = 0;
          int j=0, k=-1;
          while(j<n){
               if(k==-1 || T[j]==T[k]){
                    ++j;
                    ++k;
                    next[j] = k;
               }
               else  k = next[k];
          }
     }

3. KMP主算法：

设置比較起始下标: i=0, j=0;

循环直到 i+m>n 或者 T中全部字符都以比較完成

a. 假设 S[i]==T[j], 则继续比較S和T的下一个字符; 否则

b. 将 j=next[j]，从这位置開始继续进行比較;

c. 假设j==-1, 则将 i 和 j 分别加1, 继续比較；

假设T中全部字符均比較完成，则返回匹配的起始下标，否则返回-1；

4. KMP算法实现：

//update 2014-04-19 10:08
       int kmpmatch(const char *S, const char *T){
          if(S==NULL || T==NULL) return -1;
          int n = strlen(S);
          int m = strlen(T);
          int next[m];
          calNext(T, next);
          int i=0, j=0;
          while(i+m<=n){
               for( ; j<m&&i<n&&S[i]==T[j]; ++i, ++j) ;
               if(j==m) return i-m;
               j = next[j];
               if(j==-1){
                    ++i;
                    ++j;
               }
          }
          return -1;
     }