/**
 * File: edit_distance.kt
 * Created Time: 2024-01-25
 * Author: curtishd (1023632660@qq.com)
 */

package chapter_dynamic_programming

import java.util.*
import kotlin.math.min

/* 編輯距離：暴力搜尋 */
fun editDistanceDFS(
    s: String,
    t: String,
    i: Int,
    j: Int
): Int {
    // 若 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
    val insert = editDistanceDFS(s, t, i, j - 1)
    val delete = editDistanceDFS(s, t, i - 1, j)
    val replace = editDistanceDFS(s, t, i - 1, j - 1)
    // 返回最少編輯步數
    return (min(min(insert.toDouble(), delete.toDouble()), replace.toDouble()) + 1).toInt()
}

/* 編輯距離：記憶化搜尋 */
fun editDistanceDFSMem(
    s: String,
    t: String,
    mem: Array<IntArray>,
    i: Int,
    j: Int
): Int {
    // 若 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
    val insert = editDistanceDFSMem(s, t, mem, i, j - 1)
    val delete = editDistanceDFSMem(s, t, mem, i - 1, j)
    val replace = editDistanceDFSMem(s, t, mem, i - 1, j - 1)
    // 記錄並返回最少編輯步數
    mem[i][j] = (min(min(insert.toDouble(), delete.toDouble()), replace.toDouble()) + 1).toInt()
    return mem[i][j]
}

/* 編輯距離：動態規劃 */
fun editDistanceDP(s: String, t: String): Int {
    val n = s.length
    val m = t.length
    val dp = Array(n + 1) { IntArray(m + 1) }
    // 狀態轉移：首行首列
    for (i in 1..n) {
        dp[i][0] = i
    }
    for (j in 1..m) {
        dp[0][j] = j
    }
    // 狀態轉移：其餘行和列
    for (i in 1..n) {
        for (j in 1..m) {
            if (s[i - 1] == t[j - 1]) {
                // 若兩字元相等，則直接跳過此兩字元
                dp[i][j] = dp[i - 1][j - 1]
            } else {
                // 最少編輯步數 = 插入、刪除、替換這三種操作的最少編輯步數 + 1
                dp[i][j] =
                    (min(
                        min(dp[i][j - 1].toDouble(), dp[i - 1][j].toDouble()),
                        dp[i - 1][j - 1].toDouble()
                    ) + 1).toInt()
            }
        }
    }
    return dp[n][m]
}

/* 編輯距離：空間最佳化後的動態規劃 */
fun editDistanceDPComp(s: String, t: String): Int {
    val n = s.length
    val m = t.length
    val dp = IntArray(m + 1)
    // 狀態轉移：首行
    for (j in 1..m) {
        dp[j] = j
    }
    // 狀態轉移：其餘行
    for (i in 1..n) {
        // 狀態轉移：首列
        var leftup = dp[0] // 暫存 dp[i-1, j-1]
        dp[0] = i
        // 狀態轉移：其餘列
        for (j in 1..m) {
            val temp = dp[j]
            if (s[i - 1] == t[j - 1]) {
                // 若兩字元相等，則直接跳過此兩字元
                dp[j] = leftup
            } else {
                // 最少編輯步數 = 插入、刪除、替換這三種操作的最少編輯步數 + 1
                dp[j] = (min(min(dp[j - 1].toDouble(), dp[j].toDouble()), leftup.toDouble()) + 1).toInt()
            }
            leftup = temp // 更新為下一輪的 dp[i-1, j-1]
        }
    }
    return dp[m]
}

/* Driver Code */
fun main() {
    val s = "bag"
    val t = "pack"
    val n = s.length
    val m = t.length

    // 暴力搜尋
    var res = editDistanceDFS(s, t, n, m)
    println("將 $s 更改為 $t 最少需要編輯 $res 步")

    // 記憶化搜尋
    val mem = Array(n + 1) { IntArray(m + 1) }
    for (row in mem) Arrays.fill(row, -1)
    res = editDistanceDFSMem(s, t, mem, n, m)
    println("將 $s 更改為 $t 最少需要編輯 $res 步")

    // 動態規劃
    res = editDistanceDP(s, t)
    println("將 $s 更改為 $t 最少需要編輯 $res 步")

    // 空間最佳化後的動態規劃
    res = editDistanceDPComp(s, t)
    println("將 $s 更改為 $t 最少需要編輯 $res 步")
}