C++刷题 -- KMP算法

C++刷题 – KMP算法

---

https://leetcode.cn/problems/find-the-index-of-the-first-occurrence-in-a-string/description/

1.算法讲解

KMP算法是一种字符串匹配算法，当出现字符串不匹配时，可以记录一部分之前已经匹配的文本内容，利用这些信息避免从头再去做匹配；

前缀表（prefix table）：
前缀表是用来回退的，它记录了模式串与主串(文本串)不匹配的时候，模式串应该从哪里开始重新匹配。

• 举一个例子：
  要在文本串：aabaabaafa 中查找是否出现过一个模式串：aabaaf。
  可以看出，文本串中第六个字符b 和 模式串的第六个字符f，不匹配了。如果暴力匹配，发现不匹配，此时就要从头匹配了。
  但如果使用前缀表，就不会从头匹配，而是从上次已经匹配的内容开始匹配，找到了模式串中第三个字符b继续开始匹配。

• 前缀表是如何记录的：
  首先要知道前缀表的任务是当前位置匹配失败，找到之前已经匹配上的位置，再重新匹配，此也意味着在某个字符失配时，前缀表会告诉你下一步匹配中，模式串应该跳到哪个位置。
  前缀表：记录下标i之前（包括i）的字符串中，有多大长度的相同前缀后缀

• 最长公共前后缀：
  字符串的前缀是指不包含最后一个字符的所有以第一个字符开头的连续子串
  后缀是指不包含第一个字符的所有以最后一个字符结尾的连续子串
  前缀表要求的就是相同前后缀的长度；
  字符串a的最长相等前后缀为0。 字符串aa的最长相等前后缀为1。 字符串aaa的最长相等前后缀为2。 等等…

• 为什么一定要用前缀表
  刚刚匹配的过程在下标5的地方遇到不匹配，模式串是指向f，如图：
  然后就找到了下标2，指向b，继续匹配：如图
  
  下标5之前这部分的字符串（也就是字符串aabaa）的最长相等的前缀 和 后缀字符串是 子字符串aa ，因为找到了最长相等的前缀和后缀，匹配失败的位置是后缀子串的后面，那么我们找到与其相同的前缀的后面重新匹配就可以了。
  所以前缀表具有告诉我们当前位置匹配失败，跳到之前已经匹配过的地方的能力。

• 如何计算前缀表
  如图：
  
  长度为前1个字符的子串a，最长相同前后缀的长度为0。（注意字符串的前缀是指不包含最后一个字符的所有以第一个字符开头的连续子串；后缀是指不包含第一个字符的所有以最后一个字符结尾的连续子串。）
  
  长度为前2个字符的子串aa，最长相同前后缀的长度为1；
  
  长度为前3个字符的子串aab，最长相同前后缀的长度为0
  以此类推： 长度为前4个字符的子串aaba，最长相同前后缀的长度为1。 长度为前5个字符的子串aabaa，最长相同前后缀的长度为2。 长度为前6个字符的子串aabaaf，最长相同前后缀的长度为0。
  那么把求得的最长相同前后缀的长度就是对应前缀表的元素，如图：
  
  可以看出模式串与前缀表对应位置的数字表示的就是：下标i之前（包括i）的字符串中，有多大长度的相同前缀后缀。
  再来看一下如何利用 前缀表找到 当字符不匹配的时候应该指针应该移动的位置。如动画所示：
  
  找到的不匹配的位置， 那么此时我们要看它的前一个字符的前缀表的数值是多少，因为要找前面字符串的最长相同的前缀和后缀。
  前一个字符的前缀表的数值是2， 所以把下标移动到下标2的位置继续比配。
  最后就在文本串中找到了和模式串匹配的子串了。

• 前缀表与next数组
  很多KMP算法的实现都是使用next数组来做回退操作，那么next数组与前缀表有什么关系呢？
  next数组就可以是前缀表，但是很多实现都是把前缀表统一减一（右移一位，初始位置为-1）之后作为next数组。
  其实这并不涉及到KMP的原理，而是具体实现，next数组既可以就是前缀表，也可以是前缀表统一减一（右移一位，初始位置为-1）。

2.算法实现

• 使用next数组来匹配
  以下我们以前缀表统一减一之后的next数组来做演示。
  有了next数组，就可以根据next数组来 匹配文本串s，和模式串t了。
  注意next数组是新前缀表（旧前缀表统一减一了）。
  匹配过程动画如下

• 时间复杂度分析
  其中n为文本串长度，m为模式串长度，因为在匹配的过程中，根据前缀表不断调整匹配的位置，可以看出匹配的过程是O(n)，之前还要单独生成next数组，时间复杂度是O(m)。所以整个KMP算法的时间复杂度是O(n+m)的。
  暴力的解法显而易见是O(n × m)，所以KMP在字符串匹配中极大地提高了搜索的效率。
  为了和力扣题目28.实现strStr保持一致，方便大家理解，以下文章统称haystack为文本串, needle为模式串。

• 构造next数组
  定义一个函数getNext来构建next数组，函数参数为指向next数组的指针，和一个字符串。
  
  构造next数组其实就是计算模式串s，前缀表的过程。 主要有如下三步：

  1. 初始化：
     定义两个指针i和j，j指向前缀末尾位置，i指向后缀末尾位置。
     然后还要对next数组进行初始化赋值，如下：
     
     j 为什么要初始化为 -1呢，因为之前说过 前缀表要统一减一的操作仅仅是其中的一种实现，我们这里选择j初始化为-1，下文我还会给出j不初始化为-1的实现代码。
     next[i] 表示 i（包括i）之前最长相等的前后缀长度（其实就是j）
     所以初始化next[0] = j 。

  2. 处理前后缀不相同的情况:
     因为j初始化为-1，那么i就从1开始，进行s[i] 与 s[j+1]的比较
     所以遍历模式串s的循环下标i 要从 1开始，代码如下：
     
     如果 s[i] 与 s[j+1]不相同，也就是遇到 前后缀末尾不相同的情况，就要向前回退。
     怎么回退呢？
     next[j]就是记录着j（包括j）之前的子串的相同前后缀的长度。
     那么 s[i] 与 s[j+1] 不相同，就要找 j+1前一个元素在next数组里的值（就是next[j]）。
     
     每次求的都是当前字符串开头到i之间的子串的最大公共前后缀；
     因此j指向子串的一个前缀末尾位置，i指向子串的一个后缀末尾位置
     例如：
     
     上图中判断的是aab这个子串的前后缀，显然没有相同的前后缀，因此next[i]为0；
     
     上图判断的是aabaa这个子串的前后缀，j+1指向的是当前的前缀aa的末尾，i指向的是当前的后缀aa的末尾，显然前后缀末尾是一致的，因此next[i] = j = 1；

  3. 处理前后缀相同的情况:
     如果 s[i] 与 s[j + 1] 相同，那么就同时向后移动i 和j 说明找到了相同的前后缀，同时还要将j（前缀的长度）赋给next[i], 因为next[i]要记录相同前后缀的长度。
     

构建next数组的逻辑：

最后整体构建next数组的函数代码如下：

void getNext(int* next, const string& s){
    int j = -1;
    next[0] = j;
    for(int i = 1; i < s.size(); i++) { // 注意i从1开始
        while (j >= 0 && s[i] != s[j + 1]) { // 前后缀不相同了
            j = next[j]; // 向前回退
        }
        if (s[i] == s[j + 1]) { // 找到相同的前后缀
            j++;
        }
        next[i] = j; // 将j（前缀的长度）赋给next[i]
    }
}

• 使用next数组来做匹配
  在文本串s里 找是否出现过模式串t。
  定义两个下标j 指向模式串起始位置，i指向文本串起始位置。
  那么j初始值依然为-1，为什么呢？ 依然因为next数组里记录的起始位置为-1。
  i就从0开始，遍历文本串，代码如下：
  
  接下来就是 s[i] 与 t[j + 1] （因为j从-1开始的） 进行比较。
  如果 s[i] 与 t[j + 1] 不相同，j就要从next数组里寻找下一个匹配的位置。
  代码如下：
  
  如果 s[i] 与 t[j + 1] 相同，那么i 和 j 同时向后移动， 代码如下：
  
  如何判断在文本串s里出现了模式串t呢，如果j指向了模式串t的末尾，那么就说明模式串t完全匹配文本串s里的某个子串了。
  本题要在文本串字符串中找出模式串出现的第一个位置 (从0开始)，所以返回当前在文本串匹配模式串的位置i 减去 模式串的长度，就是文本串字符串中出现模式串的第一个位置。
  那么使用next数组，用模式串匹配文本串的整体代码如下：

int j = -1; // 因为next数组里记录的起始位置为-1
for (int i = 0; i < s.size(); i++) { // 注意i就从0开始
    while(j >= 0 && s[i] != t[j + 1]) { // 不匹配
        j = next[j]; // j 寻找之前匹配的位置
    }
    if (s[i] == t[j + 1]) { // 匹配，j和i同时向后移动
        j++; // i的增加在for循环里
    }
    if (j == (t.size() - 1) ) { // 文本串s里出现了模式串t
        return (i - t.size() + 1);
    }
}

KMP算法的整体代码：

class Solution {
public:
    void getNext(int* next, const string& s) {
        int j = -1;
        next[0] = j;
        for(int i = 1; i < s.size(); i++) { // 注意i从1开始
            while (j >= 0 && s[i] != s[j + 1]) { // 前后缀不相同了
                j = next[j]; // 向前回退
            }
            if (s[i] == s[j + 1]) { // 找到相同的前后缀
                j++;
            }
            next[i] = j; // 将j（前缀的长度）赋给next[i]
        }
    }
    int strStr(string haystack, string needle) {
        if (needle.size() == 0) {
            return 0;
        }
        int next[needle.size()];
        getNext(next, needle);
        int j = -1; // // 因为next数组里记录的起始位置为-1
        for (int i = 0; i < haystack.size(); i++) { // 注意i就从0开始
            while(j >= 0 && haystack[i] != needle[j + 1]) { // 不匹配
                j = next[j]; // j 寻找之前匹配的位置
            }
            if (haystack[i] == needle[j + 1]) { // 匹配，j和i同时向后移动
                j++; // i的增加在for循环里
            }
            if (j == (needle.size() - 1) ) { // 文本串s里出现了模式串t
                return (i - needle.size() + 1);
            }
        }
        return -1;
    }
};

• 时间复杂度: O(n + m)
• 空间复杂度: O(m), 只需要保存字符串needle的前缀表

前缀表不减1

class Solution {
public:
    void getNext(int* next, const string& s)
    {
        int j = 0;//字符串当前检查的前缀的末尾
        next[0] = j;//初始化

        for(int i = 1; i < s.size(); i++)
        {
            //前后缀两个子串都是从大小为1的子串开始检查的
            //如果前缀后缀末尾不一致，就对j进行回退，看回退后能否匹配
            while(j > 0 && s[i] != s[j])
            {
                j = next[j - 1];//next中记录着已经检查到的相同前后缀的长度
            }
            //如果前后缀末尾一致，就对j++
            if(s[i] == s[j])
            {
                j++;
            }
            //当出现一直的前后缀末尾，就需要在next表中记录当前的相同前后缀长度
            next[i] = j;
        }
    }

    int strStr(string haystack, string needle) {
        if(needle.size() == 0)
        {
            return -1;
        }
        int next[needle.size()];
        getNext(next, needle);
        int j = 0; // 指向needle的指针
        for(int i = 0; i < haystack.size(); i++)
        {
            //判断当前两指针所指字符是否相等，不相等就回退j
            while(j > 0 && haystack[i] != needle[j])
            {
                j = next[j - 1];
            }
            //如果一致，就对j++
            if(haystack[i] == needle[j])
            {
                j++;
            }
            //如果j遍历到了needle的末尾，就说明存在子串
            if(j == needle.size())
            {
                return i - needle.size() + 1;
            }
        }

        return -1;
    }
};