posts - 200, comments - 8, trackbacks - 0, articles - 0

程序员编程艺术-----第二十八 ~ 二十九章-----最大连续乘积子串、字符串编辑距离

Posted on 2013-05-28 22:18 鑫龙 阅读(317) 评论(0)  编辑 收藏 引用 所属分类: JULY_程序员编程艺术

               第二十八~二十九章：最大连续乘积子串、字符串编辑距离


前言

    时间转瞬即逝，一转眼，又有4个多月没来更新blog了，过去4个月都在干啥呢？对的，今2013年元旦和朋友利用业余时间一起搭了个方便朋友们找工作的编程面试算法论坛：为学论坛http://www.51weixue.com/。最近则开始负责一款在线编程挑战平台：英雄会http://hero.pongo.cn/，包括其产品运营，出题审题，写代码测试，制定比赛规则等等。

    前几天跟百度的几个朋友线下闲聊，听他们说，百度校招群内的不少朋友在找工作的时候都看过我的blog，一听当即便激起了自己重写此blog的欲望，恰巧眼下阳春三月(虽说已是3月，奇妙的是，前两天北京还下了一场大雪)，又是找工作的季节(相对于每年的9月份来说，3月则是一个小高潮)，那就从继续更新专为IT人员找工作时准备笔试面试的程序员编程艺术系列开始吧。

    再者从去年4月份上传的编程艺术前27章的PDF文档的1.3万下载量来看http://download.csdn.net/detail/v_july_v/4256339，此系列确确实实帮助了成千上万的人。Yeah，本文讲两个问题，

  • 第二十八章、最大连续乘积子串，
  • 第二十九章、字符串编辑距离，

    这两个问题皆是各大IT公司最喜欢出的笔试面试题，比如说前者是小米2013年校招笔试原题，而后者则更是反复出现，如去年9月26日百度一二面试题，10月9日腾讯面试题第1小题，10月13日百度2013校招北京站笔试题第二 大道题第3小题，及去年10月15日2013年Google校招笔试最后一道大题皆是考察的这个字符串编辑距离问题。

    OK，欢迎朋友们在本文下参与讨论，如果在线编译自己的代码（编程语言任选C/C++/Java/C#），可以上英雄会提交你的代码，有任何问题，欢迎随时不吝批评或指正，感谢。


第二十八章、最大连续乘积子串

题目描述：给一个浮点数序列，取最大乘积连续子串的值，例如 -2.5，4，0，3，0.5，8，-1，则取出的最大乘积连续子串为3，0.5，8。也就是说，上述数组中，3 0.5 8这3个数的乘积3*0.5*8=12是最大的，而且是连续的。

提醒：此最大乘积连续子串与最大乘积子序列不同，请勿混淆，前者子串要求连续，后者子序列不要求连续。也就是说：最长公共子串（Longest CommonSubstring）和最长公共子序列（LongestCommon Subsequence，LCS）的区别：

  • 子串（Substring）是串的一个连续的部分，
  • 子序列（Subsequence）则是从不改变序列的顺序，而从序列中去掉任意的元素而获得的新序列；
    更简略地说，前者（子串）的字符的位置必须连续，后者（子序列LCS）则不必。比如字符串“ acdfg ”同“ akdfc ”的最长公共子串为“ df ”，而它们的最长公共子序列LCS是“ adf ”，LCS可以使用动态规划法解决。

解答：

    解法一、穷举，所有的计算组合：
    或许，读者初看此题，自然会想到最大乘积子序列问题类似于最大子数组和问题：http://blog.csdn.net/v_JULY_v/article/details/6444021，可能立马会想到用最简单粗暴的方式：两个for循环直接轮询。

  1. double max=0;  
  2. double start=0;  
  3. double end=0;  
  4. for (int i=0;i<num;i++) {  
  5.     double x=arrs[i];  
  6.     for (int j = i+1; j < num; j++) {  
  7.         x*=arrs[j];  
  8.         if(x>max){  
  9.             max=x;  
  10.             start=arrs[i];  
  11.             end=arrs[j];  
  12.         }  
  13.     }  

     解法二、虽说类似于最大子数组和问题，但实际上具体处理起来诸多不同。为什么呢，因为乘积子序列中有正有负也还可能有0。我们可以把问题简化成这样：数组中找一个子序列，使得它的乘积最大；同时找一个子序列，使得它的乘积最小（负数的情况）。因为虽然我们只要一个最大积，但由于负数的存在，我们同时找这两个乘积做起来反而方便。也就是说，不但记录最大乘积，也要记录最小乘积。So，我们让

  • maxCurrent表示当前最大乘积的candidate，
  • minCurrent反之，表示当前最小乘积的candidate，
  • 而maxProduct则记录到目前为止所有最大乘积candidates的最大值。
（以上用candidate这个词是因为只是可能成为新一轮的最大/最小乘积）
    由于空集的乘积定义为1，在搜索数组前，maxCurrent，minCurrent，maxProduct都赋为1。
假设在任何时刻你已经有了maxCurrent和minCurrent这两个最大/最小乘积的candidates，新读入数组的元素x(i)后，新的最大乘积candidate只可能是maxCurrent或者minCurrent与x(i)的乘积中的较大者，如果x(i)<0导致maxCurrent<minCurrent，需要交换这两个candidates的值。
    当任何时候maxCurrent<1，由于1（空集）是比maxCurrent更好的candidate，所以更新maxCurrent为1，类似的可以更新minCurrent。任何时候maxCurrent如果比最好的maxProduct大，更新maxProduct。
    代码一：
  1. template <typename Comparable>    
  2. Comparable maxprod( const vector<Comparable>&v)    
  3. {    
  4.     int i;    
  5.     Comparable maxProduct = 1;    
  6.     Comparable minProduct = 1;    
  7.     Comparable maxCurrent = 1;    
  8.     Comparable minCurrent = 1;    
  9.     //Comparable t;    
  10.     
  11.     for( i=0; i< v.size() ;i++)    
  12.     {    
  13.         maxCurrent *= v[i];    
  14.         minCurrent *= v[i];    
  15.         if(maxCurrent > maxProduct)     
  16.             maxProduct = maxCurrent;    
  17.         if(minCurrent > maxProduct)    
  18.             maxProduct = minCurrent;    
  19.         if(maxCurrent < minProduct)    
  20.             minProduct = maxCurrent;    
  21.         if(minCurrent < minProduct)    
  22.             minProduct = minCurrent;    
  23.         if(minCurrent > maxCurrent)    
  24.             swap(maxCurrent,minCurrent);    
  25.         if(maxCurrent<1)    
  26.             maxCurrent = 1;    
  27.         //if(minCurrent>1)    
  28.         //    minCurrent =1;    
  29.     }    
  30.     return maxProduct;     
  31. }    

    代码二：思路，记录以第i个结尾的最大乘积M和最小乘积m，并且记录这两个区间的起点（终点都是i），不断更新，来源http://www.51weixue.com/thread-246-1-1.html：

  1. pair<int,int> maxproduct(double *f,int n) { //返回最大乘积的起点终点  
  2. int R = 0, r = 0;   //最大最小区间的 起点  
  3. pair<int,int> ret = make_pair(0, 0);   //最大 最小的区间下标  
  4. double M = f[0], m = f[0], answer = f[0];     // 最大 最小值  
  5.     for (int i = 1; i < n; ++i) {  
  6.         double t0 = f[i] * M, t1 = f[i] * m;  
  7.         if (t0 > t1) {  
  8.             M = t0;  
  9.             m = t1;  
  10.         }  
  11.         else {  
  12.             int t = R;  
  13.             R = r;  
  14.             r = t;  
  15.             M = t1;  
  16.             m = t0;  
  17.         }  
  18.         if (M < f[i]) {  
  19.             M = f[i];  
  20.             R = i;  
  21.         }  
  22.         if (m > f[i]) {  
  23.             m = f[i];  
  24.             r = i;  
  25.         }  
  26.         if (answer < M) {  
  27.             answer = M;  
  28.             ret = make_pair(R, i);  
  29.         }  
  30.            
  31.            
  32.     }  
  33.     return ret;  
  34. }  

    解法三、
    本题除了上述类似最大子数组和的解法，也可以直接用动态规划求解（其实，上述的解法一本质上也是动态规划，只是解题所表现出来的具体形式与接下来的解法二不同罢了。这个不同就在于下面的解法二会写出动态规划问题中经典常见的DP方程，而解法一是直接求解）。具体解法如下：
    假设数组为a[]，直接利用动归来求解，考虑到可能存在负数的情况，我们用Max来表示以a结尾的最大连续子串的乘积值，用Min表示以a结尾的最小的子串的乘积值，那么状态转移方程为：
       Max=max{a, Max[i-1]*a, Min[i-1]*a};
       Min=min{a, Max[i-1]*a, Min[i-1]*a};
    初始状态为Max[1]=Min[1]=a[1]。

    C/C++代码一，很简洁的一小段代码：

  1. double func(double *a,const int n)  
  2. {  
  3.     double *maxA = new double[n];  
  4.     double *minA = new double[n];  
  5.     maxA[0] = minA[0] =a[0];  
  6.     double value = maxA[0];  
  7.     for(int i = 1 ; i < n ; ++i)  
  8.     {  
  9.         maxA[i] = max(max(a[i],maxA[i-1]*a[i]),minA[i-1]*a[i]);  
  10.         minA[i] = min(min(a[i],maxA[i-1]*a[i]),minA[i-1]*a[i]);  
  11.         value=max(value,maxA[i]);  
  12.     }  
  13.     return value;  
  14. }  

    C/C++代码二：

  1. /*  
  2.  给定一个浮点数数组，有正有负数，0，正数组成,数组下标从1算起  
  3.  求最大连续子序列乘积，并输出这个序列，如果最大子序列乘积为负数，那么就输出-1  
  4.  用Max[i]表示以a[i]结尾乘积最大的连续子序列  
  5.  用Min[i]表示以a[i]结尾乘积最小的连续子序列  因为有复数，所以保存这个是必须的  
  6. */    
  7. void longest_multiple(double *a,int n){    
  8.  double *Min=new double[n+1]();    
  9.  double *Max=new double[n+1]();    
  10.  double *p=new double[n+1]();    
  11.  //初始化    
  12.  for(int i=0;i<=n;i++){    
  13.   p[i]=-1;    
  14.  }    
  15.  Min[1]=a[1];    
  16.  Max[1]=a[1];    
  17.  double max_val=Max[1];    
  18.  for(int i=2;i<=n;i++){    
  19.   Max[i]=max(Max[i-1]*a[i],Min[i-1]*a[i],a[i]);    
  20.   Min[i]=min(Max[i-1]*a[i],Min[i-1]*a[i],a[i]);    
  21.   if(max_val<Max[i])    
  22.    max_val=Max[i];    
  23.  }    
  24.  if(max_val<0)    
  25.   printf("%d",-1);    
  26.  else    
  27.   printf("%d",max_val);    
  28. //内存释放    
  29.  delete [] Max;    
  30.  delete [] Min;    
  31. }  

    C#版完整代码(代码来自参加英雄会在线编程挑战之1019、最大乘积连续子串：http://hero.pongo.cn/Question/Details?ID=19&ExamID=19的在线提交代码的用户)：

  1. //答题英雄：danielqkj  
  2. using System;  
  3. public class Test   
  4. {  
  5.     void Max(double a, double b, double c)  
  6.     {  
  7.         double d = (a>b)?a:b;  
  8.         return (d>c)?d:c;      
  9.     }  
  10.       
  11.     void Min(double a, double b, double c)  
  12.     {  
  13.         double d = (a>b)?b:a;  
  14.         return (d>c)?c:d;  
  15.     }  
  16.       
  17.       
  18.     public static void Main()  
  19.     {  
  20.         int n = Int32.parse(Console.readline());  
  21.         double[] a = new double[n];  
  22.         double maxvalue = a[0];  
  23.         double[] max = new double[n];  
  24.         double[] min = new double[n];  
  25.         double start, end;  
  26.           
  27.         String[] s = Console.readline().split(' ');  
  28.         for (int i = 0; i < n; i++)  
  29.         {  
  30.             a[i] = Double.parse(s[i])  
  31.         }  
  32.           
  33.         max[0] = a[0];  
  34.         min[0] = a[0];  
  35.         start = 0, end = 0;  
  36.           
  37.         for (int i = 1; i < n; i++)  
  38.         {  
  39.             max[i]=Max(a[i], a[i]*max[i-1], a[i]*min[i-1]);  
  40.             min[i]=Min(a[i], a[i]*max[i-1], a[i]*min[i-1]);  
  41.               
  42.             if (max[i] > maxvalue)  
  43.             {  
  44.                 maxvalue = max[i];  
  45.                 end = i;  
  46.             }  
  47.         }  
  48.           
  49.         double mmm = maxvalue;  
  50.         while ( (mmm - 0.0) > 0.00001 )  
  51.         {  
  52.             start = end;  
  53.             mmm = mmm / a[start];  
  54.         }  
  55.           
  56.         Console.Writeline(a[start] + " " + a[end] + " " + maxvalue);  
  57.           
  58.     }  
  59. }  

变种

  此外，此题还有另外的一个变种形式，即给定一个长度为N的整数数组，只允许用乘法，不能用除法，计算任意（N-1）个数的组合中乘积最大的一组，并写出算法的时间复杂度。

  我们可以把所有可能的（N-1）个数的组合找出来，分别计算它们的乘积，并比较大小。由于总共有N个（N-1）个数的组合，总的时间复杂度为O（N2），显然这不是最好的解法。
  OK，以下解答来自编程之美
    解法1

 解法2
  此外，还可以通过分析，进一步减少解答问题的计算量。假设N个整数的乘积为P，针对P的正负性进行如下分析（其中，AN-1表示N-1个数的组合，PN-1表示N-1个数的组合的乘积）。
1.P为0
          那么，数组中至少包含有一个0。假设除去一个0之外，其他N-1个数的乘积为Q，根据Q的正负性进行讨论：
Q为0
说明数组中至少有两个0，那么N-1个数的乘积只能为0，返回0；
Q为正数
返回Q，因为如果以0替换此时AN-1中的任一个数，所得到的PN-1为0，必然小于Q；
Q为负数
如果以0替换此时AN-1中的任一个数，所得到的PN-1为0，大于Q，乘积最大值为0。

     2.    P为负数

根据“负负得正”的乘法性质，自然想到从N个整数中去掉一个负数，使得PN-1为一个正数。而要使这个正数最大，这个被去掉的负数的绝对值必须是数组中最小的。我们只需要扫描一遍数组，把绝对值最小的负数给去掉就可以了。

      3.    P为正数

    类似地，如果数组中存在正数值，那么应该去掉最小的正数值，否则去掉绝对值最大的负数值。
    上面的解法采用了直接求N个整数的乘积P，进而判断P的正负性的办法，但是直接求乘积在编译环境下往往会有溢出的危险（这也就是本题要求不使用除法的潜在用意），事实上可做一个小的转变，不需要直接求乘积，而是求出数组中正数（+）、负数（-）和0的个数，从而判断P的正负性，其余部分与以上面的解法相同。

    在时间复杂度方面，由于只需要遍历数组一次，在遍历数组的同时就可得到数组中正数（+）、负数（-）和0的个数，以及数组中绝对值最小的正数和负数，时间复杂度为O（N）。


第二十九章、字符串编辑距离

题目描述：给定一个源串和目标串，能够对源串进行如下操作：
   1.在给定位置上插入一个字符
   2.替换任意字符
   3.删除任意字符
写一个程序，返回最小操作数，使得对源串进行这些操作后等于目标串，源串和目标串的长度都小于2000。

提醒：上文前言中已经说过了，此题反复出现，最近考的最多的是百度和Google的笔试面试经常考察。下图则是2013年Google的校招试题原景重现：


解答：

    解法一、    此题跟上面的最大连续乘积子串类似，常见的思路是动态规划，下面是简单的DP状态方程：

  1. //动态规划：    
  2.     
  3. //f[i,j]表示s[0...i]与t[0...j]的最小编辑距离。    
  4. f[i,j] = min { f[i-1,j]+1,  f[i,j-1]+1,  f[i-1,j-1]+(s[i]==t[j]?0:1) }    
  5.     
  6. //分别表示：添加1个，删除1个，替换1个（相同就不用替换）。   

    解法二、本解法来自为学论坛：http://www.51weixue.com/thread-482-1-1.html。

     编辑距离的定义和计算方法如下：
Given two strings A and B, edit A to B with the minimum number of edit operations:

  • a) .Replace a letter with another letter
  • b) .Insert a letter
  • c) .Delete a letter
    E.g.
A = interestingly    _i__nterestingly
B = bioinformatics   bioinformatics__
                     1011011011001111
Edit distance = 11
    Instead of minimizing the number of edge operations, we can associate a cost function to the
operations and minimize the total cost. Such cost is called edit distance. Instead of using string edit, in computational biology, people like to use string alignment.We use similarity function, instead of cost function, to evaluate the goodness of the alignment.
    E.g. of similarity function: match – 2, mismatch, insert, delete – -1.
Consider two strings ACAATCC and AGCATGC.
One of their alignment is

1.jpg

In the above alignment, space (‘_’) is introduced to both strings. There are 5 matches, 1
mismatch, 1 insert, and 1 delete.The alignment has similarity score 7.
    A_CAATCC
    AGCA_TGC
    Note that the above alignment has the maximum score.Such alignment is called optimal
alignment.String alignment problem tries to find the alignment with the maximum similarity
score!String alignment problem is also called global alignment problem.
Needleman-Wunsch algorithm
    Consider two strings S[1..n] and T[1..m].Define V(i, j) be the score of the optimal alignment
between S[1..i] and T[1..j].
Basis:
V(0, 0) = 0
V(0, j) = V(0, j-1) + d(_, T[j])：Insert j times
V(i, 0) = V(i-1, 0) + d(S, _)：Delete i times
that is:

2.jpg

Example :

3.jpg

4.jpg

下面是代码，测试数据比较少，若有问题请指正：

  1. //copyright@ peng_weida  
  2. //实现代码如下：  
  3. //头文件StrEditDistance.h  
  4. #pragma once  
  5. #include <string>  
  6. class CStrEditDistance  
  7. {  
  8. public:  
  9.     CStrEditDistance(std::string& vStrRow, std::string& vStrColumn);  
  10.     ~CStrEditDistance(void);  
  11.     int  getScore()    { return m_Score;   }  
  12.     int  getEditDis()  { return m_EditDis; }  
  13.     void setEditDis(int vDis) { m_EditDis = vDis; }  
  14.     void setScore(int vScore) { m_Score = vScore; }  
  15. private:  
  16.     void process(const std::string& vStrRow, const std::string& vStrColumn);  
  17.     int getMaxValue(int a, int b, int c)  
  18.     {  
  19.         if (a < b){ if (b < c) return c; return b; }  
  20.         else { if (b > c) return a; return a < c ? c : a; }  
  21.     }  
  22. private:  
  23.     int   m_EditDis;  
  24.     int   m_Score;  
  25. };  
  26. //源文件StrEditDistance.cpp  
  27. #include "StrEditDistance.h"  
  28. #include <iostream>  
  29. #include <iomanip>  
  30. #define MATCH        2  
  31. #define MISS_MATCH   -1  
  32. #define INSERT       -1  
  33. #define DELETE       -1  
  34. CStrEditDistance::CStrEditDistance(std::string& vStrRow, std::string& vStrColumn)  
  35. {  
  36.     process(vStrRow, vStrColumn);  
  37. }  
  38. CStrEditDistance::~CStrEditDistance(void)  
  39. {  
  40. }  
  41. //FUNCTION:  
  42. void CStrEditDistance::process(const std::string& vStrRow, const std::string& vStrColumn)  
  43. {  
  44.     int editDis = 0;     //编辑距离  
  45.     int row = vStrColumn.length();    
  46.     int column = vStrRow.length();  
  47.     const int sizeR = row + 1;  
  48.     const int sizeC = column + 1;  
  49.    
  50.     int **pScore = new int*[sizeR];  //二维指针  
  51.     for (int i = 0; i <= row; i++)  
  52.     pScore = new int[sizeC];  
  53.    
  54.     //初始化第一行和第一列  
  55.     for (int c = 0; c <= column; c++)  
  56.         pScore[0][c] = 0 - c;  
  57.     for (int r = 0; r <= row; r++)  
  58.         pScore[r][0] = 0 - r;  
  59.    
  60.     //从v(1,1)开始每列计算  
  61.     for (int c = 1; c <= column; c++)  
  62.     {  
  63.         for (int r = 1; r <= row; r++)  
  64.         {  
  65.           //计算v(i,j)  
  66.           int valueMatch;  
  67.           if (vStrColumn[r-1] == vStrRow[c-1])  
  68.               valueMatch = MATCH;  
  69.           else  
  70.               valueMatch = MISS_MATCH;    
  71.           int A = pScore[r-1][c] + INSERT;  
  72.           int B = pScore[r][c-1] + DELETE;  
  73.           int C = pScore[r-1][c-1] + valueMatch;  
  74.           pScore[r][c] = getMaxValue(A, B, C);  
  75.         }  
  76.     }  
  77.    
  78.     //计算编辑距离  
  79.     int r = row, c = column;  
  80.     while(r > 0 && c > 0)  
  81.     {  
  82.         if (pScore[r][c]+1 == pScore[r-1][c])      { editDis++; r--; continue; }  
  83.         else if (pScore[r][c]+1 == pScore[r][c-1]) { editDis++; c--; continue; }  
  84.         else if (pScore[r][c]+1 == pScore[r-1][c-1]){ editDis++; r--; c--; continue; }  
  85.         else { r--; c--; }  
  86.     }  
  87.     if (r > 0 && c == 0)  editDis += r;  
  88.     else if (c > 0 && r == 0) editDis += c;  
  89.    
  90.     std::cout << std::endl;  
  91.     //----------------DEBUG-------------------//  
  92.     //打印数据  
  93.     for (int i = 0; i <= row; i++)  
  94.     {  
  95.          for (int j = 0; j <= column; j++)  
  96.          std::cout << std::setw(2) << pScore[j] << "  ";  
  97.          std::cout << std::endl;  
  98.     }  
  99.     std::cout << std::endl;  
  100.    
  101.     //设置编辑距离和得分  
  102.     setEditDis(editDis);  
  103.     setScore(pScore[row][column]);  
  104.    
  105.     for (int i = 0; i <= row; i++)   //释放内存  
  106.     {  
  107.         delete pScore;  
  108.         pScore = NULL;  
  109.     }  
  110.     delete[] pScore;  
  111. }  

类似

    上述问题类似于编程之美上的下述一题「以下内容摘自编程之美第3.3节」：

许多程序会大量使用字符串。对于不同的字符串，我们希望能够有办法判断其相似程度。我们定义了一套操作方法来把两个不相同的字符串变得相同，具体的操作方法为：

  1. 修改一个字符（如把“a”替换为“b”）；
  2. 增加一个字符（如把“abdd ”变为“aebdd ”）；
  3. 删除一个字符（如把“travelling”变为“traveling”）。
    比如，对于“abcdefg”和“abcdef ”两个字符串来说，我们认为可以通过增加/减少一个“g”的方式来达到目的。上面的两种方案，都仅需要一次操作。把这个操作所需要的次数定义为两个字符串的距离，而相似度等于“距离+1”的倒数。也就是说，“abcdefg”和“abcdef”的距离为1，相似度为1 / 2 = 0.5。
给定任意两个字符串，你是否能写出一个算法来计算出它们的相似度呢？

 这样，很快就可以完成一个递归程序，如下所示：

  1. Int CalculateStringDistance(string strA, int pABegin, int pAEnd,    
  2.    string strB, int pBBegin, int pBEnd)     
  3. {    
  4.      if(pABegin > pAEnd)    
  5.      {    
  6.           if(pBBegin > pBEnd)    
  7.                return 0;     
  8.           else    
  9.      
  10.                return pBEnd – pBBegin + 1;    
  11.      }    
  12.     
  13.      if(pBBegin > pBEnd)    
  14.      {    
  15.           if(pABegin > pAEnd)    
  16.                return 0;    
  17.           else    
  18.                return pAEnd – pABegin + 1;    
  19.      }    
  20.     
  21.      if(strA[pABegin] == strB[pBBegin])    
  22.      {    
  23.           return CalculateStringDistance(strA, pABegin + 1, pAEnd,    
  24.             strB, pBBegin + 1, pBEnd);    
  25.      }    
  26.      else    
  27.      {    
  28.           int t1 = CalculateStringDistance(strA, pABegin, pAEnd, strB,     
  29.             pBBegin + 1, pBEnd);    
  30.           int t2 = CalculateStringDistance(strA, pABegin + 1, pAEnd,     
  31.             strB,pBBegin, pBEnd);    
  32.           int t3 = CalculateStringDistance(strA, pABegin + 1, pAEnd,    
  33.             strB,pBBegin + 1, pBEnd);    
  34.           return minValue(t1,t2,t3) + 1;    
  35.      }    
  36. }    
上面的递归程序，有什么地方需要改进呢？问题在于：在递归的过程中，有些数据被重复计算了。

　　 我们知道适合采用动态规划方法的最优化问题中的两个要素：最优子结构和重叠子问题。另外，还有一种方法称为备忘录（memoization），可以充分利用重叠子问题的性质。

　　下面简述一下动态规划的基本思想。和分治法一样，动态规划是通过组合子问题的解而解决整个问题的。我们知道，分治算法是指将问题划分 成一睦独立的子问题，递归 地求解各子问题，然后合并子问题的解而得到原问题的解。与此不同，动态规划适用于子问题不是独立 的情况，也就是各子问题包含公共的子子问题。在这种情况 下，若用分治法则会做许多不必要的工作，即重复地求解公共的子子问题。动态规划算法对每个子子问题只求解一次，将其结果保存在一张表中，从而避免每次遇到 各个子问题时重新计算答案。

动态规划通常应用于最优化问题。此类问题可能有很多种可行解，每个解有一个值，而我们希望找出一个具有最优（最大或最小）值的解。称这样的解为该问题的“一个”最优解（而不是“确定的”最优解），因为可能存在多个取最优值的解。

　　动态规划算法的设计可以分为如下4个步骤：

　　1）描述最优解的结构。

　　2）递归定义最优解的值。

　　3）按自底向上的方式计算最优解的值。

　　4）由计算出的结果构造一个最优解。

　　第1~3步构成问题的动态规划解的基础。第4步在只要求计算最优解的值时可以略去。如果的确做了第4步，则有时要在第3步的计算中记录一些附加信息，使构造一个最优解变得容易。

　　该问题明显完全符合动态规划的两个要素，即最优子结构和重叠子问题特性。该问题的最优指的是两个字符串的最短距离，子问题的重叠性可以从原书中的那个递归算法中看出。

　　下面再来详细说说什么是重叠子问题。适用于动态规划求解的最优化问题必须具有的第二个要素是子问题的空间要“很小”，也就是用来解原问题的递归 算法可以反复地解同样的子问题，而不是总在产生新的子问题。典型地，不同的子问题数是输入规模的一个多项式。当一个递归算法不断地调用同一问题时，我们说 该最优问题包含重叠子问题。相反地，适合用分治法解决的问题只往往在递归的每一步都产生全新的问题。动态规划算法总是充分利用重叠子问题，即通过每个子问 题只解一次，把解保存在一个需要时就可以查看的表中，而每次查表的时间为常数。

根据以上的分析，我写了如下的动态规划算法：

复制代码
/*DP Algorithm
　 * A loop method using dynamic programming.
　 * Calculate from bottom to top.
　 */
 int calculateStringDistance(string strA, string strB)
 {
　　　　 int lenA = (int)strA.length();
　　　　 int lenB = (int)strB.length();
　　　　 int c[lenA+1][lenB+1]; 
        // Record the distance of all begin points of each string
        //初始化方式与背包问题有点不同
　　　　 for(int i = 0; i < lenA; i++) c[i][lenB] = lenA - i;
　　　　 for(int j = 0; j < lenB; j++) c[lenA][j] = lenB - j;
　　　　 c[lenA][lenB] = 0;
　　　　 for(int i = lenA-1; i >= 0; i--)
　　　　　　　　 for(int j = lenB-1; j >= 0; j--)
　　　　　　　　 {
　　　　　　　　　　　　 if(strB[j] == strA[i])
　　　　　　　　　　　　　　　　 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;
　　　　　　　　 }
 
　　　　 return c[0][0];
 }

深入

  1. 详细读者朋友们也已经看到了，百度/Google经常喜欢出这个字符串编辑距离，实际上，关于这个“编辑距离”问题在搜索引擎中有着重要的作用，如搜索引擎关键字查询中拼写错误的提示，如下图所示，当你输入“Jult”后，因为没有这个单词“Jult”，所以搜索引擎猜测你可能是输入错误，进而会提示你是不是找“July”：
    但这个拼写错误检查的原理是什么呢？Google是基于贝叶斯统计推断的方法，相关原理详情可以看下Google的研发总监Peter Norvig写的这篇文章：http://norvig.com/spell-correct.html，以及fuanyif写的这篇：http://www.ruanyifeng.com/blog/2012/10/spelling_corrector.html。
  2. 关于什么是“编辑距离”：一个快速、高效的Levenshtein算法实现，这个是计算两个字符串的算法，Levenshtein距离又称为“编辑距离”，是指两个字符串之间，由一个转换成另一个所需的最少编辑操作次数。当然，次数越小越相似。这里有一个BT树的数据结构，挺有意思的：http://blog.notdot.net/2007/4/Damn-Cool-Algorithms-Part-1-BK-Trees；
  3. 最后，Lucene中也有这个算法的实现(我想，一般的搜索引擎一般都应该会有此项拼写错误检查功能的实现)，下面是lucene的源码(并没有太多优化，与实际工程中java注重实用性的原则并不背离)：
    1. public final class LevensteinDistance {  
    2.    
    3.     public LevensteinDistance () {  
    4.     }  
    5.       
    6. // Compute Levenshtein distance:   
    7. // see org.apache.commons.lang.StringUtils#getLevenshteinDistance(String, String)  
    8.       
    9.     public float getDistance (String target, String other) {  
    10.       char[] sa;  
    11.       int n;  
    12.       int p[];   
    13. //'previous' cost array, horizontally  
    14.       int d[];   
    15. // cost array, horizontally  
    16.       int _d[];   
    17. //placeholder to assist in swapping p and d  
    18.    
    19.         sa = target.toCharArray();  
    20.         n = sa.length;  
    21.         p = new int[n+1];   
    22.         d = new int[n+1];   
    23.          
    24.         final int m = other.length();  
    25.         if (n == 0 || m == 0) {  
    26.           if (n == m) {  
    27.             return 1;  
    28.           }  
    29.           else {  
    30.             return 0;  
    31.           }  
    32.         }   
    33.           
    34. // indexes into strings s and t  
    35.         int i;   
    36. // iterates through s  
    37.         int j;   
    38. // iterates through t  
    39.    
    40.         char t_j;   
    41. // jth character of t  
    42.    
    43.         int cost;   
    44. // cost  
    45.    
    46.         for (i = 0; i<=n; i++) {  
    47.             p[i] = i;  
    48.         }  
    49.    
    50.         for (j = 1; j<=m; j++) {  
    51.             t_j = other.charAt(j-1);  
    52.             d[0] = j;  
    53.    
    54.             for (i=1; i<=n; i++) {  
    55.                 cost = sa[i-1]==t_j ? 0 : 1;  
    56.                   
    57. // minimum of cell to the left+1, to the top+1, diagonally left and up +cost  
    58.                 d[i] = Math.min(Math.min(d[i-1]+1, p[i]+1),  p[i-1]+cost);  
    59.             }  
    60.    
    61.               
    62. // copy current distance counts to 'previous row' distance counts  
    63.             _d = p;  
    64.             p = d;  
    65.             d = _d;  
    66.         }  
    67.    
    68.           
    69. // our last action in the above loop was to switch d and p, so p now  
    70.           
    71. // actually has the most recent cost counts  
    72.         return 1.0f - ((float) p[n] / Math.max(other.length(), sa.length));  
    73.     }  
    74.    
    75. }  

扩展

    当然，面试官还可以继续问下去，如请问，如何设计一个比较两篇文章相似性的算法？这个问题的讨论可以看看这里：http://t.cn/zl82CAH。OK，字符串编辑距离这个问题实用性很强，限于篇幅，详情读者自己深入吧。



只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   博问   Chat2DB   管理