动态规划(DP)(1)

• 动态规划的两个要求：
  • 最优子结构：大问题的（最优）解可以由小问题的（最优）解推出。注意在问题拆解的过程中不能无限递归。
  • 无后效性：未来与过去无关，一旦得到了一个小问题的解，如何得到它的解的过程不影响大问题的求解。
• 动态规划的两个元素：
  • 状态：求解过程进行到了哪一步，可以理解为一个子问题；
  • 转移：从一个状态（小问题）的（最优）解推导出另一个状态（大问题）的（最优）解的过程。

例题

最短路

给定 n 个点 m 条边的有向图，每条边有个边权，代表经过这条边需要花费的时间，我们只能从编号小的点走到编号大的点，问从 1 号点走到 n号点最少需要花费多少时间？
输入格式
第一行两个整数 n,m。 接下来 m行，每行三个整数 u,v,w，表示存在一条从 u 到 v 的边权为 w的有向边。 保证存在一条1到 n的路径。

输出格式
输出一个数，表示答案。

数据规模 对于所有数据，
1≤n≤103,1≤m≤103,1≤u<v≤n,1≤w≤103

#include<bits/stdc++.h>

using namespace std;

int n, a[1001][1001], f[1001], m;

int main(){
    scanf("%d%d", &n, &m);
    memset(a, 127, sizeof(a));
    for(int i = 1; i <= m; i++){
        int x, y, z;
        scanf("%d%d%d", &x, &y, &z);
        a[x][y] = min(a[x][y], z);
    }
    memset(f, 127, sizeof(f));
    f[1] = 0;
    for(int i = 2; i <= n; i++)
        for(int j = 1; j <i; j++)
            if(a[j][i] < 1 << 30 && f[j] < 1 << 30)
                f[i] = min(f[i], f[j] + a[j][i]);
    printf("%d", f[n]);
}

最长上升子序列

给定一个长度为 n 的数组 a1,a2,…,an，问其中的最长上升子序列的长度。也就是说，我们要找到最大的 m 以及数组 p1,p2,…,pm，满足 1≤p1<p2<⋯<pm≤n 并且 ap1<ap2<⋯<apm。
输入格式
第一行一个数字 n。
接下来一行 n个整数 a1,a2,…,an。
输出格式
一个数，表示答案。
数据规模
对于所有数据，保证 1≤n≤1000,1≤ai≤109。

#include<bits/stdc++.h>

using namespace std;

int n, a[1001], f[1001];

int main(){
    scanf("%d", &n);
    for(int i = 1; i <= n; i++)
        scanf("%d", &a[i]);
    for(int i = 1; i <= n; i++){
        f[i] = 1;
        for(int j = 1; j < i; j++)
            if(a[i] > a[j])
                f[i] = max(f[i], f[j] + 1);
    }
    int ans = 0;
    for(int i = 1; i <= n; i++)
        ans = max(ans, f[i]);
    printf("%d", ans);
}

最长公共子序列

给定一个长度为 n 的数组 a1,a2,…,an 以及一个长度为 m 的数组 b1,b2,…,bm，问 a 和 b的最长公共子序列的长度。
也就是说，我们要找到最大的 k以及数组 p1,p2,…,pk，数组 l1,l2,…,lk 满足 1≤p1<p2<⋯<pk≤n 并且1≤l1<l2<⋯<lk≤m 并且对于所有的 i(1≤i≤k) ，api=bli。
输入格式
第一行两个整数 n,m。
接下来一行 n个整数，a1,a2,…,an。
接下来一行 m个整数，b1,b2,…,bm。
输出格式
输出一个整数，表示答案。
数据规模
对于所有数据，保证 1≤n``,m≤1000,1≤ai,bi≤103。

#include<bits/stdc++.h>

using namespace std;

int n, m, a[1001], b[1001], f[1001][1001];

int main(){
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; i++)
        scanf("%d", &a[i]);
    for(int i = 1; i <= m; i++)
        scanf("%d", &b[i]);
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++){
            f[i][j] = max(f[i-1][j], f[i][j-1]);
            if(a[i] == b[j])
                f[i][j] = max(f[i][j], f[i-1][j-1] + 1);
        }
    printf("%d", f[n][m]);
}