多重背包的问题是,给定物品,且有物品的数量限制,问方法种类或者最大的容量等,这里的限制是限制在物品上;
另外一种约束是限制在背包上,即:背包最多可以装k种物品,每种物品都是无限的(完全)的,问有多少种背包方式。或者说,这里的背包有两重的约束,即背包的容量以及背包最多能装的物品的个数都有限制,不妨成为“双约束背包”。
这种背包问题的思考方式是这样的:用两维背包,package[a][b],其中a是背包的最大容量,b是背包最多能装的数量,两个都进行约束之后,背包策略如下:
对于每一种物品,进行完全背包,对于背包容量为j的情况,考察它由b,b-1,b-2,…,1个物品构成,其递归式子为package[j][k]+=package[j-wt[i]][k-1];这里需要注意到是,背包的顺序,是从b倒序到1,也可以从1升序到b,为什么两者没有区别呢?
- #include <algorithm>
- #include <numeric>
- #include <cmath>
- #include <iostream>
- #include <cstdio>
- #include <iomanip>
- #include <vector>
- using namespace std;
- int primes[187],n=0;
- bool isPrime[1120];
- int package[1121][15];
- void getPrime()
- {
- int t;
- memset(isPrime,-1,sizeof isPrime);
- for (int i=2;i<1120;++i)
- {
- if(isPrime[i]) primes[n++]=i;
- for (int j=0;j<n && (t=i*primes[j])<1120;++j)
- {
- isPrime[t]=false;
- if(i%primes[j]==0) break;
- }
- }
- }
- int main()
- {
- getPrime();
- int a,b;
- while (scanf("%d%d",&a,&b))
- {
- if(!a && !b) break;
- memset(package,0,sizeof package);
- package[0][0]=1;
- for (int i=0;i<187 && primes[i]<=a;++i)
- for (int j=a;j>=primes[i];--j)
- for (int k=1;k<=b;++k)
- package[j][k]+=package[j-primes[i]][k-1];
- printf("%d/n",package[a][b]);
- }
- return 0;
- }