配合y总的视频和代码理解比较好
- 最最最需要注意的就是对
ne[i]
的理解,他表示p[1---i]
中最长公共前后缀的长度,并不是一般教材上写的失配时跳转的位置
时间复杂度
$O(m+n)$
参考文献
个人理解对照代码
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int N = 100010;
const int M = 1000010;
char p[N],s[M];
int n,m;
int ne[N];//ne[i]:表示p[1~i]中最长公共前后缀的长度
void getNext() {
//因为ne[1]=0 直接从第二位开始即可
for(int i=2,j=0;i<=n;++i) {
while(j&&p[i]!=p[j+1]) j=ne[j];
//匹配成功,ne[i]=j+1,同时自增继续和下一位p[i+1]比较
if(p[i]==p[j+1]) ne[i]=++j;
//匹配失败,说明完全对不上,ne[i]=0,p[i+1]从头开始比较,故j不用自增1
if(j==0) ne[i]=0;
}
}
void kmp() {
getNext();
//主串从1开始
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;
//如果完全匹配则进行输出,同时没有p[n+1],下一步必定失配,直接j=ne[j]进行一步失配操作,
//同时也可预防越界
if(j==n) printf("%d ",i-n),j=ne[j];
}
//匹配失败 j不用自增 从头开始比较 s[i+1]与p[1]比较
if(j==0) {
;
}
}
}
int main() {
#ifdef ONLINE_JUDGE
#else
freopen("in.txt","r",stdin);
freopen("out.txt","w",stdout);
#endif
scanf("%d",&n);
scanf("%s",p+1);
scanf("%d",&m);
scanf("%s",s+1);
kmp();
return 0;
}