Knuth-Morris-Pratt在纯C中实现

时间:2012-06-29 07:17:43

标签: c algorithm implementation knuth-morris-pratt

我有下一个KMP实施:

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

int kmp(char substr[], char str[])
{
   int i, j, N, M;

   N = strlen(str);
   M = strlen(substr);

   int *d = (int*)malloc(M * sizeof(int));
   d[0] = 0;

   for(i = 0, j = 0; i < M; i++)
   {
      while(j > 0 && substr[j] != substr[i])
      {
         j = d[j - 1];
      }

      if(substr[j] == substr[i])
      {
         j++;
         d[i] = j;
      }
   }

   for(i = 0, j = 0; i < N; i++)
   {
      while(j > 0 && substr[j] != str[i])
      {
         j = d[j - 1];
      }

      if(substr[j] == str[i])
      {
         j++;
      }

      if(j == M)
      {
         free(d);
         return i - j + 1;
      }
   }

   free(d);

   return -1;
}

int main(void)
{
   char substr[] = "World",
      str[] = "Hello World!";

   int pos = kmp(substr, str);

   printf("position starts at: %i\r\n", pos);

   return 0;
}

您可以在此处进行测试:http://liveworkspace.org/code/d2e7b3be72083c72ed768720f4716f80

它适用于小字符串,我用大循环测试它,这样一切都很好。

但是,如果我更改我正在搜索的子字符串以及完整的字符串:

char substr[] = "%end%",
str[] = "<h1>The result is: <%lua% oleg = { x = 0xa }
         table.insert(oleg, y) oleg.y = 5 print(oleg.y) %end%></h1>";

只有在第一次尝试后,此实施才会失败......

拜托,您能帮我修复KMP的实现,使算法能够在字符串中使用这些数据......

1 个答案:

答案 0 :(得分:2)

在某个地方偏离了你的来源,来源有

while(j>0 && p[j]!=p[i]) j = d[j-1];
    if(p[j]==p[i])
        j++;
        d[i]=j;

,而你有

while(j > 0 && substr[j] != substr[i])
{
    j = d[j - 1];
}
if(substr[j] == substr[i])
{
    j++;
    d[i] = j;
}

被来源的缩进所欺骗。在源代码中,if()分支周围没有大括号,因此只有增量j++;if控制; d[i] = j;是无条件的。

然后,源有一个错误,可能是由于索引的异常使用。设置数组的正确方法是

int *d = (int*)malloc(M * sizeof(int));
d[0] = 0;

for(i = 1, j = 0; i < M; i++)
{
    while(j > 0 && substr[j-1] != substr[i-1])
    {
        j = d[j - 1];
    }

    if(substr[j] == substr[i])
        j++;
    d[i] = j;
}

但令人困惑,因为此处的设置使用索引i-1j-1以及ij来确定d[i]。实现它的通常方式是不同的;它的实施方式in C#。既然这是你在大多数资源中找到的形式,那么说服自己的正确性要容易得多。