You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
hello-algo/zh-hant/codes/csharp/chapter_dynamic_programming/edit_distance.cs

142 lines
5.1 KiB

/**
* File: edit_distance.cs
* Created Time: 2023-07-14
* Author: hpstory (hpstory1024@163.com)
*/
namespace hello_algo.chapter_dynamic_programming;
public class edit_distance {
/* 編輯距離:暴力搜尋 */
int EditDistanceDFS(string s, string t, int i, int j) {
// 若 s 和 t 都為空,則返回 0
if (i == 0 && j == 0)
return 0;
// 若 s 為空,則返回 t 長度
if (i == 0)
return j;
// 若 t 為空,則返回 s 長度
if (j == 0)
return i;
// 若兩字元相等,則直接跳過此兩字元
if (s[i - 1] == t[j - 1])
return EditDistanceDFS(s, t, i - 1, j - 1);
// 最少編輯步數 = 插入、刪除、替換這三種操作的最少編輯步數 + 1
int insert = EditDistanceDFS(s, t, i, j - 1);
int delete = EditDistanceDFS(s, t, i - 1, j);
int replace = EditDistanceDFS(s, t, i - 1, j - 1);
// 返回最少編輯步數
return Math.Min(Math.Min(insert, delete), replace) + 1;
}
/* 編輯距離:記憶化搜尋 */
int EditDistanceDFSMem(string s, string t, int[][] mem, int i, int j) {
// 若 s 和 t 都為空,則返回 0
if (i == 0 && j == 0)
return 0;
// 若 s 為空,則返回 t 長度
if (i == 0)
return j;
// 若 t 為空,則返回 s 長度
if (j == 0)
return i;
// 若已有記錄,則直接返回之
if (mem[i][j] != -1)
return mem[i][j];
// 若兩字元相等,則直接跳過此兩字元
if (s[i - 1] == t[j - 1])
return EditDistanceDFSMem(s, t, mem, i - 1, j - 1);
// 最少編輯步數 = 插入、刪除、替換這三種操作的最少編輯步數 + 1
int insert = EditDistanceDFSMem(s, t, mem, i, j - 1);
int delete = EditDistanceDFSMem(s, t, mem, i - 1, j);
int replace = EditDistanceDFSMem(s, t, mem, i - 1, j - 1);
// 記錄並返回最少編輯步數
mem[i][j] = Math.Min(Math.Min(insert, delete), replace) + 1;
return mem[i][j];
}
/* 編輯距離:動態規劃 */
int EditDistanceDP(string s, string t) {
int n = s.Length, m = t.Length;
int[,] dp = new int[n + 1, m + 1];
// 狀態轉移:首行首列
for (int i = 1; i <= n; i++) {
dp[i, 0] = i;
}
for (int j = 1; j <= m; j++) {
dp[0, j] = j;
}
// 狀態轉移:其餘行和列
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= m; j++) {
if (s[i - 1] == t[j - 1]) {
// 若兩字元相等,則直接跳過此兩字元
dp[i, j] = dp[i - 1, j - 1];
} else {
// 最少編輯步數 = 插入、刪除、替換這三種操作的最少編輯步數 + 1
dp[i, j] = Math.Min(Math.Min(dp[i, j - 1], dp[i - 1, j]), dp[i - 1, j - 1]) + 1;
}
}
}
return dp[n, m];
}
/* 編輯距離:空間最佳化後的動態規劃 */
int EditDistanceDPComp(string s, string t) {
int n = s.Length, m = t.Length;
int[] dp = new int[m + 1];
// 狀態轉移:首行
for (int j = 1; j <= m; j++) {
dp[j] = j;
}
// 狀態轉移:其餘行
for (int i = 1; i <= n; i++) {
// 狀態轉移:首列
int leftup = dp[0]; // 暫存 dp[i-1, j-1]
dp[0] = i;
// 狀態轉移:其餘列
for (int j = 1; j <= m; j++) {
int temp = dp[j];
if (s[i - 1] == t[j - 1]) {
// 若兩字元相等,則直接跳過此兩字元
dp[j] = leftup;
} else {
// 最少編輯步數 = 插入、刪除、替換這三種操作的最少編輯步數 + 1
dp[j] = Math.Min(Math.Min(dp[j - 1], dp[j]), leftup) + 1;
}
leftup = temp; // 更新為下一輪的 dp[i-1, j-1]
}
}
return dp[m];
}
[Test]
public void Test() {
string s = "bag";
string t = "pack";
int n = s.Length, m = t.Length;
// 暴力搜尋
int res = EditDistanceDFS(s, t, n, m);
Console.WriteLine("將 " + s + " 更改為 " + t + " 最少需要編輯 " + res + " 步");
// 記憶化搜尋
int[][] mem = new int[n + 1][];
for (int i = 0; i <= n; i++) {
mem[i] = new int[m + 1];
Array.Fill(mem[i], -1);
}
res = EditDistanceDFSMem(s, t, mem, n, m);
Console.WriteLine("將 " + s + " 更改為 " + t + " 最少需要編輯 " + res + " 步");
// 動態規劃
res = EditDistanceDP(s, t);
Console.WriteLine("將 " + s + " 更改為 " + t + " 最少需要編輯 " + res + " 步");
// 空間最佳化後的動態規劃
res = EditDistanceDPComp(s, t);
Console.WriteLine("將 " + s + " 更改為 " + t + " 最少需要編輯 " + res + " 步");
}
}