埃德蒙顿华人社区-Edmonton China

 找回密码
 注册
查看: 7420|回复: 2

动态规划算法求解硬币找零问题

[复制链接]
鲜花(0) 鸡蛋(0)
发表于 2012-7-25 01:28 | 显示全部楼层 |阅读模式
老杨团队,追求完美;客户至上,服务到位!
动态规划的基本思想是将待求解问题分解成若干个子问题,先求解子问题,并将这些子问题的解保存起来,如果以后在求解较大子问题的时候需要用到这些子问题的解,就可以直接取出这些已经计算过的解而免去重复运算。保存子问题的解可以使用填表方式,例如保存在数组中。. T, E  U. l- ~/ }0 U7 d3 S
    用一个实际例子来体现动态规划的算法思想——硬币找零问题。
8 h  l  [* i6 `4 i4 w    硬币找零问题描述:现存在一堆面值为 V1、V2、V3 … 个单位的硬币,问最少需要多少个硬币才能找出总值为 T 个单位的零钱?假设这一堆面值分别为 1、2、5、21、25 元,需要找出总值 T 为 63 元的零钱。- s7 W; _( N6 D# C
    很明显,只要拿出 3 个 21 元的硬币就凑够了 63 元了。4 Z5 c7 R: X$ }' z$ N' X8 d
    基于上述动态规划的思想,我们可以从 1 元开始计算出最少需要几个硬币,然后再求 2 元、3元…每一次求得的结果都保存在一个数组中,以后需要用到时则直接取出即可。那么我们什么时候需要这些子问题的解呢?如何体现出由子问题的解得到较大问题的解呢?
" g1 V1 W! {5 v. N. {6 R. o! Z    其实,在我们从 1 元开始依次找零时,可以尝试一下当前要找零的面值(这里指 1 元)是否能够被分解成另一个已求解的面值的找零需要的硬币个数再加上这一堆硬币中的某个面值之和,如果这样分解之后最终的硬币数是最少的,那么问题就得到答案了。5 O: X" w+ N7 K' G3 {8 Z
    单是上面的文字描述太抽象,先假定以下变量:- u! [: K" X% E1 s  K, w& b
    values[] : 保存每一种硬币的币值的数组
* V* s$ Q0 o: M( o' R0 q# \    valueKinds :币值不同的硬币种类数量,即values[]数组的大小# b/ _4 w  f0 ?  R' V9 a
    money : 需要找零的面值, u% V  P# F% K9 ]8 T
    coiUsed[] : 保存面值为 i 的纸币找零所需的最小硬币数
, r: a7 C: r' S" e; [* ^    算法描述:6 w7 e0 G3 X' `! P6 |( O
    当求解总面值为 i 的找零最少硬币数 coiUsed[ i ] 时,将其分解成求解 coiUsed[ i – cents]和一个面值为 cents 元的硬币,由于 i – cents < i , 其解 coiUsed[ i – cents] 已经存在,如果面值为 cents 的硬币满足题意,那么最终解 coiUsed[ i ] 则等于 coiUsed[ i – cents] 再加上 1(即面值为 cents)的这一个硬币。
& Q2 m) z7 _+ ^5 m2 ?    下面用代码实现并测试一下:
8 P( [  h- |* M; z& Y1 N5 `# @. G    public class CoiChange {. O3 P, K1 b( p3 _0 ^2 Y. |9 y
    /**$ e. S8 \( s+ ?% `/ Y
    * 硬币找零:动态规划算法: Z( e( z# w$ V: C0 O+ f
    *
. V- n: G& J. Q    * @param values
$ b! k4 C, I2 N- b6 o. n    *  :保存每一种硬币的币值的数组  V0 ~$ |5 l5 p1 z
    * @param valueKinds
% ~- X+ @7 U3 F/ D9 `- j- `, ]7 H    *  :币值不同的硬币种类数量,即coinValue[]数组的大小. W- J( ~8 ?9 E6 L& d9 F
    * @param money4 ~- a- e  _. x
    *  :需要找零的面值
1 ?; L5 Z' _8 k: V* f! J. U0 h    * @param coiUsed+ P. v  E& G0 g  i8 C: p9 A6 N& G: D
    *  :保存面值为i的纸币找零所需的最小硬币数8 O+ b) A4 B% L( s4 ^# v. n
    */
6 x" `" `! t( }    public static void makeChange(int[] values, int valueKinds, int money,
# ]: |! V, U6 p% L- z    int[] coiUsed) {0 N  T+ i: i9 `  Q' t
    coiUsed[0] = 0;
, x' W2 a2 F4 I; [" k* U9 E* A1 V    // 对每一分钱都找零,即保存子问题的解以备用,即填表  C& j- S* o) B+ j. U
    for (int cents = 1; cents <= money; cents++) {3 j0 P: Y5 L7 a' P) ?9 X+ [8 b1 x
    // 当用最小币值的硬币找零时,所需硬币数量最多2 G: N: R6 R; V8 c) I4 x# x
    int minCoi = cents;
. @5 C2 R$ G! A8 c    // 遍历每一种面值的硬币,看是否可作为找零的其中之一. R4 [8 T7 g9 t) h# b, F$ q  M
    for (int kind = 0; kind < valueKinds; kind++) {) H' N5 V# q' I  s* m9 Z# d
    // 若当前面值的硬币小于当前的cents则分解问题并查表
( @7 R$ n5 n$ i9 @    if (values[kind] <= cents) {
/ H1 c) |; Y! X6 h    int temp = coiUsed[cents - values[kind]] + 1;6 l, {; T4 e0 o5 W) ^2 o/ o8 @* s
    if (temp < minCoi) {
2 X5 j7 O8 H. l; `0 U/ b8 \/ s     minCoi = temp;( N# P+ r: S( a3 r
    }, f7 Y- }( G! R0 u9 P  g
    }
1 A! p: r7 N) [/ ^8 [7 m2 d    }# n0 p3 m. y! u9 D$ u  U3 r* p0 l
    // 保存最小硬币数( ^+ P4 C$ t7 D9 V& v7 k
    coiUsed[cents] = minCoi;
2 m9 t1 W6 H: V* H! k0 U    System.out.println("面值为 " + (cents) + " 的最小硬币数 : "
3 s5 k  W4 ?9 s" t$ Y7 l& l- e  L    + coiUsed[cents]);
' v% e( I8 a3 V    }
+ e: K" \' D6 U: W  O! I5 P3 k    }5 f2 I, P6 z; f6 f; e7 h) n' v. V
    public static void main(String[] args) {! N9 |' P' U# d& c5 D& o
    // 硬币面值预先已经按降序排列
6 r2 t  _5 `9 p% _' p8 U( ~    int[] coinValue = new int[] { 25, 21, 10, 5, 1 };
2 {$ {9 g5 m/ v* P: h    // 需要找零的面值
8 b4 H" p  U5 I    int money = 63;
& D; O$ e" q# H' B9 x5 M    // 保存每一个面值找零所需的最小硬币数,0号单元舍弃不用,所以要多加1* x5 T5 ~, D; T6 B
    int[] coiUsed = new int[money + 1];
5 R" w  u- M, G1 `    makeChange(coinValue, coinValue.length, money, coiUsed);
- W0 i4 ]4 M" _% A* r- d    }0 D0 t) t! h4 K- k0 V
    }
  m$ Q1 s& G$ H" y4 |$ e1 u    测试结果:
, l: w  H+ K- d. ?( q    面值为 1 的最小硬币数 : 1  面值为 2 的最小硬币数 : 2  面值为 3 的最小硬币数 : 3  面值为 4 的最小硬币数 : 4  面值为 5 的最小硬币数 : 1  面值为 6 的最小硬币数 : 2  ...  ...  面值为 60 的最小硬币数 : 3  面值为 61 的最小硬币数 : 4  面值为 62 的最小硬币数 : 4  面值为 63 的最小硬币数 : 33 V, m0 E, E' H. S6 U6 w$ M- |# h
     上面的代码并没有给出具体应该是哪几个面值的硬币,这个可以再使用一些数组保存而打印出来。
6 c; s( E; C9 M  r# `  H* ^    参考lovebeyond19831993.blog.163.com/blog/static/131811444201101901819833/
) i7 ]. |; K. M8 S    问题:
: x' r) ?' `, O+ S( A    假设有n种面值不同的硬币,个个面值存于数组T[1:n]中,现在用这些硬币来找钱,各种硬币的使用个数不限;1 t! x0 D- I% q: {5 c, g; m: R
    求对于给定的钱数N,最少可以由几枚硬币组成,并输出硬币序列。
% N; Y2 ~4 R% W2 K1 W- [) C0 D* ?% ]    分析:
0 O3 e- ?, f' i2 _% O5 U    假设对于i = 1...N-1, 所需最少的硬币数Count(i) 已知, 那么对于N,所需的硬币数为Min( Count(i) + Count(N-i)) , i=1...N-1;; C+ c" t* t1 |  X) E% _
    于是一个直观的方法是用递归计算。& v0 R  [4 V/ R& t
    但是,递归过程中,每次计算Count(i),都会重复计算 Count(1)....Count(i-1); 这样时间复杂度就是O(N2);7 Z6 S* M7 I5 Z6 w7 D
    事实上,这是一个典型的动态规划问题,我们可以从1开始记录下每个钱数所需的硬币枚数,避免重复计算,为了能够输出硬币序列,我们还需要记录下每次新加入的硬币。
  g! o( h( g) T    下面给出用动态规划解决此问题的递推式:1 T3 _3 H7 [0 Y# Y6 i
    参数说明: 当只用面值为T[1],T[2],…T[n]来找出钱j时,所用的硬币的最小个数记为C(i,j),则C(i,j)的递推方程为:1 @3 X: e  O' I% Q9 l8 e; C
    算法设计:' r" D' ^; g1 {+ _( t4 ^/ h
    #include <iostream># w2 m! O/ C$ g. U
    using namespace std;' p) z1 K. i) ]- k+ C
    int min_number_of_stamps(cot int* array, size_t array_size, int request) {
5 o" X# ?0 w/ q" C/ U" P3 g' h" n    if (array == NULL || array_size <= 0 || request < 1) {) b7 z5 Q9 y: ?/ Q1 ~  _% a, d
    return -1;
! o7 G7 @3 s8 ~5 a% @. p; S0 ~    }# l6 `- `0 e5 _+ \1 ?: p# y* W
    int* numofStamps = new int[request + 1];& H0 W* ^2 A7 j3 v
    numofStamps[0] = 0;/ G8 _5 J5 ?) E0 h3 y
    for(int i = 1; i <= request; i++){
" ?0 {% g/ g: b0 H+ K/ }    int cn = i;
' C( Z) Y- A; s2 X9 @& s    for(size_t j = 0; j < array_size; j++){
* `, E6 \( W, u# z! g/ h5 s: c     if(array[j] <= i){
9 h# k0 i' E3 |- V2 x  T8 L: \4 y    int tmp = numofStamps[i - array[j]] + 1;& N) ~3 k) i1 t: d$ J
    if(tmp < cn)" T) `" r# c# @# A# e* F$ o% y  p; I  V
    cn = tmp;& f9 o& A9 `$ o% z* w
     }" w* k: t. U" }# Q
    }
' W6 D) l# ~6 [, F  ~    numofStamps[i] = cn;
: A$ ?% X* d4 U. J    }" V( S0 A& Q# R) n6 K/ M
     int t = numofStamps[request];
8 \/ q# l/ ?5 l9 k    delete[] numofStamps;: z  T/ m' M. p& w- D% J& n
    return t;
! ]" C0 i) x9 i. y    }
4 p3 ~' F8 h" u  L0 u    int main() {
0 f5 y- G  M9 X( A' }) x, F8 a2 n$ z2 O6 J    int array[] = { 90, 30, 24, 15, 12, 10, 5, 3, 2, 1 };  b0 X6 z: A0 H+ \# @
    int array_size = sizeof(array) / sizeof(int);  s( x* }; ~- E6 y" U
    int t = min_number_of_stamps(array, array_size, 32);$ X; }" I9 ]; A7 ]" A0 N
    cout << t << endl;' F. I( i2 ?1 c
    }6 E9 b( g9 K" ]
   
鲜花(3) 鸡蛋(0)
发表于 2012-7-31 09:48 | 显示全部楼层
理袁律师事务所
鲜花(3) 鸡蛋(0)
发表于 2012-7-31 09:52 | 显示全部楼层
您需要登录后才可以回帖 登录 | 注册

本版积分规则

联系我们|小黑屋|手机版|Archiver|埃德蒙顿中文网

GMT-7, 2026-9-26 01:17 , Processed in 0.337628 second(s), 13 queries , Gzip On, APC On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表