0-1背包刚接触的时候感觉很好理解,毕竟只有取或不取两种情况而已,但是后来终于开始接受0-1背包的其他解法,这里也用到了一些其他的算法,好好整理一下咯问题
问题描述:
有N件物品和一个容量为V的背包。第i件物品的费用是c[i],价值是w[i]。求解将哪些物品装入背包可使这些物品的费用总和不超过背包容量,且价值总和最大。
这个问题的特点是:每种物品只有一件,可以选择放或者不放。
dp解法
第一种就是dp
子问题:f[i][v]为放入前i件物品恰好放入体积为v的背包所获得的最大值;
dp状态转移方程为:f[i][v]=max(f[i-1][v] , f[i-1][v-c[i]]+w[i]);方程解释:
是否放入第i件物品,1)如果不放的话,和前i-1件物品放入体积为v的背包所获得的最大值是一样的,则需要求出前i-1件物品恰好放入体积为v的背包所获得的最大值 ; 2)如果放的话,它将占用c[i]体积,背包的价值将增加w[i],现在就会产生一个问题,需要求出前i-1件物品恰好放入体积为v-c[i]的背包所获得的最大值。f[i][v]也就是以上两种想法的最大值。优化空间复杂度
以上方法的空间和时间复杂度都是f(n*v),时间复杂度已经不能再小了,但是空间复杂度或许还可以减小。
其实f[i][v]只是和i-1时有关
参考: [http://www.cnblogs.com/fly1988happy/archive/2011/12/13/2285377.html]
#include<iostream>
#include<string.h>
using namespace std;
int dp[1000][1000];
int main()
{
int t,n,v,i,j;
int va[1000],vo[1000];
cin>>t;
while(t--)
{
cin>>n>>v;
for(i=1;i<=n;i++)
cin>>va[i];
for(i=1;i<=n;i++)
cin>>vo[i];
memset(dp,0,sizeof(dp));
for(i=1;i<=n;i++)
{
for(j=0;j<=v;j++)
{
if(vo[i]<=j)
dp[i][j]=max(dp[i-1][j],dp[i-1][j-vo[i]]+va[i]);
else
dp[i][j]=dp[i-1][j];
}
}
cout<<dp[n][v]<<endl;
}
return 0;
}
2.回溯
0-1背包问题的解空间可用子集树表示,在搜索解空间树时,只要其左儿子结点是一个可行结点,搜索就进入其左子树。当右子树中有可能包含最优解时才进入右子树搜索,否则将右子树剪去。设r是当前剩余物品价值总和,cp是当前价值,bestp是当前最优价值。当cp+r<=bestp时,可剪去右子树。计算右子树中解的上届的更好方法时将剩余物品依其单位质量价值排序,然后依次装入物品,直至装不下时,再装入该物品的一部分而装满背包。由此得到的价值时右子树中解的上界。
#include<iostream>
#include<algorithm>
using namespace std;
int n,C;
struct node1
{
double weight;
double value;
};
node1 a[105];
int cv,cw,bestv=0;
bool cmp(node1 a,node1 b)//按单位质量价值排序
{
return a.value/a.weight>b.value/b.weight;
}
int bound(int i)//计算右子树上界
{
int cleft=C-cw;//剩余容量
int b=cv;
while(i<=n&&a[i].weight<=cleft)
{
cleft-=a[i].weight;
b+=a[i].value;
i++;
}
if(i<=n) b+=a[i].value*cleft/a[i].weight;//若有一个物品不能整一个都放进去,就计算填满背包时最大价值
return b;
}
void backpack(int i)
{
if(i>n&&cv>bestv) //已经找到最佳方案
{
bestv=cv;
return ;
}
if(cw+a[i].weight<=C) //若可以进入左子树
{
cv+=a[i].value;
cw+=a[i].weight;
backpack(i+1);
cv-=a[i].value;
cw-=a[i].weight;
}
if(bound(i+1)>bestv)//通过右子树的上界判断是否进入右子树,否则剪去
backpack(i+1);
}
int main()
{
cin>>n>>C;
for(int i=0;i<n;i++)
{
cin>>a[i].weight>>a[i].value;
}
sort(a,a+n,cmp);
backpack(0);
cout<<bestv<<endl;
return 0;
}