KPM算法求字符串的最小周期证明
先给出公式 ans = n - LPS[n-1]
其中ans为最小周期,n为给出的由假设的周期字符串中提取出的子串长度,LPS为前缀函数,n-1为字符串最后的位置下标
证明如下
证明ans = n - LPS[n-1],思路:
(1) 证明特殊情况,即先对完整周期字符串进行证明,这时候的字符串组成是 [1][2][3][4] ,即4个周期拼接,所以由前缀函数的定义,有
[1][2][3] = [2][3][4],所以LPS[n-1] = 3*T,即三个周期,则ans = n(即4*T) - LPS[n-1] = 4*T - 3*T = T;
对于完整周期子串显然成立.
(2) 证明一般情况,即证明非完整周期字符串,假设给出的字符串是 [末部分][1][2][3][前部分] ,即中间有完整的周期,两边是不确定长度的,可能为0.
(为了方便,此处取[末部分] = [e],[前部分] = [b])
1,当len([e]) = len([b]) = 0时,即(1)的情况,显然成立.
2,当len([e]) = 0,len([b]) != 0时此时字符串为 [1][2][3][b],由于[b]是周期的一部分,则[3]中包含[b],有
[1][2][b] = [2][3][b],此时 LPS[n-1] = 2*T + len([b]), n = 3*T + len([b]),显然有
ans = n - LPS[n-1] == 3*T + len([b]) - 2*T - len([b]) = T,成立.
3,当len([b]) = 0,len([e]) != 0时同理.
4,当len([e]) != 0,且len([b]) != 0时,此时字符串为 [e][1][2][3][b],根据2,3,有
[e][1][2][b] = [e][2][3][b],符合前缀函数定义,此时LPS[n-1] = 2*T + len([b+e]),n = 3*T + len([b+e])
显然有ans = n - LPS[n-1] = T,得证
5,当[e][b]内部没有完整的周期时,显然[e][b]可以自己组成最小周期,此时的LPS[n-1] = 0,ans = len([e+b]),为自己,得证
附上例题加以理解:
1,KMP算法模板https://www.luogu.com.cn/problem/P3375
2,字符串周期模板题https://www.luogu.com.cn/problem/P4391
例题1代码解析:
点击查看代码
//背景:刚学LPS函数及其应用KMP,先来道模板题练习,结果发现细节多
//注意:以后我的所有关于KMP的算法的字符串下标均是从0开始(便于与OI-wiki统一,好记忆)
//原理:KMP匹配字符串
//时间复杂度:o(n+m)
#include <bits/stdc++.h>
using namespace std;
void Prefixion(int LPS[],string s)//求生成字符串的前缀函数
{
LPS[0] = 0;//初始化前缀
for (int i = 1,big = s.size();i<big;i++)//遍历字符串
{
int j = LPS[i-1];//获取上一个位置的LPS
while(j>0&&s[j] != s[i]) j = LPS[j-1];//寻找第一个使得当前位置i的前缀性质仍满足的j
if(s[j] == s[i]) j++;//如果是因为相等退出循环的,j是位置,那么长度为j+1
LPS[i] = j;//记录当前位置的LPS
}
}
void KMP()//KMP算法
{
string text,pattern;
cin>>text>>pattern;//输入文本字符串及模板字符串(待查找的字符串)
string cur = pattern + '#' + text;//生成新的字符串
int s1 = pattern.size();//计算模板字符串的长度
int s2 = text.size();//计算文本字符串的长度
int LPS[s2+s1+1];//注意这里LPS的数组大小,如果开小了退不出函数(我也不知道为什么)
Prefixion(LPS,cur);//求生成字符串的前缀函数
vector<int> occurrence;//记录文本字符串匹配成功的起始位置
for (int i = s1+1;i<=s2+s1;i++)//从生成字符串的文本位置开始遍历
{
if(LPS[i] == s1) occurrence.push_back(i-2*s1);//如果满足当前情况,即记录下成功匹配的位置(相对文本):当前位置 - 2*模板长度
}
for (auto v:occurrence) cout<<v+1<<"\n";//依次输出位置(相对文本的)
for (int i = 0;i<s1;i++) cout<<LPS[i]<<" ";//输出每个前缀的函数值
cout<<"\n";
}
int main()
{
ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
int t = 1;
while(t--)
{
KMP();
}
return 0;
}
例题2代码解析:
点击查看代码
//背景:字符串的周期性问题,一开始不太理解,后来看了题解自己画图推导出来了
//本题题意:本题给出的是重复拼接字符串子串,求最小周期,也就是说假定了原字符串是周期函数,给分析证明提供了思路
//公式:ans = n - LPS[n-1]
#include <bits/stdc++.h>
using namespace std;
void Prefixion(int LPS[],string s)//前缀函数
{
LPS[0] = 0;//不要忘了初始化
for (int i = 1;s[i] != '\0';i++)
{
int j = LPS[i-1];
while(j > 0&&s[j] != s[i]) j = LPS[j-1];
if(s[j] == s[i]) j++;
LPS[i] = j;
}
}
void Solve()
{
int n;
string s;
cin>>n;//输入子串长度
cin>>s;//输入字符串
int LPS[n];//定义前缀函数
Prefixion(LPS,s);//求前缀函数
cout<<n - LPS[n-1]<<"\n";//输出最小周期
}
int main()
{
ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
int t = 1;
while(t--)
{
Solve();
}
return 0;
}
/*证明ans = n - LPS[n-1],思路:
(1)证明特殊情况,即先对完整周期字符串进行证明,这时候的字符串组成是 [1][2][3][4] ,即4个周期拼接,所以由前缀函数的定义,有
[1][2][3] = [2][3][4],所以LPS[n-1] = 3*T,即三个周期,则ans = n(即4*T) - LPS[n-1] = 4*T - 3*T = T;
对于完整周期子串显然成立.
(2)证明一般情况,即证明非完整周期字符串,假设给出的字符串是 [末部分][1][2][3][前部分] ,即中间有完整的周期,两边是不确定长度的,可能为0.
(为了方便,此处取[末部分] = [e],[前部分] = [b])
1,当len([e]) = len([b]) = 0时,即(1)的情况,显然成立.
2,当len([e]) = 0,len([b]) != 0时此时字符串为 [1][2][3][b],由于[b]是周期的一部分,则[3]中包含[b],有
[1][2][b] = [2][3][b],此时 LPS[n-1] = 2*T + len([b]), n = 3*T + len([b]),显然有
ans = n - LPS[n-1] == 3*T + len([b]) - 2*T - len([b]) = T,成立.
3,当len([b]) = 0,len([e]) != 0时同理.
4,当len([e]) != 0,且len([b]) != 0时,此时字符串为 [e][1][2][3][b],根据2,3,有
[e][1][2][b] = [e][2][3][b],符合前缀函数定义,此时LPS[n-1] = 2*T + len([b+e]),n = 3*T + len([b+e])
显然有ans = n - LPS[n-1] = T,得证
5,当[e][b]内部没有完整的周期时,显然[e][b]可以自己组成最小周期,此时的LPS[n-1] = 0,ans = len([e+b]),为自己,得证
码字不易,多多支持!