KMP匹配的模式算法—，源码+原理+手写框架

4)当j=4时,j由1到j-1的串是“aaa”，由于前缀字符“aa”与后缀“aa”相

等，所以next[4]=3;

5. ………

6)当j=9时,j由1到j-1的串是“aaaaaaaa”，由于前缀字符“aaaaaaa”与后缀“aaaaaaa”相等,所以next[9]=8。

3，算法实现

---

/* 通过计算返回子串T的next数组。 */

void get_next(String T, int *next)

{

int i,k;

i=1;

k=0;

next[1]=0;

while (i<T[0]) /* 此处T[0]表示串T的长度 */

{

if(k0 || T[i] T[k])

{

++i;

++k;

next[i] = k;

}

else

k= next[k]; /* 若字符不相同，则k值回溯 */

}

}

/* 返回子串T在主串S中第pos个字符之后的位置。若不存在，则函数返回值为0。 */

/* T非空，1≤pos≤StrLength(S)。 */

int Index_KMP(String S, String T, int pos)

{

int i = pos; /* i用于主串S中当前位置下标值，若pos不为1，则从pos位置开始匹配 */

int j = 1; /* j用于子串T中当前位置下标值 */

int next[255]; /* 定义一next数组 */

get_next(T, next); /* 对串T作分析，得到next数组 */

while (i <= S[0] && j <= T[0]) /* 若i小于S的长度并且j小于T的长度时，循环继续 */

{

if (j==0 || S[i] == T[j]) /* 两字母相等则继续，与朴素算法增加了j=0判断 */

{

++i;

++j;

}

else /* 指针后退重新开始匹配 */

j = next[j];/* j退回合适的位置，i值不变 */

}

if (j > T[0])

return i-T[0];

else

4，算法改进

---

对于get_next函数来说，若T的长度为m，因只涉及到简单的单循环，其时间复杂度为0(m),而由于主值的不回溯，使得 index_KMP 算法效率得到了提高，while循环的时间复杂度为O(n)。因此，整个算法的时间复杂度为O(n+m)。相较于朴素模式匹配算法的O((n-m+1)*m)来说,是要好一些。

这里也需要强调，KMP算法仅当模式与主串之间存在许多“部分匹配”的情况下才体现出它的优势,否则两者差异并不明显。

后来有人发现，KMP还是有缺陷的。比如，如果我们的主串 S="aaaabcde”，子串T=“aaaaax””，其next 数组值分别为012345，在开始时，当 i=5、j=5时，我们发现b”与“a”不相等，如图的①，因此j=next[5]=4，如图中的②，此时“b”与第4位置的“a”依然不等，j=next[4]=3，如图中的③，后依次是④⑥，直到j=next[1]=0时,根据算法，此时i++、j++，得到 i=6、j=1，如图中的⑥。

我们发现，当中的②③④⑤步骤，其实是多余的判断。由于T串的第二、三、四、五位置的字符都与首位的“a”相等，那么可以用首位 next[1]的值去取代与它相等的字符后续next[i]的值,这是个很好的办法。因此我们对求next函数进行了改良。

假设取代的数组为nextval，增加了加粗部分，代码如下:

/* 求模式串T的next函数修正值并存入数组nextval */

void get_nextval(String T, int *nextval)

{

int i,k;

i=1;

k=0;

nextval[1]=0;

while (i<T[0]) /* 此处T[0]表示串T的长度 */

{

if(k0 || T[i] T[k]) /* T[i]表示后缀的单个字符，T[k]表示前缀的单个字符 */

{

++i;

++k;

if (T[i]!=T[k]) /* 若当前字符与前缀字符不同 */

nextval[i] = k; /* 则当前的j为nextval在i位置的值 */

else

nextval[i] = nextval[k]; /* 如果与前缀字符相同，则将前缀字符的 */

/* nextval值赋值给nextval在i位置的值 */

}

else

k= nextval[k]; /* 若字符不相同，则k值回溯 */

}

}

四，全部代码

==========

#include “string.h”

#include “stdio.h”

#include “stdlib.h”

#include “math.h”

#include “time.h”

#define OK 1

#define ERROR 0

#define TRUE 1

#define FALSE 0

#define MAXSIZE 100 /* 存储空间初始分配量 */

typedef int Status; /* Status是函数的类型,其值是函数结果状态代码，如OK等 */

typedef int ElemType; /* ElemType类型根据实际情况而定，这里假设为int */

typedef char String[MAXSIZE+1]; /* 0号单元存放串的长度 */

/* 生成一个其值等于chars的串T */

Status StrAssign(String T,char *chars)

{

int i;

if(strlen(chars)>MAXSIZE)

return ERROR;

else

{

T[0]=strlen(chars);

for(i=1;i<=T[0];i++)

T[i]=*(chars+i-1);

return OK;

}

}

Status ClearString(String S)

{

S[0]=0;/* 令串长为零 */

return OK;

}

/* 输出字符串T。 */

void StrPrint(String T)

{

int i;

for(i=1;i<=T[0];i++)

printf(“%c”,T[i]);

printf(“\n”);

}

/* 输出Next数组值。 */

void NextPrint(int next[],int length)

{

int i;

for(i=1;i<=length;i++)

printf(“%d”,next[i]);

printf(“\n”);

}

/* 返回串的元素个数 */

int StrLength(String S)

{

return S[0];

}

/* 朴素的模式匹配法 */

int Index(String S, String T, int pos)

{

int i = pos; /* i用于主串S中当前位置下标值，若pos不为1，则从pos位置开始匹配 */

int j = 1; /* j用于子串T中当前位置下标值 */

while (i <= S[0] && j <= T[0]) /* 若i小于S的长度并且j小于T的长度时，循环继续 */

{

if (S[i] == T[j]) /* 两字母相等则继续 */

{

++i;

++j;

}

else /* 指针后退重新开始匹配 */

{

i = i-j+2; /* i退回到上次匹配首位的下一位 */

j = 1; /* j退回到子串T的首位 */

}

}

if (j > T[0])

return i-T[0];

else

return 0;

}

/* 通过计算返回子串T的next数组。 */

void get_next(String T, int *next)

{

int i,k;

i=1;

k=0;

next[1]=0;

while (i<T[0]) /* 此处T[0]表示串T的长度 */

{

if(k0 || T[i] T[k])

{

++i;

++k;

next[i] = k;

}

else

k= next[k]; /* 若字符不相同，则k值回溯 */

}

}

/* 返回子串T在主串S中第pos个字符之后的位置。若不存在，则函数返回值为0。 */

/* T非空，1≤pos≤StrLength(S)。 */

int Index_KMP(String S, String T, int pos)

{

int i = pos; /* i用于主串S中当前位置下标值，若pos不为1，则从pos位置开始匹配 */

int j = 1; /* j用于子串T中当前位置下标值 */

int next[255]; /* 定义一next数组 */

get_next(T, next); /* 对串T作分析，得到next数组 */

while (i <= S[0] && j <= T[0]) /* 若i小于S的长度并且j小于T的长度时，循环继续 */

{

if (j==0 || S[i] == T[j]) /* 两字母相等则继续，与朴素算法增加了j=0判断 */

{

++i;

++j;

}

else /* 指针后退重新开始匹配 */

j = next[j];/* j退回合适的位置，i值不变 */

}

if (j > T[0])

return i-T[0];

else

return 0;

}

/* 求模式串T的next函数修正值并存入数组nextval */

void get_nextval(String T, int *nextval)

{

int i,k;

i=1;

k=0;

nextval[1]=0;

while (i<T[0]) /* 此处T[0]表示串T的长度 */

{

if(k0 || T[i] T[k]) /* T[i]表示后缀的单个字符，T[k]表示前缀的单个字符 */

{

++i;

++k;

if (T[i]!=T[k]) /* 若当前字符与前缀字符不同 */

nextval[i] = k; /* 则当前的j为nextval在i位置的值 */

else

nextval[i] = nextval[k]; /* 如果与前缀字符相同，则将前缀字符的 */

/* nextval值赋值给nextval在i位置的值 */

}

else

k= nextval[k]; /* 若字符不相同，则k值回溯 */

}

}

int Index_KMP1(String S, String T, int pos)

{

int i = pos; /* i用于主串S中当前位置下标值，若pos不为1，则从pos位置开始匹配 */

int j = 1; /* j用于子串T中当前位置下标值 */

int next[255]; /* 定义一next数组 */

get_nextval(T, next); /* 对串T作分析，得到next数组 */

while (i <= S[0] && j <= T[0]) /* 若i小于S的长度并且j小于T的长度时，循环继续 */

{

if (j==0 || S[i] == T[j]) /* 两字母相等则继续，与朴素算法增加了j=0判断 */

{

++i;

++j;

}

else /* 指针后退重新开始匹配 */

j = next[j];/* j退回合适的位置，i值不变 */

}

if (j > T[0])

return i-T[0];

else

return 0;

}

int main()

{

int i,*p;

String s1,s2;

StrAssign(s1,“abcdex”);

printf("子串为: ");

StrPrint(s1);

i=StrLength(s1);

p=(int*)malloc((i+1)*sizeof(int));

get_next(s1,p);

printf("Next为: ");

NextPrint(p,StrLength(s1));

printf(“\n”);

StrAssign(s1,“abcabx”);

printf("子串为: ");

StrPrint(s1);

i=StrLength(s1);

p=(int*)malloc((i+1)*sizeof(int));

get_next(s1,p);

printf("Next为: ");

NextPrint(p,StrLength(s1));

printf(“\n”);

小编13年上海交大毕业，曾经在小公司待过，也去过华为、OPPO等大厂，18年进入阿里一直到现在。

深知大多数初中级Java工程师，想要提升技能，往往是自己摸索成长，但自己不成体系的自学效果低效又漫长，而且极易碰到天花板技术停滞不前！

因此收集整理了一份《2024年最新Java开发全套学习资料》送给大家，初衷也很简单，就是希望能够帮助到想自学提升又不知道该从何学起的朋友，同时减轻大家的负担。

由于文件比较大，这里只是将部分目录截图出来，每个节点里面都包含大厂面经、学习笔记、源码讲义、实战项目、讲解视频

如果你觉得这些内容对你有帮助，可以添加下面V无偿领取！（备注Java）

)malloc((i+1)*sizeof(int));

get_next(s1,p);

printf("Next为: ");

NextPrint(p,StrLength(s1));

printf(“\n”);

小编13年上海交大毕业，曾经在小公司待过，也去过华为、OPPO等大厂，18年进入阿里一直到现在。

深知大多数初中级Java工程师，想要提升技能，往往是自己摸索成长，但自己不成体系的自学效果低效又漫长，而且极易碰到天花板技术停滞不前！

因此收集整理了一份《2024年最新Java开发全套学习资料》送给大家，初衷也很简单，就是希望能够帮助到想自学提升又不知道该从何学起的朋友，同时减轻大家的负担。
[外链图片转存中…(img-vZcB2IFY-1710830144145)]
[外链图片转存中…(img-7unDJRA4-1710830144147)]
[外链图片转存中…(img-pM8w9AKt-1710830144147)]

由于文件比较大，这里只是将部分目录截图出来，每个节点里面都包含大厂面经、学习笔记、源码讲义、实战项目、讲解视频

如果你觉得这些内容对你有帮助，可以添加下面V无偿领取！（备注Java）
[外链图片转存中…(img-LlJRYOXY-1710830144148)]