【总结】字符串匹配: KMP 和拓展KMP

比起ac自动机,kmp就一个next数组,理解了如何初始化next后就可以搞一些模板题了,下面是还不错的学习资料,清晰易懂,自己用的模板也来自它:

　http://chaoswork.com/blog/2011/06/14/kmp%E7%AE%97%E6%B3%95%E5%B0%8F%E7%BB%93/

kmp模板

next[0]=-1;j=-1;
for(i=0;i<m;)
{while(j>=0 && p[i]!=p[j])j=next[j];i++;j++;next[i]=j;
}for(i=0,flag=0,j=0;i<n;)
{while(j>=0 && p[j]!=t[i])j=next[j];i++;j++;//已经匹配了模式串的多少个if(j==m)//匹配成功
}

接下来需要更加深入地了解next数组,许多题目需要用到它的定义来预处理字符串:

xdu 1154

显然,用2次kmp处理处前缀和后缀在各点的匹配情况,用dp记录符合要求的子串,有几个要注意的地方:(调了2个小时T_T)

1.前缀的其实位置不但要在后缀的前面,终止位置不能超过后缀的终止位置,也就是说前缀不能包含后缀.

View Code

for(i=0,flag=0,j=0;i<n;)
{while(j>=0 && p[j]!=t[i])j=next[j];i++;j++;

　　if(j==m)//若在此处记录则会出现bug,因为此时匹配完成点是i-1
}

CF也有一道与上面类似的题,需要用kmp预处理最左端前缀和最右端后缀

http://codeforces.com/contest/149/problem/E

循环节问题

它还能用来求周期字符串的循环节HDU1358:

性质:

当且仅当len%(len-next[len])==0时,str[next[len]~len-1]为最小循环节

要证明它需要说明3点:

1.一个字符串str是周期串,假设s1为它的循环节,则str=s1 s1...s1(n个)

推导出=>len%(len-next[len])==0成立.

2.next[len]~len-1为s1 ,len%(len-next[len])==0时

推导出=>str为周期串,s1为最小循环节

3.如何保证是最小的.

证明:

1.由next的性质知道,s[1~next[len]-1]与s[1~len-1]有最长的相等的前缀和后缀s,很显然s就是n-1个s1了.

2.设s1的长度为L,由于len%L==0 , str可以分解成若干个长度为L的小串,设它们从左到右依次为

a₁, a₂... a_n

根据匹配关系得

a₁=a₂;

a₂=a₃;

a_n-1=a_n;

因此a₁=a₂=a₃...=a_n;

3.next[len]保证了前缀与后缀最大化,如果循环节s1存在而s1内还有循环节s1',则next[len]可以向后移动,与定义矛盾.

相等的循环同构问题

hdu3374 string problem

分析:

　　最大和最小的循环同构字符串可以分开处理,求位置可以利用最小最大表示法(03wc论文周源)

　　求出现的次数可以用kmp扫一遍,但利用循环节可以更快地得到答案.

性质:

　　字符串str由k个最小循环节s1组成,则它的相等循环同构数为k.

证明:

　　设相等循环同构数为p,我们可以利用循环节s1构造出k个相等的循环同构,于是p>=k;

　　下面证明p<=k:

　　　　如果p>k

　　　　假设在移动完s1到尾部前,出现了相同的循环同构串,不妨设此时移动的串为s2,

　　　　则利用循环节性质2第1种情况的推理方法可以知道s2为字符串的一个循环节

　　　　且s2的长度<s1,与s1为最小循环节的条件矛盾,因此p<=k.

　　于是p只能等于k.

拓展KMP

　　有很长一段时间单纯地以为拓展kmp只是kmp倒过来跑,后来发现很多问题其实无法转换成kmp解决,于是怒学了一下拓展kmp.

　　首先比较一下kmp和拓展kmp解决的问题:

　　kmp解决了求所有主串的前缀pre[i] (0<=i<n)的后缀与模式串前缀的最大匹配长度问题;

　　拓展kmp解决了所有主串的后缀suf[i]的前缀,与模式串前缀的最大匹配长度问题.

 1 void getNext() {
 2     int l = 1, r = -1, i, j;
 3 
 4     for (next[0] = 0, i = 1; p[i]; ++i) {
 5         if (i + next[i - l] - 1 < r) next[i] = next[i - l];
 6         else {
 7             for (j = max(r - i + 1, 0); p[i + j] && p[i + j] == p[j]; ++j);
 8             next[i] = j; l = i; r = i + j - 1;
 9         }
10     }
11     next[0] = i;
12 }