KMP算法

      今天学了KMP算法，一开始看代码看不大懂，看了一上午教学视频总算是大致了解了，难度确实并查集和字典树大。

      KMP算法主要应用于求出ex[i]数组，ex[i]数组表示主串S中以第i个字符(也即S[i-1]位置)为尾字符的后缀与模式串T的前缀的最长公共部分。

     拿一道经典的kmp模版题来介绍模版：

    POJ 3461 Oulipo(KMP:统计一个串出现的次数)      http://poj.org/problem?id=3461
题意:
        给你一个T串和P串,问你P串在T串中出现了多少次.
分析:
        直接套用KMP模板即可.

     模版1：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn = 1e6+10;
int n,Next[maxn];
char mo[maxn],str[maxn];
int main(){
    scanf("%d",&n);
    while(n--){
        memset(Next,-1,sizeof(Next));
        scanf("%s%s",mo,str);
      //next数组构建
        int i = 0,j = -1,len = strlen(mo);
        while(i < len){
            if(j == -1||mo[i] == mo[j]) Next[++i] = ++j;
            else j = Next[j];
        }
       //查找
        int ans = 0,len1 = strlen(str);
        i = j = 0;
        while(i < len1){
            if(j == -1||str[i] == mo[j]) ++i,++j;
            else j = Next[j];
            if(j == len) ans++;
        }
        printf("%d\n",ans);
    }
    return 0;
}

  模版2：//此模版应用比较广泛，刘汝佳的算法竞赛训练指南也介绍的这个模版

  #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];
        f[i+1] = P[i]==P[j]?j+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;
}