数据结构(十二)动态规划

1 递归函数建模

动态规划一般用于全局问题,在构造递归的时候,一般采用自顶向下分解的方法,先把全局问题分解成更小的子问题求解。下面举两个例子

例子1有一座高度是10阶的楼梯,从下往上走,每跨一步可以是一级或两级台阶。要求用程序求出一共一共有多少种走法。

问题分析建模:首先总共有10步,假设只剩最后一步就到达第10阶,这个时候会有两种情况:第一种是从第九阶到第十阶,第二种是从第八阶到第十阶,然后两种情况里面如何从第一个阶梯走到第九阶,从第一个阶梯走到第八阶。也就是我们可以构造出最后一步的函数:

f(10)=f(9)+f(8)

这个时候我们就可以构造出递归函数:

f(n)=f(n-1)+f(n-2)

写成通用一点的公式,也就是用for循环遍历最后一次的情况:

f(n)=sum(f(n-j)) j=1,2

因此我们只要知道初始值f(1)f(2)就可以求解f(10)。这种我们称之为一维动态规划。

例子2:给定一个m*n矩阵A,从左上角开始每次只能向右走或者向下走,最后达到右下角的位置,路径中所有数字累加起来就是路径和,返回所有路径的最小路径和,如果给定的矩阵如下,那么路径1,3,1,0,6,1,0就是最小路径和,返回12

1 3 5 9

8 1 3 4

5 0 6 1

8 8 4 0

问题分析建模:对于这个题目,假设mmn列的矩阵,那么我们用dist[m][n]来抽象这个问题,dist[i][j]表示的是从原点到i,j位置的最短路径和。我们要求解的是dist[m][n]的最大值,然而其实最后一步dist[m][n]又可以分解成两种情况,从dist[m-1][n]往下走;从dist[m][n-1]往右走,取这两种情况的最小值,那么我们就可以列出方程:

dist[m][n]=min(dist[m-1][n]+A[m][n],dist[m][n-1]+A[m][n])

抽取出A[m][n],该方程可以进一步写成:

dist[m][n]=A[m][n]+min(dist[m-1][n],dist[m][n-1])

因此接着我们只需要根据初始条件,就可以求解,初始条件dist[0][0]=1这个例子又称之为二维动态规划。

2动态规划

(1)自顶向下递归

其主要原理是加入备忘录,用于记录之前已经记录过的函数值,比如在写例子1递归函数:

f(n)=f(n-1)+f(n-2)

代码的时候,用一个数组记录一下f(i),然后在写递归函数的时候,遇到f(i)的值如果已经被计算过一次了,那么就不用继续计算了;如果还没有被记录,那么就保存一下f(i)

在写例子2递归函数的时候,用二维的数组f[m][n]记录一下已经计算过的函数值

(2)自底向上迭代

0开始往上迭代,比如对于例子1,先求:

f[2]=f[1]+f[0]

f[3]=f[2]+f[1]

f[4]=f[3]+f[3]

 直到迭代到f[10]为止。

示例源码

#include <istream>;
#include <vector>;
#include <math.h>
using namespace std;
int A[4][4] = {1,3,5,9,8,1,3,4,5,0,6,1,8,8,4,0};
int f[4][4]={0};//备忘录
int count=0;//统一一下采用备忘录的时候,减少的计算次数
int function_top2bottom(int m,int n){//自顶向下递归
    if (f[m][n]>0)
        return f[m][n];
    count++;

    int min_dist=0;
    if(m==0&&n==0)
    {
        min_dist=0;
    }
    else if (m-1>=0&&n==0)
        min_dist=function_top2bottom(m-1,n);
    else if (m==0&&n-1>=0) {
        min_dist=function_top2bottom(m,n-1);
    }
    else {
        min_dist=min(function_top2bottom(m,n-1),function_top2bottom(m-1,n));
    }
    f[m][n]=min_dist+A[m][n];

    return f[m][n];
}

int function_bottom2top(int m,int n){//自底向下迭代
    for (int i = 0; i < m; ++i) {
        for (int j = 0; j <n ; ++j) {
            int min_dist=0;
            if(m==0&&n==0)
            {
                min_dist=0;
            }
            else if (m-1>=0&&n==0)
                min_dist=f[m-1][n];
            else if (m==0&&n-1>=0) {
                min_dist=f[m][n-1];
            }
            else {
                min_dist=min(f[m-1][n],f[m][n-1]);
            }
            f[m][n]=min_dist+A[m][n];
        }

    }


    return f[m][n];
}

void dynamic_programming()
{

    std::cout<<function_top2bottom(3,3)<<std::endl;
    std::cout<<count<<std::endl;//可以测试一下,自顶向下递归的时候,少用备忘录的count值

    std::cout<<function_bottom2top(3,3)<<std::endl;

}
    原文作者:动态规划
    原文地址: https://blog.csdn.net/hjimce/article/details/79311314
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞