热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

KMP算法实现原理

KMP算法,是由Knuth,Morris,Pratt共同提出的模式匹配算法,其对于任何模式和目标序列,都可以在线性时间内完成匹配查找,而不会发生退化,是一个非常优秀的模式匹配算法。其算法的主要功能就是

KMP算法,是由Knuth,Morris,Pratt共同提出的模式匹配算法,其对于任何模式和目标序列,都可以在线性时间内完成匹配查找,而不会发生退化,是一个非常优秀的模

式匹配算法。其算法的主要功能就是寻找在给定的母串中寻找是否含有一个给定的连续字符串。下面举个例子,如图一所示:


图一:


我们需要在上面的目标串中寻找是否存在一个ABCDABD的字串,在这里我们将图一的长串称为目标串,短串叫做搜索串。图一示例了一般的寻找过程:我们从第一个字符

串开始寻找,如果第一个字符串不匹配,我们就继续检查下一个字符串是否匹配,如图2所示:

图二:


但是假设我们如果在目标串中找到了部分匹配字符,如下图所示:

图三:


从上图可以看出只有D字符与上面的目标串不一样,遇到这样的情况,如果我们还是将比较串向后移动一位的话,那么我们之前的ABCDAB就会重复搜索,这样效率就会降

低。当我们确定D不匹配时,前面6个字符是匹配的,那么我们可不可以利用这个信息来避免将搜索串一位一位向右移动以期待寻找在目标串中和搜索串第一个字符匹配的位

置,因为要确定搜索串是否存在于目标串中的第一步就是要找到两个字符串的公共字符,在上述的例子中就是字符A。那么为了解决这个问题,Knuth,Morris,Pratt利用一张

表来存放一个叫部分匹配值的数组,如下:

图四:


这个表的得到方法我们暂且不讨论,上表是基于“前缀”和“后缀”这两个概念,“前缀”是除了最后一个字符串意外,从第一个字符串开始长度依次递增1的字符串序列;“后

缀”则刚好相反;一下是前缀串和后缀串集合的产生过程:

图五:


那么知道了部分匹配数组是怎么来的,它的作用究竟是什么?在写博客之前,看过不少的人写文章来分析这个算法,但是觉得写得好的可能只有

http://www.ruanyifeng.com/blog/2013/05/Knuth–Morris–Pratt_algorithm.html,这篇博客很详细的介绍了算法的原理和部分匹配值数组的来历。下面就结合这篇博客和个人理解

来说说部分匹配值数组的作用。对于这种问题,我们一般的解决方案是将搜索串ABCDABD相对于目标串向右一位一位地移动,如果找在目标串中找到一个和搜索串的第一个字

符A相等的位置的话,我们会暂停将搜索串继续向右移动,而是继续比较它们的第二个字符串B,如果搜索串从头到尾都存在于目标串中,那么算法结束,否则继续将搜索串向右

移动。但在图三所示的情况,如果继续将搜索串向右一位一位移动以期待寻找在目标串中和搜索串第一个字符匹配的位置(因为要确定搜索串是否存在于目标串中的第一步就是

要找到两个字符串的公共字符)的话效率很低。为了确定搜索串的第一个字符在目标串的下一次出现位置(在常规的算法中,我们向右一位一位移动搜索串就是认为它们的公共字

符A在下一个位置出现),我们利用部分匹配值表来标记。还是上面的例子,我么匹配到D的时候失败了,但ABCDAB这些是匹配的,那么从部分匹配值表中可以知道,ABCDAB

这个符串中最后匹配位B的值为2,然后搜索串向后移动的位数=当前匹配的位数(6)-最后匹配字符的部分匹配值(B,2),这是什么意思????结合前面的前缀和后缀来看,2代表在

ABCDAB中,从第一个字符A开始与包含最后一个字符B为止的所有字符串集合中,它们所形成的最大公共子符串长度为2.移动的位数(4)=当前匹配的位数(6)-最后匹配字符的部

分匹配值(B,2)这个公式就代表将搜索串向后移动4个位置后又可以使目标串和搜索串的第一个字符相等(这样我们就可以后续的比较)。下面就用一张图来解释:


我们可以看到左右方框的字符长度相等且内容相等,那么将左边方框子符串向右移动(黑色区域的总长-方框字符长度)个单位后,就会重合。重合的位置的起点刚好是搜索串

的起点,这个起点的位置也是下一次搜索串的首字符和目标串当前位置相等的字符。所以现在关键的任务是求出部分匹配值表!!!!!

算法实现(c++):

vectorget_partial_match_table(string &s){
int size = s.size(),pre;
vectortab(size,0);
if (size == 1)return tab;
if (s[0] == s[1])tab[1] = 1;
else tab[1] = 0;
for (int i = 2; i if (s[i] == s[tab[i - 1]])tab[i] = tab[i - 1] + 1;
else{
if (s[i] == s[0])tab[i] = 1;
}

}
return tab;
}

算法利用动态规划方法。一开始我么将表的所有元素都置为0,我们还是用ABCDABD为例子:

index 0 1 2 3 4 5 6
char A B C D A B D
value 0 0 0 0 1 2 0

可以看到假设我们i=5,那么tab[i-1]=1,那么截止到i=4它的前缀集合是:
A,AB,ABC,ABCD
它的后缀集合是:
A,DA,CDA,BCDA,因为有一个A是一样的,所以tab[4]=1;
如果这时候我们把B(i=5)加进来,那么对于前缀来说就是增加了一项ABCDA,那么现在的前缀变成:
A,AB,ABC,ABCD,ABCDA;
后缀变成:
B,AB,DAB,CDAB,BCDAB;
那么他们的公共最长字符串就是AB,所以tab[5]=2;
总结来说就是:如果第i-1项的值为n,可以知道截止到第i-1位为止,它们的公共字符串长度为n,那么在新形成的后缀集
合中所有的字串的后面都加第i个字符,且增加一个以第i个字符单个形成的子集(以上面的字符B为例子),在新形成的前缀结

合中增加一个ABCDA即可。

下面是完整代码:

#include
#include
#include
#include
#include
#include
#include
using namespace std;


vectorget_partial_match_table(string &s){
int size = s.size(),pre;
vectortab(size,0);
if (size == 1)return tab;
if (s[0] == s[1])tab[1] = 1;
else tab[1] = 0;
for (int i = 2; i if (s[i] == s[tab[i - 1]])tab[i] = tab[i - 1] + 1;
else{
  if (s[i] == s[0])tab[i] = 1;
}

}
return tab;
}

int main(){
ifstream fin("C:\\Users\\Dell\\Desktop\\data.txt");
string tar, search;
vectortab;
while (fin >> tar >> search){
int size1 = tar.size(), size2 = search.size(), count = 0, begin = -1, pre_i;
int start = 0;//从下标start开始比较
tab = get_partial_match_table(search);
for (int i = 0; i count = 0;
start = 0;
pre_i = i;
while (start <= size2&&istart++;
count++;
i++;
}
if (count == size2){
begin = pre_i;
break;
}
else{//没有完全匹配
if (count>0)i = pre_i + count - tab[count - 1];
else i = pre_i + 1;
}
}
if (begin == -1)cout <<"Not Exist!" <else cout <<"The index of the first match char is: " <tab.clear();
}
return 0;
}

测试用例:

BBBBBBBBBBBB
AAAA
BBBBBB
BB
BBCABCDABABCDABCDABDE
ABCDABD
测试结果:





推荐阅读
  • 本文提供了一个详尽的前端开发资源列表,涵盖了从基础入门到高级应用的各个方面,包括HTML5、CSS3、JavaScript框架及库、移动开发、API接口、工具与插件等。 ... [详细]
  • 本文回顾了作者在求职阿里和腾讯实习生过程中,从最初的迷茫到最后成功获得Offer的心路历程。文中不仅分享了个人的面试经历,还提供了宝贵的面试准备建议和技巧。 ... [详细]
  • 尽管在WPF中工作了一段时间,但在菜单控件的样式设置上遇到了一些基础问题,特别是关于如何正确配置前景色和背景色。 ... [详细]
  • Hanks博士是一位著名的生物技术专家,他的儿子Hankson对数学有着浓厚的兴趣。最近,Hankson遇到了一个有趣的数学问题,涉及求解特定条件下的正整数x,而不使用传统的辗转相除法。 ... [详细]
  • 网络流24题——试题库问题
    题目描述:假设一个试题库中有n道试题。每道试题都标明了所属类别。同一道题可能有多个类别属性。现要从题库中抽取m道题组成试卷。并要求试卷包含指定类型的试题。试设计一个满足要求的组卷算 ... [详细]
  • H5技术实现经典游戏《贪吃蛇》
    本文将分享一个使用HTML5技术实现的经典小游戏——《贪吃蛇》。通过H5技术,我们将探讨如何构建这款游戏的两种主要玩法:积分闯关和无尽模式。 ... [详细]
  • 本文针对HDU 1042 N! 问题提供详细的解析和代码实现。题目要求计算给定整数N(0 ≤ N ≤ 10000)的阶乘N!。文章不仅提供了算法思路,还附上了C++语言的具体实现。 ... [详细]
  • 本文介绍了使用Python和C语言编写程序来计算一个给定数值的平方根的方法。通过迭代算法,我们能够精确地得到所需的结果。 ... [详细]
  • 本文提供了一个关于AC自动机(Aho-Corasick Algorithm)的详细解析与实现方法,特别针对P3796题目进行了深入探讨。文章不仅涵盖了AC自动机的基本概念,还重点讲解了如何通过构建失败指针(fail pointer)来提高字符串匹配效率。 ... [详细]
  • 本文详细介绍如何在SSM(Spring + Spring MVC + MyBatis)框架中实现分页功能。包括分页的基本概念、数据准备、前端分页栏的设计与实现、后端分页逻辑的编写以及最终的测试步骤。 ... [详细]
  • 笔记说明重学前端是程劭非(winter)【前手机淘宝前端负责人】在极客时间开的一个专栏,每天10分钟,重构你的前端知识体系& ... [详细]
  • 【MySQL】frm文件解析
    官网说明:http:dev.mysql.comdocinternalsenfrm-file-format.htmlfrm是MySQL表结构定义文件,通常frm文件是不会损坏的,但是如果 ... [详细]
  • 探讨了在HTML表单中使用元素代替进行表单提交的方法。 ... [详细]
  • 本文详细介绍了如何在最新版本的Xcode中重命名iOS项目,包括项目名称、应用名称及相关的文件夹和配置文件。通过本文,开发者可以轻松完成项目的重命名工作。 ... [详细]
  • 本文探讨了如何在PHP与MySQL环境中实现高效的分页查询,包括基本的分页实现、性能优化技巧以及高级的分页策略。 ... [详细]
author-avatar
mobiledu2502863723
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有