KMP算法

```c
//顺序串基本运算的算法
#include <stdio.h>
#define MaxSize 100
typedef struct
{	
	char data[MaxSize];		//串中字符
	int length;				//串长
} SqString;					//声明顺序串类型
void StrAssign(SqString &s,char cstr[])	//字符串常量赋给串s
{
	int i;
	for (i=0;cstr[i]!='\0';i++)
		s.data[i]=cstr[i];
	s.length=i;
}
void DestroyStr(SqString &s)		//销毁串
{  }

void StrCopy(SqString &s,SqString t)	//串复制
{
	for (int i=0;i<t.length;i++)
		s.data[i]=t.data[i];
	s.length=t.length;
}
bool StrEqual(SqString s,SqString t)	//判串相等
{
	bool same=true;
	if (s.length!=t.length)				//长度不相等时返回0
		same=false;
	else 
		for (int i=0;i<s.length;i++)
			if (s.data[i]!=t.data[i])	//有一个对应字符不相同时返回假
			{	same=false;
				break;
			}
	return same;
}
int StrLength(SqString s)	//求串长
{
	return s.length;
}
SqString Concat(SqString s,SqString t)	//串连接
{
	SqString str;
	int i;
	str.length=s.length+t.length;
	for (i=0;i<s.length;i++)	//s.data[0..s.length-1]→str
		str.data[i]=s.data[i];
	for (i=0;i<t.length;i++)	//t.data[0..t.length-1]→str
		str.data[s.length+i]=t.data[i];
	return str;
}
SqString SubStr(SqString s,int i,int j)	//求子串
{
	SqString str;
	int k;
	str.length=0;
	if (i<=0 || i>s.length || j<0 || i+j-1>s.length)
		return str;					//参数不正确时返回空串
	for (k=i-1;k<i+j-1;k++)  		//s.data[i..i+j]→str
		str.data[k-i+1]=s.data[k];
	str.length=j;
	return str;
} 
SqString InsStr(SqString s1,int i,SqString s2)	//插入串
{
	int j;
	SqString str;
	str.length=0;
	if (i<=0 || i>s1.length+1)		//参数不正确时返回空串
		return str;
	for (j=0;j<i-1;j++)      		//s1.data[0..i-2]→str
		str.data[j]=s1.data[j];
	for (j=0;j<s2.length;j++)		//s2.data[0..s2.length-1]→str
		str.data[i+j-1]=s2.data[j];
	for (j=i-1;j<s1.length;j++)		//s1.data[i-1..s1.length-1]→str
		str.data[s2.length+j]=s1.data[j];
	str.length=s1.length+s2.length;
	return str;
}
SqString DelStr(SqString s,int i,int j)		//串删去
{
	int k;
	SqString str;
	str.length=0;
	if (i<=0 || i>s.length || i+j>s.length+1) //参数不正确时返回空串
		return str;
	for (k=0;k<i-1;k++)       		//s.data[0..i-2]→str
		str.data[k]=s.data[k];
	for (k=i+j-1;k<s.length;k++)	//s.data[i+j-1..s.length-1]→str
		str.data[k-j]=s.data[k];
	str.length=s.length-j;
	return str;
}
SqString RepStr(SqString s,int i,int j,SqString t)	//子串替换
{
	int k;
	SqString str;
	str.length=0;
	if (i<=0 || i>s.length || i+j-1>s.length) //参数不正确时返回空串
		return str;
	for (k=0;k<i-1;k++)				//s.data[0..i-2]→str
		str.data[k]=s.data[k];
	for (k=0;k<t.length;k++)   		//t.data[0..t.length-1]→str
		str.data[i+k-1]=t.data[k];
	for (k=i+j-1;k<s.length;k++)	//s.data[i+j-1..s.length-1]→str
		str.data[t.length+k-j]=s.data[k];
	str.length=s.length-j+t.length;
	return str;
}
void DispStr(SqString s)	//输出串s
{
	if (s.length>0)
	{	for (int i=0;i<s.length;i++)
			printf("%3c",s.data[i]);
		printf("\n");
	}
}

int Index(SqString s,SqString t)	//简单匹配算法
{
	int i=0,j=0;
	while (i<s.length && j<t.length) 
	{	if (s.data[i]==t.data[j])	//继续匹配下一个字符
		{	i++;					//主串和子串依次匹配下一个字符
			j++;
		}
		else						//主串、子串指针回溯重新开始下一次匹配
		{	i=i-j+1;				//主串从下一个位置开始匹配
			j=0; 					//子串从头开始匹配
		}
	}
	if (j>=t.length)
		return(i-t.length);			//返回匹配的第一个字符的下标
	else
		return(-1);					//模式匹配不成功
}
void GetNext(SqString t,int next[])	//由模式串t求出next值
{	int j,k;
	j=0;k=-1;next[0]=-1;
	while (j<t.length-1)
	{	if (k==-1 || t.data[j]==t.data[k]) 	//k为-1或比较的字符相等时
		{	j++;k++;
			next[j]=k;
		}
		else  k=next[k];
	}
}
int KMPIndex(SqString s,SqString t)	//KMP算法
{
	int next[MaxSize],i=0,j=0;
	GetNext(t,next);
	while (i<s.length && j<t.length) 
	{	if (j==-1 || s.data[i]==t.data[j]) 
		{	i++;
			j++;				//i,j各增1
		}
		else j=next[j]; 		//i不变,j后退
	}
	if (j>=t.length)
		return(i-t.length);		//返回匹配模式串的首字符下标
	else
		return(-1);				//返回不匹配标志
}
void GetNextval(SqString t,int nextval[])  //由模式串t求出nextval值
{
	int j=0,k=-1;
	nextval[0]=-1;
	while (j<t.length)
	{	if (k==-1 || t.data[j]==t.data[k])
		{	j++;k++;
			if (t.data[j]!=t.data[k])
				nextval[j]=k;
			else
				nextval[j]=nextval[k];
		}
		else
			k=nextval[k];
	}
}
int KMPIndex1(SqString s,SqString t)	//修正的KMP算法
{
	int nextval[MaxSize],i=0,j=0;
	GetNextval(t,nextval);
	while (i<s.length && j<t.length) 
	{	if (j==-1 || s.data[i]==t.data[j]) 
		{	i++;
			j++;
		}
		else
			j=nextval[j];
	}
	if (j>=t.length)
		return(i-t.length);
	else
		return(-1);
}

int main()
{
	int j;
	int next[MaxSize],nextval[MaxSize];
	SqString s,t;
	StrAssign(s,"abcabcdabcdeabcdefabcdefg");
	StrAssign(t,"abcabcdabcdeabcdefabcdef");
	printf("串s:");DispStr(s);
	printf("串t:");DispStr(t);
	printf("简单匹配算法:\n");
	printf("  t在s中的位置=%d\n",Index(s,t));
	GetNext(t,next);			//由模式串t求出next值
	GetNextval(t,nextval);		//由模式串t求出nextval值
	printf("    j   ");
	for (j=0;j<t.length;j++)
		printf("%4d",j);
	printf("\n");
	printf(" t[j]   ");
	for (j=0;j<t.length;j++)
		printf("%4c",t.data[j]);
	printf("\n");
	printf(" next   ");
	for (j=0;j<t.length;j++)
		printf("%4d",next[j]);
	printf("\n");
	printf(" nextval");
	for (j=0;j<t.length;j++)
		printf("%4d",nextval[j]);
	printf("\n");
	printf("KMP算法:\n");
	printf("  t在s中的位置=%d\n",KMPIndex(s,t));
	printf("改进的KMP算法:\n");
	printf("  t在s中的位置=%d\n",KMPIndex1(s,t));
	DestroyStr(s); DestroyStr(t);
	return 1;
}