字符串（2）KMP算法

给你两个字符串a(len[a]=n),b(len[b]=m)，问b是否是a的子串，并且统计b在a中的出现次数，如果我们枚举a从什么位置与匹配，并且验证是否匹配，那么时间复杂度O（nm），

而n和m的范围为10^5，这样做显然超时，因此我们就要用到神奇的KMP算法，在O（n）的时间内解决这一类的问题。

首先给出两个字符串

A：abababaababacb

B：ababacb

首先我们思考朴素算法，我们枚举A串的每一位作为开始与B串匹配的位置，然后一位一位进行检验，如果匹配失败则会从B串的开始重新匹配，这也是朴素算法为什么会超时的原因所在，那么我们可不可以不从B串的开始与A串的下一位匹配，而是直接找到一个可以与A串枚举到的那一位匹配的B串的一个位置，这样我们只需要扫描一遍字符串A，然后在更新可以匹配到B串的哪一个位置。

举个例子：

对于字符串A，B，他们的前五位是一样的，可以匹配。

但是第6位就无法匹配了，所以我们需要调整B串的位置重新匹配，通过观察我们发现B串的前3位和A串的3到5位可以匹配，所以我们可以直接将j的值改为3，然后继续和A串匹配。

这样一来我们就又可以继续向后匹配了。

但是i=8时又无法匹配了，我们需要继续调整B串的位置以便于进行匹配。

然而这样仍然不能匹配，那就继续移动，直至可以匹配为止。

最终我们成功将A串与B串匹配成功。

我们根据B串的移动规律可以构造出这样的一个p数组，使得p[j]表示B[1…j]=B[j-k+1…j]的k的最大值（即B串前j个字符最长相同的前缀和后缀的长度），这样我们就可以直接将B串进行上述的移动操作了。

 void kmp()
 {
     int j=;
     for(int i=;i<n;i++)
     {
         while(j>&&b[j+]!=a[i+]) j=p[j];//如果下一位不能匹配就移动B串
         if(b[j+]==a[i+]) j++;//如果可以匹配就继续匹配下一位
         if(j==m)//这里表示完全匹配，也就是A串和B串成功匹配了
         {
             printf("%d\n",i-m+);
             j=p[j];//这样的目的是为了不遗漏重叠的匹配
         }
     }
 }

kmp匹配

接下来我们思考如何求这个p数组，假如我们知道了p[1…j-1]的值，那我们如何求p[j]的值呢？如果B[i+1]==B[j+1],那么显然p[j]=p[j-1]+1,因为这相当于前缀后缀都加了一个相同的字符，总长都加上1，那么若果B[i+1]!=B[j+1]我们可以考虑将j向后退一步，也就是减小j的值，再进行匹配。

 void pre()
 {
     p[]=;
     int j=;
     for(int i=;i<m;i++)
     {
         while(j>&&b[i+]!=b[j+]) j=p[j];
         if(b[i+]==b[j+]) j++;
         p[i+]=j;
     }
 }

求p数组

这样我们就成功完成了kmp算法