【模板】kmp以及扩展kmp

kmp算法：
最佳博客： http://www.cnblogs.com/dolphin0520/archive/2011/08/24/2151846.html  ；
这里面kmp算法讲的比较详细，我主要在这写下自己的kmp模板以及需要注意的地方。

void getnext(int t)
{
int i=0,j=-1;
while(i!=t)
   {
    if(j==-1||a[i]==a[j])
    {
     i++；
     j++;
     next[i]=j;
    }
    else j=next[j];
   }
} 

这是kmp算法中求next数组的算法，需要注意的是这个next是指的下一位的比较，而不是本位的比较。
例如：       aabaac中，c的next位是2，也就是说若在c处比较不相等，那么直接回溯到b处再进行比较
                 012345
                 -101012

扩展kmp：
可以在网上看看刘雅琼的ppt。
扩展kmp是求两个数组最大相同前缀数目的节约时间的算法，其时间是线性的，也就是a+b类的，而不是a*b。
扩展kmp模板：

#include<stdio.h>
#include<string.h>
char a[110010],b[110010];
int next[110010],ex[110010];
void getnext(int t)                //这个子函数是求扩展kmp中模板数组的next数组，和kmp算法中的next数组完全不同。
{                                          //扩展kmp算法中next数组是指从这一位开始最大的与全字符串相同的前缀。
 int i,k,L,j,p;
 next[0]=t;
 for(i=0;i<t;i++)
  if(b[i]!=b[i+1]) break;         //在getnext函数中，i是循环节的标记符，p是最大相同前缀的最长长度， 
 next[1]=i;                            // 例如：acabcaca，当循环到第二个a时，p为3，当循环到第三个a时，p为8.
 k=1;p=k+next[k]-1;
 for(i=2;i<t;i++)
 {
  L=next[i-k];
  if(i+L<=p) next[i]=L;
  else
  {
   j=p-i+1;
   if(j<0) j=0;
   while(i+j<t&&b[i+j]==b[j]) j++;    //注意添加限制条件
   next[i]=j;k=i;p=k+next[k]-1;
  }
 }
}
void getex(int t)                            //ex数组是指a串与b串最大相同前缀。
{
 int i,k,p,L,j;
 for(i=0;i<t;i++)
  if(a[i]!=b[i]) break;
 k=0;ex[0]=i;
 for(i=1;i<t;i++)
 {
  p=k+ex[k]-1;
  if(p<0) p=0;
  L=next[i-k];
  if(i+L<=p) next[i]=L;
  else
  {
   j=p-i+1;
   if(j<0) j=0;
   while(i+j<t&&a[i+j]==b[j]) j++;
   ex[i]=j;k=j;
  }
 }
}
int main()
{
 int len,i,j,max;
 while(scanf("%s",a)!=EOF)
 {
  getchar();
  len=strlen(a);
  j=len-1;
  max=0;
  for(i=0;i<len;i++)
  {
   b[j]=a[i];
   j--;
  }
  for(i=0;i<len;i++)
   printf("%c ",b[i]);
  printf("\n");
  getnext(len);
  for(i=0;i<len;i++)
   printf("%d ",next[i]);
  printf("\n");
  getex(len);
  for(i=0;i<len;i++)
   printf("%d ",ex[i]);
  printf("\n");
  /*for(i=0;i<len;i++)
   if(max<ex[i]) max=ex[i];
  printf("%d\n",max);*/
 }
}