【洛谷】P1048 [NOIP2005普及组]采药

题目描述

辰辰是个天资聪颖的孩子，他的梦想是成为世界上最伟大的医师。为此，他想拜附近最有威望的医师为师。医师为了判断他的资质，给他出了一个难题。医师把他带到一个到处都是草药的山洞里对他说：“孩子，这个山洞里有一些不同的草药，采每一株都需要一些时间，每一株也有它自身的价值。我会给你一段时间，在这段时间里，你可以采到一些草药。如果你是一个聪明的孩子，你应该可以让采到的草药的总价值最大。”

如果你是辰辰，你能完成这个任务吗？

输入格式

第一行有 2 个整数 T（1≤T≤1000）和 M（1≤M≤100），用一个空格隔开，T 代表总共能够用来采药的时间，M 代表山洞里的草药的数目。

接下来的 M 行每行包括两个在 1 到 100 之间（包括 1 和 100）的整数，分别表示采摘某株草药的时间和这株草药的价值。

输出格式

输出在规定的时间内可以采到的草药的最大总价值。

输入输出样例

输入 #1

70 3
71 100
69 1
1 2

输出 #1

3

说明/提示

【数据范围】

• 对于 30% 的数据，M≤10；
• 对于全部的数据，M≤100。

算法知识——动态规划(DP)

一.动态规划算法基本思想：

如果各个子问题不是独立的(即重复的)，且不同子问题的个数是有限的，那么如果我们可以保存已经解决的子问题的答案，而在需要时再找出已经求得的答案，这样就可以避免大量的重复计算，以此减少时间复杂度。

二.可以使用动态规划求解的问题需要满足以下基本特征：

1.最优子结构：即大问题的最优解包含子问题的最优解。

2.重叠子问题：即子问题的解会被多次用到。

3.无后效性：即求解大问题时不会改变子问题的最优解。

4.复杂问题线性化

三.实现思路：

1.定义数组，并明确其含义。

2.找出子问题和大问题之间的递推关系式(状态转移方程)，并且用编程语言表达。

3.找出初始值，即无法使用第二步中关系式求解得出的值，在程序中需要手动赋值。

4.找出边界情况，在程序中需要特殊讨论。

01背包问题思路：

1.定义数组bag[i][j],含义为前i件物品放入大小为j的背包能获得的最大价值。

2.找出状态转移方程，定义v[i]为第i个物品的价值，t[i]为第i个物品的背包消耗量：

• 若不将第i物品放入背包，那最大价值不变，即bag[ i ][ j ] = bag[ i-1 ][ j ]。
• 若将第i物品放入背包，此时背包剩下空间为j-t[ i ]。那么将前i件物品放入大小为j的背包能获得的最大价值，即为第i物品的价值v[ i ]加上前i - 1件物品放入空间为j - t[ i ]的背包能获得的最大价值，即bag[ i ][ j ] = bag[ i-1 ][ j - t[ i ] ] + v[ i ]。
• 以上两种情况求最大值。

3.无。

4.找出边界情况，即背包剩余容量不足以将第i件物品装入，此时bag[i][j] = bag[i-1][j]。

题目思路

本题属于动态规划类问题当中的01背包问题，题目特征为在给定的条件（如总耗时、总重量）之内，使背包中装入的物品价值最大。

AC代码

法一:

#include<stdio.h>
#include<algorithm>
using namespace std;

int main(){
	int t,m,herb[102][2],bag[102][1002];
	scanf("%d %d",&t,&m);
	for(int m1 = 1; m1 <= m; m1++){
		scanf("%d %d",&herb[m1][0],&herb[m1][1]);
	}
	//01背包求解
	for(int i = 1; i <= m; i++){
        for(int j = 0; j <= t; j++){
            if(j >= herb[i][0]){
                bag[i][j] = max(bag[i-1][j], bag[i-1][j-herb[i][0]] + herb[i][1]);
            }
            else{ //边界情况讨论
                bag[i][j] = bag[i-1][j];
            }
        }
    }
    printf("%d",bag[m][t]);
	return 0;
}

法二：

一维数组滚动优化，空间复杂度由O(m*n)减少至O(m)。

特别注意应当用正序物品遍历逆序背包，以此避免数据污染。

#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;

int main(){
	int t,m,herb[102][2],bag[1002];
	memset(bag,0,sizeof(bag));
	scanf("%d %d",&t,&m);
	for(int m1 = 0; m1 < m; m1++){
		scanf("%d %d",&herb[m1][0],&herb[m1][1]);
	}
	//01背包求解
	for(int i = 0; i < m; i++){
        for(int j = t; j >= 0; j--){
            if(j >= herb[i][0]){
                bag[j] = max(bag[j], bag[j - herb[i][0]] + herb[i][1]);
            }
            else{ //边界情况讨论
                bag[j] = bag[j];
            }
        }
    }
    printf("%d",bag[t]);
	return 0;
}