KMP算法的实现

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

步骤：

算法实现：

实现KMP算法需要两个步骤，首先要求出next数组的值，然后利用next数组的值与主串进行比较。next函数的编写是整个算法的核心

next函数的定义如下：

 

程序中设计了两个函数：

函数GetNext()用来求next值。

求模式串t的next函数值并存放在数组next中。

函数IndexKmp()用来实现模式匹配算法。

子串中的每个字符依次和主串中的一个连续的字符串序列相等，则称为匹配成功，反之称为不匹配成功。程序中函数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函数：

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函数。

概要设计：

GetNext()函数	求next的值
IndexKmp()函数	实现模式匹配算法

运行结果：

 

 

 代码如下：
 

#include <stdio.h>
#include <string.h>
#define MAXNUM  100
typedef char DataType;
typedef struct
{
	DataType data[MAXNUM];
	int len;
}SString;
void GetNext(DataType*t,int*next,int tlength)//求模式串t的next函数值并存入数组next 
{
	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];
	}
}
int IndexKmp(DataType *s,DataType *t,int pos,int tlength,int slength,int *next)//利用模式串t的next函数求t在串s中第pos个字符后的位置 
{
	int i=pos,j=1;//i,j分别指向串s，t中要比较的字符 
	while(i<=slength&&j<=tlength)
	{
		if(j==0||s[i]==t[j])
		{
			++j;
			++i;
		}
		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;
 }