之前在本科上数据结构与算法的课上,老师说KMP算法不考,所以在学习的时候就没怎么好好看,所以在以后的应用中找子串都是蛮力查找的办法=。=!直到在POJ做题的时候才由于始终TLE才开始看KMP算法,想来也够悲催的。于是乎,自己写了个KMP的算法实现,以告诫自己不能得过且过。(若有错误之处请大家斧正)。
// KMP.cpp : 定义控制台应用程序的入口点。
//
#include "stdafx.h"
#include "string.h"
//如果next[j]=k,表示在p[0,1...j-1]子序列中,p[0,1...k-1]=p[j-k...j-1]
//即k个前缀子串与k个后缀子串相同
void getNext(char *p,int *next)
{
int j,k;
next[0]=-1;
j=0;
k=-1;
while(j<strlen(p)-1)
{
if(k==-1||p[j]==p[k]) //匹配的情况下,p[j]==p[k]
{
j++;
k++;
next[j]=k;
}
else //p[j]!=p[k]
k=next[k];
}
for(j=0;j<=strlen(p);j++)
printf("%d ",next[j]);
printf("\n");
}
//在KMP算法下,时间复杂度由之前的O(mn)将为O(m+n)
int KMPMatch(char *s,char *p)
{
int next[100];
int i,j;
i=0;
j=0;
getNext(p,next);
while(i<strlen(s))
{
if(j==-1||s[i]==p[j])
{
i++;
j++;
}
else
{
j=next[j]; //消除了指针i的回溯
}
if(j==strlen(p))
return i-strlen(p);
}
return -1;
}
int _tmain(int argc, _TCHAR* argv[])
{
int offset;
char S[]="ababcababc";
char p[]="ababa";
offset = KMPMatch(S,p);
printf("%d\n",offset);
return 0;
}