1 CUT-ROD(p,n) 2 if n == 0 3 return 0 4 q=负无穷 5 for i = 1 to n 6 q=max(q,p[i]+CUT-ROD(p,n-i)) 7 return q
上面只用了分治策略,这个算法的性能是不好的T(n)=2n,在子问题的求解中不少都是重复的。算法
动态规划就是避免这些重复。通常有两个思路1.记录中间解 性能
1 MEM-CUT-ROD 2 let r[0..n] be a new array 3 for i = 0 to n 4 r[i]=负无穷 5 return MEM-CUT-ROD-AUX(p,n,r) 6 7 MEM-CUT-ROD-AUX(p,n,r) 8 if r[n]>=0 9 return r[n] 10 if n==0 11 q=0 12 else q=负无穷 13 for i=1 to n 14 q=max(q,p[i]+MEM-CUT-ROD-AUX(p,n-i,r)) 15 r[n]=q 16 return q
2.经过对问题求解顺序的合理安排,达到避免重复优化
1 BOTTOM-UP-CUT-ROD(p,n) 2 let r[0..n] be a new array 3 r[0]=0 4 for j=1 to n 5 q=负无穷 6 for i=1 to j 7 q=max(q,p[i]+r[j-i]) 8 r[j]=q 9 return r[n]
1 EXTEND-BOTTOM-UP-CUT-ROD(p,n) 2 let r[0..n] and s[0..n] be new arrays 3 r[0]=0 4 for j = 1 to n 5 q=负无穷 6 for i =1 to j 7 if q < p[i]+r[j-i] 8 q=p[i]+r[j-i] 9 s[j]=i 10 r[j]=q 11 return r and s 12 13 PRINT-CUT-ROD-SOLUTION(p,n) 14 (r,s)=EXTEND-BOTTOM-UP-CUT-ROD(p,n) 15 while n >0 16 print s[n] 17 n=n-s[n]
动态规划真的是个好东西,书中还讲了动态规划在最长公共子序列以及最优二叉搜索树上的应用,此处不细讲。spa
使用动态规划方法求解的最优化问题应该具有两个要素:最优子结构和子问题重叠。code