KMP算法

初看kmp算法的时候有点模糊，第一次就根本没明白过。
仔细的推敲。找相关类似的问题。现在把源程序贴出来供大家参考。
关键一点就是要了解next函数的构造，以及为什么要这么做。
在数据结构81页中的next推倒，不过不是很好理解。
其中next是按1开始。伪代码的写法

演示下推导的思路(按书上的从1开始的next求法)：

j             |  1  2  3  4  5  6  7   8
---------------------------------------------
模式串   |  a  b  a  a  b  c  a   c
next[j]    |  0  1  1  2  2  3  1   2
其中j=1时,是第一个字符'a' 而a前面没有匹配的字符所以next[1]=0;
  当 j=2时,是第二个字符'b',而b前面只有一个字符,'a',a又是第一个字符不能和自己匹配,所以next[2]=1;
 当j=3时,是第三个字符'a', 前面的b和第一个字符a不匹配,所以next[3]=1;
当j=4时,是第四个字符'a',它前面的字符是a,和第一个字符a是匹配的,所以+1
得到next[4]=2;
同理,j=5时,前面是a.和第一个字符还是匹配的.next[5]=2 ,//(一个字符长度+向前面移动一个单位)
j=6时,前面字符是abaab可以取到最长的子串是ab与第一个字符开始匹配,匹配成功next[6]=3 //最长子串长2+向前移动1个单位得到3
结合这,把83页get_next()函数仔细的把数带进去推一次就不会那么模糊了

 

#include < stdio.h >
#include < string .h >
#include < conio.h >
#define   MAX  250
void  get_next( const   char   * T, int   * next) // next
{
     int  i = 0 ,j =- 1 ;
    next[ 0 ]  =   - 1 ;  // next[0] = -1 跟数组下表相配套，从0开始
     while (i  <  strlen(T))
    {
         if (j  ==   - 1   ||  T[i]  ==  T[j])
        {
             ++ i;
             ++ j;
            next[i]  =  j; 
        }
         else
        {
            j = next[j];
        }
    }
}
int  Index_KMP( const   char   * S,  const   char   * T) // KMP
{
     int  i = 0 , j = 0 ;
     int  next[MAX];
    get_next(T,next);
     while (i  <  strlen(S)  &&  j  <  strlen(T))
    {
         if (j  ==   0   ||  S[i]  ==  T[j])
        {
             ++ i;
             ++ j;
        }
         else  
        {
            j  =  next[j];
        }
    }
     if (j  >=  strlen(T))  return  i - strlen(T);        
     else   return   0 ;
}
void  main()
{
     int  pos;
     char  str1[MAX]; // 主串
     char  str2[MAX]; // 模式串
    printf( " input first string: " );
    gets(str1);
    printf( " input second string: " );
    gets(str2);
    pos = Index_KMP(str1,str2);
    printf( " 匹配位置在主串第%d个字符之后的位置 " ,pos);
    getch();
}