KPM算法求字符串的最小周期证明

转载 作者：撒哈拉 更新时间：2024-05-26 13:06:11 58 4

 [1][2][3] = [2][3][4],所以LPS[n-1] = 3*T,即三个周期,则ans = n(即4*T) - LPS[n-1] = 4*T - 3*T = T;

 对于完整周期子串显然成立.

 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]),为自己,得证

//背景:刚学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;
}

//背景:字符串的周期性问题,一开始不太理解,后来看了题解自己画图推导出来了
//本题题意:本题给出的是重复拼接字符串子串,求最小周期,也就是说假定了原字符串是周期函数,给分析证明提供了思路
//公式: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]),为自己,得证

58 4 0

• .net - 在sharepoint中使用SPListCollection.Add方法(字符串，字符串，字符串，字符串，Int32，字符串，SPListTemplate.QuickLaunchOptions)

  如何使用 SPListCollection.Add(String, String, String, String, Int32, String, SPListTemplate.QuickLaunchO

• C++ 字符串 != 字符串

  我刚刚开始使用 C++ 并且对 C# 有一些经验，所以我有一些一般的编程经验。然而，似乎我马上就被击落了。我试过在谷歌上寻找，以免浪费任何人的时间，但没有结果。 int main(int argc,

• Java 8 : Converting Map>到映射<字符串，字符串[]>

  这个问题已经有答案了: In Java 8 how do I transform a Map to another Map using a lambda? (8 个回答) Convert a Map>

• node.js - "Type ' 字符串 |字符串[] ' is not assignable to type ' 字符串'

  我正在使用 node + typescript 和集成的 swagger 进行 API 调用。我 Swagger 提出以下要求 http://localhost:3033/employees/sear

• C++ 映射<字符串， vector <对<字符串，字符串>>> : adding a mapping to an empty vector?

  我是 C++ 容器模板的新手。我收集了一些记录。每条记录都有一个唯一的名称，以及一个字段/值对列表。将按名称访问记录。字段/值对的顺序很重要。因此我设计如下: typedef string

• java - 谁能帮我创建方法？ mystring.replacefirst(字符串,字符串);并替换(自，直到，字符串)；对于j2me，请

  我需要这两种方法，但j2me没有，我找到了一个replaceall();但这是 replaceall(string,string,string); 第二个方法是SringBuffer但在j2me中它没

• .net - 字符串 vs 字符串 - 区分大小写的联合

  If string is an alias of String in the .net framework为什么会发生这种情况，我应该如何解释它: type JustAString = string

• python - 考虑顺序如何检查列表(字符串)是否包含另一个列表(字符串)

  我有两个列表(或字符串):一个大，另一个小。 我想检查较大的(A)是否包含小的(B)。 我的期望如下: 案例 1. B 是 A 的子集 A = [1,2,3] B = [1,2] contains(A

• javascript - Jquery 字符串 + 对象 + 字符串

  我有一个似乎无法解决的小问题。 这里...我有一个像这样创建的输入... var input = $(''); 如果我这样做......一切都很好 $(this).append(input); 如果我

• c# - ienumerable <字符串>到列表<字符串>

  我有以下代码片段 string[] lines = objects.Split(new string[] { "\r\n", "\n" }, StringSplitOptions.No

• Java printf(字符串、Val、字符串)

  这可能真的很简单，但我已经坚持了一段时间了。 我正在尝试输出一个字符串，然后输出一个带有两位小数的 double ，后跟另一个字符串，这是我的代码。 System.out.printf("成本:%.2

• firebase - Cloud Firestore 字符串 >= 字符串

  以下是 Cloud Firestore 列表查询中的示例之一 citiesRef.where("state", ">=", "CA").where("state", "= 字符串，我们在Stack O

• javascript - if(字符串.包含(字符串))。那可能吗？

  我正在尝试检查一个字符串是否包含在另一个字符串中。后面的代码非常简单。我怎样才能在 jquery 中做到这一点？ function deleteRow(locName, locID) { if

• C++ 字符串 (int) + 字符串 (int)

  这个问题在这里已经有了答案: How to implement big int in C++ (14 个答案) 关闭 9 年前。 我有 2 个字符串，都只包含数字。这些数字大于 uint64_t 的

• java - 带有自定义转换器的推土机双向映射(字符串，字符串)不可能吗？

  我有一个带有自定义转换器的 Dozer 映射: com.xyz.Customer com.xyz.CustomerDAO customerName

• java - 字符串 a == 字符串 b 的规则

  这个问题在这里已经有了答案: How do I compare strings in Java? (23 个回答) 关闭 6 年前。 我想了解字符串池的工作原理以及一个字符串等于另一个字符串的规则是

• Swift 字符串 vs. 字符串!与字符串？

  我已阅读 this问题和其他一些问题。但它们与我的问题有些无关 对于 UILabel 如果你不指定 ? 或 ! 你会得到这样的错误: @IBOutlet property has non-option

• c - 字符串 [x] 与 *字符串++

  这两种方法中哪一种在理论上更快，为什么？ (指向字符串的指针必须是常量。) destination[count] 和 *destination++ 之间的确切区别是什么？ destination[co

• .net - String.Format与“字符串” +“字符串”还是StringBuilder？

  This question already has answers here: Closed 11 years ago. Possible Duplicates: Is String.Format a

• java - 流<字符串> 到映射<字符串、整数>

  我有一个Stream一个文件的，现在我想将相同的单词组合成 Map这很重要，这个词在 Stream 中出现的频率. 我知道我必须使用 collect(Collectors.groupingBy(..)