理解KMP算法

KMP算法理解

KMP算法是什么

引用来自维基百科的一段话

在计算机科学中,Knuth-Morris-Pratt字符串查找算法(简称为KMP算法)可在一个字符串S内查找一个词W的出现位置。一个词在不匹配时本身就包含足够的信息来确定下一个匹配可能的开始位置,此算法利用这一特性以避免重新检查先前配对的字符。这个算法由高德纳和沃恩·普拉特在1974年构思,同年詹姆斯·H·莫里斯也独立地设计出该算法,最终三人于1977年联合发表。

现在用我的话来解释

KMP算法原理就是说就是减少重复的比对项,也就是对字符串匹配暴力算法的一种优化。

暴力算法

一般的字符串匹配解法是将 2 个串的字符进行挨个比较,相当于是把每个字符都比较了一遍,这样是一定能得到结果的,不过显然这样操作导致的时间复杂度是O(m∗n),也就是 2 个字符串的长度之积,俗称暴力解法。

用图像表示暴力解法的过程如下:

一开始,i 和 j 分别指向主字符串和子字符串,

一般的字符串匹配解法是将 2 个串的字符进行挨个比较,相当于是把每个字符都比较了一遍,这样是一定能得到结果的,不过显然这样操作导致的时间复杂度是O(m∗n),也就是 2 个字符串的长度之积,俗称暴力解法。

用图像表示暴力解法的过程如下:

一开始,i 和 j 分别指向主字符串和子字符串,

这个时候开始比较 t[i]p[j],如果匹配则 i 和 j 同时向右移动一格,这样 i 和 j 同时移动到了位置 3,

此时 i 和 j 不匹配了,需要对 i 和 j 进行回退,i 回退到了它最开始的下一位置,即位置 1,而 j 回退到位置 0

此时 i 和 j 仍然不匹配,i 继续向右移动一格,即位置 2,j 依然保持位置 0

如此反复,每当不匹配时,i 都会回到开始匹配的位置同时移动一格,而 j 回到位置 0,当匹配时,i 和 j 同时向右移动,直到 j 到达子字符串的末尾。

代码展示

arduino 复制代码
public static int Befalgorithm(String text, String pattern) {

        for (int i = 0; i <= text.length() - pattern.length(); i++) {
            int j = 0;
            for (; j < pattern.length(); j++) {
                if (text.charAt(i+j) != pattern.charAt(j)) {
                    break;
                }
            }
            if (j == pattern.length()) {
                return i;
            }

        }
        return -1;

    }

KMP算法原理

KMP 算法用了一种聪明的办法,当发现字符串不匹配的时候,并不会从头开始比较,因为之前已经匹配成功的字符可以给我们提供一些有用的信息,利用这个信息我们可以将子串移动到某个位置,并从这个位置直接开始比较,它的时间复杂度降到了O(m+n),也就是 2 个字符串的长度之和。

字符串的前缀和后缀

首先我们需要知道字符串的前缀和后缀:

对于字符串 ababc 来说,它的前缀有 [a,ab,aba,abab],也就是以字符串第一个字符作为开头,同时不包括最后一个字符的所有子串,同理它的后缀有 [c,bc,abc,babc],也就是以字符串最后一个字符作为结尾,同时不包括第一个字符的所有字串。

字符串的最长公共前后缀

了解了这个,我们再来说什么是字符串的最长公共前后缀 ,说白了,也就是前缀和后缀这 2 个集合中的相同部分,同时取最长的那个,就是这个字符串的最长公共前后缀 。显然,在这个例子中,ababc 是没有公共前后缀的。但是对于 abab,它的前缀和后缀分别是 [a,ab,aba][b,ab,bab],那么它的最长公共前后缀 就是 ab

现在,我们的目标就是取得 ababc 所有子串 [a,ab,aba,abab,ababc]最长公共前后缀 的长度,分别保存在 next 数组中,我们只要知道最大长度即可,并不用关心串具体是什么,而我们目前通过观察即可得出 next 数组这里是 [0,0,1,2,0],我们先知道这个结果即可,此计算方法后续会说明。

KMP 的思路分析

得到这个有什么用,回到刚刚的例子,当我们第一次碰到不匹配时(i 和 j 都等于 3 的时候):

此时 i 和 j 之前的 3 个字符 aba 必然完全匹配的,因为只有前面匹配成功了 j 才能走到 3,而我们又知道 aba最长公共前后缀a,这可以给我们一个很好的提示,主串中的 aba 的后缀和字串中的 aba 前缀有最长的公共部分 a,这样一来,我们就没必要重新比较了,直接将相同部分对齐就好了,也就是说让 j 退回到位置 1 就可以了,而 i 保持不变。

分析一下,为什么 j 知道应该退回到位置 1:

  1. 我们知道 j 是在位置 3 不匹配的,那么 j 前面的字符串我们可以知道是 aba
  2. 前面我们已经得到了 next 数组,知道 aba 的最长公共前后缀长度是 1
  3. 也就是说,我们可以知道 i 之前 1 个位置(主串的后缀)和子串开头之后 1 个位置(子串的前缀)是相同的
  4. 进而我们可以让相同部分对齐,也就是让 j=next[j-1]

接下来,我们发现 i 和 j 仍然不匹配,而 j 之前的字符 a 最长公共前后缀是 0,此时 j 退到位置 0,i 仍然保持不变。

i 和 j 匹配,同时右移一格

此时 i 和 j 不匹配,j=next[j-1] 回退到 0,然后我们发现 i 和 j 仍然不匹配,同时 j 已经是 0 了,那么我们就让 i 往右移动一格。

可以看到,相比于暴力解法,i 始终在前进,并没有后退(顶多保持不变),然后我们通过 next 数组来改变 j 的值。

ini 复制代码
public static int KMP(String text, String pattern) {
        int[] next = getNext(pattern);
        int n1 = text.length();
        int n2 = pattern.length();
        if (n2 == 0) {
            return 0;
        }
        int i = 0;
        int j = 0;
        while (i < n1) {
            if (text.charAt(i) == pattern.charAt(j)) {
                i++;
                j++;
            } else {
                if (j > 0) {
                    j = next[j - 1];
                } else {
                    i++;
                }
            }
            if (j == n2) {
                return i - j;
            }
        }

       return -1;
    }

求解字符串最长公共前后缀

最后,我们还剩下一个问题,如何求 next 数组,也就是求模式字符串各个子串的最长公共前后缀长度。

我们先初始化一个和模式字符串长度相等的 next 数组,在 next 数组中,第 1 位默认为 0,因为我们规定只有一个字符的字符串没有前缀和后缀,自然公共前后缀长度只能是 0。

我们依然设定 2 个指针 i 和 j,j 指向位置 0,i 从位置 1 开始依次为 next 数组进行赋值

我们可以把这个过程依然看作是 2 个字符串的比较,j 指向的是模式字符串的前缀,而 i 指向的是模式字符串的后缀

和上面的字符串匹配一样,我们执行同样的处理:

  1. 当 i 和 j 匹配时,i 和 j 同时右移一格
  2. 当 i 和 j 不匹配时,如果 j 不在字符串开头(位置 0),就回退到上一个能匹配到的位置
  3. 当 i 和 j 不匹配时,如果 j 在字符串开头(位置 0),那么 i 就右移一格

对 next [1] 赋值:i 和 j 不匹配,同时 j 已经是字符串开头,赋值 0

对 next [2] 赋值,i 和 j 匹配,此时 j 为 0,代表只有 1 个字符匹配 (j+1),赋值 1

对 next [3] 赋值,i 和 j 匹配,此时 j 为 1,代表有 2 个字符匹配 (j+1),赋值 2

对 next [4] 赋值,i 和 j 不匹配,此时 j 为 2,可以得知 j 前面的字符是 ab,而 ab 的最长公共前后缀长度就是 next[1],这个值我们前面已经求得了结果是 0,所以 j 回退到位置 0,用代码表示就是 j=next[j-1]

此时 i 和 j 仍然不匹配,但是 j 已为 0,无法继续回退,所以直接对 next [4] 赋值为 0

实际上,我们在求解模式字符串 ababc 的最长公共前后缀长度的时候,不可避免的会对它的子串进行求解,这样可以方便在不匹配时进行回退,这也是动态规划的思想,要求的结果可以用先前已经求得的结果得出。而我们本身就是要对模式字符串的各个子串进行求解,所以可谓一举两得。

ini 复制代码
public static int[] getNext(String p) {
        int[] next = new int[p.length()];
        int n = p.length();
        int i = 1;
        int j = 0;
        while (i < n) {
            if (p.charAt(i) == p.charAt(j)) {
                next[i] = j + 1;
                i++;
                j++;
            } else {
                if (j > 0) {
                    j = next[j - 1];
                } else {
                    next[i] = 0;
                    i++;
                }
            }
        }
        return next;

    }

对比一下上一段的字符串匹配代码,可以发现,2 段代码的执行逻辑几乎一模一样,不同之处只是我们在比较过程中插入了对 next 数组进行赋值的代码。

可以这样理解,在字符串匹配中,i 指向的是主串,j 指向的子串,比较的是主串和子串。而在求公共前后缀的时候,我们相当于把模式字符串拆成了前后两部分,同样是对这 2 部分进行字符串匹配,相当于自己比自己。

以上就是 KMP 算法的核心原理及实现,当然,实现方式并不是只有我这一种,KMP 的实现算法多种多样,包括生成的 next 数组都会不一样,而我这里的代码只是希望能够把 KMP 的思路表述清楚:

  1. KMP 通过计算模式字符串自身得到 next 数组;
  2. 不匹配时通过 next 数组信息进行回退,而不用再从头开始比较;
  3. 原理就是利用模式字符串的公共前后缀信息,所以算法是奏效的。
相关推荐
秋夫人5 分钟前
B+树(B+TREE)索引
数据结构·算法
梦想科研社41 分钟前
【无人机设计与控制】四旋翼无人机俯仰姿态保持模糊PID控制(带说明报告)
开发语言·算法·数学建模·matlab·无人机
Milo_K42 分钟前
今日 leetCode 15.三数之和
算法·leetcode
Darling_001 小时前
LeetCode_sql_day28(1767.寻找没有被执行的任务对)
sql·算法·leetcode
AlexMercer10121 小时前
【C++】二、数据类型 (同C)
c语言·开发语言·数据结构·c++·笔记·算法
Greyplayground1 小时前
【算法基础实验】图论-BellmanFord最短路径
算法·图论·最短路径
Adolf_19931 小时前
Flask-JWT-Extended登录验证, 不用自定义
后端·python·flask
蓑 羽1 小时前
力扣438 找到字符串中所有字母异位词 Java版本
java·算法·leetcode
源代码:趴菜1 小时前
LeetCode63:不同路径II
算法·leetcode·职场和发展
叫我:松哥1 小时前
基于Python flask的医院管理学院,医生能够增加/删除/修改/删除病人的数据信息,有可视化分析
javascript·后端·python·mysql·信息可视化·flask·bootstrap