编程之美 3.3-- 计算字符串的相似度

许多程序会大量使用字符串。对于不同的字符串，我们希望能够有办法判断其相似程度。我们定义了一套操作方法来把两个不相同的字符串变得相同，具体的操作方法为：
1.修改一个字符（如把“a”替换为“b”）。
2.增加一个字符（如把“abdd”变为“aebdd”）。
3.删除一个字符（如把“travelling”变为“traveling”）。
比如，对于“abcdefg”和“abcdef”两个字符串来说，我们认为可以通过增加/减少一个“g“的方式来达到目的。上面的两种方案，都仅需要一次操作。把这个操作所需要的次数定义为两个字符串的距离，而相似度等于“距离+1”的倒数。也就是说，“abcdefg”和“abcdef”的距离为1，相似度 为1/2=0.5。
　　给定任意两个字符串，你是否能写出一个算法来计算它们的相似度呢？

　　原文的分析与解法　　

　　不难看出，两个字符串的距离肯定不超过它们的长度之和（我们可以通过删除操作把两个串都转化为空串）。虽然这个结论对结果没有帮助，但至少可以知道，任意两个字符串的距离都是有限的。

　　我们还是就住集中考虑如何才能把这个问题转化成规模较小的同样的子问题。如果有两个串A=xabcdae和B=xfdfa，它们的第一个字符是 相同的，只要计算A[2,...,7]=abcdae和B[2,...,5]=fdfa的距离就可以了。但是如果两个串的第一个字符不相同，那么可以进行 如下的操作（lenA和lenB分别是A串和B串的长度）。

　1.删除A串的第一个字符，然后计算A[2,...,lenA]和B[1,...,lenB]的距离。

　　2.删除B串的第一个字符，然后计算A[1,...,lenA]和B[2,...,lenB]的距离。

　　3.修改A串的第一个字符为B串的第一个字符，然后计算A[2,...,lenA]和B[2,...,lenB]的距离。

　　4.修改B串的第一个字符为A串的第一个字符，然后计算A[2,...,lenA]和B[2,...,lenB]的距离。

　　5.增加B串的第一个字符到A串的第一个字符之前，然后计算A[1,...,lenA]和B[2,...,lenB]的距离。

　　6.增加A串的第一个字符到B串的第一个字符之前，然后计算A[2,...,lenA]和B[1,...,lenB]的距离。

　　在这个题目中，我们并不在乎两个字符串变得相等之后的字符串是怎样的。所以，可以将上面的6个操作合并为：

　　1.一步操作之后，再将A[2,...,lenA]和B[1,...,lenB]变成相字符串。

　　2.一步操作之后，再将A[2,...,lenA]和B[2,...,lenB]变成相字符串。

　　3.一步操作之后，再将A[1,...,lenA]和B[2,...,lenB]变成相字符串。

这样，很快就可以完成一个递归程序。

int calStringDis(string strA, int pABegin,int pAEnd,string strB, int pBBegin,int pBEnd)
{  
    if (pABegin > pAEnd)  
    {  
        if (pBBegin > pBEnd)  
            return 0;   
        else  
            return pBEnd - pBBegin + 1;  
    }  
    if (pBBegin > pBEnd)  
    {  
        if(pABegin > pAEnd)  
            return 0;  
        else  
            return pAEnd - pABegin + 1;  
    }  
    if (strA[pABegin] == strB[pBBegin])  
    {  
        return calStringDis(strA,pABegin+1,pAEnd,strB,pBBegin+1,pBEnd);  
    }  
    else  
    {  
        int t1 = calStringDis(strA,pABegin+1,pAEnd,strB,pBBegin+2,pBEnd);  
        int t2 = calStringDis(strA,pABegin+2,pAEnd,strB,pBBegin+1,pBEnd);  
        int t3 = calStringDis(strA,pABegin+2,pAEnd,strB,pBBegin+2,pBEnd);  
  
        return minValue(t1,t2,t3)+1;  
    }  
}

上面的递归程序，有什么地方需要改进呢？问题在于：在递归的过程中，有些数据被重复计算了。


以上为编程之美的摘取。
我们知道适合采用动态规划方法的最优化问题中的两个要素：最优子结构和重叠子问题。另外，还有一种方法称为备忘录（memoization），可以充分利用重叠子问题的性质。
　　下面简述一下动态规划的基本思想。和分治法一样，动态规划是通过组合子问题的解而解决整个问题的。我们知道，分治算法是指将问题划分 成一睦独立的子问题，递归 地求解各子问题，然后合并子问题的解而得到原问题的解。与此不同，动态规划适用于子问题不是独立 的情况，也就是各子问题包含公共的子子问题。在这种情况 下，若用分治法则会做许多不必要的工作，即重复地求解公共的子子问题。动态规划算法对每个子子问题只求解一次，将其结果保存在一张表中，从而避免每次遇到 各个子问题时重新计算答案。
动态规划通常应用于最优化问题。此类问题可能有很多种可行解，每个解有一个值，而我们希望找出一个具有最优（最大或最小）值的解。称这样的解为该问题的“一个”最优解（而不是“确定的”最优解），因为可能存在多个取最优值的解。
　　动态规划算法的设计可以分为如下4个步骤：
　　1）描述最优解的结构。
　　2）递归定义最优解的值。
　　3）按自底向上的方式计算最优解的值。
　　4）由计算出的结果构造一个最优解。
　　第1~3步构成问题的动态规划解的基础。第4步在只要求计算最优解的值时可以略去。如果的确做了第4步，则有时要在第3步的计算中记录一些附加信息，使构造一个最优解变得容易。
　　该问题明显完全符合动态规划的两个要素，即最优子结构和重叠子问题特性。该问题的最优指的是两个字符串的最短距离，子问题的重叠性可以从原书中的那个递归算法中看出。
　　下面再来详细说说什么是重叠子问题。适用于动态规划求解的最优化问题必须具有的第二个要素是子问题的空间要“很小”，也就是用来解原问题的递归 算法可以反复地解同样的子问题，而不是总在产生新的子问题。典型地，不同的子问题数是输入规模的一个多项式。当一个递归算法不断地调用同一问题时，我们说 该最优问题包含重叠子问题。相反地，适合用分治法解决的问题只往往在递归的每一步都产生全新的问题。动态规划算法总是充分利用重叠子问题，即通过每个子问 题只解一次，把解保存在一个需要时就可以查看的表中，而每次查表的时间为常数。
代码如下：

int calculateStringDistance(string strA, string strB)
{
  int lenA = (int)strA.length()+1;
  int lenB = (int)strB.length()+1;
  int **c = new int*[lenA];
  for(int i = 0; i < lenA; i++)
    c[i] = new int[lenB];
  // Record the distance of all begin points of each string
  //初始化方式与背包问题有点不同
  for(int i = 0; i < lenA; i++) c[i][0] = i;
  for(int j = 0; j < lenB; j++) c[0][j] = j;
  c[0][0] = 0;
  for(int i = 1; i < lenA; i++)
  {
  for(int j = 1; j < lenB; j++)
  {
    if(strB[j-1] == strA[i-1])
      c[i][j] = c[i-1][j-1];
    else
      c[i][j] = minValue(c[i][j-1], c[i-1][j], c[i-1][j-1]) + 1;
  }
  }
  int ret =  c[lenA-1][lenB-1];
  for(int i = 0; i < lenA; i++)
    delete [] c[i];
  delete []c;

  return ret;
}

那如何理解代码和实际解决的问题之间的联系呢？
引用下网上的例子并做了相应修改：
两个字符串为：

 

红色字体标出的方法：自尾部最小的数，倒着推它产生的路径。
如（了，了）是5，它是由斜对角的（假，假）推出，（假，假）由（放，你）推出，（放，你）由（放，放）推出......
这里不理解的话可以看看http://wdhdmx.iteye.com/blog/1343856前面的“abc和abe”的分析


我们可以从前面的分析：
1.删除A串的第一个字符，然后计算A[2,...,lenA]和B[1,...,lenB]的距离。
2.删除B串的第一个字符，然后计算A[1,...,lenA]和B[2,...,lenB]的距离。
3.修改A串的第一个字符为B串的第一个字符，然后计算A[2,...,lenA]和B[2,...,lenB]的距离。
4.修改B串的第一个字符为A串的第一个字符，然后计算A[2,...,lenA]和B[2,...,lenB]的距离。
5.增加B串的第一个字符到A串的第一个字符之前，然后计算A[1,...,lenA]和B[2,...,lenB]的距离。
6.增加A串的第一个字符到B串的第一个字符之前，然后计算A[2,...,lenA]和B[1,...,lenB]的距离。
看出：
 
除了首次的次数增加，以后的每次增加，路径上如果由斜对角的格子推出的当前格子，那么其产生的动作是更改A或B；如果由上面的格子推出当前，则删除A或者增加A；如果由左边的格子推出，则删除B或增加B。
由此，备忘录的意思这样实现的。上面的例子产生的动作是：删除B（你），删除A（马上就），删除B（你）。
 
这是我的理解，请指正。


参考：http://blog.csdn.net/luxiaoxun/article/details/7493852
http://wdhdmx.iteye.com/blog/1343856
http://hi.baidu.com/zealot886/item/17acaf3463669a392f20c471