转自:http://blog.csdn.net/ddupd/article/details/19899263

KMP算法详解

KMP算法简介:

KMP算法是一种高效的字符串匹配算法,关于字符串匹配最简单的就是BF算法。BF算法是用两个游标分别指向母串S,模式串T,从开头向后面依次比较字符是否相等,如果相等继续同时向后滑动两个游标,不相等的话,T的游标回溯至开头,S的游标回溯至起初游标的下一位,这种算法原理非常简单,小学生都可以想的到。

KMP算法是在BF算法的基础上加以改进的,它的特点是在遇到字符不匹配时候维持母串T的游标不动,而把模式串向右移动,具体移动到哪一个元素下标,这就是算法的核心思想之处了。

假如母串的i处和模式串的j处不匹配,那么就令k=next(j),表示的意思就是:模式串在j处出现不匹配现象,此时应该将模式串向后一定到下标为k的游标处,在此与之前不匹配的元素进行比较。

Kmp算法的本质:

如图所示:

在下标j处出现不匹配,则k = next(j),表示此时应该把下标k移动到原本j对应的位置处,用T[k]跟s[i]进行对比。如果满足这样的条件,则有T[0],T[1],…T[k-1] = S[i-k],S[i-k+1],…S[i-1]

又因为j之前的字符串跟S都匹配,所以又有T[j-k],T[j-k+1],…T[j-1] = S[i-k],S[i-k+1],…S[i-1].所以得出  T[0],T[1],…T[k-1] = T[j-k],T[j-k+1],…T[j-1]。也就是说图中被标记出来前后两个区域的字符串相等,KMP算法的本质就是找出最大的这样一个k值满足T[0],T[1],…T[k-1] = T[j-k],T[j-k+1],…T[j-1]。

K值的求取方法:

K值的求取用到了数学中的递推的思想,求取K值只跟模式串T自身有关,跟母串S半毛钱关系都没有。先假设已经有 next(j) = k,接下来我们就去求next(j+1)的值。这个要分情况讨论:

如果T[k] = T[j]那么就很容易得到 next(j+1) = k+1 = next(j) + 1;

如果T[k] != T[j],那么此时可以将T[0],T[1],…T[k-1],T[k]看做一个模式串,T[j-k],T[j-k+1],…T[j-1],T[j]看做一个母串,此时模式串在k处出现不匹配现象,那么我们获取next(k)= k1的值,判断T[k1]跟T[j]的值是否相等,如果相等那么next(j+1) = k1+1;如果不相等再往下求新的kn的值,直到T[kn]= T[j],或者遍历到了模式串的开头都不想的话,此时就要把i向后移动一个位置,同时用模式串的开头指向i,或者抽象一点就是把模式串开头的前一位下标(-1)指向i。因为下标(-1)是没有意义的,所以此时等效于下标(0)指向母串的i+1。

算法的实现:

这里一共列出了三个版本的kmp算法,其中第一个是本人根据对算法的理解写的,也是最丑的一个,剩下的两个是改编严蔚敏版的《数据结构与算法》一书中的。

 
 //Algorithms.cpp
#include "Algorithms.h"
#include <vector>
#include <iostream>
using namespace std; Algorithms::Algorithms(void)
{
}
Algorithms::~Algorithms(void)
{
} int Algorithms::kmp1(string strS,string strT)
{
int nSize = strT.size();
vector<int> vecNext(nSize,-);
if (nSize >= )
vecNext[] =;
for (int i=;i<nSize;i++)
{
int k = vecNext[i-];
while (k>= && strT[k] != strT[i-] )
k = vecNext[k];
if(k>= && strT[i-] == strT[k])
vecNext[i] = k + ;
else
vecNext[i] = ;
}
for(int i=;i<nSize;i++)
cout<<"the vector is:"<<i<<": "<<vecNext[i]<<endl; int nLength = strS.size();
int nCount = ;
int nPoss = ;
int nPost = ; while(nPoss < nLength)
{
if ( strS[nPoss] == strT[nPost] )
{
nPoss++;
nPost++;
}
else
{
nPost = vecNext[nPost];
if (nPost == -)
{
nPoss++;
nPost++;
}
} if (nPost == nSize )
{
nCount++;
nPost = ;
}
}
return nCount;
} int Algorithms::kmp2(string strS,string strT)
{
int nSize = strT.size();
vector<int> vecNext(nSize);
int i = ;
vecNext[] = -;
int j = -;
while(i<nSize-)
{
if (j==- || strT[i]==strT[j])
{
i++;
j++;
vecNext[i] = j;
}
else
j = vecNext[j];
}
for(int i=;i<nSize;i++)
cout<<"the vector is:"<<i<<": "<<vecNext[i]<<endl; int nLength = strS.size();
int nCount = ;
int nPoss = ;
int nPost = ; while(nPoss < nLength)
{
if ( strS[nPoss] == strT[nPost] )
{
nPoss++;
nPost++;
}
else
{
nPost = vecNext[nPost];
if (nPost == -)
{
nPoss++;
nPost++;
}
} if (nPost == nSize )
{
nCount++;
nPost = ;
}
}
return nCount;
} int Algorithms::kmp3(string strS,string strT)
{
int nSize = strT.size();
vector<int> vecNext(nSize);
int i = ;
vecNext[] = -;
int j = -;
while(i<nSize-)
{
if (j==- || strT[i]==strT[j])
{
i++;
j++;
if (strT[i] == strT[j])
vecNext[i] =vecNext[j];
else
vecNext[i] = j;
}
else
j = vecNext[j];
}
for(int i=;i<nSize;i++)
cout<<"the vector is:"<<i<<": "<<vecNext[i]<<endl; int nLength = strS.size();
int nCount = ;
int nPoss = ;
int nPost = ; while(nPoss < nLength)
{
if ( strS[nPoss] == strT[nPost] )
{
nPoss++;
nPost++;
}
else
{
nPost = vecNext[nPost];
if (nPost == -)
{
nPoss++;
nPost++;
}
} if (nPost == nSize )
{
nCount++;
nPost = ;
}
}
return nCount;
}
 //main.cpp
#include <iostream>
#include <string>
#include <vector>
#include "Algorithms.h"
using namespace std; void main()
{
string str1,str2;
cout<<"please input str1:"<<endl;
cin>>str1;
cout<<"please input str2:"<<endl;
cin>>str2;
//cout<<"the number of substr in str1 is:"<<Algorithms::kmp1(str1,str2)<<endl;
//cout<<"the number of substr in str1 is:"<<Algorithms::kmp2(str1,str2)<<endl;
cout<<"the number of substr in str1 is:"<<Algorithms::kmp3(str1,str2)<<endl;
system("pause");
}

算法评析:

1:Kmp1

先说下我自己写的吧,代码没有书上的简洁,再说说几个为什么。为什么循环要从i=2开始?

因为要用到int k = vecNext[i-1],以及strT[k] != strT[i-1],如果从i=1开始的话,k的起始值=-1,这样就会出现越界的情况,所以就从i=2开始;另外

next(0)=-1,next(1)=0,这都是毫无疑问的东西,所以可以在前两者已知的情况下,向后求解。

2:Kmp2

Kmp2的巧妙之处在于用了while循环,在需要i变化时才变化,否则已知变换j的值(j表示的是next(i))

3:Kmp3

Kmp3是对kmp2的改进,它考虑到了这样的一种情况:首先在T(i) = T(j)的情况下,如果按照kmp2那么next(i+1)= j+1。但是如果T(i+1)=T(j+1)的话,那么此时就无需在拿T(j+1) 跟母串的比较,因为T(i+1)已经比较过了,并且他们不相等,所以不需要再比较T(j+1),

只需要令:       next(i+1) = next(j+1)。

按照kmp2其实是这样的:  next(i+1) = j+1,j+1 = next(j+1),

所以中间省略的一步就是:next(i+1)= j+1。

总结:

关于kmp网上的讲解很多,但是坑爹的错误程序也不少,在验证一个kmp算法程序是否坑爹,只需要把next数组的内容打印出来,将它和自己口算得到的结果比对。百度百科的C++实现的kmp算法就是错误的,坑爹的。

kmp算法详解的更多相关文章

  1. [转] KMP算法详解

    转载自:http://www.matrix67.com/blog/archives/115 KMP算法详解 如果机房马上要关门了,或者你急着要和MM约会,请直接跳到第六个自然段.    我们这里说的K ...

  2. KMP算法详解(转自中学生OI写的。。ORZ!)

    KMP算法详解 如果机房马上要关门了,或者你急着要和MM约会,请直接跳到第六个自然段. 我们这里说的KMP不是拿来放电影的(虽然我很喜欢这个软件),而是一种算法.KMP算法是拿来处理字符串匹配的.换句 ...

  3. 算法进阶面试题01——KMP算法详解、输出含两次原子串的最短串、判断T1是否包含T2子树、Manacher算法详解、使字符串成为最短回文串

    1.KMP算法详解与应用 子序列:可以连续可以不连续. 子数组/串:要连续 暴力方法:逐个位置比对. KMP:让前面的,指导后面. 概念建设: d的最长前缀与最长后缀的匹配长度为3.(前缀不能到最后一 ...

  4. 数据结构4.3_字符串模式匹配——KMP算法详解

    next数组表示字符串前后缀匹配的最大长度.是KMP算法的精髓所在.可以起到决定模式字符串右移多少长度以达到跳跃式匹配的高效模式. 以下是对next数组的解释: 如何求next数组: 相关链接:按顺序 ...

  5. KMP算法详解&&P3375 【模板】KMP字符串匹配题解

    KMP算法详解: KMP算法是一种改进的字符串匹配算法,由D.E.Knuth,J.H.Morris和V.R.Pratt(雾)提出的. 对于字符串匹配问题(such as 问你在abababb中有多少个 ...

  6. 字符串匹配KMP算法详解

    1. 引言 以前看过很多次KMP算法,一直觉得很有用,但都没有搞明白,一方面是网上很少有比较详细的通俗易懂的讲解,另一方面也怪自己没有沉下心来研究.最近在leetcode上又遇见字符串匹配的题目,以此 ...

  7. KMP算法详解-彻底清楚了(转载+部分原创)

    引言 KMP算法指的是字符串模式匹配算法,问题是:在主串T中找到第一次出现完整子串P时的起始位置.该算法是三位大牛:D.E.Knuth.J.H.Morris和V.R.Pratt同时发现的,以其名字首字 ...

  8. KMP算法详解 --- 彻头彻尾理解KMP算法

    前言 之前对kmp算法虽然了解它的原理,即求出P0···Pi的最大相同前后缀长度k. 但是问题在于如何求出这个最大前后缀长度呢? 我觉得网上很多帖子都说的不是很清楚,总感觉没有把那层纸戳破, 后来翻看 ...

  9. 字符串匹配的KMP算法详解及C#实现

    字符串匹配是计算机的基本任务之一. 举例来说,有一个字符串"BBC ABCDAB ABCDABCDABDE",我想知道,里面是否包含另一个字符串"ABCDABD" ...

随机推荐

  1. Visual Studio Online Integrations-Sync and migration

      原文:http://www.visualstudio.com/zh-cn/explore/vso-integrations-directory-vs

  2. Zhulina 的高分子刷理论

    高分子刷的解析平均场理论有两种表述方式.一个是MWC理论(Macromolecules 1988, 21, 2610-2619),另外一个就是Zhulina和Birshtein这两位俄罗斯老太太的理论 ...

  3. Dual Palindromes

    Dual PalindromesMario Cruz (Colombia) & Hugo Rickeboer (Argentina) A number that reads the same ...

  4. UIImage 和 iOS 图片压缩UIImage / UIImageVIew

    UIImageView 制作气泡 stretchableImageWithLeftCapWidth http://blog.csdn.net/justinjing0612/article/detail ...

  5. 序列化.to_sym

    <%= link_to t('links.list'), showable_paths[@subchannel_item.showable_type.to_sym], target: '_bla ...

  6. Linux / UNIX create soft link with ln command

    How to: Linux / UNIX create soft link with ln command by NIXCRAFT on SEPTEMBER 25, 2007 · 42 COMMENT ...

  7. 自动化运维之puppet的学习(如何找到你需要的模块)

    https://forge.puppetlabs.com/  puppet 模块下载 http://kisspuppet.com/2014/01/14/puppet_forge_modules/ pu ...

  8. Linux 的shell 字符串截取很有用。有八种方法。

    一 Linux 的字符串截取很有用.有八种方法. 假设有变量 var=http://www.linuxidc.com/123.htm 1  # 号截取,删除左边字符,保留右边字符. echo ${va ...

  9. ethtool 命令详解

    命令描述: ethtool 是用于查询及设置网卡参数的命令. 使用概要:ethtool ethx       //查询ethx网口基本设置,其中 x 是对应网卡的编号,如eth0.eth1等等 转自: ...

  10. ENGINE=InnoDB

    最开始用MySQL Administrator建数据库的时候,表缺省是InnoDB类型,也就没有在意.后来用Access2MySQL导数据的时候发现只能导成 MyISAM类型的表 区别如下原来是MyI ...