计算两个字符串的编辑距离的快速算法

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

/* 求两个数字的最小者 */
#define min(val1, val2) (((val1) < (val2)) ? (val1) : (val2)) 
/* 求三个数字的最小者 */
#define min3operators(val1, val2, val3) (min(val1, val2) < min(val1, val3)) ? min(val1, val2) : min(val1, val3) 

/*****************************************************************
*********************计算两个字符串的编辑距离*********************
******************************************************************
快速算法步骤：
1.设置n为字符串s的长度。("GUMBO") 
  设置m为字符串t的长度。("GAMBOL") 
  如果n等于0，返回m并退出。
  如果m等于0，返回n并退出。
  构造两个向量v0[m+1] 和v1[m+1]，串联0..m之间所有的元素。 
2.初始化 v0 to 0..m。 
3.检查 s (i from 1 to n) 中的每个字符。 
4.检查 t (j from 1 to m) 中的每个字符，j等于1时设置v1[0]=i。
5.如果 s[i] 等于 t[j]，则编辑代价为 0；
  如果 s[i] 不等于 t[j]，则编辑代价为1。 
6.设置单元v1[j]为下面的最小值之一：
  a.紧邻该单元上方+1：v1[j-1] + 1
  b.紧邻该单元左侧+1：v0[j] + 1
  c.该单元对角线上方和左侧+cost：v0[j-1] + cost 
7.在完成迭代 (3, 4, 5, 6) 之后，v1[m]便是编辑距离的值。 
*****************************************************************/
int levenshtein(const char *string1, const char *string2)
{
int n = 0;
int m = 0;
int i = 0;
int j = 0;
int cost = 0;
int *v0 = NULL;
int *v1 = NULL;


if(NULL == string1 || NULL == string2)
{
return -1;
}


n = strlen(string1);
m = strlen(string2);
if(0 == n)
{
return m;
}
if(0 == m)
{
return n;
}

v0 = (int *)malloc(sizeof(int) * (m + 1));
v1 = (int *)malloc(sizeof(int) * (m + 1));

for(i=0; i<= m; i++)
{
v0[i] = i;
}


for(i=1; i<= n; i++)
{
for(j=1; j<= m; j++)
{
if(1 == j)
{
v1[0] = i;
}


if(string1[i-1] == string2[j-1])
{
cost = 0;
}
else
{
cost = 1;
}


v1[j]=min3operators((v1[j-1]+1), (v0[j]+1), (v0[j-1])+cost);
}


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


return v1[m];
}


int main()
{
printf("levenshtein:%d\n", levenshtein("GUMBOLUabcdefgg", "GAMBOLabcedefdgh"));
printf("min3operators:%d\n", min3operators(9, 8, 7));


return 0;
}