KMP算法详解及各种应用

分享一下我老师大神的人工智能教程!零基础,通俗易懂!http://blog.csdn.net/jiangjunshow

也欢迎大家转载本篇文章。分享知识,造福人民,实现我们中华民族伟大复兴!

               

KMP算法详解:
KMP算法之所以叫做KMP算法是因为这个算法是由三个人共同提出来的,就取三个人名字的首字母作为该算法的名字。其实KMP算法与BF算法的区别就在于KMP算法巧妙的消除了指针i的回溯问题,只需确定下次匹配j的位置即可,使得问题的复杂度由O(mn)下降到O(m+n)。
在KMP算法中,为了确定在匹配不成功时,下次匹配时j的位置,引入了next[]数组,next[j]的值表示P[0...j-1]中最长后缀的长度等于相同字符序列的前缀。
对于next[]数组的定义如下:
1) next[j]=-1  j=0
2) next[j]=max k:0<k<j P[0...k-1]=P[j-k,j-1]
3) next[j]=0  其他
如:
P      a    b   a    b   a
j       0   1    2   3   4
next -1  0    0   1   2
即next[j]=k>0时,表示P[0...k-1]=P[j-k,j-1]
因此KMP算法的思想就是:在匹配过程称,若发生不匹配的情况,如果next[j]>=0,则目标串的指针i不变,将模式串的指针j移动到next[j]的位置继续进行匹配;若next[j]=-1,则将i右移1位,并将j置0,继续进行比较。
代码实现如下:

int KMPMatch(char *s,char *p){    int next[100];    int i , j;    i = 0;    j = 0;    getNext(p , next);    while(i < strlen(s))    {        if(j == -1 || s[i] == p[j])        {            i++;            j++;        }        else        {            j = next[j];       //消除了指针i的回溯        }        if(j == strlen(p))            return i - strlen(p);    }    return -1;}
因此KMP算法的关键在于求算next[]数组的值,即求算模式串每个位置处的最长后缀与前缀相同的长度, 而求算next[]数组的值有两种思路,第一种思路是用递推的思想去求算,还有一种就是直接去求解。
  1、按照递推的思想:
   根据定义next[0]=-1,假设next[j]=k, 即P[0...k-1]==P[j-k,j-1]
   1)若P[j]==P[k],则有P[0..k]==P[j-k,j],很显然,next[j+1]=next[j]+1=k+1;
   2)若P[j]!=P[k],则可以把其看做模式匹配的问题,即匹配失败的时候,k值如何移动,显然k=next[k]。
   因此可以这样去实现:
void getNext(char *p,int *next){    int j,k;    next[0] = -1;    j = 0;    k = -1;    while(j < strlen(p) - 1)    {        if(k == -1 || p[j] == p[k])    //匹配的情况下,p[j]==p[k]        {            j++;            k++;            next[j] = k;        }        else                   //p[j]!=p[k]            k = next[k];    }}
   2、直接求解方法
void getNext(char *p,int *next){    int i , j , temp;    for(i = 0 ; i < strlen(p) ; ++i)    {        if(i == 0)        {            next[i] = -1;     //next[0]=-1        }        else if(i == 1)        {            next[i] = 0;      //next[1]=0        }        else        {            temp = i - 1;            for(j = temp ; j > 0 ; --j)            {                if( equals(p , i , j) )                {                    next[i] = j;   //找到最大的k值                    break;                }            }            if(j == 0)                next[i] = 0;        }    }} bool equals(char *p,int i,int j)     //判断p[0...j-1]与p[i-j...i-1]是否相等 {    int k = 0;    int s = i - j;    for( ; k <= j - 1 && s <= i - 1 ; k++ , s++)    {        if(p[k] != p[s])            return false;    }    return true;}
http://poj.org/problem?id=2406

给定一个字符串,问最多是多少个相同子串不重叠连接构成。

KMP的next数组应用。这里主要是如何判断是否有这样的子串,和子串的个数。

若为abababa,显然除其本身外,没有子串满足条件。而分析其next数组,next[7] = 5,next[5] = 3,next[3] = 1,即str[2..7]可由ba子串连接构成,那怎么否定这样的情况呢?很简单,若该子串满足条件,则len%sublen必为0。sunlen可由上面的分析得到为len-next[len]。

因为子串是首尾相接,len/sublen即为substr的个数。

若L%(L-next[L])==0,n = L/(L-next[L]),else n = 1

扫描二维码关注公众号,回复: 4077784 查看本文章
#include<iostream>#include<cstdio>using namespace std;char pattern[1000002];int next[1000002];/*kmp算法,需要首先求出模式串的next函数值next[j] = k,说明 p0pk-1 == pj-kpj-1,也就是说k为其前面相等串的长度*/void get_nextval(const char* pattern)int i=0,j=-1; next[0]= -1while(pattern[i] != '\0') {  if(j== -1 || pattern[i]== pattern[j] )     //pattern[i]表示后缀的单个字符,pattern[j]表示前缀的单个字符  {   ++i;   ++j;   if(pattern[i] != pattern[j])    next[i]=j;   else    next[i]=next[j];  }  else   j=next[j];    //若j值不相同,则j值回溯 }}//get_nextvalint main(void)int len; while(scanf("%s",pattern)!=EOF) {  if(pattern[0]=='.')   break;  len=strlen(pattern);  get_nextval(pattern);       if(len%(len-next[len])==0)   printf("%d\n",len/(len-next[len]));  else   printf("1\n");    } return 0;}

http://poj.org/problem?id=1961

大意:
   定义字符串A,若A最多由n个相同字串s连接而成,则A=s^n,如"aaa" = "a"^3,"abab" = "ab"^2
   "ababa" = "ababa"^1
  
给出一个字符串A,求该字符串的所有前缀中有多少个前缀SA= s^n(n>1)
输出符合条件的前缀长度及其对应的n

  如aaa
  前缀aa的长度为2,由2个'a'组成
  前缀aaa的长度为3,由3个"a"组成

  分析:KMP
  若某一长度L的前缀符合上诉条件,则
    1.next[L]!=0(next[L]=0时字串为原串,不符合条件)
 2.L%(L-next[L])==0(此时字串的长度为L/next[L])

 对于2:有str[0]....str[next[L]-1]=str[L-next[L]-1]...str[L-1]
        =》str[L-next[L]-1] = str[L-next[L]-1+L-next[L]-1] = str[2*(L-next[L]-1)];
  假设S = L-next[L]-1;则有str[0]=str[s]=str[2*s]=str[3*s]...str[k*s],对于所有i%s==0,均有s[i]=s[0]
  同理,str[1]=str[s+1]=str[2*s+1]....
        str[j]=str[s+j]=str[2*s+j]....
  综上,若L%S==0,则可得L为str[0]...str[s-1]的相同字串组成,
  总长度为L,其中字串长度SL = s-0+1=L-next[L],循环次数为L/SL
   故对于所有大于1的前缀,只要其符合上述条件,即为答案之一

#include<iostream>#include<cstdio>#include<cstring>using namespace std;char pattern[1000002];int next[1000002];/*kmp算法,需要首先求出模式串的next函数值next[j] = k,说明 p0pk-1 == pj-kpj-1,也就是说k为其前面相等串的长度*/void get_nextval(const char* pattern)int i=0,j=-1; next[0]= -1while(pattern[i] != '\0') {  if(j== -1 || pattern[i]== pattern[j] )  {   ++i;   ++j;   next[i]=j;  }  else   j=next[j]; }}//get_nextvalint main(void)int i,len,n,j=1while(scanf("%d",&n)!=EOF) {  if(!n)   break;  scanf("%s",pattern);  len=strlen(pattern);  get_nextval(pattern);   printf("Test case #%d\n",j++);  for(i=2;i<=len;i++)  {   if(i%(i-next[i])==0 && i/(i-next[i])>1)    printf("%d %d\n",i,i/(i-next[i]));  }  printf("\n");    } return 0;}

http://poj.org/problem?id=2752
 大意:
给出一个字符串A,求A有多少个前缀同时也是后缀,从小到大输出这些前缀的长度。

分析:KMP
对于长度为len的字符串,由next的定义知:
A[0]A[1]...A[next[len]-1]=A[len-next[len]]...A[len-1]此时A[0]A[1]...A[next[len]-1]为一个符合条件的前缀
有A[0]A[1]....A[next[next[len]]-1] = A[len-next[next[len] - next[next[len]]]...A[next[len]-1],故A[0]A[1]....A[next[next[len]]-1]也是一个符合条件的前缀
故从len=>next[len]=>next[next[len]] ....=>直到某个next[]为0均为合法答案,注意当首位单词相同时,也为答案。

#include<iostream>#include<cstdio>#include<cstring>#include<vector>using namespace std;char pattern[400002];int next[400002];/*kmp算法,需要首先求出模式串的next函数值next[j] = k,说明 p0pk-1 == pj-kpj-1,也就是说k为其前面相等串的长度*/void get_nextval(const char* pattern)int i=0,j=-1; next[0]= -1while(pattern[i] != '\0') {  if(j== -1 || pattern[i]== pattern[j] )  {   ++i;   ++j;   next[i]=j;  }  else   j=next[j]; }}//get_nextvalint main(void)int i,len,n; vector<int>ans; while(scanf("%s",pattern)!=EOF) {  ans.clear();  len=strlen(pattern);  get_nextval(pattern);  n=len;  while(n)  {   ans.push_back(n);   n=next[n];  }  if(pattern[0]==pattern[n-1])   //首部、尾部字符相同   ans.push_back(1);  i=ans.size()-1;  for(;i>0;i--)   printf("%d ",ans[i]);  printf("%d\n",ans[0]); } return 0;}




 

           

给我老师的人工智能教程打call!http://blog.csdn.net/jiangjunshow

这里写图片描述
你好! 这是你第一次使用 **Markdown编辑器** 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章,了解一下Markdown的基本语法知识。

新的改变

我们对Markdown编辑器进行了一些功能拓展与语法支持,除了标准的Markdown编辑器功能,我们增加了如下几点新功能,帮助你用它写博客:

  1. 全新的界面设计 ,将会带来全新的写作体验;
  2. 在创作中心设置你喜爱的代码高亮样式,Markdown 将代码片显示选择的高亮样式 进行展示;
  3. 增加了 图片拖拽 功能,你可以将本地的图片直接拖拽到编辑区域直接展示;
  4. 全新的 KaTeX数学公式 语法;
  5. 增加了支持甘特图的mermaid语法1 功能;
  6. 增加了 多屏幕编辑 Markdown文章功能;
  7. 增加了 焦点写作模式、预览模式、简洁写作模式、左右区域同步滚轮设置 等功能,功能按钮位于编辑区域与预览区域中间;
  8. 增加了 检查列表 功能。

功能快捷键

撤销:Ctrl/Command + Z
重做:Ctrl/Command + Y
加粗:Ctrl/Command + B
斜体:Ctrl/Command + I
标题:Ctrl/Command + Shift + H
无序列表:Ctrl/Command + Shift + U
有序列表:Ctrl/Command + Shift + O
检查列表:Ctrl/Command + Shift + C
插入代码:Ctrl/Command + Shift + K
插入链接:Ctrl/Command + Shift + L
插入图片:Ctrl/Command + Shift + G

合理的创建标题,有助于目录的生成

直接输入1次#,并按下space后,将生成1级标题。
输入2次#,并按下space后,将生成2级标题。
以此类推,我们支持6级标题。有助于使用TOC语法后生成一个完美的目录。

如何改变文本的样式

强调文本 强调文本

加粗文本 加粗文本

标记文本

删除文本

引用文本

H2O is是液体。

210 运算结果是 1024.

插入链接与图片

链接: link.

图片: Alt

带尺寸的图片: Alt

当然,我们为了让用户更加便捷,我们增加了图片拖拽功能。

如何插入一段漂亮的代码片

博客设置页面,选择一款你喜欢的代码片高亮样式,下面展示同样高亮的 代码片.

// An highlighted block var foo = 'bar'; 

生成一个适合你的列表

  • 项目
    • 项目
      • 项目
  1. 项目1
  2. 项目2
  3. 项目3
  • 计划任务
  • 完成任务

创建一个表格

一个简单的表格是这么创建的:

项目 Value
电脑 $1600
手机 $12
导管 $1

设定内容居中、居左、居右

使用:---------:居中
使用:----------居左
使用----------:居右

第一列 第二列 第三列
第一列文本居中 第二列文本居右 第三列文本居左

SmartyPants

SmartyPants将ASCII标点字符转换为“智能”印刷标点HTML实体。例如:

TYPE ASCII HTML
Single backticks 'Isn't this fun?' ‘Isn’t this fun?’
Quotes "Isn't this fun?" “Isn’t this fun?”
Dashes -- is en-dash, --- is em-dash – is en-dash, — is em-dash

创建一个自定义列表

Markdown
Text-to- HTML conversion tool
Authors
John
Luke

如何创建一个注脚

一个具有注脚的文本。2

注释也是必不可少的

Markdown将文本转换为 HTML

KaTeX数学公式

您可以使用渲染LaTeX数学表达式 KaTeX:

Gamma公式展示 Γ ( n ) = ( n 1 ) ! n N \Gamma(n) = (n-1)!\quad\forall n\in\mathbb N 是通过欧拉积分

Γ ( z ) = 0 t z 1 e t d t &ThinSpace; . \Gamma(z) = \int_0^\infty t^{z-1}e^{-t}dt\,.

你可以找到更多关于的信息 LaTeX 数学表达式here.

新的甘特图功能,丰富你的文章

gantt
        dateFormat  YYYY-MM-DD
        title Adding GANTT diagram functionality to mermaid
        section 现有任务
        已完成               :done,    des1, 2014-01-06,2014-01-08
        进行中               :active,  des2, 2014-01-09, 3d
        计划一               :         des3, after des2, 5d
        计划二               :         des4, after des3, 5d
  • 关于 甘特图 语法,参考 这儿,

UML 图表

可以使用UML图表进行渲染。 Mermaid. 例如下面产生的一个序列图::

张三 李四 王五 你好!李四, 最近怎么样? 你最近怎么样,王五? 我很好,谢谢! 我很好,谢谢! 李四想了很长时间, 文字太长了 不适合放在一行. 打量着王五... 很好... 王五, 你怎么样? 张三 李四 王五

这将产生一个流程图。:

链接
长方形
圆角长方形
菱形
  • 关于 Mermaid 语法,参考 这儿,

FLowchart流程图

我们依旧会支持flowchart的流程图:

  • 关于 Flowchart流程图 语法,参考 这儿.

导出与导入

导出

如果你想尝试使用此编辑器, 你可以在此篇文章任意编辑。当你完成了一篇文章的写作, 在上方工具栏找到 文章导出 ,生成一个.md文件或者.html文件进行本地保存。

导入

如果你想加载一篇你写过的.md文件或者.html文件,在上方工具栏可以选择导入功能进行对应扩展名的文件导入,
继续你的创作。


  1. mermaid语法说明 ↩︎

  2. 注脚的解释 ↩︎

猜你喜欢

转载自blog.csdn.net/fguihbfg/article/details/84096790
今日推荐