程序员的自我修养
Home » Java语言, 算法拾遗 » Levenshtein距离

Levenshtein距离

0条评论2,232次浏览

Levenshtein距离,又称为编辑距离,或者成为字符串相似度,意思是一个字符串转换为另一个字符串所需要的最小编辑次数。编辑包含:插入、删除和替换三种操作。例如字符串:yurnom和yunum的编辑距离为2,因为需要插入一个r和将u替换为o两次操作。

Levenshtein算法由俄罗斯科学家Vladimir Levenshtein在1965年提出,该算法现在广泛用于DNA匹配、字符串纠错、实体共指等领域。

算法简介

将2个字符串转换为一个矩阵,为了简单起见,采用字符串heqi和hoi。转换后进行初始化,如下所示:

h e q i
0 1 2 3 4
h 1 A
o 2
i 3

其中数字部分为初始化的值,现计算绿色格子A处的值,其算法为:

  1. 若A所对应的两个字符相同,则左上角蓝色格子的值+0,否则+1,结果记为R1
  2. 选取上方和左方紫色格子中较小的值,并+1,结果记为R2
  3. 选取R1和R2中的较小值为绿色格子A处的值

经过计算可以得出A处的值为0,当A处计算完毕后,其右方和下方的值可以接着进行计算了,如此整个矩阵的值将可全部计算出来。最后矩阵右下方红色格子的值即为两个字符串的编辑距离。

至于为何左方和上方的值加1、左上角的值根据不同的情况加0或1的原因(也就是这么算的原理)比较简单,在此不再详述。

Java实现

运行结果:

不足之处与改造方法

原算法采用矩阵来记录中间的结果,矩阵的大小为str1.length()*str2.length()。当我们计算两个超长字符串的相似度时(例如计算两篇文章的相似度),所需要的空间则会非常的大,例如:str1、str2长度为1W,矩阵大小则为100M,若每个int为64字节,则该矩阵占用内存为8*100M=800MB。

然而实际计算过程中,该矩阵记录的绝大部分值都为中间值,这些中间值在完成后续的计算后就毫无用处了。通过观察Levenshtein算法的具体过程,可以发现,实际上只需要2“行”或者2“列”值来记录中间结果即可满足计算的要求。

算法过程

为了简单起见,现展示计算字符串heqi、hoi的过程:

  1. 初始化rowB,用rowB来计算rowA
    h e q i
    rowB 0 1 2 3 4
    rowA h 1 0 1 2 3
    o 2
    i 3
  2. rowA此时为历史数据,用于计算rowB
    h e q i
    0 1 2 3 4
    rowA h 1 0 1 2 3
    rowB o 2 1 1 2 3
    i 3
  3. 再次交换rowA和rowB的角色,利用rowB来计算rowA
    h e q i
    0 1 2 3 4
    h 1 0 1 2 3
    rowB o 2 1 1 2 3
    rowA i 3 2 2 2 2

代码实现:

最后

可以看到改进后的方法占用空间在数量级上下降了,时间复杂度保持不变,依然是O(mn),m和n为两个字符串的长度。对于此方法,可以进一步改进:将较长的字符串“竖着”放,较短的字符串横着放(参考上述算法过程中的三个表格理解),这样在时间复杂度不变的情况下,所占用的空间最小化了,只不过一般情况下该方法意义不大,只有当两个字符串长度相差悬殊的情况下有意义,例如一个字符串100W长度,一个字符串1W长度。

(转载本站文章请注明作者和出处 程序员的自我修养 – SelfUp.cn ,请勿用于任何商业用途)
标签:,
发表评论


profile
  • 文章总数:79篇
  • 评论总数:280条
  • 分类总数:31个
  • 标签总数:44个
  • 运行时间:1130天

大家好,欢迎来到selfup.cn。

这不是一个只谈技术的博客,这里记录我成长的点点滴滴,coding、riding and everthing!

最新评论
  • 张瑞昌: 有很多,比较常见的是Jacob迭代法,一次迭代O(n^3), 迭代次数不清楚。 如果是手动算的话按照定义求就可以了
  • Anonymous: :mrgreen:
  • lc277: 你好 我想问下一般删除节点要多久,要删除的datanode大概用了 1t,解除授权已经30多小时还没完成,请问是出现什么问题了吗 麻烦告诉下谢谢 qq1844554123
  • Anonymous: 你好 我想问下一般删除节点要多久,要删除的datanode大概用了 1t,解除授权已经30多小时还没完成,请问是出现什么问题了吗
  • Anonymous: :smile: :grin: :eek:
  • 李雪璇: 想要完整代码,可以帮忙发给我吗
  • Anonymous: 请问一下,那个 user的推荐结果楼主查看了么? 为什么输入数据 最高是五分,输出结果都是7分8分啥的?怎么设置输出的分数的最 大值?
  • Anonymous: 那个 user的推荐结果楼主查看了么? 为什么输入数据 最高是五分,输出结果都是7分8分啥的?
  • Anonymous: stopGracefullyOnShutdown在yarn- client模式下我测试的无效,你的呢
  • Anonymous: 另外,import的lib包能否发个列表.
  • Anonymous: 部分程序已经无法使用, 另外 你的import代码 也应该放上来哈
  • wzq: 赞
  • 增达网: 受教了!呵呵!
  • Anonymous: :!: :smile: :oops: :grin: :eek: :shock:
  • 27: :razz: dsa会报错,rsa不会
  • Anonymous: 看错了 忽略我
  • Anonymous: UserSideCF这个类在哪里
  • 晴子: 在主节点初始化CM5数据库的时候报错误:Verifying that we can write to /opt/cm-5.9.0/etc/cloudera-scm -server log4j:ERROR Could not...
  • zhangnew: 就4题 :?:
  • linxh: “ 但要是遇到预先并不知道数组的长度而又需要获取正确的(或者称之 为原始的)split长度时,该如何处理呢。。? ” 印象中可以split函数参数传-1?