⭐算法入门⭐《动态规划 - 串匹配》困难02 —— LeetCode 72. 编辑距离

🙉饭不食，水不饮，题必须刷🙉
C语言免费动漫教程，和我一起打卡！ 🌞《光天化日学C语言》🌞
LeetCode 太难？先看简单题！ 🧡《C语言入门100例》🧡
数据结构难？不存在的！ 🌳《画解数据结构》🌳
LeetCode 太简单？算法学起来！ 🌌《夜深人静写算法》🌌

一、题目

1、题目描述

给定两个长度不超过 500 的单词 $A$ 和 $B$，请计算出将 $A$ 转换成 $B$ 所使用的最少操作数 。可以对一个单词进行如下三种操作：
  1）插入一个字符；
  2）删除一个字符；
  3）替换一个字符；

  样例输入： A = "intentiona"、B = "executiona"
  样例输出： 5。

2、基础框架

• c++ 版本给出的基础框架代码如下：

class Solution {
public:
    int minDistance(string word1, string word2) {
    }
};

• word1就是源串 A；
• word2就是目标串 B；

3、原题链接

LeetCode 72. 编辑距离

二、解题报告

1、思路分析

• 令 $f(i,j)$ 表示 源串前缀 $a_0{a}_1...a_i$ 变成 目标串前缀 $b_0{b}_1...b_j$ 的最小操作数。
• 令 $I$、$D$、$R$ 分别为插入一个字符、删除一个字符、替换一个字符的操作数，本题中值均为 1。

1）插入

• 假设 $a_0{a}_1...a_i$ 变成 $b_0{b}_1...b_{j-1}$ 的最少操作数已经求出，等于 $f(i,j-1)$，则需要在 $a[i]$ 的后面插入一个字符 $b_j$，那么最少操作数为：
• $$f(i,j-1)+I$$
• 如图所示，源字符串为 “AGTA”，目标字符串为 “GATCGT” 的情况下，将源字符串变成 "“GATCG” 的最小消耗为 $f(i,j-1)$，那么只要在源字符串最后再插入一个 ‘T’，就可以把源字符串变成目标字符串 “GATCGT”；

2）删除

• 假设 $a_0{a}_1...a_{i-1}$ 变成 $b_0{b}_1...b_j$ 的最少操作数已经求出，等于 $f(i-1,j)$，则需要把 $a_i$ 个删掉，那么产生的最少操作数为：
• $$f(i-1,j)+D$$
• 如图所示，源字符串为 “AGTA”，目标字符串为 “GATCGT” 的情况下，将 “AGT” 变成目标字符串的最小消耗为 $f(i-1,j)$，那么只要把源字符串最后一个’A’删掉，就可以把源字符串变成目标字符串；

3）替换

• 假设 $a_0{a}_1...a_{i-1}$ 变成 $b_0{b}_1...b_{j-1}$ 的最少消耗已经求出，等于 $f(i-1,j-1)$，则将 $a_i$ 替换成 $b_j$， $a_0,a_1,...,a_i$ 就可以变成 $b_0,b_1,...b_j$。替换时需要考虑 $a_i=b_j$ 和 $a_i\ne b_j$ 的情况，所以替换产生的最小操作数为：
• $$f(i-1,j-1)+\{\begin{array}{ll}0& a_i=b_j\\ R& a_i\ne b_j\end{array}$$
• 如图所示，源字符串为 “AGTA”，目标字符串为 “GATCGT” 的情况下，将 “AGT” 变成 “GATCGT” 的最小操作数为 $f(i-1,j-1)$，那么只要将 源字符串 的最后一个字符 替换为 目标字符串 的最后一个字符 ，就可以把源字符串变成目标字符串；替换时根据 源字符串 和 目标字符串 原本是否相等来决定操作数；

4）边界处理

• 边界情况主要考虑以下几种：

a. 空串变成目标串

• 即 $f(-1,j)$，总共需要插入 $j+1$ 个字符，所以 $f(-1,j)=f(-1,j-1)+I$；

b. 源字符串变成空串

• 即 $f(i,-1)$，总共需要删除 $i+1$ 个字符，所以 $f(i,-1)=f(i-1,-1)+D$；

c. 空串变成空串

• 即 $f(-1,-1)=0$

2、状态转移方程

• 将上述所有状态进行一个整合，得到状态转移方程如下：
• $$f(i,j)=\{\begin{array}{ll}0& i=-1,j=-1\\ f(i,j-1)+I& i=-1,j\ge 0\\ f(i-1,j)+D& i\ge 0,j=-1\\ {\min }_{i\ge 0,j\ge 0}\{\begin{array}{l}f(i,j-1)+I\\ f(i-1,j)+D\\ f(i-1,j-1)+R_{a_i\ne b_j}\end{array}& \end{array}$$
• 如果源串长度为 $n$，目标串长度为 $m$，则通过这个状态矩阵，最后计算得到 $f(n-1,m-1)$ 就是该题所求 "源字符串 $a_0,a_1,...,a_{n-1}$，经过 插入、删除、替换 变成目标字符串 $b_0,b_1,...b_{m-1}$" 的最少操作数了。
• 有关编辑距离的更多内容，可以参考：夜深人静写算法（二十二）- 最小编辑距离。

3、时间复杂度

• 源串的长度为 $n$，目标串的长度为 $m$。
• 状态数：$O(nm)$
• 状态转移：$O(1)$
• 时间复杂度：$O(nm)$

4、代码详解

#define I 1
#define D 1
#define R 1
const int maxn = 505;
int dp[maxn][maxn];

class Solution {
    void setdp(int r, int c, int val) {
        dp[r+1][c+1] = val;                              // (1)
    }
    int getdp(int r, int c) {
        if(r == -1 && c == -1) {
            return 0;                                    // (2)
        }
        return dp[r+1][c+1];                             // (3)
    }
public:
    int minDistance(string word1, string word2) {
        int n = word1.size();
        int m = word2.size();
        memset(dp, 0x7f7f7f7f, sizeof(dp));              // (4)
        for(int i = 0; i < word1.size(); ++i) {
            setdp(i, -1, (i+1) * D);                     // (5)
        }
        for(int i = 0; i < word2.size(); ++i) {
            setdp(-1, i, (i+1) * I);                     // (6)
        }
        for(int i = 0; i < word1.size(); ++i) {
            for(int j = 0; j < word2.size(); ++j) {
                int ICost = getdp(i, j-1) + I;
                int DCost = getdp(i-1, j) + D;
                int RCost = getdp(i-1, j-1) + ( (word1[i] == word2[j]) ? 0 : R);
                int ans = min( min(ICost, DCost), RCost );   // (7)
                setdp(i, j, ans);
            }
        }
        return getdp(word1.size()-1, word2.size()-1);
    }
};

• $(1)$ 设置状态数组的值，偏移+1用来避免数组下标越界；
• $(2)$ 空串 变成 空串 的情况；
• $(3)$ 获取状态数组的值；
• $(4)$ 初始化将状态数组中的所有状态变成一个非常大的值；
• $(5)$ 初始状态：源串 变 空串；
• $(6)$ 初始状态：空串 变 目标串；
• $(7)$ 参考上文阐述的状态转移方程；

---

三、本题小知识

  动态规划问题需要用到数组来进行数据缓存，但是有时候会发现状态存在负数，这时候就需要加上一个偏移量，来放置数组下标越界，本文中的偏移量为1。

---