1. 程式人生 > >動態規劃---最短編輯距離

動態規劃---最短編輯距離

描述:
設A和B是2個字串。要用最少的字元操作將字串A轉換為字串B。這裡所說的字元操作包括:
(1)刪除一個字元;
(2)插入一個字元;
(3)將一個字元改為另一個字元。
將字串A變換為字串B所用的最少字元運算元稱為字串A到B的編輯距離,記為d(A,B)。試設計一個有效演算法,對任給的2個字串A和B,計算出它們的編輯距離d(A,B)。

要求:
輸入:第1行是字串A,第2行是字串B。
輸出:字串A和B的編輯距離d(A,B)

思路:
開一個二維陣列d[i][j]來記錄a0-ai與b0-bj之間的編輯距離,要遞推時,需要考慮對其中一個字串的刪除操作、插入操作和替換操作分別花費的開銷,從中找出一個最小的開銷即為所求
具體演算法:
首先給定第一行和第一列,然後,每個值d[i,j]這樣計算:d[i][j] = min(d[i-1][j]+1, d[i][j-1]+1, d[i-1][j-1]+(s1[i]==s2[j]?0:1));
 最後一行,最後一列的那個值就是最小編輯距離

二、演算法分析

dp[i][j]表示a的前i個和b的前j個相同後的最短距離

dp[i][j]來自於三種狀態{

    1、刪除,dp[i-1][j]+1;

    2、插入,dp[i][j-1]+1;

    3、替換,if(a[i]==b[j]) dp[i][j]=dp[i-1][j-1],else dp[i][j]=dp[i-1][j-1]+1;

三、程式碼:

  1. #include<stdio.h>
  2. #include<string.h>
  3. #include<algorithm>
  4. usingnamespace std;  
  5. #define N 10
  6. int dp[N][N];//dp[i][j]表示a的前i個和b的前j個相同後的最短距離
  7. int main()  
  8. {  
  9.     char a[N];  
  10.     char b[N];  
  11.     freopen("dInput.txt","r",stdin);  
  12.     freopen("doutput.txt","w",stdout);  
  13.     scanf("%s%s",a,b);  
  14.     int la=strlen(a);  
  15.     int lb=strlen(b);  
  16.     for(int i=1;i<=la;i++)  
  17.     dp[i][0]=i;  
  18.     for
    (int j=1;j<=lb;j++)  
  19.     dp[0][j]=j;  
  20.     for(int i=1;i<=la;i++)  
  21.     {  
  22.         for(int j=1;j<=lb;j++)  
  23.         {  
  24.             int insert=dp[i][j-1]+1;  
  25.             int deletes=dp[i-1][j]+1;  
  26.             int replace;  
  27.             if(a[i-1]==b[j-1])  
  28.             replace=dp[i-1][j-1];  
  29.             else
  30.             replace=dp[i-1][j-1]+1;  
  31.             dp[i][j]=min(min(insert,deletes),replace);  
  32.         }  
  33.     }  
  34.     printf("%d\n",dp[la][lb]);  
  35.     return 0;  
  36. }