3

假设我想为Levensthein 距离（编辑距离）实现通常的动态规划算法。很容易想出递归：

editDistance [] ys = length ys
editDistance xs [] = length xs
editDistance (x:xs) (y:ys) 
  | x == y = editDistance xs ys
  | otherwise = minimum [
      1 + editDistance xs (y:ys),
      1 + editDistance (x:xs) ys,
      1 + editDistance xs ys]

这受到指数运行时间的影响，因此有必要记住该函数。我想通过使用 Data.Memocombinators 来做到这一点，并且我尝试了几种方法。这是我目前的尝试：

import qualified Data.MemoCombinators as M

memLength str = 
  M.wrap 
    (\i -> drop i str) 
    (\xs -> length str - length xs)
    (M.arrayRange (0,length str))

elm xs ys = (M.memo2 memListx memListy editDistance) xs ys
  where
    memListx = memLength xs
    memListy = memLength ys

然而，记忆化似乎没有任何效果，尽管我希望任何记忆化都能显着改善运行时间，因为它至少是多项式的。我的实施有什么问题？如何在尽可能多地保留编辑距离的高级定义的同时获得良好的运行时间？

4

1 回答 1

3

如果您发布的代码实际上是您正在做的事情，那么您做错了！:-)。如果你要记忆一个递归函数，你需要将递归版本的调用回调到记忆的版本中。所以例如：

editDistance (x:xs) (y:ys)
  | x == y = elm xs ys
  | ...

否则，您将执行完整的递归计算并仅存储最终结果。您需要存储中间结果。

这里还有另一个问题。elm 的备忘录表不应该依赖于它的参数（理想情况下你甚至不应该提到参数，所以你不依赖编译器足够聪明来找出依赖关系）。如果备忘录表依赖于参数，那么您必须为每个不同的参数构建一个新表，并且我们需要为所有参数共享一个表。你可以尝试一些愚蠢的事情，比如记住论点的整个结构：

elm = M.memo2 (M.list M.char) (M.list M.char)

看看这是否加快了速度（结合前一个技巧）。然后，您可以继续尝试仅使用长度而不是整个列表结构来获得额外的提升。

希望有帮助。

于 2011-03-27T04:07:20.493 回答