c++数据结构关于KMP模式匹配next函数

Posted 2023-04-29

大家帮我解决下这个问题，可以不？
求这个模式串的next函数
'abaaabcabcbcac'

我自己是算出来了，但又不敢肯定结果，大家帮我看下，好不？

我算出来的是：
0,1,1,
2,2,2,
4,3,2,
7,3,5,
3,2.

这是我算出来的上面那个模式串1～14的next[]

有没有算对呀？

参考技术A 下标从1开始, next[i]表示已匹配上i个字符, 第i+1个字符不匹配, 余下部分所能匹配的最长前缀
next[1-14] = 0, 0, 1, 1, 1, 2, 0, 1, 2, 0, 0, 0, 1, 0

参考程序, pat为字符串指针, 首字母位于pat[1], p[]即为你所述的next[]
void compute_next(char *pat)
p[1] = 0;
int j = 0;
for(int i = 2; pat[i]; ++i)
while(j > 0 && pat[j + 1] != pat[i]) j = p[j];
if(pat[j + 1] == pat[i]) j += 1;
p[i] = j;

数据结构——KMP算法

算法介绍

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

next数组

　　我们记主串为字符串S，模式串为字符串P。

　　我们用next[j]表示以字符Pj结尾的子串的长度相等的前缀字符串与后缀字符串长度的最大值。

　　为了方便起见，我们将字符串从下标1开始匹配。如此，next数组所表示的长度就与下标数值相等了。

算法思路

　　我们从左到右依次枚举S的每一个字符Si，对于当前待匹配字符Si，我们假设当前P字符串中已匹配到Pj。

　　那么我们只需判断Si和Pj+1，若两者相同，则继续匹配。

　　若两者不相同，那么我们使j=next[j]，即可最大限度的减少匹配次数。因为S字符串的从某位置开始到前i-1的部分与P字符串的前j个字符已匹配(即完全相同)，而P1到Pnext[j]部分是长度最大的与以Pj结尾的后缀完全相同的前缀，而该以Pj结尾的后缀则必定与S中一段以Si-1结尾的子串完全相同，因而保证了上述操作的正确性。

　　接下去只需重复上述操作即可。

　　而对于next数组的预处理，也同上述操作类似，我们只需要以字符串P来匹配字符串P即可。

模板呈现

　　模板题链接：KMP字符串

　　代码如下：

#include <iostream>
#include <algorithm>
#include <cstdio>
using namespace std;
const int M = 1e5+10;
int n,m;
int ne[M];
char s[M],p[M];
int main()

    cin>>n>>p+1;
    cin>>m>>s+1;

    for(int i=2,j=0;i<=n;i++)
    
        while(j && p[i]!=p[j+1])j=ne[j];
        if(p[i]==p[j+1])j++;
        ne[i]=j;
    
    for(int i=1,j=0;i<=m;i++)
    
        while(j && s[i]!=p[j+1])j=ne[j];
        if(s[i]==p[j+1])j++;
        if(j==n)
        
            printf("%d ",i-n+1-1);
            j=ne[j];  //可有可无，好习惯要加上。若为string，不加会出错。
        
    
    printf("\n");
    return 0;