KMP总结
什么是KMP?
KMP算法,又称为模式匹配算法,能够在线性时间内判定字符串 (A[1)~(N]) 是否为字符串 (B[1)~(M]) 的子串,并求出字符串 (A) 在字符串 (B) 中各次出现的位置。(from 李煜东《算法竞赛进阶指南》)
如何进行KMP?
第一步: (A)串进行自我匹配,构造出(next)数组
第二步: 通过(next)数组和(B)串进行匹配,构造出(f)数组
什么是(next)数组?
(next[i])表示(A[1..i])的前缀和后缀最大匹配的长度(且最大匹配不能为(A)串本身),如:(下标从(1)开始)
对于串 (ababcab)来说,(next[7]=2);((A[1..2])和(A[6..7])均为(ab))
如何构造(next)数组?
(配合代码食用风味更佳)
代码如下:
int next[maxn+5];
void init(char A[],int n){
next[1]=0;
for(int i=2,j=0;i<=n;i++){
while(j>0&&A[i]!=A[j+1]) j=next[j];
if(A[i]==A[j+1]) j++;
next[i]=j;
}
}
十分简短的代码(虽然看不懂 (大雾))。
那么我们来看看是什么意思哈
- 首先,(next[1]=0),这个是显然的(定义),所以(i)直接从(2)开始;
- 接下来,就是上图的过程了;
- 对于每一次循环,(j)的初始值都等于(next[j-1]);
- 在第一次(while)中,根据定义可得(A[i-next[i-1] .. i-1])与(A[1 .. next[i-1]])相等(图中蓝色区间);于是(A[i-next[next[i-1]]..i-1])与(A[1..next[next[i-1]])(新的前后缀)相等(图中的绿色区间)
- 之后在每一次(while)循环中,由于(A[i-j .. i-1])与(A[1 .. j])相等,所以只用判断(A[i])与(A[j+1])相不相等了;也只有(A[i-j .. i-1])与(A[1 .. j])相等才有可能匹配。所以,只有(j=next[i]/next[next[i]]/...)才可能符合条件。
从某种意义而言,(next)数组的构造,就是模式串(A)自我匹配的一个过程
什么是(f)数组?
(f[i])表示(A)串的前缀和(B[1..i])的后缀最大匹配的长度,如:(下标从(1)开始)
对于串(A="abb")串 (B="ababcab")来说,(f[7]=2);((A[1..2])和(B[6..7])均为(ab))
如何构造(f)数组?
构造(f)数组的过程,和构造(next)数组的过程类似。(本质上是一样的)
话不多说,上代码:
int f[maxn+5];
void solve(char A[],int n,char B[],int m){
for(int i=1,j=0;i<=m;i++){
while(j>0&&(j==n||B[i]!=A[j+1])) j=next[j];
if(B[i]==A[j+1]) j++;
f[i]=j;
}
}
(还是看不懂)
(不如感性理解)
KMP的应用
事实上,(f)数组的构造只是其一个功能,如果将(solve)函数的一些部分进行更改,就可以用来解决一些特殊的问题。
在此不多赘述(我懒得写)