最小编辑距离 -- 解析及python实现

2401_84617246

于 2024-06-25 02:11:31 发布

阅读量413

点赞数 3

文章标签： python 开发语言

本文链接：https://blog.csdn.net/2401_84617246/article/details/139942990

版权

A[4]!=B[2]，s[4][2]=min(s[4][1],s[3][1],s[3][2])+1=3

A[4]!=B[3]，s[4][3]=min(s[4][2],s[3][2],s[3][3])+1=3

A[4]!=B[4]，s[4][4]=min(s[4][3],s[3][3],s[3][4])+1=2

A\|B		a	b	c	e
	0	1	2	3	4
a	1	0	1	2	3
c	2	1	1	1	2
e	3	2	2	2	1
f	4	3	3	3	2

所以A—>B编辑距离为2次，操作为：acef在字符ac之间插入字符b，删除字符f

三、实验代码

1. 方法一：递归


A = input("输入字符串1：")

B = input("输入字符串2：")

def recursive_edit_distance(str_a, str_b):

    if len(str_a) == 0:

        return len(str_b)

    elif len(str_b) == 0:

        return len(str_a)

    elif str_a[len(str_a)-1] == str_b[len(str_b)-1]:

        return recursive_edit_distance(str_a[0:-1], str_b[0:-1])

    else:

        return min([

            recursive_edit_distance(str_a[:-1], str_b),

            recursive_edit_distance(str_a, str_b[:-1]),

            recursive_edit_distance(str_a[:-1], str_b[:-1])

        ]) + 1

print(recursive_edit_distance(A, B))

2. 方法二：动态规划


import nltk

A = input("输入字符串1：")

B = input("输入字符串2：")

def minDistance(w1, w2):

   m, n = len(w1),len(w2)

   if(m == 0):

      return m

   if(n == 0):

      return n

   step = [[0]*(n+1)for _ in range(m + 1)]

   for i in range(1, m+1):step[i][0]=i

   for j in range(1, n+1):step[0][j]=j

   for i in range(1, m+1):

      for j in range(1, n+1):

         if w1[i-1] == w2[j-1] :

            diff=0

         else:diff=1

         step[i][j] = min(step[i-1][j-1],min(step[i-1][j],step[i][j-1]))+diff

   return step[m][n]

print(minDistance(A,B))

打印出完整变换：

A = input("输入字符串1：")

B = input("输入字符串2：")

def edit_distance_Omn(str_a, str_b):

    if str_a == str_b:

        return 0

    if len(str_a) == 0:

        return len(str_b)

    if len(str_b) == 0:

        return len(str_a)

    dp = [[0 for _ in range(len(str_a) + 1)] for _ in range(len(str_b) + 1)]

    for i in range(len(str_b) + 1):

        dp[i][0] = i

    for j in range(len(str_a) + 1):

        dp[0][j] = j

    for i in range(1, len(str_b) + 1):

        for j in range(1, len(str_a) + 1):

            dp[i][j] = dp[i-1][j-1]

            if str_a[j-1] != str_b[i-1]:

                dp[i][j] = min([dp[i-1][j-1], dp[i-1][j], dp[i][j-1]]) + 1

  #打印完整路径矩阵（这一步非必要）

    '''for i in range(len(str_b) + 1):

        for j in range(len(str_a) + 1):

            print(dp[i][j])

        # print()'''

  # 准备倒着查询编辑路径，从右下角开始

    i , j = len(str_b), len(str_a)

    op_list = []  # 记录编辑操作

    while i > 0 and j > 0:

        if dp[i][j] == dp[i-1][j-1]:

            op_list.append("keep [ {} ]".format(str_b[i-1]))

            i, j = i-1, j-1

            continue

        if dp[i][j] == dp[i-1][j]  + 1:

            op_list.append("remove [ {} ]".format(str_b[i-1]))

            i, j = i-1, j

            continue

        if dp[i][j] == dp[i-1][j-1] + 1:

            op_list.append("change [ {} ] to [ {} ]".format(str_b[i-1], str_a[j-1]))

            i, j = i-1, j-1

            continue

        if dp[i][j] == dp[i][j-1] + 1:

            op_list.append("insert [ {} ]".format(str_a[j-1]))

            i, j = i, j-1

    for i in range(len(op_list)):

        print(op_list[len(op_list)-i-1])

    return dp[len(str_b)][len(str_a)]

print(edit_distance_Omn(A, B))

四、结果