实验二 图论、动态规划求解
西华大学上机实验报告
一、实验目的
1、掌握网络图的计算机输入，求解最小树、最短路、最大流问题。
2、掌握动态规划的基本原理及求解方法，熟悉软件操作，正确完成模型求解过程及分析过程。
二、实验内容或设计思想
1、首先将欲求的网络用计算机语言表达，再用lingo计算软件求出模型问题的解。
2、最小树应用破圈方法求解
3、最大流算法用找流量可增链的方法求解
4、网络最短路运用了动态规划，函数迭代、矩阵运算等的基本原理求解
三、实验环境与工具
计算机，lingo软件，运筹学软件。
四、实验过程或实验数据
1、最短路问题
【例题】给定10个点,各点城市之间的距离和路线如下程序所示。试用LINGO软件解网络最短路。解：（1）lingo程序代码：
!最短路问题;
model:
data:
n=10;
enddata
sets:
cities/1..n/: F; !10个城市;
roads(cities,cities)/
1,2 1,3
2,4 2,5 2,6
3,4 3,5 3,6
4,7 4,8
5,7 5,8 5,9
6,8 6,9
7,10
8,10
9,10
/: D, P;
endsets
data:
D=
6 5
3 6 9
7 5 11
9 1
8 7 5
4 10
5
7
9;
enddata
F(n)=0;
@for(cities(i) | i #lt# n:
F(i)=@min(roads(i,j): D(i,j)+F(j));
);
!显然，如果P(i,j)=1,则点i到点n的最短路径的第一步是i --> j，否则就不是。
由此，我们就可方便的确定出最短路径;
@for(roads(i,j):
P(i,j)=@if(F(i) #eq# D(i,j)+F(j),1,0)
);
End
（2）计算的部分结果为：
Feasible solution found at iteration: 0
Variable Value
N 10.00000
F( 1) 17.00000
F( 2) 11.00000
F( 3) 15.00000
F( 4) 8.000000
F( 5) 13.00000
F( 6) 11.00000
F( 7) 5.000000
F( 8) 7.000000
F( 9) 9.000000
F( 10) 0.000000
P( 1, 2) 1.000000
P( 1, 3) 0.000000
P( 2, 4) 1.000000
P( 2, 5) 0.000000
P( 2, 6) 0.000000
P( 3, 4) 1.000000
P( 3, 5) 0.000000
P( 3, 6) 0.000000
P( 4, 7) 0.000000
P( 4, 8) 1.000000
P( 5, 7) 1.000000
P( 5, 8) 0.000000
P( 5, 9) 0.000000
P( 6, 8) 1.000000
P( 6, 9) 0.000000
P( 7, 10) 1.000000
P( 8, 10) 1.000000
P( 9, 10) 1.000000
（3）结论：
2、动态规划求解
【例题】书本例6－4
解：（1）lingo程序代码：
Model:
Sets:
Nodes/a,b1,b2,b3,c1,c2,c3,d1,d2,e/:d;
Arcs(nodes,nodes)/a,b1 a,b2 a,b3 b1,c1 b1,c2 b1,c3 b2,c1 b2,c2 b2,c3 b3,c1
b3,c2 b3,c3
C1,d1 c1,d2 c2,d1 c2,d2 c3,d1 c3,d2 d1,e d2,e/:w,p;
Endsets
N = @size(nodes);
D(n) = 0;
@for(nodes(i)|i#LT#n:d(i)= @min(arcs(i,j):w(i,j)+d(j)));
@for(arcs(i,j):
P(i,j)=@if(d(i)#eq#w(I,j)+d(j),1,0)
);
Data:
W=2 5 3 7 5 6 3 2 4 5 1 5 1 4 6 3 3 3 3 4;
Enddata
End
（2）求解结果：
Feasible solution found.
Total solver iterations: 0
Variable Value
N 10.00000
D( A) 11.00000
D( B1) 11.00000
D( B2) 7.000000
D( B3) 8.000000
D( C1) 4.000000
D( C2) 7.000000
D( C3) 6.000000
D( D1) 3.000000
D( D2) 4.000000
D( E) 0.000000
W( A, B1) 2.000000
W( A, B2) 5.000000
W( A, B3) 3.000000
W( B1, C1) 7.000000 W( B1, C2) 5.000000 W( B1, C3) 6.000000 W( B2, C1) 3.000000 W( B2, C2) 2.000000 W( B2, C3) 4.000000 W( B3, C1) 5.000000 W( B3, C2) 1.000000 W( B3, C3) 5.000000 W( C1, D1) 1.000000 W( C1, D2) 4.000000 W( C2, D1) 6.000000 W( C2, D2) 3.000000 W( C3, D1) 3.000000 W( C3, D2) 3.000000 W( D1, E) 3.000000
W( D2, E) 4.000000
P( A, B1) 0.000000
P( A, B2) 0.000000
P( A, B3) 1.000000
P( B1, C1) 1.000000 P( B1, C2) 0.000000 P( B1, C3) 0.000000 P( B2, C1) 1.000000 P( B2, C2) 0.000000 P( B2, C3) 0.000000 P( B3, C1) 0.000000
P( B3, C2) 1.000000
P( B3, C3) 0.000000
P( C1, D1) 1.000000
P( C1, D2) 0.000000
P( C2, D1) 0.000000
P( C2, D2) 1.000000
P( C3, D1) 1.000000
P( C3, D2) 0.000000
P( D1, E) 1.000000
P( D2, E) 1.000000
Row Slack or Surplus
1 0.000000
2 0.000000
3 0.000000
4 0.000000
5 0.000000
6 0.000000
7 0.000000
8 0.000000
9 0.000000
10 0.000000
11 0.000000
12 0.000000
13 0.000000
14 0.000000
15 0.000000
16 0.000000
17 0.000000
18 0.000000
19 0.000000
20 0.000000
21 0.000000
22 0.000000
23 0.000000
24 0.000000
25 0.000000
26 0.000000
27 0.000000
28 0.000000
29 0.000000
30 0.000000
31 0.000000
（3）结论：
五、总结
对上机实践结果进行分析，问题回答，上机的心得体会及改进意见。
动态规划实验报告
华东师范大学计算机科学技术系上机实践报告 一、 内容与设计思想 1．对于以下5 个矩阵： M 1: 2?3, M 2: 3?6, M 3: 6?4, M 4: 4?2, M 5: 2?7 , (a) 找出这5个矩阵相乘需要的最小数量乘法的次数。 (b) 请给出一个括号化表达式，使在这种次序下达到乘法的次数最少。 输入： 第一行为正整数N,表示有N 组测试数据； 每组测试数据的第一行为n,表示有n 个矩阵,2<=n<=50； 接下去的n 行,每行有两个整数x 和y,表示第ni 个矩阵是x*y 的。 输出： 对行每组数据,输出一行,每行一个整数,最小的矩阵连乘积。 我们保证输出的结果在2^64之内。 基本思想： 对于n 个矩阵的连乘积，设其不同的计算次序为P(n)。 由于每种加括号方式都可以分解为两个子矩阵的加括号问题：(A1...Ak)(Ak+1…An)可以得到关于P(n)的递推式如下： 2．定义0/1/2背包问题为：}x p max{n 1i i i ∑=。限制条件为：c x w n 1i i i ≤∑=，且 n i 1},2,1,0{x i ≤≤∈。设f(i , y)表示剩余容量为y ，剩余物品为：i ，i+1，…，n 时的最优解的值。 1．）给出f(i , y)的递推表达式； 2．）请设计求解f(i , y)的算法，并实现你的算法； 3．）设W=[10,20,15,30]，P=[6,10,15,18],c=48，请用你的算法求解。 输入： 第一行为一个正整数N ，表示有几组测试数据。 每组测试数据的第一行为两个整数n 和M ，0
动态规划练习二
动态规划练习二 1、乘积最大 [问题描述] 在一次数学智力竞赛活动中，主持人给所有参加竞赛的选手出了一到题目：设有一个长度为N的数字串，要求选手使用K个乘号将它分成K+1个部分，找出一种分法，使得这K+1个部分的成绩最大。 同时为了帮助选手能够理解题意，主持人还举了如下一个例子： 有一个数字串：312，当N=3，K=1时有两种分法： （1）3*12=36； （2）31*2=62 这时，符合题目要求的结果是：31*2=62。现在要求设计一个程序，以求得正确的答案。 输入 Input.in文件共有两行：第一行有两个自然数N，K（2<=N<=40, 1<=K<=6）；第二行是一个长度为N的数字串。 输出 一个自然数，即所求得的最大乘积。 输入输出样例 输入（input.in） 4 2 1231 输出（ans.out） 62
2、数字加法问题 [问题描述] 有一个由数字1，2，... ，9组成的数字串（长度不超过200），问如何将M(M<=20)个加号("+")插入到这个数字串中，使所形成的算术表达式的值最小。请编一个程序解决这个问题。 注意： 加号不能加在数字串的最前面或最末尾，也不应有两个或两个以上的加号相邻。M保证小于数字串的长度。 例如：数字串79846，若需要加入两个加号，则最佳方案为79+8+46，算术表达式的值133。 [输入格式］ 从键盘读入输入文件名。数字串在输入文件的第一行行首（数字串中间无空格且不折行），Ｍ的值在输入文件的第二行行首。 [输出格式］ 在屏幕上输出所求得的最小和的精确值。 [输入输出举例］ 82363983742 3 输入输出 2170
动态规划算法实验
一、实验目的 (2) 二、实验内容 (2) 三、实验步骤 (3) 四.程序调试及运行结果分析 (5) 附录：程序清单（程序过长，可附主要部分） (7)
实验四动态规划算法的应用 一、实验目的 1．掌握动态规划算法的基本思想，包括最优子结构性质和基于表格的最优值计算方法。 2．熟练掌握分阶段的和递推的最优子结构分析方法。 3．学会利用动态规划算法解决实际问题。 二、实验内容 1.问题描述： 题目一：数塔问题 给定一个数塔，其存储形式为如下所示的下三角矩阵。在此数塔中，从顶部出发，在每一节点可以选择向下走还是向右走，一直走到底层。请找出一条路径，使路径上的数值和最大。 输入样例（数塔）： 9 12 15 10 6 8 2 18 9 5 19 7 10 4 16 输出样例（最大路径和）： 59 题目二：最长单调递增子序列问题(课本184页例28) 设有由n个不相同的整数组成的数列，记为:a(1)、a(2)、……、a(n)且a(i)<>a(j) (i<>j) 若存在i1 题目三 0-1背包问题 给定n种物品和一个背包。物品i的重量是wi，其价值为vi，背包的容量为c，。问应如何选择装入背包中的物品，使得装入背包中物品的总价值最大? 在选择装入背包的物品时，对每种物品只有两个选择：装入或不装入，且不能重复装入。输入数据的第一行分别为：背包的容量c，，物品的个数n。接下来的n 行表示n个物品的重量和价值。输出为最大的总价值。 输入样例： 20 3 11 9 9 10 7 5 输出样例 19 2.数据输入：个人设定，由键盘输入。 3.要求： 1）上述题目任选一做。上机前，完成程序代码的编写 2）独立完成实验及实验报告 三、实验步骤 1.理解算法思想和问题要求； 2.编程实现题目要求； 3.上机输入和调试自己所编的程序； 4.验证分析实验结果； 5.整理出实验报告。 实验二动态规划算法棋盘覆盖 最长公共子序列问题 一、实验目的 : 1、熟悉最长公共子序列问题的算法； 2、初步掌握动态规划算法。 二、实验内容 若给定序列X={x1,x2,…,xm}，则另一序列Z={z1,z2,…,zk}，是X的子序列是指存在一个严格递增下标序列{i1,i2,…,ik}使得对于所有j=1,2,…,k有：zj=xij。例如，序列Z={B，C，D，B}是序列X={A，B，C，B，D，A，B}的子序列，相应的递增下标序列为{2，3，5，7}。给定2个序列X和Y，当另一序列Z既是X的子序列又是Y的子序列时，称Z是序列X和Y的公共子序列。给定2个序列X={x1,x2,…,xm}和Y={y1,y2,…,yn}，找出X和Y的最长公共子序列。 三、实验步骤 1、代码 // ZXL_1.cpp : Defines the entry point for the console application. // #include "stdafx.h" #include "iostream.h" #include //为每一行分配空间 for( int i= 0 ; i 实验二动态规划算法的应用 一、实验目的 1．掌握动态规划算法的基本思想，包括最优子结构性质和基于表格的最优值计算方法。 2．熟练掌握分阶段的和递推的最优子结构分析方法。 3．学会利用动态规划算法解决实际问题。 二、实验内容 1.问题描述： 题目一：数塔问题 给定一个数塔，其存储形式为如下所示的下三角矩阵。在此数塔中，从顶部出发，在每一节点可以选择向下走还是向右走，一直走到底层。请找出一条路径，使路径上的数值和最大。 输入样例（数塔）： 9 12 15 10 6 8 2 18 9 5 19 7 10 4 16 输出样例（最大路径和）： 59 三、算法设计 void main() { 申明一个5*5的二维数组； for(int i=0;i<5;i++) { for(int j=0;j<=i;j++) { 输入数组元素p[i][j]; } } for(int k=0;k<5;k++) { for(int w=0;w<=k;w++) { 输出数组元素p[k][w]; } } for(int a=4;a>0;a--) { for(int s=0;s<=a;s++) { if(p[a][s]大于p[a][s+1]) p[a-1][s]等于p[a-1][s]加p[a][s]; else p[a-1][s] 等于p[a-1][s] 加p[a][s+1]; } } 输出p[0][0] } 四.程序调试及运行结果分析 五.实验总结 虽然这个实验比较简单，但是通过这次实验使我更加了解的动态规划法的好处和、，在解决问题时要尝试使用动态规划，这样就有可能得到一种即简单复杂性有不高的算法。 2设计动态规划算法的主要步骤为： （1）找出最优解的性质，并刻划其结构特征。（2）递归地定义最优值。（3）以自底向上的方式计算出最优值。（4）根据计算最优值时得到的信息，构造最优解。 3. 分治法与动态规划法的相同点是：将待求解的问题分解成若干个子问题，先求解子问题，然后从这些子问题的解得到原问题的解。 两者的不同点是：适合于用动态规划法求解的问题，经分解得到的子问题往往不是互相独立的。而用分治法求解的问题，经分解得到的子问题往往是互相独立的。 贪心选择算法与动态规划算法的异同点：同：都要求问题具有最优子结构性质；异：动态规划算法为自底向上的方式解各子问题，贪心算法为自顶向下的方式进行，以迭代的方式作出相继的贪心选择，每做一次贪心选择问题就转换为规模更小的字问题。 6. 分治法所能解决的问题一般具有的几个特征是：（1）该问题的规模缩小到一定的程度就可以容易地解决； （2）该问题可以分解为若干个规模较小的相同问题，即该问题具有最优子结构性质; （3）利用该问题分解出的子问题的解可以合并为该问题的解； （4）原问题所分解出的各个子问题是相互独立的，即子问题之间不包含公共的子问题。 P：也即是多项式复杂程度的问题。 NP就是多项式复杂程度的非确定性问题。 NPC(NP Complete)问题 ADT 抽象数据类型 分析问题→设计算法→编写程序→上机运行和测试 算法特性1. 确定性、可实现性、输入、输出、有穷性 算法分析目的2. 分析算法占用计算机资源的 情况，对算法做出比较和评价，设计出额更好 的算法。 3. 算法的时间复杂性与问题的规模相关，是 问题大小n的函数。 算法的渐进时间复杂性的含义：当问题的规模 n趋向无穷大时，影响算法效率的重要因素是 T(n)的数量级，而其他因素仅是使时间复杂度 相差常数倍，因此可以用T(n)的数量级(阶) 评价算法。时间复杂度T(n)的数量级(阶)称为 渐进时间复杂性。 最坏情况下的时间复杂性和平均时间复杂性有什么不同？ 最坏情况下的时间复杂性和平均时间复杂性 考察的是n固定时，不同输入实例下的算法所 耗时间。最坏情况下的时间复杂性取的输入实 例中最大的时间复杂度： W(n) = max{ T(n，I) } , I∈Dn 平均时间复杂性是所有输入实例的处理时间 与各自概率的乘积和： A(n) =∑P(I)T(n，I) I∈Dn 为什么要分析最坏情况下的算法时间复杂 性？最坏情况下的时间复杂性决定算法的优 劣，并且最坏情况下的时间复杂性较平均时间 复杂性游可操作性。 1.贪心算法的基本思想？ 是一种依据最优化量度依次选择输入的分级处理方法。基本思路是：首先根据题意，选取一种量度标准；然后按这种量度标准对这n个输入排序，依次选择输入量加入部分解中。如果当前这个输入量的加入，不满足约束条件，则不把此输入加到这部分解中。 贪心选择算法与动态规划算法的异同点：同：都要求问题具有最优子结构性质；异：动态规划算法为自底向上的方式解各子问题，贪心算法为自顶向下的方式进行，以迭代的方式作出相继的贪心选择，每做一次贪心选择问题就转换为规模更小的字问题。 Farmer John's farm consists of a long row of N (1 <= N <= 100,000)fields. Each field contains a certain number of cows, 1 <= ncows <= 2000. FJ wants to build a fence around a contiguous group of these fields in order to maximize the average number of cows per field within that block. The block must contain at least F (1 <= F <= N) fields, where F given as input. Calculate the fence placement that maximizes the average, given the constraint. Input * Line 1: Two space-separated integers, N and F. * Lines 2..N+1: Each line contains a single integer, the number of cows in a field. Line 2 gives the number of cows in field 1,line 3 gives the number in field 2, and so on. Output * Line 1: A single integer that is 1000 times the maximal average.Do not perform rounding, just print the integer that is 1000*ncows/nfields. Sample Input 10 6 6 4 2 10 3 8 5 9 4 1 Sample Output 6500 实验三：动态规划法 【实验目的】 深入理解动态规划算法的算法思想，应用动态规划算法解决实际的算法问题。 【实验性质】 验证性实验。 【实验要求】 对于下列所描述的问题，给出相应的算法描述，并完成程序实现与时间复杂度的分析。该问题描述为：一般地，考虑矩阵A1，A2，…，An的连乘积，它们的维数分别为d0,d1,…,dn,即Ai的维数为di-1×di (1≤i≤n)。确定这n个矩阵的乘积结合次序，使所需的总乘法次数最少。对应于乘法次数最少的乘积结合次序为这n个矩阵的最优连乘积次序。按给定的一组测试数据对根据算法设计的程序进行调试：6个矩阵连乘积A=A1×A2×A3×A4×A5×A6，各矩阵的维数分别为：A1：10×20，A2：20×25，A3：25×15，A4：15×5，A5：5×10，A6：10×25。完成测试。 【算法思想及处理过程】 【程序代码】 printf ("\n\n矩阵连乘次数的最优值为:\n"); printf ("-----------------------------------------------\n"); print2 (0, 6-1, s); printf ("\n-----------------------------------------------\n\n"); return 0; } void MatrixChain (int p[], int m[][6], int s[][6], int n) { int i, j, k, z, t; for (i=0; i 题目：用动态规划法实现求两序列的最长公共子序列。 程序代码 #include 实验报告 （2009/2010学年第一学期） 课程名称算法分析与设计A 实验名称动态规划法 实验时间2009 年11 月20 日指导单位计算机学院软件工程系 指导教师张怡婷 学生姓名丁力琪班级学号B07030907 学院(系) 计算机学院专业软件工程 实验报告 实验名称动态规划法指导教师张怡婷实验类型验证实验学时2×2实验时间2009-11-20一、实验目的和任务 目的：加深对动态规划法的算法原理及实现过程的理解，学习用动态规划法解决实际应用中的最长公共子序列问题。 任务：用动态规划法实现求两序列的最长公共子序列，其比较结果可用于基因比较、文章比较等多个领域。 要求：掌握动态规划法的思想，及动态规划法在实际中的应用；分析最长公共子序列的问题特征，选择算法策略并设计具体算法，编程实现两输入序列的比较，并输出它们的最长公共子序列。 二、实验环境(实验设备) 硬件：计算机 软件：Visual C++ 三、实验原理及内容（包括操作过程、结果分析等） 1、最长公共子序列（LCS）问题是：给定两个字符序列X={x1,x2,……,x m}和Y={y1,y2,……,y n}，要求找出X和Y的一个最长公共子序列。 例如：X={a,b,c,b,d,a,b},Y={b,d,c,a,b,a}。它们的最长公共子序列LSC={b,c,d,a}。 通过“穷举法”列出所有X的所有子序列，检查其是否为Y的子序列并记录最长公共子序列并记录最长公共子序列的长度这种方法，求解时间为指数级别的，因此不可取。 2、分析LCS问题特征可知，如果Z={z1,z2,……，z k}为它们的最长公共子序列，则它们一定具有以下性质： （1）若x m=y n，则z k=x m=y n，且Z k-1是X m-1和Y n-1的最长公共子序列； （2）若x m≠y n且x m≠z k，则Z是X m-1和Y的最长公共子序列； （3）若x m≠y n且z k≠y n，则Z是X和Y的最长公共子序列。 这样就将求X和Y的最长公共子序列问题，分解为求解较小规模的问题： 若x m=y m，则进一步分解为求解两个（前缀）子字符序列X m-1和Y n-1的最长公共子序列问题； 如果x m≠y n，则原问题转化为求解两个子问题，即找出X m-1和Y的最长公共子序列与找出X 和Y n-1的最长公共子序列，取两者中较长者作为X和Y的最长公共子序列。 由此可见，两个序列的最长公共子序列包含了这两个序列的前缀的最长公共子序列，具有最优子结构性质。 3、令c[i][j]保存字符序列X i={x1,x2,……,x i}和Y j={y1,y2,……,y j}的最长公共子序列的长度，由上述分析可得如下递推式： 0 i=0或j=0 c[i][j]= c[i-1][j-1]+1 i,j>0且x i=y j max{c[i][j-1],c[i-1][j]} i,j>0且x i≠y j 由此可见，最长公共子序列的求解具有重叠子问题性质，如果采用递归算法实现，会得到一个指数时间算法，因此需要采用动态规划法自底向上求解，并保存子问题的解，这样可以避免重复计算子问题，在多项式时间内完成计算。 4、为了能由最优解值进一步得到最优解（即最长公共子序列），还需要一个二维数组s[][]，数组中的元素s[i][j]记录c[i][j]的值是由三个子问题c[i-1][j-1]+1,c[i][j-1]和c[i-1][j]中的哪一个计算得到，从而可以得到最优解的当前解分量（即最长公共子序列中的当前字符），最终构造出最长公共子序列自身。 本科实验报告 课程名称：算法设计与分析 实验项目：递归与分治算法 实验地点：计算机系实验楼110 专业班级：物联网1601 学号：2016002105 学生姓名：俞梦真 指导教师：郝晓丽 2018年05月04 日 实验一递归与分治算法 1.1 实验目的与要求 1．进一步熟悉C/C++语言的集成开发环境； 2．通过本实验加深对递归与分治策略的理解和运用。 1.2 实验课时 2学时 1.3 实验原理 分治（Divide-and-Conquer）的思想：一个规模为n的复杂问题的求解，可以划分成若干个规模小于n的子问题，再将子问题的解合并成原问题的解。 需要注意的是，分治法使用递归的思想。划分后的每一个子问题与原问题的性质相同，可用相同的求解方法。最后，当子问题规模足够小时，可以直接求解，然后逆求原问题的解。 1.4 实验题目 1．上机题目：格雷码构造问题 Gray码是一个长度为2n的序列。序列无相同元素，每个元素都是长度为n的串，相邻元素恰好只有一位不同。试设计一个算法对任意n构造相应的Gray码（分治、减治、变治皆可）。 对于给定的正整数n，格雷码为满足如下条件的一个编码序列。 （1）序列由2n个编码组成，每个编码都是长度为n的二进制位串。 （2）序列中无相同的编码。 （3）序列中位置相邻的两个编码恰有一位不同。 2．设计思想： 根据格雷码的性质，找到他的规律，可发现，1位是0 1。两位是00 01 11 10。三位是000 001 011 010 110 111 101 100。n位是前n-1位的2倍个。N-1个位前面加0，N-2为倒转再前面再加1。 3．代码设计： 实验二用MATLAB解决动态规划问题 问题:有一部货车每天沿着公路给四个售货店卸下6箱货物,如果各零售店出售该货物所得利润如下表所示,试求在各零售店卸下几箱货物,能使获得总利润最 解: 1)将问题按售货店分为四个阶段 2)设s k表示为分配给第k个售货店到第n个工厂的货物数, x k设为决策变量,表示为分配给第k个售货店的货物数, 状态转移方程为s k＋1＝s k－x k。 P k(x k)表示为x k箱货物分到第k个售货店所得的盈利值。 f k(s k)表示为s k箱货物分配给第k个售货店到第n个售货店的最大盈利值。 3)递推关系式: f k(s k)＝max[ P k(x k)+ f k+1(s k－x k) ] k=4,3,2,1 边界条件:f5(s5)＝0 4)从最后一个阶段开始向前逆推计算。 第四阶段: 设将s4箱货物(s4＝0,1,2,3,4,5,6)全部分配给4售货店时,最大盈利值为: f4(s4)＝max[P4(x4)] 其中x4＝s4＝0,1,2,3,4,5,6 x4*表示使得f4(s4)为最大值时的最优决策。 第三阶段: 设将s3箱货物(s3＝0,1,2,3,4,5,6)分配给3售货店与4售货店时,对每一个s3值,都有一种最优分配方案,使得最大盈利值为:f3(s3)＝max[ P3(x3)+ f4(s3－x3) ] ,x3＝ 第二阶段: 设将s2箱货物(s2＝0,1,2,3,4,5,6)分配给2售货店、3售货店与4售货店时,则最大盈利值为:f2(s2)＝max[ P2(x2)+ f3(s2－x2) ] 第一阶段: 设将s2箱货物(s1＝0,1,2,3,4,5,6)分配给1售货店、2售货店、3售货店与4售货店时,则最大盈利值为:f1(s1)＝max[ P1(x1)+ f2(s1－x1) ] 按计算表格的顺序反推,可知最优分配方案有6个: 1） x1*＝1,x2*＝1,x3*=3,x4*=1。 2） x1*＝1,x2*＝2,x3*=2,x4*=1。 3） x1*＝1,x2*＝3,x3*=1,x4*=1。 《算法设计与分析》实验报告实验二递归与分治策略 Module 1: 免费馅饼 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 59327 Accepted Submission(s): 20813 Problem Description 都说天上不会掉馅饼，但有一天gameboy正走在回家的小径上，忽然天上掉下大把大把的馅饼。说来gameboy的人品实在是太好了，这馅饼别处都不掉，就掉落在他身旁的10米范围内。馅饼如果掉在了地上当然就不能吃了，所以gameboy马上卸下身上的背包去接。但由于小径两侧都不能站人，所以他只能在小径上接。由于gameboy平时老呆在房间里玩游戏，虽然在游戏中是个身手敏捷的高手，但在现实中运动神经特别迟钝，每秒种只有在移动不超过一米的范围内接住坠落的馅饼。现在给这条小径如图标上坐标： 为了使问题简化，假设在接下来的一段时间里，馅饼都掉落在0-10这11个位置。开始时gameboy站在5这个位置，因此在第一秒，他只能接到4,5,6这三个位置中其中一个位置上的馅饼。问gameboy最多可能接到多少个馅饼？（假设他的背包可以容纳无穷多个馅饼） Input 输入数据有多组。每组数据的第一行为以正整数n(0 实验3动态规划上机 [实验目的] 1．掌握动态规划的基本思想和效率分析方法； 2．掌握使用动态规划算法的基本步骤； 3．学会利用动态规划解决实际问题。 [实验要求] 按以下实验内容完成题目，并把编译、运行过程中出现的问题以及解决方法填入实验报告中，按时上交。 [实验学时] 2学时。 [实验内容] 一、实验内容 利用动态规划算法编程求解多段图问题，要求读入多段图，考虑多段图的排序方式，求源点到汇点的最小成本路径。并请对自己的程序进行复杂性分析。 二、算法描述 先输入点的个数和路径数以及每段路径的起点、长度、终点，再计算所有路径的值大小，比较输出后最小值。 三、源程序 #define N 2147483647 #include } array2[pointnum]=N; array3[pointnum]=N; for(i=0;i XXXX大学计算机学院实验报告计算机学院2017级软件工程专业 5 班指导教师 学号姓名2019年10 月21 日成绩 实验内容、上机调试程序、程序运行结果 System.out.println("选中的物品是第"); for(int i=1;i<=n;i++){ for(int j=1;j<=maxweight;j++){ //当前最大价值等于放前一件的最大价值 maxvalue[i][j] = maxvalue[i-1][j]; //如果当前物品的重量小于总重量，可以放进去或者拿出别的东西再放进去 if(weight[i-1] <= j){ //比较（不放这个物品的价值）和（这个物品的价值放进去加上当前能放的总重量减去当前物品重量时取i-1个物品是的对应重量时候的最高价值） if(maxvalue[i-1][j-weight[i-1]] + value[i - 1] > maxvalue[i-1][j]){ maxvalue[i][j] = maxvalue[i-1][j-weight[i-1]] + value[i - 1]; } } } } return maxvalue[n][maxweight]; } public static void main(String[] args) { int weight[] = {2,3,4,5}; int value[] = {3,4,5,7}; int maxweight = 8; System.out.println(knapsack(weight,value,maxweight)); } } 完成效果： 实验二动态规划算法（2学时） 一、实验目的与要求 1、熟悉最长公共子序列问题的算法； 2、初步掌握动态规划算法； 二、实验题 若给定序列X={x1,x2,…,xm}，则另一序列Z={z1,z2,…,zk}，是X的子序列是指存在一个严格递增下标序列{i1,i2,…,ik}使得对于所有j=1,2,…,k有：zj=xij。例如，序列Z={B，C，D，B}是序列X={A，B，C，B，D，A，B}的子序列，相应的递增下标序列为{2，3，5，7}。给定2个序列X和Y，当另一序列Z既是X的子序列又是Y的子序列时，称Z是序列X 和Y的公共子序列。 给定2个序列X={x1,x2,…,xm}和Y={y1,y2,…,yn}，找出X和Y的最长公共子序列。 三、实验提示 include "stdlib.h" #include "string.h" void LCSLength(char *x ,char *y,int m,int n, int **c, int **b) { int i ,j; for (i = 1; i <= m; i++) c[i][0] = 0; for (i = 1; i <= n; i++) c[0][i] = 0; for (i = 1; i <= m; i++) for (j = 1; j <= n; j++) { if (x[i]==y[j]) { c[i][j]=c[i-1][j-1]+1; b[i][j]=1; } else if (c[i-1][j]>=c[i][j-1]) { c[i][j]=c[i-1][j]; b[i][j]=2; } else { c[i][j]=c[i][j-1]; b[i][j]=3; } } } void LCS(int i ,int j, char *x ,int **b) { if (i ==0 || j==0) return; §1 二次规划模型 数学模型： ub x lb beq x Aeq b x A x f Hx x T T x ≤≤=?≤?+21min 其中H 为二次型矩阵，A 、Aeq 分别为不等式约束与等式约束系数矩阵，f,b,beq,lb,ub,x 为向量。 求解二次规划问题函数为quadprog( ) 调用格式： X= quadprog(H,f,A,b) X= quadprog(H,f,A,b,Aeq,beq) X= quadprog(H,f,A,b,Aeq,beq,lb,ub) X= quadprog(H,f,A,b,Aeq,beq,lb,ub,x0) X= quadprog(H,f,A,b,Aeq,beq,lb,ub,x0,options) [x,fval]= quadprog(…) [x,fval,exitflag]= quadprog(…) [x,fval,exitflag,output]= quadprog(…) [x,fval,exitflag,output,lambda]= quadprog(…) 说明：输入参数中，x0为初始点；若无等式约束或无不等式约束，就将相应的矩阵和向量设置为空；options 为指定优化参数。输出参数中，x 是返回最优解；fval 是返回解所对应的目标函数值；exitflag 是描述搜索是否收敛；output 是返回包含优化信息的结构。Lambda 是返回解x 入包含拉格朗日乘子的参数。 例1：求解：二次规划问题 min f(x)= x 1-3x 2+3x 12+4x 22 -2x 1x 2 s.t 2x 1+x 2≤2 -x 1+4x 2≤3 动态规划算法实验报告 ————————————————————————————————作者: ————————————————————————————————日期: 实验标题 1、矩阵连乘 2、最长公共子序列３、最大子段和 4、凸多边形最优三角剖分 5、流水作业调度 6、0－１背包问题 7、最优二叉搜索树 实验目的掌握动态规划法的基本思想和算法设计的基本步骤。 实验内容与源码1、矩阵连乘 #inclｕde 算法设计与分析基础 ——实验报告 姓名：周建权 学号：0909122820 班级：信安1202 实验一分治 —最近点对 一．问题 Problem Have you ever played quoit in a playground? Quoit is a game in which flat rings are pitched at some toys, with all the toys encircled awarded. In the field of Cyberground, the position of each toy is fixed, and the ring is carefully designed so it can only encircle one toy at a time. On the other hand, to make the game look more attractive, the ring is designed to have the largest radius. Given a configuration of the field, you are supposed to find the radius of such a ring. Assume that all the toys are points on a plane. A point is encircled by the ring if the distance between the point and the center of the ring is strictly less than the radius of the ring. If two toys are placed at the same point, the radius of the ring is considered to be 0. Input The input consists of several test cases. For each case, the first line contains an integer N (2 <= N <= 100,000), the total number of toys in the field. Then N lines follow, each contains a pair of (x, y) which are the coordinates of a toy. The input is terminated by N = 0. Output For each test case, print in one line the radius of the ring required by the Cyberground manager, accurate up to 2 decimal places. 二．分析思路 题目是给n个点的坐标，求距离最近的一对点之间距离的一半。第一行是一个数n表示有n个点，接下来n行是n个点的x坐标和y坐标。 首先，假设点是n个，编号为1到n。找一个中间的编号mid，先求出1到mid点的最近距离设为d1，还有mid+1到n的最近距离设为d2。如果说最近点对中的两点都在1-mid 集合中，或者mid+1到n集合中，则d就是最小距离了。但是还有可能的是最近点对中的两点分属这两个集合，若存在，则把这个最近点对的距离记录下来，去更新d。这样就得到最小的距离d了。 三．源代码 #include 动态规划求解资源分配 实验目标： （1）掌握用动态规划方法求解实际问题的基本思路。 （2）进一步理解动态规划方法的实质，巩固设计动态规划算法的基本步骤。 实验任务： （1）设计动态规划算法求解资源分配问题，给出算法的非形式描述。 （2）在Windows环境下用C语言实现该算法。计算10个实例，每个实例中n=30，m=10，C i j为随机产生于范围（0，103）内的整数。记录各实例的数据及执行结果（即最优分配方案、最优分配方案的值）、运行时间。 （3）从理论上分析算法的时间和空间复杂度，并由此解释相应的实验结果。 实验设备及环境： PC；C/C++等编程语言。 实验主要步骤： (1)认真阅读实验目的与实验任务，明确本次实验的内容； (2)分析实验中要求求解的问题，根据动态规划的思想，得出优化方程； (3)从问题出发，设计出相应的动态规划算法，并根据设计编写程序实现算法； (4)设计实验数据并运行程序、记录运行的结果； (5)分析算法的时间和空间复杂度，并由此解释释相应的实验结果； 问题描述：资源分配问题 某厂根据计划安排，拟将n台相同的设备分配给m个车间，各车间获得这种设备后，可以为国家提供盈利C i j(i台设备提供给j号车间将得到的利润，1≤i≤n，1≤j≤m) 。问如何分配，才使国家得到最大的盈利？ 1.问题分析： 本问题是一简单资源分配问题，由于具有明显的最优子结构，故可以使用动态规划求解，用状态量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)。实验二 动态规划算法棋盘覆盖
动态规划算法的应用实验报告
2设计动态规划算法的主要步骤为
动态规划(2)
算法实验动态规划----矩阵连乘
算法分析与设计实验二：动态规划法
南京邮电大学算法设计实验报告——动态规划法
算法设计与分析实验报告
运筹学实验_动态规划
算法设计与分析---动态规划实验
算法实验 动态规划上机
实验报告：动态规划---0-1背包问题)
实验二 动态规划算法—李明明
二次_动态规划-图论
动态规划算法实验报告
中南大学算法实验报告
动态规划 求解资源分配 实验报告