? ? ? ?本文全章節一共一萬七千多字,詳細介紹動態規劃基礎與進階技巧,全篇以代碼為主,認真讀完理解,你對動態規劃的理解一定會有一個質的飛躍。
一、動態規劃簡介:? ? ? ??
? ? ? ?動態規劃(Dynamic Programming,簡稱DP)?是一種通過把原問題分解為相對簡單的子問題的方式求解復雜問題的方法。它的核心思想是:將復雜問題分解成子問題,保存子問題的解,避免重復計算。
動態規劃本質上是一種用空間換時間的算法思想:
- 時間優化:避免重復計算相同的子問題
- 空間代價:需要額外的存儲空間來保存子問題的解
?動態規劃 vs 其他算法思想
算法思想 | 特點 | 適用場景 | 時間復雜度 |
---|---|---|---|
暴力遞歸 | 直觀,但有重復計算 | 小規模問題 | 指數級 |
分治算法 | 分解獨立子問題 | 子問題無重疊 | 通常O(nlogn) |
貪心算法 | 局部最優選擇 | 具有貪心選擇性質 | 通常O(n)或O(nlogn) |
動態規劃 | 保存子問題解 | 有重疊子問題和最優子結構 | 通常O(n2)或O(n3) |
動態規劃的適用條件,動態規劃能夠解決的問題必須具備以下兩個特征:
1.1 最優子結構(Optimal Substructure)
問題的最優解包含其子問題的最優解。
// 示例:最短路徑問題
// 如果A到C的最短路徑是A->B->C,那么B到C的路徑也必須是B到C的最短路徑
1.2? 重疊子問題(Overlapping Subproblems)
在用遞歸算法自頂向下求解時,每次產生的子問題并不總是新問題,有些子問題會被重復計算多次。
// 斐波那契數列的遞歸實現中,f(n-1)和f(n-2)會被重復計算
int fibonacci(int n) {if (n <= 1) return n;return fibonacci(n-1) + fibonacci(n-2); // 大量重復計算
}
1.3? 動態規劃的設計步驟:
步驟1:狀態定義
確定用什么變量來表示子問題的解。
步驟2:狀態轉移方程
找出狀態之間的遞推關系。
步驟3:初始條件和邊界處理
確定初始狀態的值。
步驟4:確定計算順序
確保計算每個狀態時,所依賴的狀態都已經計算出來。
二、經典入門案例
案例1:斐波那契數列
問題描述:計算斐波那契數列的第n項,其中F(0)=0, F(1)=1, F(n)=F(n-1)+F(n-2)。
方法1:暴力遞歸(指數時間復雜度)
#include <iostream>
#include <chrono>// 暴力遞歸實現(效率極低)
long long fibonacciRecursive(int n) {if (n <= 1) return n;return fibonacciRecursive(n - 1) + fibonacciRecursive(n - 2);
}
時間復雜度:O(2^n)
空間復雜度:O(n)(遞歸棧空間)
方法2:自頂向下的記憶化搜索
#include <vector>
#include <unordered_map>class FibonacciMemo {
private:std::unordered_map<int, long long> memo;public:long long fibonacci(int n) {// 基礎情況if (n <= 1) return n;// 檢查是否已經計算過if (memo.find(n) != memo.end()) {return memo[n];}// 計算并存儲結果memo[n] = fibonacci(n - 1) + fibonacci(n - 2);return memo[n];}
};
時間復雜度:O(n)
空間復雜度:O(n)
方法3:自底向上的動態規劃
class FibonacciDP {
public:// 使用數組存儲long long fibonacciArray(int n) {if (n <= 1) return n;std::vector<long long> dp(n + 1);dp[0] = 0;dp[1] = 1;for (int i = 2; i <= n; i++) {dp[i] = dp[i - 1] + dp[i - 2];}return dp[n];}// 空間優化版本long long fibonacciOptimized(int n) {if (n <= 1) return n;long long prev2 = 0, prev1 = 1;long long current;for (int i = 2; i <= n; i++) {current = prev1 + prev2;prev2 = prev1;prev1 = current;}return current;}
};
時間復雜度:O(n)
空間復雜度:O(1)
案例2:爬樓梯問題
問題描述:假設你正在爬樓梯。需要n階你才能到達樓頂。每次你可以爬1或2個臺階。你有多少種不同的方法可以爬到樓頂呢?
分析過程:
- 狀態定義:
dp[i]
?表示爬到第i級臺階的方法數 - 狀態轉移方程:
dp[i] = dp[i-1] + dp[i-2]
- 要到達第i級臺階,可以從第i-1級臺階爬1步,或從第i-2級臺階爬2步
- 初始條件:
dp[1] = 1
,?dp[2] = 2
代碼實現:
#include <iostream>
#include <vector>class ClimbingStairs {
public:// 基礎DP版本int climbStairs(int n) {if (n <= 2) return n;std::vector<int> dp(n + 1);dp[1] = 1;dp[2] = 2;for (int i = 3; i <= n; i++) {dp[i] = dp[i - 1] + dp[i - 2];}return dp[n];}// 空間優化版本int climbStairsOptimized(int n) {if (n <= 2) return n;int prev2 = 1, prev1 = 2;int current;for (int i = 3; i <= n; i++) {current = prev1 + prev2;prev2 = prev1;prev1 = current;}return current;}// 擴展版本:每次可以爬1到k個臺階int climbStairsK(int n, int k) {std::vector<int> dp(n + 1, 0);dp[0] = 1; // 0級臺階有1種方法(不爬)for (int i = 1; i <= n; i++) {for (int j = 1; j <= k && j <= i; j++) {dp[i] += dp[i - j];}}return dp[n];}
};
案例3:背包問題系列
背包問題是動態規劃的經典應用,我們來看最重要的幾種類型。
0-1背包問題
問題描述:給定n個物品,每個物品有重量weight[i]
和價值value[i]
。背包容量為W,每個物品只能取一次,求背包能裝入物品的最大價值。
狀態定義與轉移:
- 狀態定義:
dp[i][w]
?表示前i個物品,背包容量為w時的最大價值 - 狀態轉移方程:
- 不選第i個物品:
dp[i][w] = dp[i-1][w]
- 選第i個物品:
dp[i][w] = dp[i-1][w-weight[i]] + value[i]
(前提是w >= weight[i]
) dp[i][w] = max(dp[i-1][w], dp[i-1][w-weight[i]] + value[i])
- 不選第i個物品:
代碼實現
#include <iostream>
#include <vector>
#include <algorithm>class Knapsack01 {
public:// 二維DP實現int knapsack2D(std::vector<int>& weights, std::vector<int>& values, int capacity) {int n = weights.size();// dp[i][w] 表示前i個物品,容量為w的最大價值std::vector<std::vector<int>> dp(n + 1, std::vector<int>(capacity + 1, 0));for (int i = 1; i <= n; i++) {for (int w = 0; w <= capacity; w++) {// 不選擇第i個物品dp[i][w] = dp[i - 1][w];// 如果能選擇第i個物品if (w >= weights[i - 1]) {dp[i][w] = std::max(dp[i][w], dp[i - 1][w - weights[i - 1]] + values[i - 1]);}}}return dp[n][capacity];}// 一維DP實現(空間優化)int knapsack1D(std::vector<int>& weights, std::vector<int>& values, int capacity) {std::vector<int> dp(capacity + 1, 0);for (int i = 0; i < weights.size(); i++) {// 從后往前遍歷,避免重復使用同一個物品for (int w = capacity; w >= weights[i]; w--) {dp[w] = std::max(dp[w], dp[w - weights[i]] + values[i]);}}return dp[capacity];}// 輸出具體的選擇方案std::vector<int> knapsackSolution(std::vector<int>& weights, std::vector<int>& values, int capacity) {int n = weights.size();std::vector<std::vector<int>> dp(n + 1, std::vector<int>(capacity + 1, 0));// 填充DP表for (int i = 1; i <= n; i++) {for (int w = 0; w <= capacity; w++) {dp[i][w] = dp[i - 1][w];if (w >= weights[i - 1]) {dp[i][w] = std::max(dp[i][w], dp[i - 1][w - weights[i - 1]] + values[i - 1]);}}}// 回溯找出具體選擇的物品std::vector<int> selected;int w = capacity;for (int i = n; i > 0; i--) {if (dp[i][w] != dp[i - 1][w]) {selected.push_back(i - 1); // 選擇了第i-1個物品(0-indexed)w -= weights[i - 1];}}std::reverse(selected.begin(), selected.end());return selected;}
};
完全背包問題
問題描述:與0-1背包類似,但每個物品可以選擇無限次。
class KnapsackComplete {
public:int completeKnapsack(std::vector<int>& weights, std::vector<int>& values, int capacity) {std::vector<int> dp(capacity + 1, 0);for (int i = 0; i < weights.size(); i++) {// 從前往后遍歷,允許重復使用同一個物品for (int w = weights[i]; w <= capacity; w++) {dp[w] = std::max(dp[w], dp[w - weights[i]] + values[i]);}}return dp[capacity];}
};
案例四:序列DP問題
最長遞增子序列(LIS)
問題描述:給定一個數組,找到其中最長的嚴格遞增子序列的長度。
#include <vector>
#include <algorithm>
#include <iostream>class LIS {
public:// O(n2) DP解法int lengthOfLIS_DP(std::vector<int>& nums) {if (nums.empty()) return 0;int n = nums.size();std::vector<int> dp(n, 1); // dp[i]表示以nums[i]結尾的最長遞增子序列長度int maxLength = 1;for (int i = 1; i < n; i++) {for (int j = 0; j < i; j++) {if (nums[j] < nums[i]) {dp[i] = std::max(dp[i], dp[j] + 1);}}maxLength = std::max(maxLength, dp[i]);}return maxLength;}// 獲取具體的LIS序列std::vector<int> findLIS(std::vector<int>& nums) {if (nums.empty()) return {};int n = nums.size();std::vector<int> dp(n, 1);std::vector<int> parent(n, -1); // 記錄前驅元素的索引int maxLength = 1;int maxIndex = 0;for (int i = 1; i < n; i++) {for (int j = 0; j < i; j++) {if (nums[j] < nums[i] && dp[j] + 1 > dp[i]) {dp[i] = dp[j] + 1;parent[i] = j;}}if (dp[i] > maxLength) {maxLength = dp[i];maxIndex = i;}}// 構造LIS序列std::vector<int> lis;int current = maxIndex;while (current != -1) {lis.push_back(nums[current]);current = parent[current];}std::reverse(lis.begin(), lis.end());return lis;}// O(nlogn) 二分搜索解法int lengthOfLIS_Binary(std::vector<int>& nums) {if (nums.empty()) return 0;std::vector<int> tails; // tails[i] 表示長度為i+1的遞增子序列的最小尾部元素for (int num : nums) {auto it = std::lower_bound(tails.begin(), tails.end(), num);if (it == tails.end()) {tails.push_back(num);} else {*it = num;}}return tails.size();}
};
最長公共子序列(LCS)
問題描述:給定兩個字符串,找到它們的最長公共子序列的長度。
class LCS {
public:int longestCommonSubsequence(std::string text1, std::string text2) {int m = text1.length(), n = text2.length();// dp[i][j] 表示text1[0..i-1]和text2[0..j-1]的LCS長度std::vector<std::vector<int>> dp(m + 1, std::vector<int>(n + 1, 0));for (int i = 1; i <= m; i++) {for (int j = 1; j <= n; j++) {if (text1[i - 1] == text2[j - 1]) {dp[i][j] = dp[i - 1][j - 1] + 1;} else {dp[i][j] = std::max(dp[i - 1][j], dp[i][j - 1]);}}}return dp[m][n];}// 獲取具體的LCS字符串std::string findLCS(std::string text1, std::string text2) {int m = text1.length(), n = text2.length();std::vector<std::vector<int>> dp(m + 1, std::vector<int>(n + 1, 0));// 填充DP表for (int i = 1; i <= m; i++) {for (int j = 1; j <= n; j++) {if (text1[i - 1] == text2[j - 1]) {dp[i][j] = dp[i - 1][j - 1] + 1;} else {dp[i][j] = std::max(dp[i - 1][j], dp[i][j - 1]);}}}// 回溯構造LCSstd::string lcs;int i = m, j = n;while (i > 0 && j > 0) {if (text1[i - 1] == text2[j - 1]) {lcs = text1[i - 1] + lcs;i--;j--;} else if (dp[i - 1][j] > dp[i][j - 1]) {i--;} else {j--;}}return lcs;}
};
案例五:路徑問題
最小路徑和
問題描述:給定一個包含非負整數的m×n網格,請找出一條從左上角到右下角的路徑,使得路徑上的數字總和為最小。
class PathSum {
public:int minPathSum(std::vector<std::vector<int>>& grid) {if (grid.empty() || grid[0].empty()) return 0;int m = grid.size(), n = grid[0].size();std::vector<std::vector<int>> dp(m, std::vector<int>(n));// 初始化第一個格子dp[0][0] = grid[0][0];// 初始化第一行for (int j = 1; j < n; j++) {dp[0][j] = dp[0][j - 1] + grid[0][j];}// 初始化第一列for (int i = 1; i < m; i++) {dp[i][0] = dp[i - 1][0] + grid[i][0];}// 填充其余位置for (int i = 1; i < m; i++) {for (int j = 1; j < n; j++) {dp[i][j] = std::min(dp[i - 1][j], dp[i][j - 1]) + grid[i][j];}}return dp[m - 1][n - 1];}// 空間優化版本int minPathSumOptimized(std::vector<std::vector<int>>& grid) {if (grid.empty() || grid[0].empty()) return 0;int m = grid.size(), n = grid[0].size();std::vector<int> dp(n);dp[0] = grid[0][0];for (int j = 1; j < n; j++) {dp[j] = dp[j - 1] + grid[0][j];}for (int i = 1; i < m; i++) {dp[0] += grid[i][0];for (int j = 1; j < n; j++) {dp[j] = std::min(dp[j], dp[j - 1]) + grid[i][j];}}return dp[n - 1];}// 輸出具體路徑std::vector<std::pair<int, int>> findMinPath(std::vector<std::vector<int>>& grid) {if (grid.empty() || grid[0].empty()) return {};int m = grid.size(), n = grid[0].size();std::vector<std::vector<int>> dp(m, std::vector<int>(n));// 計算DP表(同上)dp[0][0] = grid[0][0];for (int j = 1; j < n; j++) {dp[0][j] = dp[0][j - 1] + grid[0][j];}for (int i = 1; i < m; i++) {dp[i][0] = dp[i - 1][0] + grid[i][0];}for (int i = 1; i < m; i++) {for (int j = 1; j < n; j++) {dp[i][j] = std::min(dp[i - 1][j], dp[i][j - 1]) + grid[i][j];}}// 回溯找路徑std::vector<std::pair<int, int>> path;int i = m - 1, j = n - 1;while (i > 0 || j > 0) {path.push_back({i, j});if (i == 0) {j--;} else if (j == 0) {i--;} else {if (dp[i - 1][j] < dp[i][j - 1]) {i--;} else {j--;}}}path.push_back({0, 0});std::reverse(path.begin(), path.end());return path;}
};
案例六:區間DP
區間DP是在區間上進行動態規劃,通常用來解決關于區間的最優化問題。
矩陣鏈乘法
問題描述:給定n個矩陣的維度,求這些矩陣相乘時的最少標量乘法次數。
class MatrixChainMultiplication {
public:int matrixChainOrder(std::vector<int>& dims) {int n = dims.size() - 1; // n個矩陣// dp[i][j] 表示矩陣i到矩陣j相乘的最少乘法次數std::vector<std::vector<int>> dp(n, std::vector<int>(n, 0));// l是鏈的長度,從2開始for (int l = 2; l <= n; l++) {for (int i = 0; i <= n - l; i++) {int j = i + l - 1;dp[i][j] = INT_MAX;// 嘗試所有可能的分割點kfor (int k = i; k < j; k++) {int cost = dp[i][k] + dp[k + 1][j] + dims[i] * dims[k + 1] * dims[j + 1];dp[i][j] = std::min(dp[i][j], cost);}}}return dp[0][n - 1];}// 輸出最優括號化方案void printOptimalParens(std::vector<std::vector<int>>& s, int i, int j, std::string& result) {if (i == j) {result += "M" + std::to_string(i);} else {result += "(";printOptimalParens(s, i, s[i][j], result);printOptimalParens(s, s[i][j] + 1, j, result);result += ")";}}std::string getOptimalParentheses(std::vector<int>& dims) {int n = dims.size() - 1;std::vector<std::vector<int>> dp(n, std::vector<int>(n, 0));std::vector<std::vector<int>> s(n, std::vector<int>(n, 0)); // 記錄分割點for (int l = 2; l <= n; l++) {for (int i = 0; i <= n - l; i++) {int j = i + l - 1;dp[i][j] = INT_MAX;for (int k = i; k < j; k++) {int cost = dp[i][k] + dp[k + 1][j] + dims[i] * dims[k + 1] * dims[j + 1];if (cost < dp[i][j]) {dp[i][j] = cost;s[i][j] = k;}}}}std::string result;printOptimalParens(s, 0, n - 1, result);return result;}
};
?案例七:樹形DP
樹形DP是在樹結構上進行的動態規劃。
二叉樹的最大路徑和
問題描述:給定一個二叉樹,找到任意兩個節點間路徑上數值和的最大值。
#include <algorithm>
#include <climits>struct TreeNode {int val;TreeNode *left;TreeNode *right;TreeNode() : val(0), left(nullptr), right(nullptr) {}TreeNode(int x) : val(x), left(nullptr), right(nullptr) {}TreeNode(int x, TreeNode *left, TreeNode *right) : val(x), left(left), right(right) {}
};class TreeDP {
private:int maxSum = INT_MIN;// 返回以當前節點為起點的最大路徑和(只能向下)int maxPathDown(TreeNode* node) {if (!node) return 0;// 遞歸計算左右子樹的最大貢獻值// 如果子樹路徑和為負,則不選擇該路徑(用0代替)int leftMax = std::max(maxPathDown(node->left), 0);int rightMax = std::max(maxPathDown(node->right), 0);// 當前節點的最大路徑和(可能經過當前節點連接左右子樹)int currentMax = node->val + leftMax + rightMax;// 更新全局最大值maxSum = std::max(maxSum, currentMax);// 返回當前節點向上的最大貢獻值(只能選擇一邊)return node->val + std::max(leftMax, rightMax);}public:int maxPathSum(TreeNode* root) {maxSum = INT_MIN;maxPathDown(root);return maxSum;}
};
打家劫舍III(樹上版本)
問題描述:小偷在二叉樹形狀的房屋中偷盜,不能偷相鄰的房屋(父子節點不能同時偷)。
class HouseRobberTree {
public:// 返回一個pair:{不偷當前節點的最大收益, 偷當前節點的最大收益}std::pair<int, int> robHelper(TreeNode* node) {if (!node) return {0, 0};auto left = robHelper(node->left);auto right = robHelper(node->right);// 不偷當前節點:左右子樹可偷可不偷,取最大值int notRob = std::max(left.first, left.second) + std::max(right.first, right.second);// 偷當前節點:左右子樹都不能偷int rob = node->val + left.first + right.first;return {notRob, rob};}int rob(TreeNode* root) {auto result = robHelper(root);return std::max(result.first, result.second);}
};
?案例八:狀態壓縮DP
當狀態數量較少時,可以用位運算來壓縮狀態,提高效率。
旅行商問題(TSP)
問題描述:給定n個城市和它們之間的距離,找到訪問所有城市恰好一次并回到起點的最短路徑。
class TSP {
public:int tsp(std::vector<std::vector<int>>& dist) {int n = dist.size();int fullMask = (1 << n) - 1; // 所有城市都訪問過的狀態// dp[mask][i] 表示當前訪問過的城市集合為mask,當前在城市i的最小距離std::vector<std::vector<int>> dp(1 << n, std::vector<int>(n, INT_MAX));// 從城市0開始dp[1][0] = 0; // 狀態1表示只訪問了城市0for (int mask = 1; mask < (1 << n); mask++) {for (int u = 0; u < n; u++) {if (!(mask & (1 << u)) || dp[mask][u] == INT_MAX) continue;// 嘗試從城市u轉移到城市vfor (int v = 0; v < n; v++) {if (mask & (1 << v)) continue; // 城市v已經訪問過int newMask = mask | (1 << v);dp[newMask][v] = std::min(dp[newMask][v], dp[mask][u] + dist[u][v]);}}}// 找到訪問所有城市后回到起點的最小距離int result = INT_MAX;for (int i = 1; i < n; i++) {if (dp[fullMask][i] != INT_MAX) {result = std::min(result, dp[fullMask][i] + dist[i][0]);}}return result;}// 輸出具體路徑std::vector<int> getTSPPath(std::vector<std::vector<int>>& dist) {int n = dist.size();int fullMask = (1 << n) - 1;std::vector<std::vector<int>> dp(1 << n, std::vector<int>(n, INT_MAX));std::vector<std::vector<int>> parent(1 << n, std::vector<int>(n, -1));dp[1][0] = 0;for (int mask = 1; mask < (1 << n); mask++) {for (int u = 0; u < n; u++) {if (!(mask & (1 << u)) || dp[mask][u] == INT_MAX) continue;for (int v = 0; v < n; v++) {if (mask & (1 << v)) continue;int newMask = mask | (1 << v);int newDist = dp[mask][u] + dist[u][v];if (newDist < dp[newMask][v]) {dp[newMask][v] = newDist;parent[newMask][v] = u;}}}}// 找到最優的結束城市int minDist = INT_MAX;int lastCity = -1;for (int i = 1; i < n; i++) {if (dp[fullMask][i] + dist[i][0] < minDist) {minDist = dp[fullMask][i] + dist[i][0];lastCity = i;}}// 重構路徑std::vector<int> path;int currentMask = fullMask;int currentCity = lastCity;while (currentCity != -1) {path.push_back(currentCity);int prevCity = parent[currentMask][currentCity];currentMask ^= (1 << currentCity);currentCity = prevCity;}std::reverse(path.begin(), path.end());path.push_back(0); // 回到起點return path;}
};
?三、DP優化技巧
1. 滾動數組優化空間
當DP狀態只依賴于前一個或前幾個狀態時,可以使用滾動數組節省空間。
// 示例:背包問題的空間優化
class SpaceOptimization {
public:// 原始二維DPint knapsack2D(std::vector<int>& weights, std::vector<int>& values, int capacity) {int n = weights.size();std::vector<std::vector<int>> dp(n + 1, std::vector<int>(capacity + 1, 0));for (int i = 1; i <= n; i++) {for (int w = 0; w <= capacity; w++) {dp[i][w] = dp[i-1][w];if (w >= weights[i-1]) {dp[i][w] = std::max(dp[i][w], dp[i-1][w-weights[i-1]] + values[i-1]);}}}return dp[n][capacity];}// 空間優化后的一維DPint knapsack1D(std::vector<int>& weights, std::vector<int>& values, int capacity) {std::vector<int> dp(capacity + 1, 0);for (int i = 0; i < weights.size(); i++) {for (int w = capacity; w >= weights[i]; w--) {dp[w] = std::max(dp[w], dp[w-weights[i]] + values[i]);}}return dp[capacity];}
};
2. 單調隊列優化
在某些DP轉移中,可以使用單調隊列來優化時間復雜度。
#include <deque>class MonotonicQueueOptimization {
public:// 滑動窗口最大值的DP應用std::vector<int> maxSlidingWindow(std::vector<int>& nums, int k) {std::deque<int> dq; // 存儲數組下標std::vector<int> result;for (int i = 0; i < nums.size(); i++) {// 移除超出窗口范圍的元素while (!dq.empty() && dq.front() <= i - k) {dq.pop_front();}// 維護單調遞減隊列while (!dq.empty() && nums[dq.back()] <= nums[i]) {dq.pop_back();}dq.push_back(i);// 當窗口大小達到k時,記錄最大值if (i >= k - 1) {result.push_back(nums[dq.front()]);}}return result;}
};
3. 矩陣快速冪優化
對于線性遞推關系,可以使用矩陣快速冪將時間復雜度從O(n)降到O(logn)。
class MatrixFastPower {
private:using Matrix = std::vector<std::vector<long long>>;Matrix multiply(const Matrix& A, const Matrix& B) {int n = A.size();Matrix C(n, std::vector<long long>(n, 0));for (int i = 0; i < n; i++) {for (int j = 0; j < n; j++) {for (int k = 0; k < n; k++) {C[i][j] += A[i][k] * B[k][j];}}}return C;}Matrix matrixPower(Matrix base, long long exp) {int n = base.size();Matrix result(n, std::vector<long long>(n, 0));// 初始化為單位矩陣for (int i = 0; i < n; i++) {result[i][i] = 1;}while (exp > 0) {if (exp & 1) {result = multiply(result, base);}base = multiply(base, base);exp >>= 1;}return result;}public:// 使用矩陣快速冪計算斐波那契數列long long fibonacciMatrix(long long n) {if (n <= 1) return n;// 轉移矩陣: [[1,1],[1,0]]Matrix base = {{1, 1}, {1, 0}};Matrix result = matrixPower(base, n - 1);// F(n) = result[0][0] * F(1) + result[0][1] * F(0)// = result[0][0] * 1 + result[0][1] * 0// = result[0][0]return result[0][0];}
};
四、綜合練習題
練習1:最大子數組和(Kadane算法)
class MaxSubarray {
public:int maxSubArray(std::vector<int>& nums) {int maxSum = nums[0];int currentSum = nums[0];for (int i = 1; i < nums.size(); i++) {// 要么繼續之前的子數組,要么重新開始currentSum = std::max(nums[i], currentSum + nums[i]);maxSum = std::max(maxSum, currentSum);}return maxSum;}// 返回具體的最大子數組std::vector<int> findMaxSubArray(std::vector<int>& nums) {int maxSum = nums[0];int currentSum = nums[0];int start = 0, end = 0, tempStart = 0;for (int i = 1; i < nums.size(); i++) {if (currentSum < 0) {currentSum = nums[i];tempStart = i;} else {currentSum += nums[i];}if (currentSum > maxSum) {maxSum = currentSum;start = tempStart;end = i;}}return std::vector<int>(nums.begin() + start, nums.begin() + end + 1);}
};
練習2:編輯距離
class EditDistance {
public:int minDistance(std::string word1, std::string word2) {int m = word1.length(), n = word2.length();// dp[i][j] 表示word1[0..i-1]轉換為word2[0..j-1]的最小操作數std::vector<std::vector<int>> dp(m + 1, std::vector<int>(n + 1));// 初始化邊界條件for (int i = 0; i <= m; i++) dp[i][0] = i;for (int j = 0; j <= n; j++) dp[0][j] = j;for (int i = 1; i <= m; i++) {for (int j = 1; j <= n; j++) {if (word1[i-1] == word2[j-1]) {dp[i][j] = dp[i-1][j-1]; // 不需要操作} else {dp[i][j] = 1 + std::min({dp[i-1][j], // 刪除dp[i][j-1], // 插入dp[i-1][j-1] // 替換});}}}return dp[m][n];}// 輸出具體的編輯操作序列std::vector<std::string> getEditOperations(std::string word1, std::string word2) {int m = word1.length(), n = word2.length();std::vector<std::vector<int>> dp(m + 1, std::vector<int>(n + 1));// 填充DP表(同上)for (int i = 0; i <= m; i++) dp[i][0] = i;for (int j = 0; j <= n; j++) dp[0][j] = j;for (int i = 1; i <= m; i++) {for (int j = 1; j <= n; j++) {if (word1[i-1] == word2[j-1]) {dp[i][j] = dp[i-1][j-1];} else {dp[i][j] = 1 + std::min({dp[i-1][j], dp[i][j-1], dp[i-1][j-1]});}}}// 回溯操作序列std::vector<std::string> operations;int i = m, j = n;while (i > 0 || j > 0) {if (i > 0 && j > 0 && word1[i-1] == word2[j-1]) {i--; j--;} else if (i > 0 && j > 0 && dp[i][j] == dp[i-1][j-1] + 1) {operations.push_back("Replace " + std::string(1, word1[i-1]) + " with " + std::string(1, word2[j-1]));i--; j--;} else if (i > 0 && dp[i][j] == dp[i-1][j] + 1) {operations.push_back("Delete " + std::string(1, word1[i-1]));i--;} else {operations.push_back("Insert " + std::string(1, word2[j-1]));j--;}}std::reverse(operations.begin(), operations.end());return operations;}
};
?總結與進階建議
DP解題模板:
// 通用DP解題模板
class DPTemplate {
public:int solveProblem(/* 輸入參數 */) {// 第1步:確定狀態定義// dp[i][j][...] 的含義是什么?// 第2步:找出狀態轉移方程// dp[i][j] = f(dp[...], dp[...], ...)// 第3步:確定初始條件和邊界// dp的初始值是什么?// 第4步:確定計算順序// 保證計算dp[i][j]時,所依賴的狀態都已計算// 第5步:優化空間復雜度(如果可能)// 使用滾動數組等技巧return 0; // 返回最終答案}
};
DP問題分類總結
類型 | 特征 | 常見問題 | 時間復雜度 |
---|---|---|---|
線性DP | 狀態轉移呈線性關系 | 爬樓梯、最大子序列和 | O(n) |
區間DP | 在區間上進行DP | 矩陣鏈乘法、回文子串 | O(n3) |
背包DP | 選擇物品的最優化問題 | 0-1背包、完全背包 | O(nW) |
樹形DP | 在樹結構上的DP | 樹的直徑、打家劫舍III | O(n) |
狀態壓縮DP | 用位運算壓縮狀態 | TSP、集合覆蓋 | O(n22?) |
數位DP | 按數位進行的DP | 統計特定數字個數 | O(logn) |
進階學習路徑
- 掌握基礎問題:從斐波那契、爬樓梯開始
- 熟練背包問題:0-1背包、完全背包、多重背包
- 掌握序列DP:LIS、LCS、編輯距離
- 學習區間DP:矩陣鏈乘法、石子合并
- 深入樹形DP:樹的直徑、重心分解
- 挑戰狀態壓縮:TSP、狀態空間搜索
- 學習優化技巧:單調隊列、斜率優化、矩陣快速冪
調試技巧
// DP調試輔助函數
class DPDebugger {
public:template<typename T>void printDP(const std::vector<std::vector<T>>& dp, const std::string& name) {std::cout << "=== " << name << " ===" << std::endl;for (int i = 0; i < dp.size(); i++) {for (int j = 0; j < dp[i].size(); j++) {std::cout << std::setw(6) << dp[i][j] << " ";}std::cout << std::endl;}std::cout << std::endl;}template<typename T>void print1D(const std::vector<T>& arr, const std::string& name) {std::cout << name << ": ";for (const auto& x : arr) {std::cout << x << " ";}std::cout << std::endl;}
};
🏆 最后的話
? ? ? ?動態規劃是算法設計中的一顆明珠,它教會我們如何化繁為簡,化重復為存儲。掌握DP不僅能幫你解決編程競賽中的難題,更重要的是能培養你分解復雜問題的思維能力。
記住DP的精髓:
- 最優子結構:大問題的最優解包含小問題的最優解
- 重疊子問題:避免重復計算,用空間換時間
- 狀態轉移:找出問題內在的遞推關系
多練習,多思考,相信你一定能掌握這個強大的算法思想!
練習建議:每天至少做1-2道不同類型的DP題目,并嘗試優化空間復雜度。熟能生巧,DP的精髓就在于大量的練習和思考。