4个背包问题 动态规划 折半枚举

01背包问题

//有n个重量和价值分别为wi,vi的物品，从中选出总重量不超过W的物品，求最大总价值的挑选方案
// 1<=n<=100
// 1<=wi,vi<=100
// 1<=W<=10000
//思路：dp[i][j]从前i个物品开始挑选总重量小于j时，总价值的最大值
//dp[0][j]=0
//dp[i+1][j] = | dp[i][j] (j<w[i])
//             | max(dp[i][j],dp[i][j-w[i]]+v[i])

int dp[MAX_N+1][MAX_W+1];
void solve(){
    //不要忘记初始化
    memset(dp,0,sizeof(dp[0]));
    for(int i=0;i<n;i++){
        for(int j=0;j<=W;j++){
            if(j<w[i])dp[i+1][j]=dp[i][j];
            else dp[i+1][j] = max(dp[i][j],dp[i][j-w[i]]+v[i]);
        }
    }
    printf("%d\n",dp[n][W]);
}

完全背包问题

//有n种重量和价值分别为wi,vi的物品，从中挑选总重量不超过W的物品，求出最大物品总价值。注意：每种物品可以挑选多件
// 1<=n<=100
// 1<=wi,vi<=100
// 1<=W<=10000

//令dp[i][j]=从前i种物品中挑选重量不超过j时总价值最大值。
//dp[0][j]=0
//dp[i+1][j]=max{dp[i][j-k*w[i]]+k*v[i] | 0<=k}
//这样子需要通过一个三重循环，复杂度o(nW^2)，并不是最优
//改进为：dp[i+1][j]=max(dp[i][j],dp[i+1][j-w[i]]+v[i])

void solve(){
    for(int i=0;i<n;i++){
        for(int j=0;j<=W;j++){
            if(j<w[i]){dp[i+1][j]=dp[i][j];}
            else dp[i+1][j]=max(dp[i][j],dp[i+1][j-w[i]]+v[i]);
        }
    }
    printf("%d\n",dp[n][W]);
}

01背包问题升级版(改变限制条件)

//有n个重量和价值分别为wi,vi的物品，从中选出总重量不超过W的物品，求最大总价值的挑选方案
// 1<=n<=100
// 1<=wi<=1e7
// 1<=vi<=100
// 1<=W<=1e9
//之前是针对不同重量直接计算最大的价值，但是因为改变了限制条件，我们可以通过dp针对不同的价值计算最小的重量。
// dp[i+1][j]=前i个物品中挑选出来总价值是j时重量的最小值（不存在用INF表示）
// dp[0][0]=0;dp[0][j]=INF
// dp[i+1][j]=min(dp[i][j],dp[i][j-v[i]]+w[i])

int dp[MAX_N+1][MAX_N*MAX_V+1];

void solve(){
    fill(dp[0],dp[0]+MAX_N*MAX_V+1,INF);
    dp[0][0]=0;
    for(int i=0;i<n;i++){
        for(int j=0;j<=MAX_N*MAX_V;j++){
            if(j<v[i]){
                dp[i+1][j]=dp[i][j];
            }else{
                dp[i+1][j]=min(dp[i][j],dp[i][j-v[i]]+w[i]);
            }
        }
    }
    int res=0;
    for(int i=0;i<=MAX_N*MAX_V;i++)if(dp[n][i]<=W)res=i;
    printf("%d\n",res);
}

超大背包问题（再次升级限制条件）

//有n个重量和价值分别为wi,vi的物品，从中选出总重量不超过W的物品，求最大总价值的挑选方案
// 1<=n<=40
// 1<=wi<=1e15
// 1<=vi<=1e15
// 1<=W<=1e15

//这种情况下是不能用动态规划了。可以使用折半枚举法。挑选方法一共有2^n种，拆成两半。前半部分选取重量价值对应总和记为w1,v1，后半部分总重量w2<=W-w1时v2最大就完事了。

typedef long long ll;

//输入
int n;
ll w[MAX_N],v[MAX_N];
ll W;

pair<ll,ll> ps[1<<(MAX_N/2)]; //(重量，价值)

void solve(){
    //枚举前半部分
    int n2=n>>1;
    for(int i=0;i<1<<n2;i++){
        ll sw=0,sv=0;
        for(int j=0;j<n2;j++){
            if(i>>j &1){
                sw+=w[j];
                sv+=v[j];
            }
        }
    }
    //去除多余元素(w1[i]<=w2[j]&&v2[i]>v2[j]的元素)
    sort(ps,ps+(1<<n2));
    int m=1;
    for(int i=1;i<1<<n2;i++){
        if(ps[m-1].second<ps[i].second)ps[m++]=ps[i];
    }

    //枚举后半部分
    ll res=0;
    for(int i=0;j< 1<<(n-n2);i++){
        ll sw=0,sv=0;
        for(int j=0;j<n-n2;j++){
            if(i>>j & 1){
                sw+=w[n2+j];
                sv+=v[n2+j];
            }
        }
        if(sw<=W){
            ll tv=(lower_bound(ps,ps+m,make_pair(W-sw,INF))-1)->second;
            res=max(res,sv+tv);
        }
        printf("%lld\n",res);
    }
}