原题链接http://bailian.openjudge.cn/practice/4120/
本题参考了网上大佬的代码,此贴目的是记录本人做题时的思考
思路就是枚举每一个硬币,然后dp,看剩下的硬币能否构成目标数字(sum),不能构成的话,这枚硬币就是目标硬币。
思路虽然是这样,但每次从头枚举后dp会超时。
这里的处理方法是先dp(dp)所有硬币(使用了滚动数组),另开一个数组存枚举硬币时的dp(ans)(也是滚动数组) 从原始dp到ans的状态转移方程是ans【j】=dp【j】-ans【j-a[i]】。
解释一下这个方程为什么成立:
1.默认的公理(?):0 1背包问题的物品无论怎么排序,最后的状态一定是一样的(可以自己用纸比划一下)。
2.在1.的前提下,我们可以知道最终的dp滚动数组可以代表任何一个硬币的状态转移过程(把要枚举的硬币放最后)
3.我们知道dp的状态转移方程是dp[i][j]=dp[i-1][j]+dp[i-1][j-a[i]],其中前面的表示不选择a[i],后面的表示选择a[i],那么去除a[i]的情况数量就是dp[i-1][j]=dp[i][j]-dp[i-1][j-a[i]],及dp[i-1][j]为此时的ans[i-1(-1是去除了a【i】,物品总数减一)][j],变成滚动数组就是ans[j]=dp[j]-ans[j-a[i]]。 所以可以证明状态转移方程的正确性 然后就可以愉快的敲代码了
#include<iostream>
#include<cstdio>
#include<string.h>
using namespace std;
int a[210],dp[10010],ans[10010],res[200];
int main(){
int n,sum;
cin>>n>>sum;
for(int i=1;i<=n;i++)scanf("%d",a+i);
dp[0]=1;//初始化
for(int i=1;i<=n;i++)
for(int j=sum;j>=a[i];j--){
dp[j]+=dp[j-a[i]];//先求出所有硬币的dp状态
}
int x=0;
for(int i=1;i<=n;i++){//枚举每一个要去除的硬币
memset(ans,0,sizeof(ans));
for(int j=0;j<=sum;j++){
if(j<a[i])ans[j]=dp[j];
else ans[j]=dp[j]-ans[j-a[i]];
}
if(ans[sum]==0)res[x++]=a[i];/*ans【sum】等于0,满足条件,存答案*/
}
cout<<x<<endl;
for(int i=0;i<x;i++)cout<<res[i]<<" ";
}