算法导论学习-动态规划之记忆化搜索

一. 动态规划

动态规划(dynamic programming),与“分治思想”有些相似,都是利用将问题分 为子问题,并通过合并子问题的解来获得整个问题的解。于“分治”的不同之处在 于,对于一个相同的子问题动态规划算法不会计算第二次,其实现原理是将每一个计算过的子问题的值保存在一个表中。

二. 记忆化搜索

我们常见的动态规划问题,比如流水线调度问题,矩阵链乘问题等等都是“一步接着一步解决的”,即规模为 i 的问题需要基于规模 i-1 的问题进行最优解选择,通常的递归模式为DP(i)=optimal{DP(i-1)}。而记忆化搜索本质上也是DP思想,当子问题A和子问题B存在子子问题C时,如果子子问题C的最优解已经被求出,那么子问题A或者是B只需要“查表”获得C的解,而不需要再算一遍C。记忆化搜索的DP模式比普通模式要“随意一些”,通常为DP(i)=optimal(DP(j)), j < i。

三. 滑雪问题

上图显示为R*C的雪场,R是行数,C是列数。圆圈内的数字表示的是雪场的海拔高度h,根据常识我们知道,滑雪只有从上往下滑行才可能滑的动,现在我们想要求出能够滑行的最长距离,上面的例子我们可以很直接判断出25-24-......-1这个顺时针方向螺旋的滑雪方式可以滑的最远。

那么这个问题如何用编程来实现呢?我们发现这是一个典型的递推,DP(i, j)表示从坐标(i,j)出发所能滑行的最大长度,且有:DP(i, j)=optimal{DP(i±1, j±1)}+1。下面貼上源代碼。

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 using namespace std;
 6 const int max_size=110;
 7 int R,C;
 8 int dir[4][2]={{-1,0},{0,1},{1,0},{0,-1}};
 9 int h[max_size][max_size],dp[max_size][max_size];
10 int inMap(int x,int y){
11     if(x>=0&&x<=R-1&&y>=0&&y<=C-1) return 1;
12     return 0;
13 }
14 int max2(int a,int b,int c,int d){
15     return max(max(a,b),max(c,d));
16 }
17 int dfs(int i,int j){
18     int nx,ny,down=0,up=0,left=0,right=0;
19     if(dp[i][j]) return dp[i][j];
20     nx=i+dir[0][0]; ny=j+dir[0][1];
21     if(inMap(nx,ny)){
22         if(h[i][j]>h[nx][ny]) up=dfs(nx,ny);
23     }
24     nx=i+dir[1][0]; ny=j+dir[1][1];
25     if(inMap(nx,ny)){
26         if(h[i][j]>h[nx][ny]) right=dfs(nx,ny);
27     }
28     nx=i+dir[2][0]; ny=j+dir[2][1];
29     if(inMap(nx,ny)){
30         if(h[i][j]>h[nx][ny]) down=dfs(nx,ny);
31     }
32     nx=i+dir[3][0]; ny=j+dir[3][1];
33     if(inMap(nx,ny)){
34         if(h[i][j]>h[nx][ny]) left=dfs(nx,ny);
35     }
36     dp[i][j]=max2(up,down,left,right)+1;
37     return dp[i][j];
38 }
39 int main(){
40     scanf("%d%d",&R,&C);
41     memset(h,0,sizeof(h));
42     memset(dp,0,sizeof(dp));
43     for(int i=0;i<R;i++){
44         for(int j=0;j<C;j++){
45             scanf("%d",&h[i][j]);
46         }
47     }
48     int ans=-1;
49     for(int i=0;i<R;i++){
50         for(int j=0;j<C;j++){
51             ans=max(ans,dfs(i,j));
52         }
53     }
54     printf("%d\n",ans);
55 }

四. 切棒子问题

给你一根长n英尺的棒子和一份关于该棒子的价目表如下(其中 i = 1,2,3,…,n),请问如何将这根棒子卖出最高的价格,可以对棒子进行切割。

这个题同样是可以利用DP记忆化搜索来实现的,递推公式为DP(n)=optimal{max{price(i)+DP(n-i)|1≤i≤n}}。实现代码如下:

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 using namespace std;
 5 const int max_size=50;
 6 const int inf=1<<30;
 7 int price[max_size],dp[max_size];
 8 int n;
 9 int dfs(int n){
10     if(dp[n]) return dp[n];
11     if(n==0) return 0;
12     int mmax=-inf;
13     for(int i=1;i<=n;i++){
14         mmax=max(mmax,price[i]+dfs(n-i));
15     }
16     dp[n]=mmax;
17     return dp[n];
18 }
19 int main(){
20     while(scanf("%d",&n)!=EOF){
21         for(int i=1;i<=n;i++) scanf("%d",&price[i]);
22         printf("%d\n",dfs(n));
23     }
24 }

五. 总结

通过前两个例子分析,我们可以得出DP记忆化搜索的算法模板(自己DIY的,大家可以选择参考)

1 dfs(problem a){
2     if(a has been solved)
3         then: consult the record.
4     else//get the optimal solution of problem a.
5         divide the problem a into several sub-problems(a1,a2,...,ak)
6         get the solution of problem a by dfs(a1),dfs(a2),...,dfs(ak).
7     finally write the optimal solution into record.
8 }
时间: 2024-10-13 16:23:04

算法导论学习-动态规划之记忆化搜索的相关文章

动态规划、记忆化搜索、Dijkstra算法的总结

动态规划 动态规划算法通常用于求解具有某种最优性质的问题.在这类问题中,可能会有许多可行解.每一个解都对应于一个值,我们希望找到具有最优值的解.动态规划算法与分治法类似,其基本思想也是将待求解问题分解成若干个子问题,先求解子问题,然后从这些子问题的解得到原问题的解.与分治法不同的是,适合于用动态规划求解的问题,经分解得到子问题往往不是互相独立的.若用分治法来解这类问题,则分解得到的子问题数目太多,有些子问题被重复计算了很多次.如果我们能够保存已解决的子问题的答案,而在需要时再找出已求得的答案,这

hdu 1978 How many ways (动态规划、记忆化搜索)

How many ways Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 2957    Accepted Submission(s): 1733 Problem Description 这是一个简单的生存游戏,你控制一个机器人从一个棋盘的起始点(1,1)走到棋盘的终点(n,m).游戏的规则描述如下: 1.机器人一开始在棋盘的起始点并

区间型动态规划的记忆化搜索实现与环形动态规划的循环数组实现

区间型动态规划的典型例题是石子归并,同时使用记忆化搜索实现区间动归是一种比较容易实现的方式,避免了循环数组实现的时候一些边界的判断 n堆石子排列成一条线,我们可以将相邻的两堆石子进行合并,合并之后需要消耗的代价为这两堆石子的质量之和,问最小的合并代价 状态转移方程很容易给出: f[i][j]=min(f[i][j],f[i][k]+f[k+1][j]+sum[i][j]) 因为要计算区间和,考虑前缀和进行预处理 然后我们给出用记忆化搜索形式实现的代码,这里的记忆化搜索形式可以作为后续问题的一个模

leetcode_1048. Longest String Chain_[DP,动态规划,记忆化搜索]

https://leetcode.com/problems/longest-string-chain/ Let's say word1 is a predecessor of word2 if and only if we can add exactly one letter anywhere in word1 to make it equal to word2.  For example, "abc" is a predecessor of "abac". A w

FZU 2092 bfs+记忆化搜索

晚上团队训练赛的题 和普通bfs不同的是 这是同时操纵人与影子两个单位进行的bfs 由于可能发生人和影子同时接触水晶 所以不可以分开操作 当时使用node记录人和影子的位置 然后进行两重for循环来分别改变位置 结果超内存 分析了一下应该是队列超了内存 毕竟如果每个点都存入的话一个点最多可以衍生出25个node 然后t最大为200s 一定会超 之间还发生了一些并不能理解的bug 被逼到最后重构才拿到了一个超内存 名次也不好 急需一个ac来赶上去 简直要烧起来了 侧面反映心理素质还是差一些 当时未

【算法学习笔记】55.DFS 记忆化搜索 SJTU OJ 1063 小M爱滑雪

Description 小M超级喜欢滑雪~~ 滑雪的确很刺激.可是为了获得速度,滑的区域必须向下倾斜,而且当小M滑到坡底,便不得不再次走上坡或者等待升降机来载你.小M想知道滑雪场中最长底的滑坡.滑雪场由一个二维数组给出.数组的每个数字代表点距离水平面的相对距离.下面是一个例子 1 2 3 4 5 16 17 18 19 6 15 24 25 20 7 14 23 22 21 8 13 12 11 10 9 小M可以从某个点滑向上下左右相邻四个点之一,当且仅当高度减小.在上面的例子中,一条可滑行的

动态规划 数字三角形(递归,递推,记忆化搜索)

题目要求: 7 3 8 8 1 0 2 7 4 4 4 5 2 6 5 在上面的数字三角形中寻找在上面的数字三角形中寻找一条从顶部到底边的路径,使得路径上所经过的数字之和最大.路径上的每一步都只能往左下或右下走.只需要求出这个最大和即可,不必给出具体路径. 三角形的行数大于1小于等于100,数字为 0 - 99 输入格式: 5 //三角形行数.下面是三角形 7 3 8 8 1 0 2 7 4 4 4 5 2 6 5 解题思路: 用二维数组存放数字三角形 D[r][j] //表示第i行第j个元素的

sicily 1176. Two Ends (Top-down 动态规划+记忆化搜索 v.s. Bottom-up 动态规划)

DescriptionIn the two-player game "Two Ends", an even number of cards is laid out in a row. On each card, face up, is written a positive integer. Players take turns removing a card from either end of the row and placing the card in their pile. T

Codevs_1017_乘积最大_(划分型动态规划/记忆化搜索)

描述 http://codevs.cn/problem/1017/ 给出一个n位数,在数字中间添加k个乘号,使得最终的乘积最大. 1017 乘积最大 2000年NOIP全国联赛普及组NOIP全国联赛提高组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 黄金 Gold 题目描述 Description 今年是国际数学联盟确定的“2000——世界数学年”,又恰逢我国著名数学家华罗庚先生诞辰90周年.在华罗庚先生的家乡江苏金坛,组织了一场别开生面的数学智力竞赛的活动,你的一个好朋友