设为首页 加入收藏

TOP

VC++2012编程演练数据结构 KMP算法(一)
2014-11-23 17:37:49 】 浏览:362
Tags:2012 编程 演练 数据结构 KMP 算法

KMP算法是通过分析子串,预先计算每个位置发生不匹配的时候,所需GOTO的下一个比较位置,整理出来一个next数组,然后在上面的算法中使用。
Knuth-Morris-Pratt Algorithm,简称KMP算法。
 一种由Knuth(D.E.Knuth)、Morris(J.H.Morris)和Pratt(V.R.Pratt)三人设计的线性时间字符串匹配算法。这个算法不用计算变迁函数δ,匹配时间为Θ(n),只用到辅助函数π[1,m],它是在Θ(m)时间内,根据模式预先计算出来的。数组π使得我们可以按需要,“现场”有效的计算(在平摊意义上来说)变迁函数δ。粗略地说,对任意状态q=0,1,…,m和任意字符a∈Σ,π[q]的值包含了与a无关但在计算δ(q,a)时需要的信息。由于数组π只有m个元素,而δ有Θ(m Σ )个值,所以通过预先计算π而不是δ,使得时间减少了一个Σ因子。[1]
  KMP算法是通过分析子串,预先计算每个位置发生不匹配的时候,所需GOTO的下一个比较位置,整理出来一个next数组,然后在上面的算法中使用。
 当我们分析一个子串时,例如:abcabcddes. 需要分析一下,每个字符x前面最多有多少个连续的字符和字符串从初始位置开始的字符匹配。然后+1就行了(别忘了,我们的字符串都是从索引1开始的)当然,不要相同位置自己匹配,默认第一个字符的匹配数是0。
 设字符串为 x1x2x3...xn,其中x1,x2,x3,... xi,... xn均是字符,设ai为字符xi对应的整数。则a=m,当且仅当满足如下条件:字符串x1x2...xm equals 字符串x(i-m+1)...xi-1 xi 并且x1x2...xm x(m+1) unequals x(i-m) x(i-m+1)...xi-1 xi。
 abcabcddes
  0111234111
  |----------------------默认是0
  --| | |-----------------不能自己在相同位置进行字符匹配,所以这里认为没有匹配字符串,所以0+1 =1,继续从1开始匹配
  ------| | |-----------前面的字符和开始位置的字符相同,所以是2,3,4
  -----------| | | |-------不匹配只能取1。
  希望能明白的是,如果开始字符是 Ch1的话,那么我们就是要在串中第2个Ch1后面的位置开始自己和自己匹配,计算最大的吻合度。
我们创建一个工程

KMP实现如下
[cpp]
/字符串的模式匹配
#include
#include
#include
#include
#include
#define MAXSTRLEN 64
using namespace std;
void GetNext(char T[MAXSTRLEN],int (&next)[64])
{int i,j;i=1;next[1]=j=0;
while(i<(int)T[0])
if(j==0||T[i]==T[j])
{++i;++j;next[i]=j;}
else j=next[j];
for(j=1;j<(int)T[0];j++)
{printf("next[%d]=%-3d",j,next[j]);
if((j)%5==0) printf("\n");}
cout< }
void GetNext(char T[MAXSTRLEN],int (&next)[64],int m)
{int i,j;i=0;next[0]=-1;
for(j=1;j {i=next[j-1];
while(T[j]!=T[i+1]&&i>=0) i=next[i];
if(T[j]==T[i+1]) next[j]=i+1;
else next[j]=-1;}
for(j=0;j<=m;j++)
{printf("next[%d]=%-3d",j,next[j]);
if((j+1)%5==0) printf("\n");}
cout< }
void GetNextVal(char T[MAXSTRLEN],int (&next)[64])
{int i,j;i=1;next[1]=j=0;
while(i<(int)T[0])
if(j==0||T[i]==T[j])
{++i;++j;
if(T[i]!=T[j]) next[i]=j;
else next[i]=next[j];}
else j=next[j];
for(i=1;i<(int)T[0];i++)
{printf("next[%d]=%-3d",i,next[i]);
if(i%5==0) cout< cout< }
int IndexKMP(char S[MAXSTRLEN],char T[MAXSTRLEN],int (&next)[64])
{int i,j,n,m;i=j=1;
n=(int)S[0];m=(int)T[0];
while((i if(j==0||S[i]==T[j]) {++i;++j;}
else j=next[j];
if(j>=m) return i+1-m;
else return 0;
}
int IndexKMP(char S[MAXSTRLEN],char T[MAXSTRLEN],int (&next)[64],int pos)
{int i,j;i=pos;j=1;
while((i<(int)S[0]||j<(int)T[0])&&T[j]!='\0'&&S[i]!='\0')
if(j==0||S[i]==T[j]) {++i;++j;}
else j=next[j];
if(j>=(int)T[0]) return i+1-(int)T[0];
else return 0;
}
int IndexKMP(char S[MAXSTRLEN],char T[MAXSTRLEN],int (&next)[64],int n,int m)
{int i,j;i=j=0;
while(i if(S[i]==T[j]) {++i;++j;}
else if(j==0) i++;
else j=next[j-1]+1;
if(j else return i-m+1;
}
int IndexBF(char S[MAXSTRLEN],char T[MAXSTRLEN],int pos)
{int i,j;i=pos;j=1;
while(i<=S[0]&&j<=T[0])
if(S[i]==T[j]) {++i;++j;}
else {i=i-j+2;j=1;}
if(j>T[0]) return i-T[0];
else return 0;
}


调用如下
[cpp]
void main()
{printf("Findstr.cpp运行结果:\n");
int Index,N,M,next[64]={0};
char s[MAXSTRLEN],t[MAXSTRLEN];
printf("GetNext-IndexKMP的结果:\n");
printf("输入主串s:");gets(s);
printf("输入模式串t:");gets(t);
N=strlen(s);M=strlen(t);
printf("主串s长=%d\n",N);
printf(" 模式串t长=%d\n",M);
GetNext(t,next,M);
Index=IndexKMP(s,t,next,N,M);
if(Index!=-1)
printf("模式串在主串的位置从第%d个字符开始\n",Index);
else print

首页 上一页 1 2 下一页 尾页 1/2/2
】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
上一篇VC++2012编程演练数据结构稀疏矩阵 下一篇VC++2012编程演练数据结构散列文件

最新文章

热门文章

Hot 文章

Python

C 语言

C++基础

大数据基础

linux编程基础

C/C++面试题目