python数据结构与算法 26 动态规划

动态规划

很多程序的目的在数值优化。象计算两地之间的最短距离,一个点集的最佳拟合曲线,符合某些条件的最小集合等。这些问题有很多的求解方法,本节的目标就是向你展示不同的问题解决策略。动态规划就是优化问题的解决方法之一。

找零问题是数值优化的经典问题之一,假设你是一个自动售货机厂商的程序员,公司决定,为每次找交易的找零,计算需要的最小的硬币数量,从而实现流程化。如果一个顾客投入了1美元,买了37美分的东西,需要最小的硬币是多少?答案是6个:225美分(译者:不知道作者为什么不使用50美分的硬币),110美分,31美分。答案怎么来的?我们从最大的硬币开始,尽可能多地使用大额硬币,直到余额不够这个面值,然后再尽可能多地使用第二面额的硬币,以此类推。这种方法叫做“贪婪方法”,因为我们总是对正在计算的问题,作出当前看来最好的结果。

贪婪方法在美国硬币问题上回答正确,但如果你的公司要把自动售货机出口到Lower Elbonia(译注:Elbonia 是系列漫画 Dilbert 中一个杜撰的原东欧共产主义国家。Lower Elbonia 指这国家的南部。),那里除了平时所用的151025美分硬币,还有一种21美分的。这时贪婪算法就失灵了,因为它仍然会算出6个硬币,但正确答案是321美分的硬币。

我们在研究找出一种确保正确答案的方法。既然本章是在谈递归,你也许想到用递归解决。先说递归的基点开始。

如果我们要找的零钱,正好是某个硬币的面值,那就简单了,一个硬币。

如果找零和手上的硬币都对不上号,那么答案是以下四种情况的最小值:

从要找的零钱中拿掉一个1分币,计算余额最少需要多少硬币,然后加上这个硬币

从要找的零钱中拿掉一个5分币,计算余额最少需要多少硬币,然后加上这个硬币

从要找的零钱中拿掉一个10币,计算余额最少需要多少硬币,然后加上这个硬币

从要找的零钱中拿掉一个25分币,计算余额最少需要多少硬币,然后加上这个硬币。

也就是按以下公式计算:

python数据结构与算法 26 动态规划

以上算法的实现如下。第三行是查检基点,就是说看一下找钱是否就是某个硬币的面值,如果不是,就把找钱的值减去每种硬币的面值,然后依次递归计算。第6行加了一个过滤,就是只比较面币比找钱小的硬币,递归过程也一直在降低找钱的数值。第7行进行递归调用,注意这一行我们对硬币数量有一个加1的的过程,因为后面的递归是减少一个硬币的递归,趋向基点的过程。

Listing 7

def recMC(coinValueList,change):

   minCoins = change

   if change incoinValueList:

     return 1

   else:

      for i in [c for c incoinValueList if c <= change]:

         numCoins = 1 +recMC(coinValueList,change-i)

         if numCoins <minCoins:

            minCoins = numCoins

   return minCoins

 

print(recMC([1,5,10,25],63))

可惜的是,上面算法效率非常之差,实际上,找到63分钱的硬币数量,进行了67,716,925次递归。为了帮助理解,图5显示了一个较小的为了找到26的硬币377次函数调用的过程。

图中,每个节点对一次recMC函数的调用。节点上的标签表示当前求解的找钱数量,箭头上的标签表示减掉的那个硬币。跟踪图形可以得到任意的硬币组合。主要问题在于重复计算太多了,例如图中15分的数值至少有3次计算,每次计算包括52次函数调用,很明显我们很多时间浪费在重复计算上。

python数据结构与算法 26 动态规划

减少计算量的关键是记住已经得到的中间结果。一个简单方法是把已经计算出最小值的找钱保存在表里,每次计算之前,先查一下表看看这个结果是否已知。如果已知就直接套用不再计算。下面是优化的算法:

defrecDC(coinValueList,change,knownResults):

  minCoins = change

  if change in coinValueList:

     knownResults[change] = 1

     return 1

  elif knownResults[change] > 0:

     return knownResults[change]

  else:

      for i in [c for c in coinValueList if c <= change]:

        numCoins = 1 + recDC(coinValueList, change-i,

                              knownResults)

        if numCoins < minCoins:

            minCoins = numCoins

            knownResults[change] = minCoins

  return minCoins

 

print(recDC([1,5,10,25],63,[0]*64))

注意第6行增加一个测试表中是否包括了计算结果,如果没有,计算并保存到表中。这样优化后,计算63分钱的最小硬币数量,只用了221次调用。

虽然这个结果是正确的,但似乎有点黑客的味道。并且如果查看一下knownResults列表,会发现有太多的洞(0值),事实上上面的算法不能称之为动态规划,只是利用了记忆或缓冲的做法优化了我们程序的性能。

真正的动态规划算法要更加系统化地逼近问题答案。解决方法从一个硬币开始,逐步靠近我们要计算的数值,在中间过程中每个找钱数字我们都得到了所需的最少硬币数量。

我们来研究一下,对于11分的找钱计算最小需要的硬币数量,怎样填表。如图4所示的过程,从1分开始,唯一可能是1个硬币。下一行显示的是1分和2分的情况,当然,2分也只有一个答案就是2个硬币。从第5行开始事情变得有意思起来,我们有2个选择,是51分币或15分币。哪个方案好?我们查表发现对要找4分钱的情况是4个硬币,加上1个就5个硬币,或者是05分币加上15分币为1个硬币。因为最小值是1,我们在表中填上1。继续向前到表的末端考虑11分的情况,图5所示为我们要考虑的3个选项。

  1. 1个1分币加上10分找钱的最小硬币数量(1个)
  2. 1个5分币加上6分找钱的最小硬币数量(2个)
  3. 1个10分币加上1分找钱的最小硬币数量(1个)

13都给出了最小值2的答案。

python数据结构与算法 26 动态规划

python数据结构与算法 26 动态规划

下面是找零的程序代码。dpMakeChange3个参数:可用硬币的列表,找钱的数量,每个找零值所需硬币的最小数量列表。函数计算完成后,minCoin将包括小于找钱的所有找零值所需要的最小硬币数量。

Listing 8

defdpMakeChange(coinValueList,change,minCoins):

   for cents inrange(change+1):

     coinCount = cents

      for j in [c for c in coinValueList if c <= cents]:

           if minCoins[cents-j] +1< coinCount:

              coinCount = minCoins[cents-j]+1

     minCoins[cents] = coinCount

   return minCoins[change]

注意dpMakeChange不是一个递归函数,虽然我们是从递归开始的。你要明白,虽然你可以写一个递归解法,并不意味着递归就是最好或最有效率的。这个函数的大部分工作是从第4行开始的循环中完成的,这个循环里,我们通过cents值考虑了所有可能的硬币找法,象在上面我们找11分的例子一样,把小于change值的所有最小找法存放在minCoins列表里。


虽然上面的找零算法很好地给出了最小的数量,但是不能用于找零操作,因为没有跟踪用到的零钱。我们可以很容易扩展算法,通过在minCoin表的入口记住最后加入的每个硬币来得到实际找的零钱。只要知道了最后加入的硬币,就可以减去这个硬币值得到上一个表的入口,依此类推直到最开始的第一个硬币。

下面代码显示了dpMakeChange的改进算法,加入一个函数printCoins用来倒查表格来打印每种硬币的使用数量。这里用到了我们解决Lower Elbonia问题时的做法。Main函数前两行给出了找零的数量和硬币种类列表,下面两行创建两个后面要存储结果的列表,coinsUsed是存储找零用到的硬币的列表,coinCount是找零最小硬币数量的列表,这个表里找零值就是表索引。

注意我们打印硬币的时候,是从coinsUsed数组里直接输出的,例如查表63的位置,得到21,然后63-21=42,查找42的位置,再次得到21,最后在21的位置也查到了21,也就是说是321分的硬币。

defdpMakeChange(coinValueList,change,minCoins,coinsUsed):

  for cents in range(change+1):

     coinCount = cents

     newCoin = 1

     for j in [c for c in coinValueList if c <= cents]:

            if minCoins[cents-j] + 1 <coinCount:

               coinCount = minCoins[cents-j]+1

               newCoin = j

     minCoins[cents] = coinCount

     coinsUsed[cents] = newCoin

  return minCoins[change]

 

def printCoins(coinsUsed,change):

  coin = change

  while coin > 0:

     thisCoin = coinsUsed[coin]

     print(thisCoin)

     coin = coin - thisCoin

 

def main():

   amnt = 63

   clist = [1,5,10,21,25]

   coinsUsed = [0]*(amnt+1)

   coinCount = [0]*(amnt+1)

 

   print("Making change for",amnt,"requires")

   print(dpMakeChange(clist,amnt,coinCount,coinsUsed),"coins")

   print("They are:")

   printCoins(coinsUsed,amnt)

   print("The used list is as follows:")

   print(coinsUsed)

 

main()



python数据结构与算法 26 动态规划,布布扣,bubuko.com

python数据结构与算法 26 动态规划

上一篇:java 基本知识


下一篇:php:关于微信一直无法通过php token验证, ehcho echo $echoStr; 总是提示不对的问题.