kmp算法大全，历史详解

三种next数组推导：

(1)Poj-2752       前边字串等于最后字串

#include <stdio.h>

#include <string.h>

#include<iostream>

using namespace std;

int next[500000],a[500000];

char p[500000];

int i,j,k,len,t;

void naxt(int *next)

{

    int i,j,k;

    next[0]=-1;

    j=0;

    k=-1;

 

    while(j<len)

    {

        if(k==-1||p[j]==p[k])

        {

            j++;

            k++;

            next[j]=k;

        }

        else

           k=next[k];

    }

}

int main()

{

    while(scanf("%s",p)!=EOF)

    {

        memset(next,0,sizeof(next));

        len=strlen(p);

        naxt(next);

        for(i=0; len!=0; i++)

        {

            a[i]=len;

            len=next[len];

        }

        i--;

        for(; i>=0; i--)

            printf("%d ",a[i]);

        printf("\n");

    }

    return 0;

}

(2)Poj-3461       字串在母串出现次数

#include <iostream>

#include<cstdio>

#include<algorithm>

#include<cstring>

#include<cmath>

using namespace std;

char T[1000000+100],P[10000+100];

int f[10000+100];

int cnt;//统计出现的次数

void find(char *T,char *P,int *f)

{

    int n=strlen(T);

    int m=strlen(P);

    int j=0;

    for(int i=0; i<n; i++)

    {

        if(j && T[i]!=P[j])

            j=f[j];

        if(T[i]==P[j])

            j++;

        if(j==m)

            cnt++;

    }

}

void getFail(char *P,int *f)

{

    int m =strlen(P);

    f[0]=f[1]=0;

    for(int i=1; i<m; i++)

    {

        int j=f[i];

        while(j && P[i]!=P[j])

            j=f[j];

        if(P[i]==P[j])

 

        f[i+1] = j+1;

        else

            f[i+1]= 0;

    }

}

int main()

{

    int K;

    scanf("%d",&K);

    while(K--)

    {

        cnt=0;

        scanf("%s %s",P,T);

        getFail(P,f);

        find(T,P,f);

        printf("%d\n",cnt);

    }

    return 0;

}

(3)HDU-1711      子序列在母序列位置

#include<iostream>

#include<stdlib.h>

#include<stdio.h>

#include<string.h>

using namespace std;

const int p=1000005;

int ne[p],a[p];

int b[p],u,n,g,m;

void naxt()

{

    ne[1]=0;

    int i=1, j=0;

    while(i<m)

    {

        if(j==0||b[i]==b[j])

        {

            ++j;

            ++i;

            if(b[i]!=b[j])

                ne[i]=j;

            else

                ne[i]=ne[j];

        }

        else

            j=ne[j];

    }

}

int kmp()

{

    int i=1;

    int j=1;

    while(i<=n&&j<=m)

    {

        if(j==0||a[i]==b[j])

        {

            ++i;

            ++j;

        }

        else

        {

            j=ne[j];

        }

    }

    if(j>m)

        return i-m;

    return 0 ;

}

int main()

{

   scanf("%d",&u);

        while(u--)

        {

            memset(ne,0,sizeof(ne));

            memset(a,0,sizeof(a));

            memset(b,0,sizeof(b));

            scanf("%d%d",&n,&m);

            for(int i=1; i<=n; i++)

                scanf("%d",&a[i]);

            for(int j=1; j<=m; j++)

                scanf("%d",&b[j]);

            naxt();

            int o;

            o=kmp();

            if(o==0)

                o=-1;

            printf("%d\n",o);

        }

    return 0;

}