PAT A1068 Find More Coins
- 这题花了挺长时间 ——虽然还是抄的书上的,,刚开始想到应该是排序,然后每个位置有选和不选两种,复杂度O(2^n),然后应该是利用动态规划优化,想了半天睡着了。。。
- 看了书上发现是叫背包问题,以前还以为背包是个贪心。。。——真是to young
- 看了书上的解析和套路,理解了好久,睡着了好几次,终于大概归纳如下:首先是一个二维的描述,状态转换方程为dp[i][j] = max(dp[i - 1][j], dp[i - 1][j - w[i]] + v[i]);这里dp的第一维表示已经考虑过的物品范围,第二维表示到达的重量,元素值本身表示此种情况下能获得的最大价值。则dp[i - 1][j]就表示在前i-1个物品中挑挑拣拣能组合出来重量为j的情况下,所能获得的最大价值。
- 所以状态方程描述的是在面对第i件物品之时,有要和不要两种选择:if不要,则此处的最大价值dp[i][j]只需继承前面的dp[i - 1][j];if要,因为此时要满足的重量是j,当前物品重量是w[i],则之前的重量应该是j - w[i],所以dp[i][j]就取前面i-1个物品挑挑拣拣达到重量j - w[i]时的最大价值 + 当前物品的价值,即dp[i - 1][j - w[i]] + v[i]
- 下面紧接着其一维描述:dp[j] = max(dp[j], dp[j - w[i]] + v[i]);直接把上面的二维复制下来删去第一个维度就是了 ——因为这两个东西在意义上是完全等价的。我最开始分开看怎么也联系不到一起,书上关于二维到一维转换的解释对我来说就像是某些数学题解析中的“显然有~”,当然我更不明白为啥内圈循环必须要逆序
- 看了一些文章后在https://blog.csdn.net/wumuzi520/article/details/7014559里面终于看懂了。因为对于dp[i][j]来说,他的取值只需要考虑dp[i - 1][j]和dp[i - 1][j - w[i]],这两个数的第一维都是i-1(即i的前一个状态),and第二维都<=j(即j的左边) ——根据第一维可知,我们只需要存dp[i - 1]这个一维数组就行了;根据第二维可知,对于重量j来说我们只需要用到dp[i - 1]这个一维数组下标从0到j的部分(其实只有j和j-w[i])。
- 所以我们可以用一个以重量为下标的全0数组为初始状态(i=0),内层循环j从后向前依次赋值,每次赋值只会关系到当前下标及其之前的部分,而与他后面的部分无关(后面都是本状态下已经赋值好了的),一趟结束后,本状态赋值完毕,进入下一状态(i++),j同样从后向前更新赋值
- so内圈循环必须从后往前,如果从前往后的话,当你计算下标j上值的时候,0~j上的值已经不是前一个状态的了(当然这些值也都是错的)。
- 对应到状态方程dp[j] = max(dp[j], dp[j - w[i]] + v[i]);dp就是这样的数组,在此步操作之前,下标0~j上是前一个状态的值(即前i-1个物品挑挑拣拣总重量<=j下的所有结果),j后面是本状态(i个物品时)已经赋值好了的。
- 或许j应该叫做“假设背包的容量”或者“背包的子容量”更合适一点?毕竟大部分时候物品的重量不够j。也许之后有大佬看到这里能给我生动而科学的解释一下吧,不知道以上自圆其说的总结是否正确
- 后半个问题就是回溯路径了,如果用的是二维的方程,则之前的所有状态上的数据都保存了,直接逆向操作就行了(起点是dp[amount] == amount);一维方程的话没有了之前的状态,只能在赋值dp的时候顺便用二维数组记下路径,之后我们就得到了在哪些重量之下有哪些物品被选择了,然后从重量=amount开始回溯
- dp.resize(amount + 1);下标是重量莫要搞成num
- 保证输出最小序列的关键是sort,循环的顺序,以及if(dp[j] <= dp[j - val[i]] + val[i])
#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <queue>
#include <math.h>
#include <set>
#include <map>
#include <unordered_map>
#include <stack>
using namespace std;
#define MAXSIZE 10010
bool cmp(int a,int b){
return a > b;
}
vector<int> val,dp,res;
int path[MAXSIZE][MAXSIZE] = {0};
#define DEBUG 1
int main(){
#ifdef DEBUG
freopen("1.txt","r",stdin);
#endif
int num,amount;
cin >> num >> amount;
val.resize(num + 1);
dp.resize(amount + 1);
for(int i = 1; i <= num;i ++) cin >> val[i];
sort(val.begin() + 1,val.end(),cmp);
for(int i = 1; i <= num;i ++){
for(int j = amount; j >= val[i];j --){
if(dp[j] <= dp[j - val[i]] + val[i]){
dp[j] = dp[j - val[i]] + val[i];
path[i][j] = 1;
}
}
}
if(dp[amount] != amount) cout << "No Solution";
else{
int value = amount;
for(int i = num; i >= 1;i --){
if(path[i][value] == 1){
res.push_back(val[i]);
value -= val[i];
}
}
for(int i = 0; i < res.size();i ++){
cout << res[i];
if(i < res.size() - 1) cout << ' ';
}
}
return 0;
}