《王道论坛计算机考研机试指南》第七章【动态规划】

• 递推求解
• 最长递增子序列（LIS）
• 最长公共子序列
• 动态规划问题分析举例
• 背包
  
  • 0-1背包
  • 完全背包
  • 多重背包

一、递推求解

N阶楼梯上楼问题

#include<cstdio>
using namespace std;
int F[91];
int main() {
    F[1] = 1;
    F[2] = 2;
    for (int i = 3; i <= 90; i++)
        F[i] = F[i-1] + F[i-2];
    int n;
    while (scanf("%d", &n) != EOF)
        printf("%d\n", F[n]);
    return 0;
}

不容易系列之一

#include<cstdio>
using namespace std;
long long F[91];
int main() {
    F[1] = 0;
    F[2] = 1;
    for (int i = 3; i <= 20; i++)
        F[i] = (i-1)*F[i-1] + (i-1)*F[i-2];
    int n;
    while (scanf("%d", &n) != EOF)
        printf("%lld\n", F[n]);
    return 0;
}

二、最长递增子序列（LIS）

拦截导弹

#include <cstdio>
#include <algorithm>
using namespace std;
int list[26];  // 按袭击事件顺序保存各导弹高度
int dp[26]; // dp[i]保存以第i个导弹结尾的最长不增子序列长度
int main() {
    int n;
    while (scanf("%d", &n) != EOF) {
        for (int i = 1; i <=n ; i++)
            scanf("%d", &list[i]);
        for (int i = 1; i <= n; i++) {
            int tmax = 1;
            for (int j = 1; j < i; j++) {
                if (list[j] >= list[i]) {
                    tmax = max(tmax, dp[j]+1);
                }
            }
            dp[i] = tmax;
        }
        int ans = 1;
        for (int i = 1; i <= n; i++) {
            if (dp[i] > ans)
                ans = dp[i];
        }
        printf("%d\n", ans);
    }
    return 0;
}

三、最长公共子序列（LCS）

Coincidence

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int dp[101][101];
int main() {
    char s1[101], s2[101];
    while (scanf("%s%s", s1, s2) != EOF) {
        int l1 = strlen(s1);
        int l2 = strlen(s2);
        for (int i = 0; i <= l1; i++)
            dp[i][0] = 0;
        for (int j = 0; j <= l2; j++)
            dp[0][j] = 0;
        for (int i = 1; i <= l1; i++) {
            for (int j = 1; j <= l2; j++) {
                if (s1[i-1] != s2[j-1])
                    dp[i][j] = max(dp[i-1][j], dp[i][j-1]);
                else
                    dp[i][j] = dp[i-1][j-1] + 1;
            }
        }
        printf("%d\n", dp[l1][l2]);
    }
    return 0;
}

四、动态规划问题分析举例

搬寝室

#include <cstdio>
#include <algorithm>
using namespace std;
const int Max = 0x0fffffff;   // 预定义最大的int取值为无穷
int list[2001];        // 保存每个物品的重量
int dp[1001][2001];    // 保存每个状态，dp[i][j]为在前j件物品中选择i对物品时最小的疲劳度
int main() {
    int n, k;
    while (scanf("%d%d", &n, &k) != EOF) {
        for (int i = 1; i <= n; i++)
            scanf("%d", &list[i]);
        sort(list+1, list+n+1);   // 使所有物品按照重量递增排序
        for (int i = 0; i <= n; i++) {
            for (int j = 0; j <= k; j++) {
                dp[j][i] = Max;
            }
        }
        for (int i = 0; i <= n; i++)   // 初始值
            dp[0][i] = 0;
        for (int i = 1; i <= k; i++) {  // 递推求得每个状态
            for (int j = 2*i; j <= n; j++) {
                dp[i][j] = min(dp[i][j-1], dp[i-1][j-2] + (list[j]-list[j-1]) * (list[j]-list[j-1]) );
            }
        }
        printf("%d\n", dp[k][n]);
    }
    return 0;
}

五、背包

0-1背包

采药

#include <cstdio>
#include <algorithm>
using namespace std;
struct E { // 保存物体信息结构体
    int w;  // 物品的体积
    int v;  // 物品的价值
} list[101];
int dp[101][1001];    // 记录状态数组，dp[i][j]表示前i个物品的总体积不大于j的最大价值和

int main() {
    int s, n;
    while (scanf("%d%d", &s, &n) != EOF) {
        for (int i = 1; i <= n; i++)
            scanf("%d%d", &list[i].w, &list[i].v);
        for (int i = 1; i <= n; i++)   // 初始化状态
            dp[0][i] = 0;
        for (int i = 1; i <= n; i++) {  // 循环每一个物品
            for (int j = s; j >= list[i].w; j--) {
                dp[i][j] = max(dp[i-1][j], dp[i-1][j-list[i].w] + list[i].v);
            }
            for (int j = list[i].w-1; j >= 0; j--)
                dp[i][j] = dp[i-1][j];
        }
        printf("%d\n", dp[n][s]);

    }
    return 0;
}

#include <cstdio>
#include <algorithm>
using namespace std;
struct E { // 保存物体信息结构体
    int w;  // 物品的体积
    int v;  // 物品的价值
} list[101];
int dp[1001];

int main() {
    int s, n;
    while (scanf("%d%d", &s, &n) != EOF) {
        for (int i = 1; i <= n; i++)
            scanf("%d%d", &list[i].w, &list[i].v);
        for (int i = 1; i <= n; i++)   // 初始化状态
            dp[i] = 0;
        for (int i = 1; i <= n; i++) {  // 循环每一个物品
            for (int j = s; j >= list[i].w; j--) {
                dp[j] = max(dp[j], dp[j-list[i].w] + list[i].v);
            }
            for (int j = list[i].w-1; j >= 0; j--)
                dp[j] = dp[j];
        }
        printf("%d\n", dp[s]);

    }
    return 0;
}

完全背包

Piggy-Bank

#include <cstdio>
#include <algorithm>
using namespace std;
const int INF = 0x7fffffff;
struct E { // 代表钱币结构体
    int w;  // 重量
    int v;  // 价值
} list[501];
int dp[10001];

int main() {
    int T;
    scanf("%d", &T);
    while (T--) {
        int s, tmp;
        scanf("%d%d", &tmp, &s); // 输入空储蓄罐重量，和装满钱币储蓄罐重量
        s -= tmp;   // 计算钱币所占重量
        int n;
        scanf("%d", &n);
        for (int i = 1; i <= n; i++)
            scanf("%d%d", &list[i].v, &list[i].w);
        for (int i = 0; i <= s; i++)
            dp[i] = INF;
        dp[0] = 0;  // 因为要求所有物品恰好装满，所以初始时，除dp[0]外，其余dp[j]均为无穷（或者不存在）
        for (int i = 1; i <= n; i++) {   // 遍历所有物品
            for (int j = list[i].w; j <= s; j++) {   // 完全背包，顺序遍历所有可能转移的状态
                if (dp[j-list[i].w] != INF)
                       dp[j] = min(dp[j], dp[j-list[i].w] + list[i].v);
            }
        }
        if (dp[s] != INF)
            printf("The minimum amount of money in the piggy-bank is %d.\n", dp[s]);
        else
            printf("This is impossible.\n");

    }
    return 0;
}

多重背包

珍惜现在，感恩生活

#include <cstdio>
#include <algorithm>
using namespace std;
struct E { // 大米
    int v;  // 价格
    int w;  // 重量
} list[2001];
int dp[101];

int main() {
    int T;
    scanf("%d", &T);
    while (T--) {
        int s, n;
        scanf("%d%d", &s, &n);
        int cnt = 0;       // 拆分后的物品总数
        for (int i = 1; i <= n; i++) {
            int v, w, k;   // 价格，重量和袋数
            scanf("%d%d%d", &v, &w, &k);
            int c = 1;
            while (k-c > 0) {  // 对输入的数字k，拆分成1,2,4,k-(2^c-1)，其中c为使最后一项大于0的最大整数
                k -= c;
                list[++cnt].v = c*v;  // 拆分后的大米重量和价格均为组成该物品的大米的重量价格和
                list[cnt].w = c*w;
                c *= 2;
            }
            list[++cnt].v = v*k;
            list[cnt].w = w*k;
        }
        for (int i = 1; i <= s; i++)
            dp[i] = 0;
        for (int i = 1; i <= cnt; i++) {   // 对拆分后的物品进行0-1背包
            for (int j = s; j >= list[i].v; j--) {
                dp[j] = max(dp[j], dp[j-list[i].v] + list[i].w);
            }
        }
        printf("%d\n", dp[s]);      // 输出答案
    }
    return 0;
}