【模板】KMP字符串匹配
题目
P3375
ybtoj只需要输出出现次数
解析
KMP算法的名字是来源于三位发明者*Knuth(D.E.Knuth)&Morris(J.H.Morris)&Pratt(V.R.Pratt)*的名字,其特点是线性
考虑设
n
x
t
i
nxt_i
nxti为
s
1
,
i
s_{1,i}
s1,i的最大公共前后缀(即题意中的border)
显然有
n
x
t
1
=
0
nxt_1=0
nxt1=0
考虑
n
x
t
nxt
nxt递推,发现可以成为
n
x
t
i
nxt_i
nxti的备用解只有
n
x
t
i
−
1
,
n
x
t
n
x
t
i
−
1
.
.
.
nxt_{i-1},nxt_{nxt_{i-1}}...
nxti−1,nxtnxti−1...
直接暴力跳失配即可
因为每次最多让它多配对一个,又不可能跳到负数,所以是线性的
再考虑求出后的配对
仿照递推,直接改一下配对位置即可
code:
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int kmp[1000010],l1,l2,q;
char a[1000010],b[1000010];
int main()
{
cin>>a+1>>b+1;
l1=strlen(a+1),l2=strlen(b+1);
for(int i=2;i<=l2;++i)
{
while(q&&b[q+1]!=b[i])q=kmp[q];
if(b[q+1]==b[i])++q;
kmp[i]=q;
}
q=0;
for(int i=1;i<=l1;++i)
{
while(q&&b[q+1]!=a[i])q=kmp[q];
if(b[q+1]==a[i])++q;
if(q==l2)printf("%d\n",i-l2+1),q=kmp[q];
}
for(int i=1;i<=l2;++i)printf("%d ",kmp[i]);
return 0;
}