poj3260 （多重背包+完全背包）

买家的硬币是有限制的，是多重背包。
卖家每种硬币没有限制，是完全背包。
dp_buy[i]表示买家凑齐i需要的最少硬币个数。
dp_sale[i]表示专家凑齐i需要的最少硬币个数。

源码

#include<iostream>
#include<cstdio>
#define inf 1<<25
using namespace std;
int n,m,cnt[105],val[105];
int dp_buy[20005],dp_sale[20005];
void completepack(int *dp,int cost,int limit){
for(int i=cost;i<=limit;i++)
dp[i]=min(dp[i],dp[i-cost]+1);
}
void zeroonepack(int *dp,int cost,int limit,int count){
for(int i=limit;i>=cost;i--)
dp[i]=min(dp[i],dp[i-cost]+count);
}
void doublepack(int *dp,int cost,int count,int limit){
if(cost*count>=limit)
completepack(dp,cost,limit);
else{
int tmp=1;
while(tmp<=count){
zeroonepack(dp,cost*tmp,limit,tmp);
count-=tmp;
tmp<<=1;
}
if(count)
zeroonepack(dp,cost*count,limit,count);
}
}
int main(){
while(scanf("%d%d",&n,&m)!=EOF){
for(int i=0;i<n;i++)
scanf("%d",&val[i]);
for(int i=0;i<n;i++)
scanf("%d",&cnt[i]);
for(int i=1;i<=m+10000;i++)
dp_sale[i]=inf;
for(int i=1;i<=m+10000;i++)
dp_buy[i]=inf;
dp_buy[0]=dp_sale[0]=0;
for(int i=0;i<n;i++)
completepack(dp_sale,val[i],m+10000);
for(int i=0;i<n;i++)
doublepack(dp_buy,val[i],cnt[i],m+10000);
int ans=inf;
for(int i=m;i<=m+10000;i++)
if(dp_buy[i]+dp_sale[i-m]<ans)
ans=dp_buy[i]+dp_sale[i-m];
if(ans==inf)
printf("-1\n");
else
printf("%d\n",ans);
}
return 0;
}