【代码记录】 多重背包的二进制优化与单调队列优化

多重背包问题

有 N 种物品和一个容量是 V 的背包。

第 i 种物品最多有 si 件，每件体积是 vi，价值是 wi。

求解将哪些物品装入背包，可使物品体积总和不超过背包容量，且价值总和最大。求最大价值

二进制优化(例题：QAQ)

简单思路；将每件物品的数量分情况讨论，如果物品数量乘物品的体积大于背包的体积那么就可以看做一个完全背包，否则将第i种物品分成若干件物品，其中每件物品有一个系数，这件物品的费用和价值均是原来的费用和价值乘以这个系数。使这些系数分别为1,2,4,...,2^(k-1),n[i]-2^k+1，且k是满足n[i]-2^k+1>0的最大整数。

代码：

#include<bits/stdc++.h>
using namespace std;
const int M=2e4;
int number[M],weight[M],value[M];
int bag;
int dp[M+5];
void zeronepack(int weight,int value)
{
	for(int j=bag;j>=weight;j--)
	  dp[j]=max(dp[j],dp[j-weight]+value);
}
void completepack(int weight,int value)
{
	for(int j=weight;j<=bag;j++)
	 dp[j]=max(dp[j],dp[j-weight]+value);
}
void multiplepack(int weight,int value,int number)
{
	if(number*weight>=bag)
	{
	 	completepack(weight,value);
        return ; 
	}
	else
	{
		int k=1;
		while(k<=number)
		{
			zeronepack(k*weight,k*value);
			number-=k;
			k=k*2;
		}
		zeronepack(number*weight,number*value);
	}
}
int main()
{
	int n;
	cin>>n>>bag;
	for(int i=1;i<=n;i++)
	 cin>>weight[i]>>value[i]>>number[i];
	for(int i=1;i<=n;i++)
	  multiplepack(weight[i],value[i],number[i]);
	cout<<dp[bag]<<endl;
	return 0;
}

单调队列优化：

例题：0V0

在面对数据范围较大时（20000左右）二进制优化仍然会超时，这时候就需要用单调队列优化

简单思路：

   多重背包的转移方程;

         dp[i][j]=max{dp [ i - 1 ] [ j - k * v [ i ] ] + k * w [ i ] }  0<=k<=c[i]， k是物品数量。根据v[i]将（每组）背包数量分为v[i]组。每组从0到v[i]-1。

令b=j % v[i] , x=j / v[i];

则 j = b + x * v [i]

因为每组之间相互独立所以状态转移方程可以改为

   dp[i][ b + x * v [ i ] = max{ dp[ i - 1 ] [ b + ( x -k ) * v[ i ] ] + k * v[i] }

令 t= x-k;

 则原方程转为：dp[i][b+x*[vi]]=max{dp[i-1][b+t*v[i]]+k*v[i]}

代码：

#include<bits/stdc++.h>
using namespace std;
const int M=20010;
int n,m;
int dp[M],g[M],q[M];
int main()
{
   cin>>n>>m;
   ios::sync_with_stdio(false);
   cin.tie(0);
   cout.tie(0);
   for(int i=1;i<=n;i++)
   {
   	    int weight,value,number;
   	    cin>>weight>>value>>number;
   	    memcpy(g,dp,sizeof dp);
   	    //复制 dp到g 
		for(int j=0;j<weight;j++)
   	    {//枚举m%weight的余数 
   	        int hh=0,tt=-1;
			for(int k=j;k<=m;k+=weight)
			{//枚举所有的加上余数里的数 weight+j,2*weight+j..... 
			  	dp[k]=g[k];
			  	if(hh<=tt&&k-number*weight>q[hh]) hh++;
			  	if(hh<=tt) dp[k]=max(dp[k],g[q[hh]]+(k-q[hh])/weight*value); 
			  	while(hh<=tt&&g[q[tt]]-(q[tt]-j)/weight*value<=g[k]-(k-j)/weight*value) tt--;
			    q[++tt]=k;
			 //单调队列 
			}	
		}
   }
   cout<<dp[m]<<"\n";
   return 0;
}