算法——数学与简单DP（2）

DP：从集合的角度分析DP问题（一般用来求最值/个数问题）

ps：DP问题的数量级很大，用暴力会超时

1.步骤一：化零为整（把零散的化为集合）

2.步骤二：化整为零（把集合划分）

以01背包问题为例：(组合模型，所有选择问题都可以用01背包的思想)

有 N件物品和一个容量是 V的背包。每件物品只能使用一次。

第 i件物品的体积是 vi，价值是 wi。

求解将哪些物品装入背包，可使这些物品的总体积不超过背包容量，且总价值最大。
输出最大价值。

输入格式

第一行两个整数，N，V，用空格隔开，分别表示物品数量和背包容积。

接下来有 N 行，每行两个整数 vi,wi，用空格隔开，分别表示第 i 件物品的体积和价值。

输出格式

输出一个整数，表示最大价值。

数据范围

0<N,V≤1000
0<vi,wi≤1000

输入样例

4 5
1 2
2 4
3 4
4 5

输出样例：

8

 按以上方式分析：

（1）状态f[i][j]定义：前 i 个物品，背包容量 j下的最优解（最大价值）：当前的状态依赖于之前的状态，可以理解为从初始状态f[0][0] = 0开始决策，有 N件物品，则需要 N次决 策，每一次对第 i
 件物品的决策，状态f[i][j]不断由之前的状态更新而来。
（2）当前背包容量不够（j < v[i]），没得选，因此前 i 个物品最优解即为前 i−1个物品最优解：对应代码：f[i][j] = f[i - 1][j]。
（3）当前背包容量够，可以选，因此需要决策选与不选第 i个物品：

选：f[i][j] = f[i - 1][j - v[i]] + w[i]。
不选：f[i][j] = f[i - 1][j] 。
我们的决策是如何取到最大价值，因此以上两种情况取 max() 。

代码实现：

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;

const int N = 1010;
int n,m;
int v[N],w[N];//分别记录第i个物品的容量和价值
int f[N][N];//记录状态的集合，在选择前n件物品，且容量不超过m的集合

int main()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>v[i]>>w[i];
    
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            f[i][j] = f[i-1][j];//考虑左半边的情况
            if(j>=v[i]) f[i][j] = max(f[i-1][j],f[i-1][j-v[i]]+w[i]);//考虑右半边的情况，要先判断右半边是否存在
        }
        cout<<f[n][m]<<endl;
        return 0;
}

例题2：路线问题

注：边界的划分：从题目的实际含义出发，是否从0开始，还是从1开始

一般来说，如果题目中出现 i-1，j-1 则下标从1开始，否则从0开始 

按照上述的分析思路进行集合的划分

例题：摘花生

Hello Kitty想摘点花生送给她喜欢的米老鼠。

她来到一片有网格状道路的矩形花生地(如下图)，从西北角进去，东南角出来。

地里每个道路的交叉点上都有种着一株花生苗，上面有若干颗花生，经过一株花生苗就能摘走该它上面所有的花生。

Hello Kitty只能向东或向南走，不能向西或向北走。

问Hello Kitty最多能够摘到多少颗花生。

输入格式

第一行是一个整数T，代表一共有多少组数据。

接下来是T组数据。

每组数据的第一行是两个整数，分别代表花生苗的行数R和列数 C。

每组数据的接下来R行数据，从北向南依次描述每行花生苗的情况。每行数据有C个整数，按从西向东的顺序描述了该行每株花生苗上的花生数目M。

输出格式

对每组输入数据，输出一行，内容为Hello Kitty能摘到得最多的花生颗数。

数据范围

1≤T≤100,
1≤R,C≤100,
0≤M≤1000

输入样例：

2
2 2
1 1
3 4
2 3
2 3 4
1 6 5

输出样例：

8
16

代码实现：

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;

const int N = 110;
int n,m;
int f[N][N];//所有路径的最大值集合
int w[N][N];//每个节点的花生数量

int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        cin>>n>>m;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++) cin>>w[i][j];
        
        //枚举每一种情况
        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])+w[i][j];
            }
            cout<<f[n][m]<<endl;
    }
    return 0;
}

 例题3：最长上升子序列问题

每次需要判断一下，若a（k）>= a(i) 则该方法不合法 

这里是指从第一个数字到a[k]结尾的最长子序列

最长上升子序列

给定一个长度为 N的数列，求数值严格单调递增的子序列的长度最长是多少。

输入格式

第一行包含整数 N。

第二行包含 N个整数，表示完整序列。

输出格式

输出一个整数，表示最大长度。

数据范围

1≤N≤1000，
−1e9≤数列中的数≤1e9

输入样例：

7
3 1 2 1 8 5 6

输出样例：

4

代码实现：

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;

const int N = 1010;
int f[N],w[N];//f[n]为记录最长序列状态的数组，w[n]为原数组
int n;

int main()
{
    cin>>n;
    for(int i=1;i<=n;i++) cin>>w[i];
    for(int i=1;i<=n;i++)
    {
       f[i] = 1; //如果集合中只有一个数字的情况
       for(int j=1;j<i;j++)
       {
           if(w[j]<w[i]) f[i] = max(f[i],f[j]+1);
       }
    }
    //依次遍历集合，找到最大的结果
    int res = 0;
    for(int i=1;i<=n;i++)
    {
        res = max(res,f[i]);
    }
    cout<<res<<endl;
    return 0;
}