POJ 3132 双约束背包问题

    多重背包的问题是，给定物品，且有物品的数量限制，问方法种类或者最大的容量等，这里的限制是限制在物品上；

    另外一种约束是限制在背包上，即：背包最多可以装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，为什么两者没有区别呢？

1. #include <algorithm>
2. #include <numeric>
3. #include <cmath>
4. #include <iostream>
5. #include <cstdio>
6. #include <iomanip>
7. #include <vector>
8. using namespace std;
9. int primes[187],n=0;
10. bool isPrime[1120];
11. int package[1121][15];
12. void getPrime()
13. {
14.     int t;
15.     memset(isPrime,-1,sizeof isPrime);
16.     for (int i=2;i<1120;++i)
17.     {
18.         if(isPrime[i]) primes[n++]=i;
19.         for (int j=0;j<n && (t=i*primes[j])<1120;++j)
20.         {
21.             isPrime[t]=false;
22.             if(i%primes[j]==0) break;
23.         }
24.     }
25. }
26. int main()
27. {
28.     getPrime();
29.     int a,b;
30.     while (scanf("%d%d",&a,&b))
31.     {
32.         if(!a && !b) break;
33.         memset(package,0,sizeof package);
34.         package[0][0]=1;
35.         for (int i=0;i<187 && primes[i]<=a;++i)
36.             for (int j=a;j>=primes[i];--j)
37.                 for (int k=1;k<=b;++k)
38.                 package[j][k]+=package[j-primes[i]][k-1];
39.         printf("%d/n",package[a][b]);
40.     }
41.     return 0;
42. }