(01背包详解) 顺带几题完全背包问题
(update:)
2019.4.4 初稿。
2019.4.13 重改加上一些注释 顺便加几道完全背包题目。以及调整Markdown。
本文涉及到的题目
(small P1048 采药)
(small P1049 装箱问题)
(small P1060 开心的金明)
(small P1164 小A点菜)
(small P2639 [USACO09OCT]Bessie的体重问题 Bessie's We…)
(small P1794 装备运输_NOI导刊2010提高 （04）)
(small P1877 [HAOI2012]音量调节)
(small P1910 L国的战斗之间谍)
(small P2871 [USACO07DEC]手链Charm Bracelet)
(small P1455 搭配购买)
(small P1616 疯狂的采药)
(small P2722 总分 Score Inflation)
(small P2918 [USACO08NOV]买干草Buying Hay)
前言：DP 快接触半年了。 还是想起来把曾经
虐我的(01背包) 好好写写。
下面借鉴此处
其中(F[i-1][j])表示前(i-1)件物品中选取若干件物品放入剩余空间为(j)的背包中所能得到的最大价值
而(F[i-1][j-C_i]+W_i)表示前(i-1)件物品中选取若干件物品放入剩余空间为(j-C_i)的背包中所能取得的最大价值加上第(i)件物品的价值
根据第 (i) 件物品放或是不放确定遍历到第 (i) 件物品时的状态(F[i][j])
设物品件数为(N) 背包容量为(V) 第 (i)件物品体积为(C_i) 第(i)件物品价值为(W_i)
所以得出代码
for(register int i=1; i<=n; i++)
for(register int j=m; j>=c[i]; j--) dp[i][j] = max(dp[i-1][j-c[i]] + w[i] , dp[i][j]) ;
[F[i][j]只与F[i-1][j]和F[i-1][j-C_i]有关 ]即只和(i-1)时刻状态有关 所以我们只需要用一维数组(F[])来保存(i-1)时的状态(F[])。
假设(i-1)时刻的(F[])为({a[0] a[1] a[2]… a[v]})
那么(i)时刻的(F[])中第(k)个应该为(max(a_k,a[k]-C[i]+W_i))
即 (max(F_k,F[k-C_i]+W_i))
这就需要我们遍历(V)时逆序遍历
这样才能保证求 (i)时刻 (F[k]) 时 (F[k-C_i]) 是 (i-1) 时刻的值
如果正序遍历则当求(F[k])时其前面的(F_0 F_1 F_k-1)都已经改变过，里面存的都不是(i-1)时刻的值，这样求(F_k)时利用(F[K-C_i])必定是错的值。最后(F_v)即为最大价值。
下面的完全背包会仔细解释为什么是倒序遍历（(update on 4.13)）
for(register int i=1; i<=n; i++)
for(register int j=m; j>=c[i]; j--) dp[j] = max(dp[j-c[i]] + w[i] , dp[j]) ;
虽然时间复杂度不变 但是内存减少很多了。
=========================================update：大致模板
#include <bits/stdc++.h>
#define rep(i,j,n) for(register int i=j;i<=n;i++)
#define Rep(i,j,n) for(register int i=j;i>=n;i--)
#define low(x) x&(-x)
using namespace std ;
typedef long long LL ;
const int inf = INT_MAX >> 1 ;
inline LL In() { LL res(0) , f(1) ; register char c ;
#define gc c = getchar()
while(isspace(gc)) ; c == '-' ? f = - 1 , gc : 0 ;
while(res = (res << 1) + (res << 3) + (c & 15) , isdigit(gc)) ;
return res * f ;
#undef gc
}
int n , m ;
const int N = 100000 + 5 ;
int f[N] ; //数组
inline void Ot() {
n = In() , m = In() ;
for(register int i=1;i<=n;i++) {
int w = In() , c = In() ; //输入 体积&&价值
for(register int j=m;j>=w;j--)
f[j] = max(f[j] , f[j-w] + c) ; //状态转移方程
}
cout << f[m] << endl ; //f[m] 是 目标值。 所以输出。
}
signed main() {
// freopen("test.in","r",stdin) ;
return Ot() , 0 ;
}
可以根据这个代码来模拟过程。
=========================================例题。
P1048 采药
P1049 装箱问题
P1060 开心的金明
P1164 小A点菜
P2639 [USACO09OCT]Bessie的体重问题Bessie's We…
=========================================
这题就比较裸了。直接套进去板子。
#include <bits/stdc++.h>
#define rep(i,j,n) for(register int i=j;i<=n;i++)
#define Rep(i,j,n) for(register int i=j;i>=n;i--)
#define low(x) x&(-x)
using namespace std ;
typedef long long LL ;
const int inf = INT_MAX >> 1 ;
inline LL In() { LL res(0) , f(1) ; register char c ;
#define gc c = getchar()
while(isspace(gc)) ; c == '-' ? f = - 1 , gc : 0 ;
while(res = (res << 1) + (res << 3) + (c & 15) , isdigit(gc)) ;
return res * f ;
#undef gc
}
int T , M ;
const int N = 1000 + 5 ;
int dp[N] ;
inline void Ot() {
T = In() , M = In() ;
for(register int i=1;i<=M;i++) {
int x = In() , y = In() ;
for(register int j=T;j>=x;j--) dp[j] = max(dp[j-x]+y, dp[j]) ;
}
cout << dp[T] << endl ;
}
signed main() {
return Ot() , 0 ;
}
这题转移方程略有不同。 重量就是价格（应该能这么理解）
所以 (W[i] = C[i]) 所以转移方程为(dp[j] = max ( dp[j] , dp[j-w[i]] + w[i]);)
#include<bits/stdc++.h>
#define mx 25000
using namespace std;
int w[mx],n,V,dp[mx];
int main() {
ios::sync_with_stdio(false);
cin>>V>>n;
for(int i=1; i<=n; i++) {
cin>>w[i];
}
for(int i=1; i<=n; i++) {
for(int j=V; j>=w[i]; j--) {
dp[j]=max(dp[j],dp[j-w[i]]+w[i]);
}
}
int ans=0;
for(int i=1; i<=V; i++) {
ans=max(ans,dp[i]);
}
cout<<V-ans<<'
';
return 0;
}
这题就是把模板套一下。 价值先预处理出来。
#include<bits/stdc++.h>
#define f(i,j,n) for(int i=j;i<=n;i++)
#define fa(i,j,n) for(int i=j;i>=n;i--)
using namespace std;
int w[30],v[30],f[50000];
int n,m;
int main() {
cin>>m>>n;
f(i,1,n) {
cin>>v[i]>>w[i];
w[i]*=v[i];
}
f(i,1,n)
fa(j,m,v[i]) f[j]=max(f[j],f[j-v[i]]+w[i]);
cout<<f[m]<<endl;
return 0;
}
01背包的变式 (f_j = f_j + f[j-w_i];)
#include <bits/stdc++.h>
#define rep(i,j,n) for(register int i=j;i<=n;i++)
#define Rep(i,j,n) for(register int i=j;i>=n;i--)
#define low(x) x&(-x)
using namespace std ;
typedef long long LL ;
const int inf = INT_MAX >> 1 ;
inline LL In() { LL res(0) , f(1) ; register char c ;
#define gc c = getchar()
while(isspace(gc)) ; c == '-' ? f = - 1 , gc : 0 ;
while(res = (res << 1) + (res << 3) + (c & 15) , isdigit(gc)) ;
return res * f ;
#undef gc
}
int n , m ;
const int N = 100 + 5 ;
const int M = 10000 + 5 ;
int w[N] ;
int f[M] ;
inline void Ot() {
n = In() , m = In() ;
for(register int i=1;i<=n;i++) w[i] = In() ;
f[0] = 1 ;
for(register int i=1;i<=n;i++)
for(register int j=m;j>=w[i];j--) f[j] += f[j-w[i]] ;
cout << f[m] << endl ;
}
signed main() {
// freopen("test.in","r",stdin) ;
return Ot() , 0 ;
}
重量就是价格
所以 (W_i = C_i) 所以转移方程为(dp_j = max ( dp_j , dp[j-w_i] + w_i);)
但是这里赋值过了 所以就没什么关系了 直接跑(01)背包。
#include<iostream>
using namespace std;
int f[450001],w[100001],c[100001],n,m;
int main() {
ios::sync_with_stdio(false);
cin>>m>>n;
for(int i=1; i<=n; i++) {
cin>>c[i];
w[i]=c[i];
}
for(int i=1; i<=n; i++) {
for(int j=m; j>=c[i]; j--) {
if(f[j-c[i]]+w[i]>f[j])
f[j]=f[j-c[i]]+w[i]; //这里其实应该是max 当时就这么写了。
}
}
cout<<f[m];
return 0;
}
=========================================二维例题。
P1794 装备运输_NOI导刊2010提高（04）
P1877 [HAOI2012]音量调节
P1910 L国的战斗之间谍
P2871 [USACO07DEC]手链Charm Bracelet
=========================================
同上面说明
#include <bits/stdc++.h>
#define rep(i,j,n) for(register int i=j;i<=n;i++)
#define Rep(i,j,n) for(register int i=j;i>=n;i--)
#define low(x) x&(-x)
using namespace std ;
typedef long long LL ;
const int inf = INT_MAX >> 1 ;
inline LL In() { LL res(0) , f(1) ; register char c ;
#define gc c = getchar()
while(isspace(gc)) ; c == '-' ? f = - 1 , gc : 0 ;
while(res = (res << 1) + (res << 3) + (c & 15) , isdigit(gc)) ;
return res * f ;
#undef gc
}
int v , g ;
int n ;
const int N = 500 ;
int dp[N][N] ;
inline void Ot() {
v = In() ;
g = In() ;
n = In() ;
for(register int i=1;i<=n;i++) {
int c = In() ;
int x = In() , y = In() ;
for(register int j=v;j>=x;j--)
for(register int k=g;k>=y;k--) dp[j][k] = max(dp[j][k] , dp[j-x][k-y] + c) ;
}
cout << dp[v][g] << endl ;
}
signed main() {
// freopen("testdata.txt","w",stdout) ;
return Ot() , 0 ;
}
这题在某谷上写过(题解) 了。
#include <bits/stdc++.h>
#define rep(i,j,n) for(register int i=j;i<=n;i++)
using namespace std;
typedef long long LL;
inline LL read(){ LL x=0;int f(1);char ch=getchar();
while(!isdigit(ch)) { if(ch=='-') f=-1;
ch=getchar();
}
while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar(); return x*f;
}
int n;
int Begin,Max;
int a[1<<6];
int dp[1<<6][1<<12];
signed main(){
memset(dp,0,sizeof(dp));
n=read(); Begin=read(); Max=read();
dp[0][Begin]=1;
rep(i,1,n) a[i]=read();
rep(i,1,n) rep(j,0,Max) {
if (j+a[i] <= Max) dp[i][j]=dp[i][j]||dp[i-1][j+a[i]];
if (j-a[i] >= 0) dp[i][j]=dp[i][j]||dp[i-1][j-a[i]];
}
LL ans = -0x7f;
for(register int i=1;i<=Max;i++) {
if(dp[n][i]) ans = i;
}
if ( ans != -0x7f) cout << ans << endl ;
else puts("-1") ;
return 0;
}
根据上面这部分。 把几个变量整合到一块。
#include <bits/stdc++.h>
#define rep(i,j,n) for(register int i=j;i<=n;i++)
#define Rep(i,j,n) for(register int i=j;i>=n;i--)
#define low(x) x&(-x)
using namespace std ;
typedef long long LL ;
const int inf = INT_MAX >> 1 ;
inline LL In() { LL res(0) , f(1) ; register char c ;
#define gc c = getchar()
while(isspace(gc)) ; c == '-' ? f = - 1 , gc : 0 ;
while(res = (res << 1) + (res << 3) + (c & 15) , isdigit(gc)) ;
return res * f ;
#undef gc
}
int v , m ;
int n ;
const int N = 1000 + 5 ;
int dp[N][N] ;
inline void Ot() {
n = In() ;
v = In() , m = In() ;
rep(i,1,n) {
int x = In() , y = In() , z = In() ;
Rep(i,v,y)
Rep(j,m,z) dp[i][j] = max(dp[i-y][j-z] + x , dp[i][j]) ;
}
cout << dp[v][m] << endl ;
}
signed main() {
return Ot() , 0 ;
}
#include <bits/stdc++.h>
#define rep(i,j,n) for(register int i=j;i<=n;i++)
#define Rep(i,j,n) for(register int i=j;i>=n;i--)
#define low(x) x&(-x)
using namespace std ;
typedef long long LL ;
const int inf = INT_MAX >> 1 ;
inline LL In() { LL res(0) , f(1) ; register char c ;
#define gc c = getchar()
while(isspace(gc)) ; c == '-' ? f = - 1 , gc : 0 ;
while(res = (res << 1) + (res << 3) + (c & 15) , isdigit(gc)) ;
return res * f ;
#undef gc
}
int v , m ;
int n ;
const int N = 400 + 5 ;
int dp[N][N] ;
inline void Ot() {
v = In() , m = In() ;
n = In() ;
rep(i,1,n) {
int x = In() , y = In() , z = In() ;
Rep(i,v,x)
Rep(j,m,y) dp[i][j] = max(dp[i-x][j-y] + z , dp[i][j]) ;
}
cout << dp[v][m] << endl ;
}
signed main() {
return Ot() , 0 ;
}
带点[并查集]的
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
inline LL read () { LL res = 0 ;int f (1) ;char ch = getchar ();
while (!isdigit(ch)) { if (ch == '-') f = -1 ;ch = getchar();}
while (isdigit(ch)) res = (res << 1) + (res << 3) + (ch ^ 48) ,ch = getchar(); return res * f ;
}
int n,m,we,fa[10005],w[10005],c[10005],dp[10005];
inline int find(int x){
if(fa[x]!=x) fa[x]=find(fa[x]);
return fa[x];
}
inline void merge(int x,int y){
int f1=find(x),f2=find(y);
if(f1!=f2)fa[f1]=f2,c[f2]+=c[f1],w[f2]+=w[f1];
}
int main(){
n=read(),m=read(),we=read();
for(register int i=1;i<=n;i++) fa[i]=i;
for(register int i=1;i<=n;i++) c[i]=read(),w[i]=read();
for(register int i=1;i<=m;i++){
int u=read(),v=read();
merge(u,v);
}
for(register int i=1;i<=n;i++) if(fa[i]==i)
for(register int j=we;j>=c[i];j--) dp[j]=max(dp[j],dp[j-c[i]]+w[i]);
cout << dp[we] << endl ;
return 0;
}
HDU - 2026 - Bone Collector
HDU - 2546 - 饭卡
HDU - 2955 - Robberies
HDU - 1203 - I NEED A OFFER!
HDU - 1171 - Big Event in HDU
====================================================完全背包。（(update on 4.13)）
那么我们再次仔细思考。 为什么01背包要反过来？
就假设是正过来的好了吖。
正序遍历 是 (for(register int i = w[i];i<=m;i++))
根据循环条件 我们可以得知 (w_i) 可以放 (inf) 次（假设(m=∞)）
这样子的话。就是完全背包。
(∵01背包一件物品只能放一次)
(∴就只能更新一次值 所以是 from w_i to m)
=========================================例题。
P1616 疯狂的采药
P2722 总分 Score Inflation
(P2918 [USACO08NOV]买干草Buying Hay)
=========================================
内循环倒过来。 因为这是一个完全背包（也是01背包的变式
在 <采药> 中，每种草药只允许 采一次 。（所以是标准的01背包 上面也有了。
我们将采 第(i)种草药 所需的 时间 设为 (t_i) 价值 设为 (p_i)
如果有一个数组 (f[i][j]) 来表示 从前往后 到第(i)种草药 (当然前面可能有草药不采) ，花费了最多 (j) 时间 (意思是可能花费了少于(j)时间) 时 能采到草药的 最大价值
到第(i)种草药时有两种情况：采与不采。
若不采这种草药，则 时间花费没有增多 ，经过的 草药种数增加了(1) ， 采到草药价格不变 ，所以 (f[i][j]=f[i-1][j]) ；
若采这种草药 则 时间花费增加了(t_i) 种数增加(1) 采到草药价格增加了(p_i)，所以 (f[i][j]=f[i-1][j-t_i]+p_i) 。
我们当然要使 (f[i][j])尽可能大 即有 (f[i][j]=max(f[i-1][j],f[i-1][j-t_i]+p_i))
#include <bits/stdc++.h>
#define rep(i,j,n) for(register int i=j;i<=n;i++)
#define Rep(i,j,n) for(register int i=j;i>=n;i--)
#define low(x) x&(-x)
using namespace std ;
typedef long long LL ;
const int inf = INT_MAX >> 1 ;
inline LL In() { LL res(0) , f(1) ; register char c ;
#define gc c = getchar()
while(isspace(gc)) ; c == '-' ? f = - 1 , gc : 0 ;
while(res = (res << 1) + (res << 3) + (c & 15) , isdigit(gc)) ;
return res * f ;
#undef gc
}
int t , m ;
const int T = 100000 + 5 ;
LL f[T] ;
inline void Ot() {
t = In() , m = In() ;
rep(i,1,m) {
int x = In() , y = In() ;
rep(j,x,t) f[j] = max(f[j] , f[j-x] + y) ;
}
cout << f[t] << endl ;
}
signed main() {
// freopen("test.in","r",stdin) ;
return Ot() , 0 ;
}
这和P1616 疯狂的采药 基本是一样的
都是完全背包。
#include<bits/stdc++.h>
#define f(i,j,n) for(int i=j;i<=n;i++)
using namespace std;
int V, n;
int a[10001],b[10001],f[10001];
void read(int &x) {
int f=1;
x=0;
char s=getchar();
while(s<'0' or s>'9') {
if(s=='-') f=-1;
s=getchar();
}
while(s>='0' and s<='9') {
x=x*10+s-'0';
s=getchar();
}
x*=f;
}
int main() {
read(V),read(n);
f(i,1,n) read(a[i]),read(b[i]);
f(j,1,n)
f(k,b[j],V) f[k]=max(f[k],f[k-b[j]]+a[j]);
cout<<f[V]<<endl;
return 0;
}
这题就有点坑了。
坑点1 ：要用最小值。所以初值赋值最大值。 但是 (dp_0 = 0) 不然如何都是 最大值
坑点2 ：目标不是(dp_h) 而是 (max) (dp_h...dp_{h+5000}) 因为
第(i)公司卖的干草包重量 为(P_i (1<=P_i<=5000))磅，需要的开销为(C_i (1<C_i <=5000))美元
没准 多一丢丢还便宜点。（大雾。
(code)
#include <bits/stdc++.h>
#define rep(i,j,n) for(register int i=j;i<=n;i++)
#define Rep(i,j,n) for(register int i=j;i>=n;i--)
#define low(x) x&(-x)
using namespace std ;
typedef long long LL ;
const int inf = INT_MAX >> 1 ;
inline LL In() { LL res(0) , f(1) ; register char c ;
#define gc c = getchar()
while(isspace(gc)) ; c == '-' ? f = - 1 , gc : 0 ;
while(res = (res << 1) + (res << 3) + (c & 15) , isdigit(gc)) ;
return res * f ;
#undef gc
}
int n , h ;
const int N = 55000 + 5 ;
int dp[N] ;
inline void Ot() {
n = In() , h = In() ;
memset(dp,0x7f,sizeof(dp)) ;
dp[0] = 0 ;
rep(i,1,n) {
int x = In() , y = In() ;
rep(j,x,h+5000) dp[j] = min(dp[j] , dp[j-x] + y) ;
}
int ans = 0x7f7f7f7f7f ;
rep(i,h,h+5000) ans = min(ans , dp[i]) ;
cout << ans << endl ;
}
signed main() {
// freopen("test.in","r",stdin) ;
return Ot() , 0 ;
}