【算法】Levenshtein Distance算法 求字符串相似度
93°2009/08/15软件综合 IP:广东
Levenshtein Distance算法 也就是 “编辑距离算法”,是俄国科学家Levenshtein提出的

具体就是,一个字符串,通过加一个字符,删一个字符,替换一个字符得到另一个字串。我们把从字符串1转换成字符串2,前面3种操作用的最少次数称为12相似度

具体步骤如图所示

ld.png

XXXXXt实现方法

调用Levenshtein_Distance返回相似度百分比,浮点数

Levenshtein_Distance(ByVal str1 As String, ByVal str2 As String) As Double

    Public Function Levenshtein_Distance(ByVal str1 As String, ByVal str2 As String) As Double
        Dim d As Long = 0
        Dim result As Double
        If str1.Length > str2.Length Then
            d = str1.Length
        Else
            d = str2.Length
        End If

        result = 1 - (LD(str1, str2) / d)

        Return result

    End Function

    Private Function Minimum(ByVal a As Integer, ByVal b As Integer, ByVal c As Integer) As Integer
        Dim min As Integer

        min = a
        If b < min Then
            min = b
        End If
        If c < min Then
            min = c
        End If

        Minimum = min

    End Function

    Public Function LD(ByVal s As String, ByVal t As String) As Integer
        Dim dis(,) As Integer
        Dim m As Integer
        Dim n As Integer
        Dim i As Integer
        Dim j As Integer
        Dim s_i As String
        Dim t_j As String
        Dim cost As Integer

        n = Len(s)
        m = Len(t)
        If n = 0 Then
            LD = m
            Exit Function
        End If
        If m = 0 Then
            LD = n
            Exit Function
        End If

        ReDim dis(0 To n, 0 To m)

        For i = 0 To n
            dis(i, 0) = i
        Next i

        For j = 0 To m
            dis(0, j) = j
        Next j

        For i = 1 To n

            s_i = Mid$(s, i, 1)

            For j = 1 To m

                t_j = Mid$(t, j, 1)

                If s_i = t_j Then
                    cost = 0
                Else
                    cost = 1
                End If

                dis(i, j) = Minimum(dis(i - 1, j) + 1, dis(i, j - 1) + 1, dis(i - 1, j - 1) + cost)

            Next j

        Next i

        LD = dis(n, m)
        Erase dis

    End Function
+1000  科创币    novakon    2009/08/15 multilang
来自:计算机科学 / 软件综合
12
 
已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也
93° 作者
15年5个月前 IP:未同步
141029
C#的实现方法,同样是调用第一个函数

public double Levenshtein_Distance(string str1, string str2)
{
    long d = 0L;
    if (str1.Length > str2.Length)
    {
        d = str1.Length;
    }
    else
    {
        d = str2.Length;
    }
    return (1.0 - (((double) this.LD(str1, str2)) / ((double) d)));
}


public int LD(string s, string t)
{
    int i;
    int n = Strings.Len(s);
    int m = Strings.Len(t);
    if (n == 0)
    {
        return m;
    }
    if (m == 0)
    {
        return n;
    }
    int[,] dis = new int[n + 1, m + 1];
    int VB$t_i4$L0 = n;
    for (i = 0; i <= VB$t_i4$L0; i++)
    {
        dis[i, 0] = i;
    }
    int VB$t_i4$L1 = m;
    int j = 0;
    while (j <= VB$t_i4$L1)
    {
        dis[0, j] = j;
        j++;
    }
    int VB$t_i4$L2 = n;
    for (i = 1; i <= VB$t_i4$L2; i++)
    {
        string s_i = Strings.Mid(s, i, 1);
        int VB$t_i4$L3 = m;
        for (j = 1; j <= VB$t_i4$L3; j++)
        {
            int cost;
            string t_j = Strings.Mid(t, j, 1);
            if (s_i == t_j)
            {
                cost = 0;
            }
            else
            {
                cost = 1;
            }
            dis[i, j] = this.Minimum(dis[i - 1, j] + 1, dis[i, j - 1] + 1, dis[i - 1, j - 1] + cost);
        }
    }
    int LD = dis[n, m];
    dis = null;
    return LD;
}


private int Minimum(int a, int b, int c)
{
    int min = a;
    if (b < min)
    {
        min = b;
    }
    if (c < min)
    {
        min = c;
    }
    return min;
}

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
93°作者
15年5个月前 IP:未同步
141032
Java的实现方法,那个,返回百分比的函数没写 。。。自己根据上面写就可以了

public class Distance {

  private int Minimum (int a, int b, int c) {
  int mi;

    mi = a;
    if (b < mi) {
      mi = b;
    }
    if (c < mi) {
      mi = c;
    }
    return mi;

  }

  public int LD (String s, String t) {
  int d[][];
  int n;
  int m;
  int i;
  int j;
  char s_i;
  char t_j;
  int cost;


    n = s.length ();
    m = t.length ();
    if (n == 0) {
      return m;
    }
    if (m == 0) {
      return n;
    }
    d = new int[n+1][m+1];



    for (i = 0; i <= n; i++) {
      d[i][0] = i;
    }

    for (j = 0; j <= m; j++) {
      d[0][j] = j;
    }


    for (i = 1; i <= n; i++) {

      s_i = XXXXarAt (i - 1);


      for (j = 1; j <= m; j++) {

        t_j = XXXXarAt (j - 1);


        if (s_i == t_j) {
          cost = 0;
        }
        else {
          cost = 1;
        }


        d[i][j] = Minimum (d[i-1][j]+1, d[i][j-1]+1, d[i-1][j-1] + cost);

      }

    }

    return d[n][m];

  }

}

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
93°作者
15年5个月前 IP:未同步
141033
Delphi的实现方法,写了百分比的函数 。。。【注意form1,直接复制要改代码】


function Form1.Levenshtein_Distance(str1: string; str2: string): Double;
begin
    d := 0;
    if (str1.Length > str2.Length) then
        d := str1.Length
    else
        d := str2.Length;
    begin
        Result := (1 - ((self.LD(str1, str2) as Double) div (d as Double)));
        exit
    end
end;


function Form1.LD(s: string; t: string): Integer;
var
    i: Integer;
    cost: Integer;
begin
    n := Strings.Len(s);
    m := Strings.Len(t);
    if (n = 0) then
        begin
            Result := m;
            exit
        end;
    if (m = 0) then
        begin
            Result := n;
            exit
        end;
    dis := New(array[(n + 1), (m + 1)] of Integer);
    VB$t_i4$L0 := n;
    i := 0;
    while ((i <= VB$t_i4$L0)) do
    begin
        dis[i, 0] := i;
        inc(i)
    end;
    VB$t_i4$L1 := m;
    j := 0;
    while ((j <= VB$t_i4$L1)) do
    begin
        dis[0, j] := j;
        inc(j)
    end;
    VB$t_i4$L2 := n;
    i := 1;
    while ((i <= VB$t_i4$L2)) do
    begin
        s_i := Strings.Mid(s, i, 1);
        VB$t_i4$L3 := m;
        j := 1;
        while ((j <= VB$t_i4$L3)) do
        begin
            t_j := Strings.Mid(t, j, 1);
            if (s_i = t_j) then
                cost := 0
            else
                cost := 1;
            dis[i, j] := self.Minimum((dis[(i - 1), j] + 1), (dis[i, (j - 1)] + 1), (dis[(i - 1), (j - 1)] + cost));
            inc(j)
        end;
        inc(i)
    end;
    LD := dis[n, m];
    dis := nil;
    begin
        Result := LD;
        exit
    end
end;

function Form1.Minimum(a: Integer; b: Integer; c: Integer): Integer;
begin
    min := a;
    if (b < min) then
        min := b;
    if (c < min) then
        min := c;
    begin
        Result := min;
        exit
    end
end;
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
93°作者
15年5个月前 IP:未同步
141034
cpp的实现方法,cpp有点那啥,算了。 = =

记得#include别丢了。。。

public: Double __gc* Levenshtein_Distance(String __gc* str1, String __gc* str2)
{
    Int64 __gc* d = 0;
    if (str1->Length > str2->Length)
    {
        d = str1->Length;
    }
    else
    {
        d = str2->Length;
    }
    return (1 - (*static_cast<__box Double*>(this->LD(str1, str2)) / *static_cast<__box Double*>(d)));
}

public: Int32 __gc* LD(String __gc* s, String __gc* t)
{
    Int32 __gc* i;
    Int32 __gc* n = Strings:[s:10]en(s);
    Int32 __gc* m = Strings:[s:10]en(t);
    if (n == 0)
    {
        return m;
    }
    if (m == 0)
    {
        return n;
    }
    Int32 __gc* dis __gc [,] = __gc new Int32 __gc*[(n + 1), (m + 1)];
    Int32 __gc* R$t_i4$L0 = n;
    for (i = 0; (i <= R$t_i4$L0); i++)
    {
        dis[i, 0] = i;
    }
    Int32 __gc* R$t_i4$L1 = m;
    Int32 __gc* j = 0;
    while ((j <= R$t_i4$L1))
    {
        dis[0, j] = j;
        j++;
    }
    Int32 __gc* R$t_i4$L2 = n;
    for (i = 1; (i <= R$t_i4$L2); i++)
    {
        String __gc* s_i = Strings::Mid(s, i, 1);
        Int32 __gc* R$t_i4$L3 = m;
        for (j = 1; (j <= R$t_i4$L3); j++)
        {
            Int32 __gc* cost;
            String __gc* t_j = Strings::Mid(t, j, 1);
            if (s_i == t_j)
            {
                cost = 0;
            }
            else
            {
                cost = 1;
            }
            dis[i, j] = this->Minimum((dis[(i - 1), j] + 1), (dis[i, (j - 1)] + 1), (dis[(i - 1), (j - 1)] + cost));
        }
    }
    Int32 __gc* LD = dis[n, m];
    dis = 0;
    return LD;
}

private: Int32 __gc* Minimum(Int32 __gc* a, Int32 __gc* b, Int32 __gc* c)
{
    Int32 __gc* min = a;
    if (b < min)
    {
        min = b;
    }
    if (c < min)
    {
        min = c;
    }
    return min;
}
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
93°作者
15年5个月前 IP:未同步
141037
为了让代码更好看,很黄很暴力!

LD1.png
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
93°作者
15年5个月前 IP:未同步
141067
引用第6楼我说要有光于2009-08-15 23:44发表的  :
Delphi 的很囧...好多exit....

等了半天 。。。好不容易等到一个回帖,居然是这破帖。。 = = 不活了 555555555 pia
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
93°作者
15年5个月前 IP:未同步
143988
…………………………………………………………………………………………………………………………………………
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论

想参与大家的讨论?现在就 登录 或者 注册

所属专业
上级专业
同级专业
93°
学者 笔友
文章
651
回复
6032
学术分
30
2007/04/10注册,6年11个月前活动
暂无简介
主体类型:个人
所属领域:无
认证方式:邮箱
IP归属地:未同步
文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
支持的图片格式:jpg, jpeg, png
插入公式
评论控制
加载中...
文号:{{pid}}
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}