C++动态规划dp算法题

大数据分析BDA 2019-06-02

问题1:找硬币,换钱的方法

输入:

penny数组代表所有货币的面值,正数不重复
aim小于等于1000,代表要找的钱
输出:
换钱的方法总数

解法1:经典dp,空间复杂度O(n*aim)

class Exchange {
public:
    int countWays(vector<int> penny, int n, int aim) {
        if (penny.empty()||n == 0)
            return 0;
        vector<vector<int> > dp(n,vector<int>(aim+1)); //二维数组dp
        for (int i = 0;i < n;i++) {
            dp[i][0] = 1;
        }
        for (int j = 1;j < aim+1;j++) {
            dp[0][j] = j%penny[0] == 0?1:0;  //只需要算dp[0][j]
        } 
         
        for (int i = 1;i < n;i++) {
            for (int j = 1;j < aim+1;j++) {
                dp[i][j] = (j-penny[i]) >= 0?(dp[i-1][j] + dp[i][j-penny[i]]):dp[i-1][j];  //这是关键,不用管penny【i】到底使用了几次,直接减去1次使用就好               
            }
        }
        return dp[n-1][aim];
    }
};

解法2:与上面的问题一样,只不过在求dp时只使用1维数组来做;使用迭代,时间复杂度一样:

class Exchange {
public:
    int countWays(vector<int> penny, int n, int aim) {
        vector<int> dp(aim + 1);
        for (int i = 0; i <= aim; i++)
            if (i % penny[0] == 0)
                dp[i] = 1;
 
        for (int i = 1; i < n; i++)
            for (int j = 1; j <= aim; j++)
                if ( j >= penny[i]) //条件,如果不满足就直接等于上轮的结果,不用做修改
                    dp[j] += dp[j - penny[i]];
        return dp[aim];
    }
};

问题2:跳台阶问题:

其实是斐波那契问题,f(n)=f(n-1)+f(n-2)

#include <iostream>
using namespace std;
int main(){
    int step;
    while(cin>>step){
        vector<int> dp(2,1); //初始化赋值
        dp[1]=2;
        int temp;
        for(int i=3;i<=step;i++){
            temp=dp[0];
            dp[0]=dp[1];
            dp[1]=dp[1]+temp;
        }
        if(step==1) dp[1]=1;;
        cout<<dp[1]<<endl;
    }
    return 0;
}

问题3:走矩阵,求路劲最小和,或者是求整个路径

n×m的map,则 f(n,m)=min(f(n-1,m),f(n,m-1))+map[n][m];
由于这里和问题1类似,可以只用到一个一维数组求解;

class MinimumPath {
public:
    int getMin(vector<vector<int> > map, int n, int m) {
        vector<int> dp(m,0);
        dp[0] = map[0][0];
        for (int i = 1,j = 0;i < m;i++,j++) {
            dp[i] = map[0][i]+dp[j];
        }
         
        for (int i = 1;i < n;i++) {
            dp[0] += map[i][0];    //不能忘了dp[0]的更新
            for (int j = 1;j < m;j++) {
                dp[j] = min(dp[j],dp[j-1])+map[i][j]; //如果求路径,则在这里记录,需要额外存储空间
            }
        }
        return dp[m-1];
    }
};

问题4:最长上升子序列问题(LIS)

解法:O(N方)用dp数组的dp[i]记录下以A[i]结尾的递增子序列中最长的长度,计算dp[i+1]时,遍历A[0~i]找到比A[i+1]小的元素,再比较与这些元素对应的dp数组中的值,找到最大的一个再加1,赋值给dp[i+1]。

class LongestIncreasingSubsequence {
public:
    int getLIS(vector<int> A, int n) {
        if (A.empty()||n == 0)
            return 0; 
        vector<int> dp(n,0);
        dp[0] = 1;
        int resMax = 0;
        for (int i = 1;i < n;i++) {
            int tempMax = 0;
            for (int j = 0;j < i;j++) {
                if (A[i] > A[j])
                    tempMax = max(tempMax,dp[j]);
            }
            dp[i] = ++tempMax;
            resMax = max(resMax,dp[i]);  //记录最大的上升子序列长度,因为当前i可能并不在最长上升子序列中
        }
        return resMax;
    }
};

如上的实现复杂度为N方,可以增加归纳的假设,增加b[k]存储长度为k最长子序列最小结尾元素,那么可以利用二分查找,使用logn查找到插入点,对于每次比较,要么直接比较b【k】比它大直接k+1,更新b【k+1】,要么二分查找到位置,更新b【j】,所以最终复杂度为nlogn(如果数据量大的话,使用该算法较好)

C++动态规划dp算法题

问题5:最长公共子序列长度(LCS)

C++动态规划dp算法题

C++动态规划dp算法题

上图可以看出使用了斜侧的比较,所以不能再使用1维数组了

class LCS {
public:
    int findLCS(string A, int n, string B, int m) {
        if (A.empty()||n==0||B.empty()||m==0)
            return 0;
        vector<vector<int> > dp(n,vector<int>(m));
        //下面是两个for的初始化,当出现第一个相等的时,后面的都直接赋值为1;
        for (int i = 0;i < m;i++) {
            if (A[0] == B[i]) {
                for (int j = i;j < m;j++)
                    dp[0][j] = 1;
                break ;
            }             
        }
        for (int i = 0;i < n;i++) {
            if (B[0] == A[i]) {
                for (int j = i;j < n;j++)
                    dp[j][0] = 1;
                break ;
            }             
        }
        for (int i = 1;i < n;i++) {
            for (int j = 1;j < m;j++) {
                if (A[i] == B[j])
                    dp[i][j] = dp[i-1][j-1]+1;
                else
                    dp[i][j] = max(dp[i-1][j],dp[i][j-1]);
            }
        }
        return dp[n-1][m-1];
    }
};

上面的方法中初始化第一行和第一列有点麻烦,增加了额外的语句,可以增加数组一行和一列来优化代码:

class LCS {
public:
    int findLCS(string A, int n, string B, int m) {       
        vector<vector<int> > dp(n+1,vector<int>(m+1,0));
        for (int i =1;i<=n ;++i){           
            for (int j=1; j<=m; ++j){
                if (A[i-1] == B[j-1]){
                    dp[i][j]  = dp[i-1][j-1]+1; //第1行也可以照此直接初始化
                }
                else {
                    dp[i][j] = max( dp[i-1][j] ,dp[i][j-1]);
                }               
            }
        }
        return dp[n][m];
    }
};

问题6:背包

N件物品,价值记录在数组V,重量记录在数组W,背包总重量最大为cap,要求总价值最大;

class Backpack {
public:
    int maxValue(vector<int> w, vector<int> v, int n, int cap) {
        if (w.empty()||v.empty()||n==0||cap==0)
            return 0;
        vector<vector<int> > dp(n,vector<int>(cap+1));
         
        for (int j = 1;j < cap+1;j++) {
            dp[0][j] = w[0] <= j?v[0]:0;
        }
        for (int i = 0;i < n;i++) {
            dp[i][0] = 0;
        }
         
        for (int i = 1;i < n;i++) {
            for (int j = 1;j < cap+1;j++) {
                if (w[i] > j)
                    dp[i][j] = dp[i-1][j];
                else
                    dp[i][j] = max(dp[i-1][j],v[i]+dp[i-1][j-w[i]]); //由于该问题每个物品最多只能放1件,如果不限制个数的话,则在这里修改条件
            }
        }
        return dp[n-1][cap];
    }
};

由于没有用到斜侧的比较,所以可以使用1维的数组:

class Backpack {
public:
    int maxValue(vector<int> w, vector<int> v, int n, int cap) {
        if (w.empty()||v.empty()||n==0||cap==0)
            return 0;
        vector<int> dp(cap+1,0);       
        for (int i = 0;i < n;i++) {
            vector<int> last(dp);
            for (int j = 1;j < cap+1;j++) {
                dp[j] = j < w[i]?last[j]:max(last[j],v[i]+last[j-w[i]]);               
            }
        }
        return dp[cap];
    }
};

相关推荐