poj 1742 多重背包,可行性

点击打开链接

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
using namespace std;
const int M =1100;
long n,m;
int a[M],c[M];
bool dp[M*M]; // 只需判断能否凑成某个价钱k即可

// 给出物品种数和价值 每种物品数量有限 求这些物品能凑成不超过m的价钱有多少种

// 等价于 背包体积最大为m 求能填满 1~m 背包中成立的个数

// c[i]*a[i]>=m时 即物品可以看成无限件 就不需要拆分 -> 完全背包

void CompletePack(int value)
{

for(int j=value;j<=m;j++)
{
dp[j]=dp[j]||dp[j-value];
}

}
void ZeroPack(int value)
{
for(int j=m;j>=value;j--)
{
dp[j]=dp[j]||dp[j-value];
}
}
void MultiplePack(int value,int num)
{
if(value*num>=m)
CompletePack(value);

else
{
int k=1;
while(k<=num)
{
ZeroPack(value*k);
num-=k;
k=(k<<1);
}
if(num)
{
ZeroPack(value*num);
}
}

}
int main()
{

while(scanf("%d%d",&n,&m)&&(n||m))
{

for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]); //价值
}
for(int i=1;i<=n;i++)
{
scanf("%d",&c[i]); // 件数
}
dp[0]=1; // 前0件能装满体积为0的背包
for(int j=1;j<=m;j++)
{
dp[j]=0;
}
for(int i=1;i<=n;i++)
{
if(c[i])
MultiplePack(a[i],c[i]);
}
int ans=0;
for(int j=1;j<=m;j++)
{
if(dp[j])
ans++;
}
printf("%d\n",ans);
}
return 0;
}