KMP算法【代码】

废话不多说，看毛片算法的核心在于next数组。

很多地方用的都是严书上的方法来求解next数组，代码确实简洁明了，但是可能对于初学者来说不容易想到，甚至不是一下子能理解。(好了，其实我说的就是自己，别笑)

以下为严师太的代码，也是很多人用的

照着代码用手调试了一遍，确实厉害

但是就跟我看别人写的递归算法时，照着推一遍，能懂，但是自己想不到，水平低，怪我。

int Next(char *p,int * next)//课本
{
    int s=0,i=1,j=0;
    char *q=p;
    while(*q!='\0')
    {
        s++;
        q++;
    }
    s--;
    next[1]=0;
    while(i<s)
    {
        if(j==0||p[i]==p[j])
        {
            i++;
            j++;
            next[i]=j;
        }
        else j=next[j];
    }
    return 0;
}

　　

这里本人贴一段自己经过一晚上思考得出来的思路，也许不完善，也许别人已经有过类似更好的，但是我只是想分享给初学者一点思路，希望共同进步。

 int Next(char *p,int * next)//p为模式串
 {
     char *q=p;
     int i,j,k,l,s=,m=;
     while(*q!='\0')
     {
         s++;
         q++;
     }
     s--;//求出模式串的长度，这里依照书上的规定从下标为1的算起
     q=p;
     next[]=s;
     next[]=;
     next[]=;
     for(i=;i<=s;i++)//所需求解的next数组元素个数
     {
         for(l=;l<=i-;l++)//每次最多进行i-2回循环，因为模式串第i个元素之前一共有i-1个元素，而i-1个元素至多有 i-2 对相等的元素
         {
             for(j=,k=j+l;j<=i--l;j++,k++)//从相等对数最多对的情况判断
             {
                 if(q[j]!=q[k]){ m=; break;}//一旦出现不相等，直接跳出循环，进行下一次判断，并且把m置为0,表示这种情况不符合
             }
             if(m) break;//经过上面的循环，若m仍然为1，则代表符合情况，不用再判断相等对数更少的情况了，直接跳出求解下一个next数组中元素的值
             m=;
         }
         next[i]=i-l;//此处是化简后的结果，可以自己演算一下。
     }
     return ;
 }
过程如下：

KMP完整代码如下：

 #include <iostream>

 using namespace std;
 #define MAX_SIZE 10

 int Next(char *p,int * next)
 {
     char *q=p;
     int i,j,k,l,s=,m=;
     while(*q!='\0')
     {
         s++;
         q++;
     }
     s--;
     q=p;
     next[]=s;
     next[]=;
     next[]=;
     for(i=;i<=s;i++)
     {
         for(l=;l<=i-;l++)
         {
             for(j=,k=j+l;j<=i--l;j++,k++)
             {
                 if(q[j]!=q[k]){ m=; break;}
             }
             if(m) break;
             m=;
         }
         next[i]=i-l;
     }
     return ;
 }

 int Index(char *S,char *T)
 {
     int next[MAX_SIZE],i=,j=,s=,t=,pos=;
     char *p=S,*q=T;
     while(*p!='\0')
     {
         s++;
         p++;
     }
     s--;
     while(*q!='\0')
     {
         t++;
         q++;
     }
     t--;
     Next(T,next);
     while(i<=s&&j<=t)
     {
         if(j==) {i++;j++;}

　　　　　　 if(i > s)
　　　　　　　 {
　　　　　　　　　 break;
　　　　　　　 }

         if(S[i]==T[j]){i++;j++;}
         else j=next[j];
     }
     if(j>t)
     {
         pos=i-t;
         cout<<"模式串在主串中第"<<pos<<"个位置匹配"<<endl;
         return ;
     }
     if(i>s&&j<=t)
     {
         cout<<"模式串无法与主串匹配!"<<endl;
         return ;
     }
 }
 int main()
 {
     char *S1="#ababcabcacbab";
     char *T1="#abcac";
     char *S2="#www.FishC.com";
     char *T2="#ww.";
     char *S3="#bbsbbs.FishC";
     char *T3="#bbsbbc";
     char *S4="#00000000000000000001";
     char *T4="#00000001";
     Index(S1,T1);
     Index(S2,T2);
     Index(S3,T3);
     Index(S4,T4);
     return ;
 }