using System;
namespace kmp
{
/// <summary>
/// Summary description for Class1.
/// </summary>
class Class1
{
/// <summary>
/// The main entry point for the application.
/// </summary>
static int[] next=new int[20];
static int count=0;
static int kmp(char[] a, char[] b, int pos)
{
int i=pos, j=0;
int lena=a.Length,lenb=b.Length;
count=0;
while((i<=lena-1)&&(j<=lenb-1))
{
if(a[i]==b[j])
{++i;++j;}
else
{i=i-j+1; j=0;}
count++;
}
if(j>lenb-1)
return i-lenb;
else
return 0;
}
[STAThread]
static void Main(string[] args)
{
//
// TODO: Add code to start application here
//
int pos=0,reval;
Console.WriteLine("input string1");
string s1=Console.ReadLine();
Console.WriteLine("input string2");
string s2=Console.ReadLine();
char[] p1=s1.ToCharArray();
char[] p2=s2.ToCharArray();
Console.WriteLine("lena={0},lenb={1}",p1.Length,p2.Length);
reval=kmp(p1,p2,pos);
Console.WriteLine("====not kmp ====");
Console.WriteLine("value={0},count={1}",reval,count);
get_next(p2,next);
Console.WriteLine(next);
reval=index_kmp(p1,p2,pos);
Console.WriteLine("====this is kmp ====");
Console.WriteLine("value={0},count={1}",reval,count);
get_goodnext(p2,next);
reval=index_kmp(p1,p2,pos);
Console.WriteLine("====this is good kmp ====");
Console.WriteLine("value={0},count={1}",reval,count);
Console.ReadLine();
}
static int index_kmp(char[] a,char[] b, int pos)
{
int i=pos, j=0;
int lena=a.Length,lenb=b.Length;
count=0;
while((i<=lena-1)&&(j<=lenb-1))
{
if(j==-1||a[i]==b[j])
{++i;++j;}
else
{j=next[j];}
count++;
}
if(j>lenb-1)
return i-lenb;
else
return 0;
}
static void get_next(char[] a,int[] next)
{
int i=0,j=-1;
next[0]=-1;
while(i<a.Length-1)
{
if(j==-1||a[i]==a[j])
{++i;++j;next[i]=j;}
else
{j=next[j];}
}
}
static void get_goodnext(char[] a,int[] next)
{
int i=0,j=-1;
next[0]=-1;
while(i<a.Length-1)
{
if(j==-1||a[i]==a[j])
{
++i;++j;
if(a[i]!=a[j])
{
next[i]=j;
}
else
{
next[i]=next[j];
}
}
else
{j=next[j];}
}
}
}
}
前几天怎么调试都不通过,头都大了。回家休息了一下,又仔细看了看书上对next[j]函数的定义,再结合程序中数组的特性,修改了一下next函数的定义
next[j]=-1 当j=0时
next[j]=max{k/1<k<j 且 'p1....pk-1'='pj-k+1...pj-1'}当此集合不空时
next[j]=0 其他情况
和书上的区别就在于起始下标的不同和函数定义的不同。
调试的时候还有一个地方花了很多时间就是求next函数值时的循环条件while(i<a.Length-1)是没有‘=’号的,这个和主算法中的条件不同,感觉是next是递推而来的,由n-1就可以推出n了:)