4120:硬币

原题链接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]<<" ";
}