扩展KMP(记录学习过程)

首先看了几篇博客，发现还得耐下心来理解，动手画图 ，才能理解俩字符串还有数组之间神奇的关系

这篇博客图不错：https://segmentfault.com/a/1190000008663857

/*扩展kmp */
/*一开始假设求出了 next数组，然后去求extend数组，然后方法一样
next数组可以自己求自己*/
#include<stdio.h>
#include<string.h>
#define S 1000
#define T 1000
char s[S],t[T];
/* 母串  子串*/
int extend[S],nex[T],p,a;/*extend[a]=p  p最大值 */
/* extend[i]数组 ：s母串的 i开头后缀串 与 t子串的最长公共前缀
   nex[i] 数组：t子串以 i开头后缀串 与 t子串的最长公共前缀
*/
void getnext()
{
    a=0,p=0;
    int m=strlen(t);
    memset(nex,0,sizeof(nex));
    nex[0]=m;/*自己本身作为前后缀*/
    for(int i=1;i<m;i++)
    {
        if(i>=p||i+nex[i-a]>=p)
        {
            if(i>=p)
                p=i;
            while(p<m&&t[p]==t[p-i])
                p++;
            nex[i]=p-i;
            a=i;
        }
        else
            nex[i]=nex[i-a];
    }
}
void getextend()
{
    getnext();
    a=0,p=0;
    int n=strlen(s),m=strlen(t);
    for(int i=0;i<n;i++)
    {
        if(i>=p||i+nex[i-a]>=p)/*i-a的含义，多看图*/
        {
            if(i>=p)
                p=i;
            while(p<n&&p-i<m&&s[p]==t[p-i])/*可以继续往后匹配*/
                p++;
            extend[i]=p-i;
            a=i;/*记录下标*/
        }
        else
            extend[i]=p-i;
    }
//    for(int i=0;i<n;i++)
//        printf("%d ",extend[i]);
//    printf("\n");
//    for(int i=0;i<m;i++)
//        printf("%d ",nex[i]);
//    printf("\n");
    return ;
}
int main()
{
    while(~scanf("%s%s",s,t))
    {
        getextend();
    }
    return 0;
}