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，为什么两者没有区别呢？

#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;

}