devbox
正经夜光杯
这个屌丝很懒，什么也没留下！
关注作者
热门标签
热门文章
当前位置:   article > 正文

动态规划的应用_动态规划应用

作者：正经夜光杯 | 2024-06-29 21:48:14

赞

踩

动态规划应用

动态规划的基本思想：

在用分治法求解时，有些问题被重复计算了好多次。如果能够保存已解决的子问题的答案，而在需要时再找出已求得的答案，就可以避免大量重复计算。为达到这个目的，可以用一个表来记录所有已解决的子问题的答案，只要他被计算过，就将其填入表中。

矩阵连乘问题

这种计算次序可以用加括号的方式来确定。若一个矩阵连乘积的计算次序完全确定,也就是说该连乘积已完全加括号,则可以依此次序反复调用2个矩阵相乘的标准算法计算出矩阵连乘积。


 

2矩阵乘积的标准算法：

  1. public static void matrixMultiply(int [][]a,int [][]b,int [][]c,int ra,int ca,int rb,int cb){
  2.         for (int i = 0;i<ra;i++){
  3.             for (int j = 0;j<cb;j++){
  4.                 int sum = 0;
  5.                 for (int k = 0;k<ca;k++){
  6.                     sum = sum + a[i][k]*b[k][j];
  7.                 }
  8.                 c[i][j] = sum;
  9.             }
  10.         }
  11.     }

计算3个矩阵{A1,A2,A3 }连乘积的例子。设这3个矩阵的维数分别为10 X 100,100 X 5和5 X 50。若按第一种加括号方式((A1,A2)A3)计算,3个矩阵连乘积需要的数乘次数为10 X 100 X 5 + 10  X 5 X 50 = 7500。若按第二种加括号方式(A1(A2A3))计算,3个矩阵连乘积总共需要10X5X50+ 10 X 100 X 50 = 75 000次数乘。第二种加括号方式的计算量是第一种加括号方式计算量的10倍。

 数组p(n)表示矩阵的维数，上述数组可表示为

p(0)=30,p(1)=35,p(2)=15,p(3)=5,p(4)=10,p(5)=20,p(6)=25

数组m[i][j]表示最少乘的次数

s[i][j]表示m[i][j]最少乘次数时断开点k的值

获取最优值数组与记录最优断开位置的数组s 

输入参数{p0,p1,p2,p3....}存储于数组p中

特别注意m数组的构造顺序问题，不是一行一行的构造。

  1. public static void maxtrixChain(int []p,int [][]m,int [][]s) {
  2.         int len = p.length - 1;
  3.         for (int i = 1;i<=len;i++){
  4.             m[i][i] = 0;//m数组的第一次构造（对角线）
  5.         }
  6.         //对m数组的赋值从对角线开始向顶点移动，越来越少
  7.         for (int r = 2;r <= len;r++){//对m数组的第r次赋值及相应的顺序。
  8.             for (int i =1;i <= len - r + 1;i++){
  9.                 int j = i + r - 1;
  10.                 m[i][j] =m[i][i]+ m[i+1][j] + p[i-1]*p[i]*p[j];
  11.                 s[i][j] = i;
  12.                 for (int k = i + 1;k < j;k++){
  13.                     int t =m[i][k]+ m[k+1][j] + p[i-1]*p[k]*p[j];
  14.                     if (t < m[i][j]){
  15.                         m[i][j] = t;
  16.                         s[i][j] = k;
  17.                     }
  18.                 }
  19.             }
  20.         }
  21.     }

算法Trackback利用算法matrixChain计算出的断点矩阵指示的加括号方式输出计算矩阵连乘的最优计算次序。

  1. public static void Traceback(int [][]s,int i,int j){
  2.         //m[i][j]的最优断点s[i][j]，将其分为m[i][ s[i][j] ]与 m[ s[i][j]+1 ][j]
  3.         //按其规律递归，无限利用断点划分
  4.         if (i==j) return;
  5.         Traceback(s,i,s[i][j]);//m[i][ s[i][j] ]部分的递归
  6.         Traceback(s,s[i][j]+1,j);//m[ s[i][j]+1 ][j]部分的递归
  7.  
  8.         System.out.println("A"+i+","+s[i][j]+"and A"+(s[i][j]+1)+","+j);
  9.     }

单词划分问题

给定一个非空字符串 s 和一个包含非空单词的列表 wordDict，判定 s 是否可以被空格拆分为一个或多个在字典中出现的单词。

输入: s = “leetcode”, wordDict = [“leet”, “code”] 输出: true 解释: 返回 true 因为 “leetcode” 可以被拆分成 “leet code”。

将大问题分解为规模小一点的子问题：

  • 前 i 个字符的子串，能否分解成单词
  • 剩余子串，是否为单个单词。

dp[i]：长度为i的s[0:i-1]子串是否能拆分成单词。题目求:dp[s.length]

再利用j对s[0,i]进行分割，s[0,i]能否被划分为字典里存在的词取决于：

  • 前缀部分s[0,j-1]的dp[j]是否为true
  • 剩余部分s[j,i]是否能拆分为单词（原问题）

最长公共子序列

设序列X={x_{1},x_{2},x_{3},x_{4},x_{5}....x_{m}}，Y={y_{1},y_{2},y_{3},y_{4}....y_{n}}的最长公共子序列Z={z_{1},z_{2},z_{3}.....z_{k}},则

解释（1）,若Xm=Yn，即表示序列X与序列Y的公共子序列的最后一项相同，即z的最后一项也是这个，因此求Xm与Yn的公共子序列Zk变成求 

 用c[i][j]来表示序列Xi与Yi的最长公共子序列的长度。

c[i-1][j-1]+1表示Xm=Yn=Zk

max{c[i][j-1],c[i-1][j]}表示Z是Xi-1与Yi的最长子序列或者Z是Xi与Yi-1的最长子序列

b[i][j]记录c[i][j]的值是由哪一个子问题的解得到的。

计算最优值（计算最优值的数组c[i][j]）

  1. public static int lcsLength(char []x,char []y,int[][]b){
  2.         int m = x.length-1;
  3.         int n = y.length-1;
  4.         int [][]c = new int[m+1][n+1];
  5.         for (int i=1;i<=m;i++) c[i][0]=0;//y序列为空
  6.         for (int i=1;i<=n;i++) c[0][i]=0;//x序列为空
  7.         for (int i=1;i<=m;i++){
  8.             for (int j=1;j<=n;j++){
  9.                 if (x[i]==y[j]){//公共元素
  10.                     c[i][j] = c[i-1][j-1]+1;
  11.                     b[i][j]=1;//子问题1
  12.                 }else {
  13.                     if (c[i-1][j]>=c[i][j-1]){
  14.                         c[i][j] = c[i-1][j];
  15.                         b[i][j] = 2;
  16.                     }else {
  17.                         c[i][j] = c[i][j-1];
  18.                         b[i][j] = 3;
  19.                     }
  20.                 }
  21.             }
  22.         }
  23.         return c[m][n];
  24.     }

构造最长公共子序列（打印输出最长子序列）

  1. public static void lcs(int i,int j,char []x,int [][]b){
  2.         if (i==0||j==0) return;
  3.         if (b[i][j]==1){
  4.             lcs(i-1,j-1,x,b);
  5.             System.out.println(x[i]);//此时当b[i][j]==1，第i项必为公共子序列的元素
  6.         }
  7.         if (b[i][j]==2) lcs(i-1,j,x,b);
  8.         if (b[i][j]==3) lcs(i,j-j,x,b);
  9.     }

凸多边形最优三角剖分

给定凸多边形P={V0,V1,V2,,,,Vn-1}，以及定义在由多边形的边和弦组成的三角形上的权函数w。要求确定该凸多边形的三角剖分，使得该三角剖分所对应的权，即该三角剖分中诸三角形上权之和为最小。

权函数W    w(Vi,Vj,Vk) = Vj*Vk+Vi*Vj+Vk*Vj

数组t[i][j]表示凸多边形{Vi-1，Vi，，，，Vj}的最优三角剖分对应的权函数值。即最优值。

设多边形{Vi-1，Vi}的权值为0

最优子结构

实现与矩阵连乘类似

数组t[i][j]表示凸多边形{Vi-1，Vi，，，，Vj}的最优三角剖分对应的权函数值。

s[i][j]表示t[i][j]对应的k的位置

t[i][j]的构建顺序

  1. public static void minWeightTriangulation(int n,int[][]t,int [][]s){//t[][]表示
  2.         for (int i=1;i<=n;i++) t[i][i]=0;
  3.         for (int r=2;r<=2;r++){//赋值是一个三角形状
  4.             for (int i=1;i<=n-r+1;i++){
  5.                 int j = i+r-1;
  6.                 t[i][j] = t[i][i]+t[i+1][j]+w(i-1,i,j);//给一个初值便于比较大小
  7.                 s[i][j] = i;
  8.                 for (int k=i+1;k<i+r-1;k++){
  9.                     int u =t[i][k]+t[k+1][j]+w(i-1,k,j);
  10.                     if (u<t[i][j]){
  11.                         t[i][j] = u;
  12.                         s[i][j] = k;
  13.                     }
  14.                 }
  15.             }
  16.         }
  17.     }

电路布线问题

  1.     public  static void mns(int []c,int[][]size) {
  2.         int n=c.length-1;
  3.         for(int j=0;j<c[1];j++)
  4.             size[1][j]=0;
  5.         for(int j=c[1];j<=n;j++) {
  6.             size[1][j]=1;
  7.         }
  8.         for(int i=2;i<n;i++) {
  9.             for(int j=0;j<c[i];j++) {
  10.                 size[i][j]=size[i-1][j];
  11.             }
  12.             for(int j=c[i];j<=n;j++) {
  13.                 size[i][j]=Math.max(size[i-1][j], size[i-1][c[i]-1]+1);
  14.             }
  15.             size[n][n]=Math.max(size[n-1][n], size[n-1][c[n]-1]+1);
  16.  
  17.         }
  18.     }
  19.  
  20.     //构造最优解
  21.     public static int traceback(int []c,int [][]size,int []net) {
  22.         int n=c.length-1;
  23.         int j=n;
  24.         int m=0;
  25.         for(int i=n;i>1;i--) {
  26.             if(size[i][j]!=size[i-1][j]) {
  27.                 net[m++]=i;
  28.                 j=c[i]-1;
  29.             }}
  30.  
  31.         if(j>=c[1])
  32.         {
  33.             net[m++]=1;
  34.         }
  35.  
  36.  
  37.         return m;
  38.     }

0-1背包问题

  1. /*
  2.     v[]  价值
  3.     W[]  质量
  4.     c  总质量
  5.     m 存储m(i,j)的值
  6.      */
  7.     public static void knapack(int []v,int []w,int c,int [][]m){
  8.         int n = v.length-1;//物品种类数
  9.  
  10.         // i=n
  11.         int Wn = Math.min(w[n],c);//判断w(n)与c的大小
  12.         for (int j=0;j<=c;j++){
  13.             if (j<Wn) m[n][j]=0;
  14.             if (j>=Wn) m[n][j]=v[n];
  15.         }
  16.         //i=n-1,n-2,,,,,1
  17.         for (int i=n-1;i>=1;i--){
  18.             int Wi = Math.min(w[i],c);
  19.             for (int j=1;j<=c;j++){
  20.                 if (j>=Wi) m[i][j]=Math.max(m[i+1][j],m[i+1][j-w[i]]+v[i]);
  21.                 if (j<Wi) m[i][j]=m[i+1][j];
  22.  
  23.             }
  24.         }
  25.     }

m[1][c]则是背包问题的最优值

x[i]判断物品i是否放入背包，1放入，0不放入

 根据次可得，

若m( i , j ) = m( i+1 , j )可以得知i没有放入背包

若m( i , j ) = m( i + 1 , j - Wi ) + Vi可以得知将物品i放入了背包，此时背包容量为j-Wi

若m(n,j) = 0 得知物品n没有放入背包

若m(n,j) > 0 得知物品n放入背包

  1. //m[1][c]则是背包问题的最优值
  2.      /*
  3.     m 存储m(i,j)的值
  4.     W[]  质量
  5.     c  总质量
  6.     x[i]判断物品i是否放入背包，1放入，0不放入
  7.      */
  8.     public static void traceback(int[][]m,int[]w,int c,int []x){
  9.         int n = w.length-1;
  10.         for (int i=1;i<n;i++){
  11.             if (m[i][c]==m[i+1][c]) x[i]=0;//物品i没有放入背包
  12.             else {//物品i放入背包
  13.                 x[i] = 1;
  14.                 c=c-w[i];
  15.             }
  16.         }
  17.         //最后一个特殊讨论
  18.         if (m[n][c]>0) x[n]=1;
  19.         else x[n]=0;
  20.     }

  

二维背包问题

给定n种物品和一个背包。物品i的重量是wi，体积是bi，其价值是vi，背包的容量为c,容积为d。问：应该如何选择装入背包中的物品，使得装入背包中的物品的总价值最大？在选择装入背包额物品时，对每种物品i只有两种选择，即装入背包或不装入背包。不能将物品i装入背包多次，也不能只装入部分的物品i。

输入的第一行包含三个整数n, c ，d，分别表示物品的个数，背包能装重量和背包的容积。w[],v[],b[]分别表示各物品的重量、价值、体积。

与0-1问题类似，在此基础上再加上一个体积限制。

m[1][c][d]则是背包问题的最优值

x[i]判断物品i是否放入背包，1放入，0不放入

若m( i , j , k) = m( i+1 , j , k )可以得知i没有放入背包

若m( i , j , k) = m( i + 1 , j - Wi , k - Di) + Vi可以得知将物品i放入了背包，此时背包容量为j-Wi ,体积为k-Di

若m(n,j,k) = 0 得知物品n没有放入背包

若m(n,j,k) > 0 得知物品n放入背包

  1. //调用部分
  2. 输入 5 8 8
  3. v[]={2,1,4,3,5}
  4. w[]={1,4,2,3,5}
  5. b[]={2,3,5,2,4}
  6.  
  7. public static void main(String[] args) {
  8.         int c = 8;
  9.         int d = 8;
  10.         int []v = {0,2,1,4,3,5};//0用来占位，程序从下标1开始
  11.         int []w = {0,1,4,2,3,5};
  12.         int []b = {0,2,3,5,2,4};
  13.         int [][][]m = new int[v.length][c+1][d+1];
  14.         knapack(v,w,b,c,d,m);
  15.         System.out.println("装入背包中物品总价值最大为："+m[1][8][8]);
  16.         int x[] = new int[v.length];
  17.         traceback(m,w,b,c,d,x);
  18.         System.out.print("装入的物品的序号为:");
  19.         for (int i=0;i<v.length;i++){
  20.             if (x[i]==1) System.out.print(i+" ");
  21.         }
  22.     }
  23.  
  24.  
  25.  
  26.     public static void knapack(int []v,int []w,int []b,int c,int d,int [][][]m){
  27.         int n = v.length-1;
  28.         int Wn = Math.min(w[n],c);
  29.         int Dn = Math.min(b[n],d);
  30.         for (int j=0;j<=c;j++){
  31.             for (int k=0;k<=d;k++){
  32.                 if (j<Wn || k<Dn) m[n][j][k]=0;
  33.                 else m[n][j][k]=v[n];
  34.             }
  35.         }
  36.         //i=n-1,n-2,,,,,1
  37.         for (int i=n-1;i>=1;i--){
  38.             int Wi = Math.min(w[i],c);
  39.             int Di = Math.min(b[i],d);
  40.             for (int j=1;j<=c;j++){
  41.                 for (int k=1;k<=d;k++){
  42.                     if (j>=Wi&&k>=Di) m[i][j][k]=Math.max(m[i+1][j][k],m[i+1][j-w[i]][k-b[i]]+v[i]);
  43.                     if (j<Wi||k<Di) m[i][j][k]=m[i+1][j][k];
  44.                 }
  45.             }
  46.         }
  47.     }
  48.  
  49.  
  50.  
  51.     public static void traceback(int[][][]m,int[]w,int b[],int c,int d,int []x){
  52.         int n = w.length-1;
  53.         for (int i=1;i<n;i++){
  54.             if (m[i][c][d]==m[i+1][c][d]) x[i]=0;//物品i没有放入背包
  55.             else {//物品i放入背包
  56.                 x[i] = 1;
  57.                 c=c-w[i];
  58.                 d=d-b[i];
  59.             }
  60.         }
  61.         //最后一个特殊讨论
  62.         if (m[n][c][d]>0) x[n]=1;
  63.         else x[n]=0;
  64.     }

声明：本文内容由网友自发贡献，转载请注明出处：【wpsshop博客】
推荐阅读
相关标签

Copyright © 2003-2013 www.wpsshop.cn 版权所有，并保留所有权利。

  

闽ICP备14008679号