KMP算法

KMP算法

KMP算法是一种改进的字符串匹配算法，由D.E.Knuth，J.H.Morris和V.R.Pratt提出的，因此人们称它为克努特—莫里斯—普拉特操作（简称KMP算法）。KMP算法的核心是利用匹配失败后的信息，尽量减少模式串与主串的匹配次数以达到快速匹配的目的。具体实现就是通过一个next()函数实现，函数本身包含了模式串的局部匹配信息。KMP算法的时间复杂度O(m+n)。

Next数组

首先要明白什么是Next数组（以下简称next数组），
next[i]表示“T 中以 i 结尾的非前缀子串”与“T 的前缀”能够匹配的最长长度，即：next[i] = max{j}，其中j < i 并且 T[i-j+1, i] = T[i, j]。
next 数组的求法:
1.初始化 nex[1] = j = 0，假设nex[1, i-1] 已求出，下面求nex[i]。
2.不断尝试拓展匹配长度 j，如果拓展失败（下一个字符不相等），令 j 变为nex[j]，直至 j 为0（应该从头开始匹配）。
3.如果能够拓展成功，匹配长度 j 就增加1。nex[i] 的值就是 j 。

void getNex(const char *s){
 /*更新模式串s的nex数组*/
 int len = strlen(s);
 memset(nex,0,sizeof nex);
 for(int i = 2,j = 0;i < len;i++){
  while(j > 0 && s[i] != s[j+1]) j = nex[j];
  if(s[i] == s[j+1]) j++;
  nex[i] = j;
 }
}

代码

#include<stdio.h>
#include<string.h>
char w[10001],t[10001];
int next[10001];
int lenw,lent;
void getnext(int lenw)
{
    int i=0,j=-1;
    next[i]=-1;
    while(i<lenw)
    {
        if(j==-1||w[i]==w[j])
        {
            next[++i]=++j;
        }
        else
            j=next[j];
    }
}
int kmp(int pos,int lenw,int lent)
{
    int i=pos,j=0,ans=0;
    while(i<lent)
    {
        if(t[i]==w[j]||j==-1)
        {
            ++i;
            ++j;
        }
        else
            j=next[j];
        if(j==lenw)
        {
            ans++;
            j=next[j-1];
            i--;
        }
    }
    return ans;
}
int main()
{
    int n;
    scanf("%d",&n);
    while(n--)
    {
        scanf("%s%s",w,t);
        lenw=strlen(w);
        lent=strlen(t);
        getnext(lenw);
        printf("%d\n",kmp(0,lenw,lent));
    }
    return 0;
}