动态规划学习笔记 初识动态规划

写在前面

注意:此文章仅供参考,如发现有误请及时告知。

更新日期:2018/3/16,2018/12/03


动态规划介绍

动态规划,简称DP(Dynamic Programming)
简介1 简介2

动态规划十分奇妙,它可以变身为记忆化搜索,变身为递推,甚至有时可以简化成一个小小的算式。

动态规划十分灵活,例如 NOIP2018 PJ T3 摆渡车 ,写法有很多很多,但时间、内存却各有差异。

动态规划十分简单,有时候一个小小的转移方程就能解决问题。

动态规划十分深奥,有时你会死也想不出合适的转移方程,有时你会被后效性困扰,有时动态规划的同时还有许多蜜汁优化。

动态规划在NOIP中十分重要,我目前为止参加的\(NOIP_{2017 PJ} \& NOIP_{2018PJ}\)都有一道动态规划,而且都是\(T3\)。(估计普及考纲比较窄,要出难题只有DP了)


问题引入

还是这道题...... 数塔问题!!!

img

这里我们选择动态规划来解决.
我们不难理解,对于每一个元素,它到顶层的最大值是确定的,也就是说,从顶层到任何一个元素的最大值都是确定的.比如,对于第3层的第2个元素6,顶层到它的最大值只有一个(9 + 15 + 6 = 30)(但不代表路径只有一条),不会改变.

所以,我们用一个数组dp来存储从元素(i, j)到底层的最大值.

#define MAXN 100
int dp[MAXN + 5][MAXN + 5];

仔细观察分析,不难发现,对于每一个元素dp[i][j],都存在

dp[i][j] = max( dp[i + 1][j], dp[i + 1][j + 1] ) + a[i][j];

即每一个元素到(1, 1)的最大值都是上一层与它相连的两个元素中较大的一个,再加上这个元素本身的值. 最后的答案即为dp[1][1].

不过,我们自顶向下分析,但是却要自底向上实现,即从最顶层开始分析,写代码时却要注意for语句要倒过来写:

for ( int i = N; i >= 1; --i )
    for ( int j = 1; j <= i; ++j )
        dp[i][j] = max( dp[i + 1][j], dp[i + 1][j + 1] ) + a[i][j];

为什么会这样呢?其实不难分析,在算dp[i][j]时,你必须确保dp[i + 1][j]dp[i + 1][j + 1]已经完成,如果没有完成,dp[i + 1][j]dp[i + 1][j + 1]的值就是错误的,算出的dp[i][j]也是错误的,这样结果就不对了。而反过来做,你就会发现i从大的开始,在做dp[i][j]的时候dp[i + 1][1 ~ N]都已经做过了。还有,要注意,动态规划的初始化很重要,有时初始化就会决定你结果对不对。这里的初始化很简单,现在给出两种方法:

memset( dp[N + 1], 0, sizeof( dp[N + 1] ) );//即把dp[N + 1][0...]全部初始化为0.
for ( int i = 1; i <= N; ++i )
    dp[i] = a[i];
//下面这个与上面等价:
copy( a[N] + 1, a[N] + N + 1, dp[N] );// copy( 开始地址, 结束地址, 放到的数组 ); copy( a, a + n, b );即为把a数组下标为0~n按次序复制到b数组.
//当然,这样写,实现时要注意少一层循环:(下面这个是修改后的)
for ( int i = N - 1; i >= 1; --i )
    for ( int j = 1; j <= i; ++j )
        dp[i][j] = max( dp[i + 1][j], dp[i + 1][j + 1] ) + a[i][j];
//至于为什么这样,这里不再赘述,请自己思考. 

这里再完整地放一放代码,实在不会写的可以参考.

#include<bits/stdc++.h>
using namespace std;
#define MAXN 100

int C, N;
int a[MAXN + 5][MAXN + 5];
int dp[MAXN + 5][MAXN + 5];

void solve(){
    scanf( "%d", &N );
    memset( dp, 0, sizeof dp ); 
    for ( int i  = 1; i <= N; ++i )
        for ( int j = 1; j <= i; ++j )
            scanf( "%d", &a[i][j] );
    for ( int i = N; i >= 1; --i )
        for ( int j = 1; j <= i; ++j )
            dp[i][j] = max( dp[i + 1][j], dp[i + 1][j + 1] ) + a[i][j];
    printf( "%d\n", dp[1][1] );
}

int main(){
    scanf( "%d", &C );
    while( C-- ) solve();
    return 0;
}

事实上,可以做一个优化:去掉dp数组,直接用a数组来做:(节约空间,人人有责)

#include<bits/stdc++.h>
using namespace std;
#define MAXN 100

int C, N;
int a[MAXN + 5][MAXN + 5];

void solve(){
    scanf( "%d", &N ); 
    for ( int i  = 1; i <= N; ++i )
        for ( int j = 1; j <= i; ++j )
            scanf( "%d", &a[i][j] );
    for ( int i = N - 1; i >= 1; --i )
        for ( int j = 1; j <= i; ++j )
            a[i][j] += max( a[i + 1][j], a[i + 1][j + 1] );
    printf( "%d\n", a[1][1] );
}

int main(){
    scanf( "%d", &C );
    while( C-- ) solve();
    return 0;
}

至于为什么,请诸位自己理解(很好理解的,选个小一点的数据自己算一算就知道了)。


总结

怎么样,找到些感觉了吧?现在我们来学习怎么写动态规划的程序.

第一步,我们要观察题目是否可以用动态规划实现。怎么判断呢?我们要看它是否可以分成几个阶段,如上题,可以分成1~N层共N个阶段,每个阶段还可以分成1~i个元素共i个小阶段。然后,我们要看看每个阶段的答案是不是确定的,上题中,每一个元素到底层的最大值就是确定的。再看看每个阶段是不是有关联,如果有,还要确定有什么关联,是否对于每一个阶段都满足。

第二步,就是确定关联啦。怎么确定呢?我们要仔细分析题目,观察每两个阶段之间的关系。动态规划的重点也就在这里,关联确定了,动态规划基本上就可以写下来了。

第三步,确定边界条件,比如,上题就要把dp[N+1][...]全部赋值为0,否则就会出错。

除此之外,还要确定完成的顺序,要做某个阶段,它需要用到的阶段必须先做完。

当然,有时还要添加滚动数组、优化等。

这样,一个动态规划程序就完成啦。


尾声

当然,动态规划还有许多分支(背包DP、区间DP等),以上讲的都是最表皮的。那些难一点的,都只好下次再讲吧。

最好拿点题目来练一下:洛谷的DP

猜你喜欢

转载自www.cnblogs.com/louhancheng/p/10061320.html