c# kmp算法

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了：）