思路来源:
https://www.cnblogs.com/mfrank/p/10533701.html
目录
01背包 模型问题
n 件物品,物品重量 w[i],物品价值 v[i]
挑物品放入背包,背包承受最大重量 W
如何选择装入的物品,使背包中物品总价值最大?
https://www.cnblogs.com/fengziwei/p/7750849.html
1.1 思路1 - 暴力法
把所有的可能性求出来,选一个最大的
2
n
2^n
2n 种可能,太复杂!
1.2 思路2 - 动态规划
假设有 4 件物品
每个物品有 2 种可能:不装入包中(0),或 装入包中(1),这就 0 1 状态
最优性原理:
必定至少有一种选择方案,使得背包物品总价值最大(易知)
为什么是至少:
可能有两个物品的 w 和 v 相等
或有 m 个物品 与 另 n 个物品的 w+v 相等
假定 就想着其中的一个最优选择方案A
那物品就一个一个看,进行选择
若前面 (m-1) 个物品已经选择完毕,现在看第 m 个物品
若 最优方案A中第 m 个物品应该选择
可现在我是不知道该不该选择第 m 个物品的!
最优方案A中,进行至前 m 个物品,装进背包的物品有价值和M1 的
若不选择 m,背包的价值 M2 > M1(前m个物品可能和最优解A的前m个物品选择不一致)
在后面的物品选择中,都按最优方案A来做,
最后我的的选择方案是不是比最优方案A 的还好?
这与最优方案A是最优的 矛盾!
所以最优解的子解仍是子问题的最优解
怎么看子问题的最优解?
若前面 (m-1) 个物品已经选择完毕,和最优方案A是一致的,对应价值M1和体积V1
现在看第 m 个物品
体积V1若含有 物品m的话,就要把已经在背包的物品放弃一部分
对应体积是(V1 - 物品m的体积) 对应价值,加上m的价值,和为A
A和M1比较,选大的数
1.2.1 二维列表
1.2.1.1 案例特殊代码
#include<bits/stdc++.h>
#define max 100
using namespace std;
void solve(int v[], int w[], int r[][max], int m, int n)
{
int i, j;
for(i = 0 ;i <= n; i++)
r[0][i] = 0;
for(i = 0 ;i <= m; i++)
r[i][0] = 0;
for(i = 1; i <= m; i++)
{
for(j = 1; j <= n; j++)
{
if(j < w[i])
r[i][j] = r[i-1][j];
else
{
if( r[i-1][j] > (r[i-1][j-w[i]] + v[i]) )
r[i][j] = r[i-1][j];
else
r[i][j] = r[i-1][j-w[i]] + v[i];
}
}
}
}
int main()
{
int v[max] = {0,2,4,3,7};
int w[max] = {0,2,3,5,5};
int r[max][max];
int m,n;
while(cin>>m>>n)
{
solve(v,w,r,m,n);
cout << r[m][n] <<endl;
memset(v, 0, sizeof(v));
memset(w, 0, sizeof(w));
}
return 0;
}
1.2.1.2 随机数据代码
#include<bits/stdc++.h>
#define max 100
#define num 100
using namespace std;
void solve(int v[], int w[], int r[][max], int m, int n)
{
int i, j;
for(i = 0 ;i <= n; i++)
r[0][i] = 0;
for(i = 0 ;i <= m; i++)
r[i][0] = 0;
for(i = 1; i <= m; i++)
{
for(j = 1; j <= n; j++)
{
if(j < w[i])
r[i][j] = r[i-1][j];
else
{
if( r[i-1][j] > (r[i-1][j-w[i]] + v[i]) )
r[i][j] = r[i-1][j];
else
r[i][j] = r[i-1][j-w[i]] + v[i];
}
}
}
}
void myrand(int v[], int w[], int m)
{
int i;
v[0] = w[0] = 0;
for(i = 1 ; i <= m;i++)
{
v[i] = rand()%num;
w[i] = rand()%num;
}
}
int main()
{
srand(time(NULL));
int v[max], w[max], r[max][max];
int m,n;
while(cin>>m>>n)
{
myrand(v,w,m);
solve(v,w,r,m,n);
cout << r[m][n] <<endl;
memset(v, 0, sizeof(v));
memset(w, 0, sizeof(w));
}
return 0;
}
1.2.2 递归
1.2.2.1 案例特殊代码
#include<bits/stdc++.h>
#define max 100
using namespace std;
int solve(int v[], int w[], int m, int n)
{
int r = 0;
if(m==0 || n == 0)
r = 0;
else if(w[m] > n)
r = solve(v,w,m-1,n);
else
{
int t1 = solve(v,w,m-1,n);
int t2 = solve(v,w,m-1,n-w[m]) +v[m];
r = (t1>t2)?t1:t2;
}
return r;
}
int main()
{
int v[max] = {0,2,4,3,7};
int w[max] = {0,2,3,5,5};
int m,n;
while(cin>>m>>n)
{
cout << solve(v,w,m,n) <<endl;
memset(v, 0, sizeof(v));
memset(w, 0, sizeof(w));
}
return 0;
}
1.2.2.2 随机数据代码
#include<bits/stdc++.h>
#define max 100
#define num 100
using namespace std;
int solve(int v[], int w[], int m, int n)
{
int r = 0;
if(m==0 || n == 0)
r = 0;
else if(w[m] > n)
r = solve(v,w,m-1,n);
else
{
int t1 = solve(v,w,m-1,n);
int t2 = solve(v,w,m-1,n-w[m]) +v[m];
r = (t1>t2)?t1:t2;
}
return r;
}
void myrand(int v[], int w[], int m)
{
int i;
v[0] = w[0] = 0;
for(i = 1 ; i <= m;i++)
{
v[i] = rand()%num;
w[i] = rand()%num;
}
}
int main()
{
srand(time(NULL));
int v[max];
int w[max];
int m,n;
while(cin>>m>>n)
{
myrand(v,w,m);
cout << solve(v,w,m,n) <<endl;
memset(v, 0, sizeof(v));
memset(w, 0, sizeof(w));
}
return 0;
}
1.3 思路3 - 贪心法
排序!
在重量一定的情况下,价值最大 = 单位重量的价值最大
所以按照 w/v 由高到低排序
单位体积的价值高的 在最优方案里面
知道背包装满
1.3.1 特殊案例代码
#include<bits/stdc++.h>
#define max 100
#define num 100
using namespace std;
typedef struct
{
int v,w;
double p;
}vw;
int solve(vw s[], int m, int n)
{
int r = 0;
int w = 0;
int i = 0;
while(w<n && i < m)
{
r += s[i].v;
w += s[i].w;
i++;
}
return r;
}
int com(vw a, vw b)
{
return a.p > b.p;
}
int main()
{
vw s[max];
s[0].v = 2;s[1].v = 4;s[2].v = 3;s[3].v = 7;
s[0].w = 2;s[1].w = 3;s[2].w = 5;s[3].w = 5;
int m,n, i;
while(cin>>m>>n)
{
for(i = 0; i < m ; i++)
s[i].p = (1.0*s[i].v)/s[i].w;
sort(s,s+m,com);
cout << solve(s,m,n) <<endl;
memset(s, 0, sizeof(s));
}
return 0;
}
1.3.2 随机数据代码
#include<bits/stdc++.h>
#define max 100
#define num 100
using namespace std;
typedef struct
{
int v,w;
double p;
}vw;
int solve(vw s[], int m, int n)
{
int r = 0;
int w = 0;
int i = 0;
while(w<n && i < m)
{
r += s[i].w;
w += s[i].v;
i++;
}
return r;
}
int com(vw a, vw b)
{
return a.p > b.p;
}
void myrand(vw s[], int m)
{
int i;
for(i = 0 ; i < m;i++)
{
s[i].v = rand()%num;
s[i].w = rand()%num;
s[i].p = (1.0*s[i].v)/s[i].w;
}
sort(s,s+m,com);
}
int main()
{
srand(time(NULL));
vw s[max];
int m,n;
while(cin>>m>>n)
{
myrand(s, m);
cout << solve(s,m,n) <<endl;
memset(s, 0, sizeof(s));
}
return 0;
}
背包问题
完全背包
01背包的基础上,背个物品有无限个
牛客题
https://www.nowcoder.com/practice/f9c6f980eeec43ef85be20755ddbeaf4?tab=answerKey
题目1
王强有 N 元的年终奖。
购物,想买的物品分两类:主件与附件,附件从属于某个主件
下表就是一些主件与附件的例子:
主件 附件
电脑 打印机,扫描仪
书柜 图书
书桌 台灯,文具
工作椅 无
如果要买附件,必须先买该附件所属的主件。
每个主件可以有 0 个、 1 个或 2 个附件。
附件没有从属于自己的附件。
王强想买的东西很多,为了不超出预算,他把每件物品规定了一个重要度,分为 5 等:用整数 1 ~ 5 表示,第 5 等最重要。
每件物品的价格(都是 10 元的整数倍)。
他希望在不超过 N 元(可以等于 N 元)的前提下,使每件物品的价格与重要度的乘积的总和最大。
设第 j 件物品的价格为 v[j] ,重要度为 w[j] ,共选中了 k 件物品
编号依次为 j1 , j2 ,……, jk ,则所求的总和为:
v[ j1 ]*w[ j1 ]+v[ j2 ]*w[ j2 ]+ … +v[ jk ]*w[ jk ] 。(其中 *
为乘号)
请你帮助王强设计一个满足要求的购物单。
输入描述
输入的第 1 行,为两个正整数,用一个空格隔开:N m
( N ( <32000 )表示总钱数, m ( <60 )为希望购买物品的个数)
从第 2 行到第 m+1 行,第 j 行给出了编号为 j-1 的物品的基本数据,每行有 3 个非负整数 v p q
v 表示该物品的价格( v<10000 ), p 表示该物品的重要度( 1 ~ 5 ), q 表示该物品是主件还是附件。如果 q=0 ,表示该物品为主件,如果 q>0 ,表示该物品为附件, q 是所属主件的编号
输出描述
输出文件只有一个正整数,是不超过总钱数的物品的价格与重要度乘积的总和的最大值( <200000 )。
输入
1000 5
800 2 0
400 5 1
300 5 1
400 3 0
500 2 0
输出
2200
#include<bits/stdc++.h>
#define max 61
using namespace std;
int main()
{
unsigned int N, m;
unsigned int v,p,q;
unsigned int w[max][3]={0}, vm[max][3]={0};
int i,j, t;
while(cin>>N>>m)
{
unsigned int dp[max][3201] = {0};
N = N/10;
for(i = 1 ; i <= m; i++)
{
cin>> v >> p >> q;
p = p*v;
v=v/10;
if ( q == 0)
{
w[i][0] = p;
vm[i][0] = v;
}
else if ( vm[q][1] == 0)
{
w[q][1] = p;
vm[q][1]=v;
}
else if ( vm[q][2] == 0)
{
w[q][2] = p;
vm[q][2]=v;
}
}
for( i = 1; i <= m; i++)
{
for (j = 1; j<= N; j++)
{
dp[i][j] = dp[i-1][j];
if(j >= vm[i][0])
{
t = dp[i-1][j] > (dp[i-1][j-vm[i][0]]+w[i][0]) ? dp[i-1][j]:(dp[i-1][j-vm[i][0]]+w[i][0]);
if( dp[i][j] < t )
dp[i][j] = t;
}
if(j >= (vm[i][0] + vm[i][1]) )
{
t = dp[i-1][j-vm[i][0]-vm[i][1]]+w[i][0] + w[i][1];
if( dp[i][j] < t )
dp[i][j] = t;
}
if(j >= (vm[i][0] + vm[i][2]) )
{
t = dp[i-1][j-vm[i][0]-vm[i][2]]+w[i][0] + w[i][2];
if( dp[i][j] < t )
dp[i][j] = t;
}
if(j >= (vm[i][0] + vm[i][1] + vm[i][2]) )
{
t = dp[i-1][j-vm[i][0]-vm[i][1]-vm[i][2]] + w[i][0] + w[i][1] + w[i][2];
if( dp[i][j] < t )
dp[i][j] = t;
}
}
}
cout << dp[m][N] << endl;
}
return 0;