算法学习-动态规划

动态规划(Dynamic Programming,简称DP)是一种常用于优化问题的算法思想,通常用于在计算机科学、经济学和数学等领域解决最优化问题。它的核心思想是将大问题划分为小问题,通过解决小问题来解决大问题。在本文中,我们将介绍动态规划的原理、经典案例以及示例代码。

原理

动态规划的核心思想是将一个问题分解成若干子问题,逐步求解子问题的最优解,从而得到原问题的最优解。动态规划包括以下几个步骤:

  • 定义问题的状态:首先需要确定问题的状态,即描述问题的变量。状态可以是一个数值、一个数组、一个字符串或者更复杂的结构。
  • 定义状态转移方程:状态转移方程是将问题的状态转移为下一个状态的函数。通过状态转移方程,可以计算出当前状态下的最优解。
  • 确定初始状态:初始状态是问题的起点,它是无法通过状态转移方程计算出来的。因此,需要明确地给出初始状态。
  • 确定终止状态:终止状态是问题的终点,它是已知的,通常是问题的目标状态。
  • 求解最优解:通过状态转移方程和初始状态,可以递推计算出所有状态的最优解。最终,可以得到问题的最优解。

经典案例

1. Fibonacci数列

Fibonacci数列是一种经典的动态规划问题。其定义如下:
F(0) = 0
F(1) = 1
F(n) = F(n-1) + F(n-2) (n>=2)

这个问题可以使用DP来解决。首先定义一个状态数组f,其中f[i]表示第i个Fibonacci数的值。由于f[0]和f[1]已知,所以可以先将它们填充到数组中。然后,从f[2]开始,根据状态转移方程计算出每个Fibonacci数的值,最后返回f[n]即可。

C++代码示例:

int fibonacci(int n) {
    if (n <= 1) {
        return n;
    }
    int f[n+1];
    f[0] = 0;
    f[1] = 1;
    for (int i = 2; i <= n; i++) {
        f[i] = f[i-1] + f[i-2];
    }
    return f[n];
}

2.最长公共子序列

最长公共子序列(Longest Common Subsequence,简称LCS)是另一个经典的动态规划问题。给定两个字符串X和Y,求它们的最长公共子序列。其中,子序列指的是从原序列中删除一些元素后得到的序列,不一定要相邻。
例如,对于字符串"ABCDGH"和" AEDFHR",最长公共子序列为"ADH",长度为3。

我们需要创建一个二维数组dp,其中dp[i][j]表示序列1的前i个字符和序列2的前j个字符之间的LCS的长度。然后,我们可以根据下面的递归式填充dp数组:
if (sequence1[i-1] == sequence2[j-1])
dp[i][j] = dp[i-1][j-1] + 1
else
dp[i][j] = max(dp[i-1][j], dp[i][j-1])

最终结果将存储在dp[m][n]中,其中m是序列1的长度,n是序列2的长度。示例代码如下:

int lcs(string s1, string s2) {
    int m = s1.length();
    int n = s2.length();
    vector<vector<int>> dp(m+1, vector<int>(n+1, 0));
    for (int i = 1; i <= m; i++) {
        for (int j = 1; j <= n; j++) {
            if (s1[i-1] == s2[j-1])
                dp[i][j] = dp[i-1][j-1] + 1;
            else
                dp[i][j] = max(dp[i-1][j], dp[i][j-1]);
        }
    }
    return dp[m][n];
}

这段代码中,我们首先定义了两个字符串s1和s2,然后创建了一个二维向量dp。接下来,我们使用嵌套的for循环来填充dp数组,最后返回dp[m][n]作为最长公共子序列的长度。

3.背包问题

背包问题,是指有一个固定大小的背包,需要选择一些物品放入其中,每个物品有自己的重量和价值,要求选择的物品重量不超过背包容量,使得放入背包的物品总价值最大。

0/1背包问题:每个物品只能选择 0 或 1 个

定义状态 f(i,j) 表示前 i 个物品中,选择若干个放入重量不超过 j 的背包中所能获得的最大价值。

需要考虑两种情况:
不选物品 i,此时 f(i,j)=f(i-1,j)
选物品 i,此时 f(i,j)=f(i-1,j-w_i)+v_i
最终,所求解即为 f(n,W)

代码实现:

// w: 物品重量数组
// v: 物品价值数组
// c: 背包容量
int Knapsack01(vector<int>& w, vector<int>& v, int c)
{
    int n = w.size();
    // dp[i][j]表示将前i个物品放入容量为j的背包中所能获得的最大价值
    vector<vector<int>> dp(n+1, vector<int>(c+1, 0));
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= c; j++) {
            if (j < w[i-1]) {
                dp[i][j] = dp[i-1][j];
            } else {
                dp[i][j] = max(dp[i-1][j], dp[i-1][j-w[i-1]]+v[i-1]);
            }
        }
    }
    return dp[n][c];
}

int main()
{
    vector<int> w = {2, 3, 4, 5};
    vector<int> v = {3, 4, 5, 6};
    int c = 8;
    cout << "0/1 Knapsack Problem: " << Knapsack01(w, v, c) << endl;
    return 0;
}
完全背包问题:每个物品可以选择任意个

需要考虑两种情况:
不选物品 i,此时 f(i,j)=f(i-1,j)
选物品 i,此时 f(i,j)=f(i,j-w_i)+v_i
最终,所求解即为 f(n,W)

// w: 物品重量数组
// v: 物品价值数组
// c: 背包容量
int KnapsackUnbounded(vector<int>& w, vector<int>& v, int c)
{
    int n = w.size();
    // dp[i][j]表示将前i个物品放入容量为j的背包中所能获得的最大价值
    vector<vector<int>> dp(n+1, vector<int>(c+1, 0));
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= c; j++) {
            if (j < w[i-1]) {
                dp[i][j] = dp[i-1][j];
            } else {
                dp[i][j] = max(dp[i-1][j], dp[i][j-w[i-1]]+v[i-1]);
            }
        }
    }
    return dp[n][c];
}

int main()
{
    vector<int> w = {2, 3, 4, 5};
    vector<int> v = {3, 4, 5, 6};
    int c = 8;
    cout << "Unbounded Knapsack Problem: " << KnapsackUnbounded(w, v, c) ;
    return 0;
}

总结一下,动态规划是解决许多优化问题的有用工具。通过将大问题分解为子问题,我们可以使用动态规划来找到最优解。在实践中,动态规划通常用于优化问题的解决方案,以便在较短的时间内找到最优解。在使用动态规划时,我们需要确定状态转移方程,并设计适当的数据结构来存储计算过程中产生的中间结果。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 211,423评论 6 491
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,147评论 2 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 157,019评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,443评论 1 283
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,535评论 6 385
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 49,798评论 1 290
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,941评论 3 407
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,704评论 0 266
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,152评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,494评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,629评论 1 340
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,295评论 4 329
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,901评论 3 313
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,742评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,978评论 1 266
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,333评论 2 360
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,499评论 2 348

推荐阅读更多精彩内容