动态规划求解资源最优分配问题

为第i个工程计算最优时,只需将分配若干给前i-1个工程、剩下的留给当前工程得到的各种分配方案,和只分配给当前工程的方案,取其中的最大值即可;而第一个工程解就是只分配给自身的解。

Code:
  1. #include<iostream>   
  2. using namespace std;   
  3.     
  4. #define M 10    //  最多资源份量   
  5. #define N 10    //  最多工程数量   
  6.     
  7. int g[M+1][N];  //        资源分配利润函数   
  8. int solve[N];      // 最优分配方案   
  9. int m, n;              //     资源数与工程数   
  10.     
  11. int dp(){   
  12.         int f[M+1][N];  //        资源分配最多利润   
  13.         int d[M+1][N];  //        资源分配方案   
  14.         int best;                     //     资源分配最优方案可得利润best,   
  15.         int bi, bj;          //   最优时资源分配给了前bi个项目,共分配bj份资源   
  16.         int i,j,k;   
  17.     
  18.         //      先得到只有第一个工程时的利润   
  19.         for(i=0; i<=m; i++){   
  20.                 f[0][i] = g[0][i];   
  21.                 d[0][i] = i;   
  22.         }   
  23.     
  24.         //      计算n个工程的资源分配表   
  25.         for(i=1; i<n; i++){   
  26.                 f[i][0] = g[i][0] + f[i-1][0], d[i][0] = 0;     //   当为0时解直接相加   
  27.                 for(j=1; j<=m; j++){   
  28.                         f[i][j] = f[i][0], d[i][j] = 0; //       初始化为不分配给该工程资源时的解   
  29.                         for(k=0; k<=j; k++){        //  计算分配给其若干个资源时的最优解   
  30.                                 if(f[i-1][j-k] + g[i][k] > f[i][j]){   
  31.                                         f[i][j] = f[i-1][j-k] + g[i][k];   
  32.                                         d[i][j] = k;   
  33.                                 }   
  34.                         }   
  35.                 }   
  36.         }   
  37.     
  38.         //      寻找最优分配方案   
  39.         for(i=0; i<n; i++){   
  40.                 for(j=0; j<=m; j++){   
  41.                         if(f[i][j]>best){   
  42.                                 best = f[i][j];   
  43.                                 bi = i, bj = j;   
  44.                         }   
  45.                 }   
  46.         }   
  47.     
  48.         //      计算每个工程各消耗多少资源   
  49.         memset(solve, 0, sizeof(solve));   
  50.         k = bj;   
  51.         for(i=bi; i>=0; i--){   // bi以后的项目都不用分配   
  52.                 solve[i] = d[i][k];   
  53.                 k-=solve[i];   
  54.         }   
  55.     
  56.         return best;   
  57. }   
  58.     
  59. int main(){   
  60.         int i,j;   
  61.     
  62.         cin>>m>>n;   
  63.           
  64.         for(i=0; i<n; i++){   
  65.                 for(j=0; j<=m; j++){   
  66.                         cin>>g[i][j];   
  67.                 }   
  68.         }   
  69.     
  70.         cout<<"最优分配方案可得利润:"<<dp()<<endl;   
  71.     
  72.         cout<<"最优分配方案:"<<endl;   
  73.         for(i=0; i<n; i++){   
  74.                 cout<<"工程"<<(i+1)<<"分配:"<<solve[i]<<"资源"<<endl;   
  75.         }   
  76.     
  77.         return 0;   
  78. }  

 

  • 1
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
实验课程:算法分析与设计 实验名称:用动态规划求解资源分配问题 (验证型实验) 实验目标: (1)掌握用动态规划方法求解实际问题的基本思路。 (2)进一步理解动态规划方法的实质,巩固设计动态规划算法的基本步骤。 实验任务: (1)设计动态规划算法求解资源分配问题,给出算法的非形式描述。 (2) 在Windows环境下用C 语言实现该算法。计算10个实例,每个实例中n=30, m=10, Ci j为随机产生于范围(0,103)内的整数。记录各实例的数据及执行结果(即最优分配方案、最优分配方案的值)、运行时间。 (3)从理论上分析算法的时间和空间复杂度,并由此解释相应的实验结果。 实验设备及环境: PC;C/C++等编程语言。 实验主要步骤: (1) 根据实验目标,明确实验的具体任务; (2) 分析资源分配问题,获得计算其最优值的递推计算公式; (3) 设计求解问题动态规划算法,并编写程序实现算法; (4) 设计实验数据并运行程序、记录运行的结果; (5) 分析算法的时间和空间复杂度,并由此解释释相应的实验结果; 问题分析: 问题描述: 某厂根据计划安排,拟将n台相同的设备分配给m个车间,各车间获得这种设备后,可以为国家提供盈利Ci j(i台设备提供给j号车间将得到的利润,1≤i≤n,1≤j≤m) 。问如何分配,才使国家得到最大的盈利? 算法基本思想: 本问题是一简单资源分配问题,由于具有明显的最优子结构,故可以使用动态规划求解,用状态量f[i][j]表示用i台设备分配给前j个车间的最大获利,那么显然有f[i][j] = max{ f[k][j–1] + c[i-k][j] },0<=k<=i。再用p[i][j]表示获得最优解时第j号车间使用的设备数为i-p[i][j],于是从结果倒推往回求即可得到分配方案。程序实现时使用顺推,先枚举车间数,再枚举设备数,再枚举状态转移时用到的设备数,简单3重for循环语句即可完成。时间复杂度为O(n^2*m),空间复杂度为O(n*m),倘若此题只需求最大获利而不必求方案,则状态量可以减少一维,空间复杂度优化为O(n)。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值