【Ybtoj 第19章例4】硬币方案【背包问题】

---

解题思路

依旧是多重背包。。用二进制拆分优化了一下。

设$f[j]$表示前i种硬币能否组成面值j，$d[j]$表示f[j]真的时候至少需要多少种第i种硬币。
所以在$f[j-a_i]$为真 $f[j]$为假的时候，f[j]的状态由$f[j-a_i]$转移而来，这个时候$d[j]=d[j-a_i]+1;$
最后看有多少个为真的状态就有多少种面值

---

代码

#include<iostream>
#include<cstdio>
#include<iomanip>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;

int n,m,N,ans,f[100010],aa[100010],a[100010],b[100010],c[100010],d[100010];

int main() {
	while(scanf("%d%d",&n,&m)) {
		if(n==0&&m==0)return 0;
		ans=0;
		memset(f,0,sizeof(f));
		memset(c,0,sizeof(c));
		for(int i=1; i<=n; i++)
			scanf("%d",&aa[i]);
		for(int i=1; i<=n; i++)
			scanf("%d",&b[i]);
		N=0;
		for(int i=1; i<=n; i++) {
			int bb=b[i];
			for(int j=1; j<bb; j*=2) {
				N++;
				a[N]=aa[i]*j;
				c[N]++;
				bb-=j;
			}
			if(bb) {
				a[++N]=aa[i]*bb;
				c[N]++;
			}
		}
		f[0]=1;
		sort(a+1,a+N+1);
		for(int i=1; i<=N; i++) {
			memset(d,0,sizeof(d));
			for(int j=a[i]; j<=m; j++)
				if(!f[j]&&f[j-a[i]]&&d[j-a[i]]<c[i])
					d[j]=d[j-a[i]]+1,f[j]=1;

		}
		for(int i=1; i<=m; i++) {
			if(f[i]) {
				ans++;
			}
		}
		printf("%d\n",ans);
	}
}