KMP算法的实现

 

目录

问题描述

数据结构

算法设计 

算法流程图

 程序实现

运行结果 

---

问题描述

求一个字符串在另一个字符串中出现的位置。要求：利用键盘输入两个字符串，一个设定为主串，另一个设定为子串，对这两个字符串应用KMP算法，求出子串在主串中第一次出现的位置。

数据结构

本设计使用串的顺序存储结构来实现。

算法设计 

当子串中每个字符依次和主串中的一个连续的字符序列相等，则成为匹配成功，反之，成为匹配不成功。程序中函数GetNext()是求出模式串t的next函数值并存入数组next中，函数IndexKmp()为模式匹配函数，是利用模式串的next函数求t在主串s中第pos个字符之前的位置。

当某个位置匹配不成功的时候，应该从子串的下一个位置开始新的比较。将这个位置的值存放在next数组中，其中next数组中的元素满足条件next[j]=k,表示当子串中的第j+1个字符发生匹配不成功的情况下，应该从子串的第k+1个字符开始新的匹配。如果已经得到了next数组，匹配可进行如下：

1. 将指针i、j分别指向主串s和模式串t中比较字符，初值i=pos、j=1;
2. 如果si=tj,则++i、++j顺次比较后面的字符；
3. 如si<>tj,则指针i不动，指针j退到next[j]位置再比较。然后指针i和指针j所指向的字符按此种方法继续比较，直到j=m-1,即在主串s中中找到模式串t为止。

     next函数的编写是整个算法的核心。这里利用递推思想设计next函数。

1. 令next[0]=-1(next[j]=-1时，证明字符串匹配要从模式串的第0个字符开始，且第0个字符并不和主串的第i个字符相等，i指针向前移动）。
2. 假设next[j]=k,说明T[0~k-1]=T[j-k~j-1]。
3. 现在求next[j+1]

①当T[j]=T[k]时，说明T[0~k]=T[j-k~j]，这时分两种情况讨论：

a:当T[j+1]!=T[k+1]时，显然next[j+1]=k+1;

b:当T[j+1]=T[k+1]时，说明T[k+1]和T[j+1]一样，都不和主串的字符相匹配，因此m=k+1,j=next[m],直到T[m]!=T[j+1],next[j+1]=m。

②当T[j]!=T[k]时，必须在T[0~k-1]中找到next[j+1]。这时k=next[k],直到T[j]=T[k],next[j+1]=next[k]。这样就通过递推思想求得了匹配串T的next函数。

算法流程图

 

 程序实现

/********
date:2021-10-28
author:sy
version:1.0
Description:应用KMP算法，求子串在主串中第一次出现的位置 
*********/
#include<stdio.h>
#include<string.h>
#define MAXNUM 100

typedef char DataType;
typedef struct
{
	DataType data[MAXNUM];
	int len;
}SString;

/*求模式串t的next函数值并存入数组next*/
void GetNext(DataType*t,int*next,int tlength)
{
	int i=1,j=0;
	next[1]=0;
	while(i<tlength)
	{
		if(j==0||t[i]==t[j])
		{
			++i;
			++j;
			next[i]=j;
		}
		else
		    j=next[j];
	}
}

/*利用模式串t的next函数求t在主串s中第pos个字符之后的位置*/
int IndexKmp(DataType*s,DataType*t,int pos,int tlength,int slength,int*next)
{
	int i=pos,j=1;
	while(i<=slength&&j<=tlength)
	{
		if(j==0||s[i]==t[j])  /*继续比较后继字符*/
		{
			 ++i;
			 ++j;
		 } 
		 else   
		     j=next[j];  /*模式串向后移动*/
	}
	if(j>tlength)  /*匹配成功，返回匹配起始位置*/
	    return i-tlength;
	else
	    return 0;
}
int main()
{
	int locate,tlength,slength,next[256];
	DataType s[256],t[256];
	printf("请输入第一个串（母串）：");
	slength=strlen(gets(s+1));
	printf("请输入第二个串（子串）：");
	tlength=strlen(gets(t+1));
	GetNext(t,next,tlength);
	locate=IndexKmp(s,t,0,tlength,slength,next);
	printf("匹配位置：%d\n",locate);
	return 0; 
}

运行结果