KMP算法精解及其Python版的代码示例


Posted in Python onJune 01, 2016

KMP算法是经典的字符串匹配算法,解决从字符串S,查找模式字符串M的问题。算法名称来源于发明者Knuth,Morris,Pratt。
假定从字符串S中查找M,S的长度ls,M的长度lm,且(ls > lm)。

朴素的字符串查找方法
从字符串S的第一个字符开始与M进行比较,如果匹配失败。从下一字符开始,重新比较。指导第 (ls - lm) 个字符。
这种方法容易想到并且容易理解,效率不高。
问题在于每次匹配失败后,移动的步伐固定为 1,其实步子可以迈得再大一些。

KMP的字符串查找方法
假定在模式串的连续字串M[0, i] 且 i < lm,已经成功匹配字符串S。但是不巧第 i+1 个字符失败了,怎么办?移动一个字符,重头再来?当然不好,那就是朴素路线了。我们能否从跌倒的地方继续走呢?
既然字串M[0 - i]已经匹配成功,那就从这个子串上做文章。举个栗子     

S序号 j j + 1  j + 2 j + 3 j + 4 j + 5  j+6 j + 7 。。。
S串 a b c a b c d e 。。。
M串 a b c a b d
M序号 0 1 2 3 4 5

此时匹配失败在M串的第5个字符,前4个字符已经匹配成功。
如果从跌倒的地方出发,则需要存在M[0, 4]的子串M[0, k] == S[j+4-k , j+4]。
由于M[0, 4] == S[j ,  j+4] 则有 字串S[j+4-k, j+4] == M[4-k, 4]。综上有M[0, k] == M[4-k, 4]
如果这样的k不存在,那就老老实实的朴素了。
从上面的表格可以直观的看出,下一次匹配只要把M串移动到 j + 3 位置,从 j+5 开始匹配就可以。很容易看出来 在已经匹配成功的字串M[0 , 4]中有最长的子串 (M[0 , 1] == M[3 , 4]),这个就是问题的关键。
因此KMP的核心部分就是计算模式串的各个子串的 k。

实例
首先我们来看一下字符串的朴素匹配.
可以想象成把文本串s固定住,模式串p从s最左边开始对齐,如果对齐的部分完全一样,则匹配成功,失败则将模式串p整体往右移1位,继续检查对齐部分,如此反复.

#朴素匹配 
def naive_match(s, p): 
 m = len(s); n = len(p) 
 for i in range(m-n+1):#起始指针i 
  if s[i:i+n] == p: 
   return True 
 return False

关于kmp算法,讲的最好的当属阮一峰的<字符串匹配的KMP算法>.一路读下来,豁然开朗.
其实就是,对模式串p进行预处理,得到前后缀的部分匹配表,使得我们可以借助已知信息,算出可以右移多少位.即 kmp = 朴素匹配 + 移动多位.
更多细节请看阮一峰的文章,这里就不展开了.
下面给出python的代码实现.

#KMP 
def kmp_match(s, p): 
 m = len(s); n = len(p) 
 cur = 0#起始指针cur 
 table = partial_table(p) 
 while cur<=m-n: 
  for i in range(n): 
   if s[i+cur]!=p[i]: 
    cur += max(i - table[i-1], 1)#有了部分匹配表,我们不只是单纯的1位1位往右移,可以一次移动多位 
    break 
  else: 
   return True 
 return False 
 
#部分匹配表 
def partial_table(p): 
 '''''partial_table("ABCDABD") -> [0, 0, 0, 0, 1, 2, 0]''' 
 prefix = set() 
 postfix = set() 
 ret = [0] 
 for i in range(1,len(p)): 
  prefix.add(p[:i]) 
  postfix = {p[j:i+1] for j in range(1,i+1)} 
  ret.append(len((prefix&postfix or {''}).pop())) 
 return ret 
 
print naive_match("BBC ABCDAB ABCDABCDABDE", "ABCDABD") 
print partial_table("ABCDABD") 
print kmp_match("BBC ABCDAB ABCDABCDABDE", "ABCDABD")

Python 相关文章推荐
python网络编程学习笔记(四):域名系统
Jun 09 Python
python根据给定文件返回文件名和扩展名的方法
Mar 27 Python
提升Python程序运行效率的6个方法
Mar 31 Python
python执行子进程实现进程间通信的方法
Jun 02 Python
Windows下Eclipse+PyDev配置Python+PyQt4开发环境
May 17 Python
python中模块查找的原理与方法详解
Aug 11 Python
Python+tkinter使用40行代码实现计算器功能
Jan 30 Python
python邮件发送smtplib使用详解
Jun 16 Python
Python 分发包中添加额外文件的方法
Aug 16 Python
Django接收自定义http header过程详解
Aug 23 Python
python如何实现单链表的反转
Feb 10 Python
Python爬虫获取豆瓣电影并写入excel
Jul 31 Python
Python缩进和冒号详解
Jun 01 #Python
Python注释详解
Jun 01 #Python
深入理解python try异常处理机制
Jun 01 #Python
python学习 流程控制语句详解
Jun 01 #Python
python+Django+apache的配置方法详解
Jun 01 #Python
python中函数默认值使用注意点详解
Jun 01 #Python
Python中基础的socket编程实战攻略
Jun 01 #Python
You might like
php select,radio和checkbox默认选择的实现方法
2010/05/15 PHP
用php来改写404错误页让你的页面更友好
2013/01/24 PHP
php生成mysql的数据字典
2016/07/07 PHP
PHP实现按之字形顺序打印二叉树的方法
2018/01/16 PHP
php 可变函数使用小结
2018/06/12 PHP
Extjs ajax同步请求时post方式参数发送方式
2009/08/05 Javascript
jQuery EasyUI API 中文文档 - Draggable 可拖拽
2011/09/29 Javascript
js判断输入是否为数字的具体实例
2013/08/03 Javascript
javascript的数组和常用函数详解
2014/05/09 Javascript
Node.js安装教程和NPM包管理器使用详解
2014/08/16 Javascript
jquery使用remove()方法删除指定class子元素
2015/03/26 Javascript
javascript实现倒计时(精确到秒)
2015/06/26 Javascript
js实现多图左右切换功能
2016/08/04 Javascript
Javascript实现代码折叠功能
2016/08/25 Javascript
JavaScript 计算笛卡尔积实例详解
2016/12/02 Javascript
vue2.0获取自定义属性的值
2017/03/28 Javascript
node.js利用mongoose获取mongodb数据的格式化问题详解
2017/10/06 Javascript
jQuery实现动态添加节点与遍历节点功能示例
2017/11/09 jQuery
JS实现的简单下拉框联动功能示例
2018/05/11 Javascript
jQuery实现购物车的总价计算和总价传值功能
2018/11/28 jQuery
Nodejs模块的调用操作实例分析
2018/12/25 NodeJs
JS实现页面侧边栏效果探究
2021/01/08 Javascript
Python入门必须知道的11个知识点
2018/03/21 Python
mac下给python3安装requests库和scrapy库的实例
2018/06/13 Python
Sanic框架配置操作分析
2018/07/17 Python
CSS3 对过渡(transition)进行调速以及延时
2020/10/21 HTML / CSS
Omio葡萄牙:全欧洲低价大巴、火车和航班搜索和比价
2019/02/09 全球购物
Speedo速比涛法国官方网站:泳衣、泳镜、泳帽、泳裤
2019/07/30 全球购物
高中生自我鉴定范文
2013/10/30 职场文书
迎接领导欢迎词
2014/01/11 职场文书
放飞蜻蜓反思
2014/02/05 职场文书
村党的群众路线教育实践活动工作总结
2014/10/25 职场文书
幼儿园教师考核评语
2014/12/31 职场文书
故意伤害罪辩护词
2015/05/21 职场文书
大学生入党自传2015
2015/06/26 职场文书
欧元符号 €
2022/02/17 杂记