动态规划. 换钱的最少货币数和最多方法数

通过对换钱类题目的学习,我们将了解到

  • 暴力递归及优化方法
  • 记忆搜索(优化一)
  • 动态规划的基本实现方法(优化二)
  • 动态规划的空间优化(优化三)

1. 换钱的最少货币数,货币可重复使用

给定数组arr,arr中所有的值都为正数且不重复。每个值代表一种面值的货币,每种货币都可以使用任意张,再给定一个整数aim代表要找的钱数,求组成aim的最少货币数。

例如
arr=[5, 2, 3], aim=20, 返回4
arr=[5, 2, 3], aim=0, 返回0
arr=[3, 5], aim=2, 返回-1

动态规划方法的关键是构造动态规划表dp。
动态规划的实质就是在计算过程中,根据动态规划表计算下一个目标值,并更新动态规划表
在本题中,构造这样二维动态规划表,dp[i][j],其中,i j的含义如下:

  • i: 代表可以使用的货币种类为 arr[0..i]
  • j: 代表需要兑换的面值数,其取值范围为[0..aim],因此实现时,二维数组的列数应该为aim + 1

构造过程:

  1. 构造边界值。边界值是更新规划表的起始值,也是很容易犯错的地方,需要谨慎设置起始边界值。
  • dp[0..N-1][0]:规划表的第一列值,表示当需要兑换0元时,需要的货币数,显然货币数为0,直接设置为0.
  • dp[0][0..aim]:规划表的第一行值,表示只使用arr[0]一种货币兑换[0..aim]时,需要的货币数,因此,只要tmp_aim可以被arr[0]整除,返回整除后的数,即为对应的值。在下面的实现中,我们使用了更为通用的方法来得到规划表的第一行的值,可以根据下面的算法,理解一下实现中的计算方法。
  1. 更新动态规划表
  • 更新方向:逐行更新,每行从左到右更新。
  • 更新值:对于dp[i][j],更新的依据有两个值,一个是“不使用当前种类的货币时,组成总数的j的最小方法,即dp[i-1][j]”,另外一个是“使用并且仅使用一张当前种类的货币时,组成总数的j的最小方法,即 dp[i][j-arr[i]] + 1”,取这两个值中的最小值即为dp[i][j]的值。

关于依据中的第二个值,可以通过公式推导得到。可以查看文末的参考资料详细了解。

实现

class Solution
{
public:
    int ExchangeMoney(std::vector<uint32_t>& arr, uint32_t aim);
};

int Solution::ExchangeMoney(std::vector<uint32_t>& arr, uint32_t aim)
{
    // fix: precheck
    if (arr.size() == 0 || aim == 0)
    {
        return 0;
    }
    int arrSize = arr.size();
    int dp[arrSize][aim + 1];
    for (uint32_t lineIndex = 0; lineIndex < arrSize; ++lineIndex)
    {
        dp[lineIndex][0] = 0;
    }
    for (uint32_t rowIndex = 1; rowIndex < aim + 1; ++rowIndex)
    {
        // 动态初始化边界值的方法
        if (int(rowIndex - arr[0]) >= 0 && dp[0][rowIndex - arr[0]] != UINT32_MAX)
        {
            dp[0][rowIndex] = dp[0][rowIndex - arr[0]] + 1;
        }
        else
        {
            dp[0][rowIndex] = UINT32_MAX;
        }
    }
    for (uint32_t lineIndex = 1; lineIndex < arrSize; ++lineIndex)
    {
        for (uint32_t rowIndex = 1; rowIndex < aim + 1; ++rowIndex)
        {
            int subCurArrValue = rowIndex - arr[lineIndex];
            uint32_t upValue = dp[lineIndex - 1][rowIndex];
            uint32_t leftValue = UINT32_MAX;
            if (subCurArrValue >= 0 && dp[lineIndex][subCurArrValue] != UINT32_MAX)
            {
                leftValue = dp[lineIndex][subCurArrValue] + 1;
            }
            dp[lineIndex][rowIndex] = leftValue < upValue ? leftValue : upValue;
        }
    }
    return dp[arrSize - 1][aim] == UINT32_MAX ? -1 : dp[arrSize - 1][aim];
}

2. 动态规划的空间优化

上面的实现,需要的空间复杂度是O(N * aim),下面的方法,可以将空间复杂度优化到O(aim)。
优化方法:

  • 生成一个一维动态规划数组dp[aim + 1]
  • 构造初始值:参考上面的方法,初始值表示只使用arr[0]一种货币时,兑换[0..aim]的最少货币数。
  • 更新动态规划表:更新方向,从左到右更新。对于dp[j],更新的依据有两个,一个是 dp[j](old),换算成二维数组,即为不使用arr[j]货币时的最少货币数,即为dp[line-1][j]。另外一个是dp[j-arr[j]] + 1,即为“使用并且仅使用一张当前种类的货币时的最少货币数”,换算成二维数组,即为dp[line][j-arr[j]] + 1.

实现

int Solution::ExchangeMoney(std::vector<uint32_t>& arr, uint32_t aim)
{
    if (arr.size() == 0 || aim == 0)
    {
        return 0;
    }
    int arrSize = arr.size();
    int dp[aim + 1]; // updated by line
    dp[0] = 0; // first element is always 0
    for (uint32_t rowIndex = 1; rowIndex < aim + 1; ++rowIndex)
    {
        if (int(rowIndex - arr[0]) >= 0 && dp[rowIndex - arr[0]] != UINT32_MAX)
        {
            dp[rowIndex] = dp[rowIndex - arr[0]] + 1;
        }
        else
        {
            dp[rowIndex] = UINT32_MAX;
        }
    }
    for (uint32_t lineIndex = 1; lineIndex < arrSize; ++lineIndex)
    {
        for (uint32_t rowIndex = 1; rowIndex < aim + 1; ++rowIndex)
        {
            int subCurArrValue = rowIndex - arr[lineIndex];
            uint32_t upValue = dp[rowIndex];
            uint32_t leftValue = UINT32_MAX;
            if (subCurArrValue >= 0 && dp[subCurArrValue] != UINT32_MAX)
            {
                leftValue = dp[subCurArrValue] + 1;
            }
            dp[rowIndex] = leftValue < upValue ? leftValue : upValue;
        }
    }
    return dp[aim] == UINT32_MAX ? -1 : dp[aim];
}

从上面的实现,我们可以看到,我们将二维动态规划表压缩成一维动态规划表,依据是:

  • 对于表中的每个元素,我们只会在更新下一个值时,使用一次,之后便不再使用。
  • 更新下一个值,依赖“向左跳一次”和“向上跳一位”,而这两个数可以在一维动态规划表中保存。如果依赖多个“向上跳”的值,则无法使用一维动态规划表实现空间优化。

3. 还钱的最少货币数,货币不可重复使用

给定数组arr,arr中所有的值都为正数。每个值仅代表一张钱的面值,再给定一个数aim代表要找的钱数,求组成aim的最少货币数。

例如
arr = [5, 2, 3], aim = 20,无法组成,返回-1
arr = [5, 2, 5, 3], aim = 10,返回2
arr = [5, 2 ,5 ,3], aim = 15,返回4
arr = [5, 2, 5, 3], aim = 0,返回0

构造二维动态规划表dp[i][j],i表示可以使用货币arr[0..i],j表示组成的货币总数。
构造过程

  • 构造边界值
  • dp[0..N-1][0]:规划表的第一行值,表示当前需要兑换0元时需要的货币数,显然货币数为0,直接设置为0
  • dp[0][0..aim]:规划表的第一列,表示仅使用arr[0]一种货币兑换[0..aim]时,是否可以兑换,如果arr[0]就等于j,则dp[0][j] = 1,否则等于UINT32_MAX,表示无法兑换。
  • 更新动态规划表
  • 更新方向:逐行更新,每行从左到右更新。
  • 更新值:对于dp[i][j],更新的依据有两个值,一个是“不适用当前面值的货币时,组成总数j的最少方法,即dp[i-1][j]“,另一个是”使用当前面值的货币时,组成总数j的最少方法,即dp[i-1][j-arr[j]] + 1“,取这两个值中的最小值即为dp[i][j]的值。
    这里同样可以用一维动态规划表对空间优化到O(aim)。
    实现的思路基本和上面两例一样。

4 换钱的方法数,可重复使用

给定数组arr,arr中所有的值都为正数,且不重复。每个值代表一种面值的货币,每种面值的货币可以使用任意张,再给定一个整数aim代表要找的钱数,求换钱有多少种方法。

例如
arr = [5, 10, 25, 1], aim = 0,返回1
arr = [5, 10, 25, 1], aim = 15,返回6
arr = [3, 5],aim = 2,返回0

暴力递归法

递归终止条件:aim == 0
循环条件:使用当前面值货币0..k张,其中 k * arr[curIndex] <= aim
递归:对每个面值的货币,返回组成余下aim的方法数

int Solution::ExchangeMoney(std::vector<uint32_t>& arr, uint32_t aim)
{
    if (aim == 0)
    {
        return 1;
    }
    if (arr.size() == 0)
    {
        return 0;
    }
    return process(arr, 0, aim);
}

int Solution::process(std::vector<uint32_t>& arr, uint32_t index, uint32_t aim)
{
    int ret = 0;
    if (index == arr.size())
    {
        ret = 0 == aim ? 1 : 0;
    }
    else
    {
        for (uint32_t k = 0; k * arr[index] <= aim; ++k)
        {
            ret += process(arr, index + 1, aim - k * arr[index]);
        }
    }
    return ret;
}

记忆搜索法

上面的暴力递归法,存在着大量的重复计算。
优化方法是将中间结果保存下来,在下一次计算的时候,查看表中是否已经有结果。

int Solution::ExchangeMoney(std::vector<uint32_t>& arr, uint32_t aim)
{
    if (aim == 0)
    {
        return 1;
    }
    if (arr.size() == 0)
    {
        return 0;
    }
    std::map<std::pair<uint32_t,uint32_t>, uint32_t> resultMap;
    return process(arr, 0, aim, resultMap);
}
int Solution::process(std::vector<uint32_t>& arr, uint32_t index, uint32_t aim,
                std::map<std::pair<uint32_t, uint32_t>, uint32_t>& resultMap)
{
    int ret = 0;
    if (index == arr.size())
    {
        ret = 0 == aim ? 1 : 0;
    }
    else
    {
        for (uint32_t k = 0; k * arr[index] <= aim; ++k)
        {
            uint32_t nextIndex = index + 1;
            uint32_t nextAim = aim - k * arr[index];
            std::map<std::pair<uint32_t, uint32_t>, uint32_t>::iterator it = resultMap.find(std::make_pair(nextIndex, nextAim));
            if (it != resultMap.end())
            {
                ret += it->second;
            }
            else
            {
                ret += process(arr, nextIndex, nextAim, resultMap);
            }
        }
    }
    return ret;
}

动态规划

构造二维动态规划表dp[i][j]

  • 构造边界值
  • dp[0][0..j],动态规划表第一行,表示只使用arr[0]一个种类的货币时,可以兑换货币总数j的方法数,无法兑换时,则设置为0。
  • dp[0..N-1][0],动态规划表第一列,表示兑换货币总数0的方法数,设置为1.
  • 更新动态规划表
  • dp[i][j]的取值依据有两个,第一个是,"不使用当前面值的货币时,组成货币总数的方法数,dp[i-1][j]",第二个是,”分别使用1..k张当前面值的货币时,组成货币总数的方法数之和,由公式可以推导出,该值为dp[i][j-arr[i]]“,dp[i][j] = dp[i-1][j] + dp[i][j - arr[i]]。

代码实现请参考上面的几例,稍加需改即可得出。

参考资料

《程序员代码面试指南:IT名企算法与数据结构题目与最优解》

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

推荐阅读更多精彩内容

  • 动态规划(Dynamic Programming) 本文包括: 动态规划定义 状态转移方程 动态规划算法步骤 最长...
    廖少少阅读 3,256评论 0 18
  • 回溯算法 回溯法:也称为试探法,它并不考虑问题规模的大小,而是从问题的最明显的最小规模开始逐步求解出可能的答案,并...
    fredal阅读 13,623评论 0 89
  • 背景 一年多以前我在知乎上答了有关LeetCode的问题, 分享了一些自己做题目的经验。 张土汪:刷leetcod...
    土汪阅读 12,723评论 0 33
  • (欢迎转载,但请注明出处并附带链接)算法好久没复习了,今天看见一妹子在办公室刷Leetcode,顿时我也来了兴趣,...
    Nick_Zuo阅读 655评论 0 3
  • 分治方法 将问题划分成互不相交的子问题 递归地求解子问题 将子问题的解组合起来 动态规划(两个要素:最优子结构、子...
    superlj666阅读 492评论 0 0