这道题大意是有六种石子,分别价值1,2,3,4,5,6,每个可能有m个,m<=20000。问能否他们分成价值相等的两部分。
思路:多重背包,总价值为奇数直接判定不可分,总价值为偶数时,令背包容量为总价值的一半,另每个石子的cost = 价值,恰好填满说明有解。这里涉及到一个状态初始化的技巧,《背包九讲》里面有说过:比如F[]为记录总价值的数组,那么如果只要恰好填满的最优解,就另F[0] = 0,其他F[1]到F[V]都为负无穷。如果不讨论是否恰好填满,只要价值最大的解,就都初始化为0。理解起来就是:F[v]的初始化的值,就代表容量为v的背包不装东西且“恰好填满”时的最大价值,那么当然只有容量为0的背包可以再不装东西的情况下“恰好填满”。
P.S. 最近在读《背包九讲》,读完可能写一些笔记。
P.S. 这题的函数可以做背包问题的模板。
如有错误请指出,谢谢~
#include <iostream>
#include <stdio.h>
#define inf 100000000
using namespace std;
int f[220000];
int vol; // 背包容量 为弹珠总价值的一半
int max (int a, int b)
{
return a>b ? a : b;
}
// 处理一件01背包问题的物品 a为记录总价值的数组 c为物品大小 w为物品价值 v为背包容量
void ZeroOnePack (int *a, int c, int w, int v)
{
for(int i = v; i >= c; i--)
a[i] = max(a[i], a[i - c] + w);
}
//处理一件完全背包问题的物品(物品无限)
void CompletePack (int *a, int c, int w, int v)
{
for(int i = c; i <= v; i++)
a[i] = max(a[i], a[i-c] + w);
}
//处理一件多重背包问题的物品 m为该物品数量
void MultiPack (int *a, int c, int w, int v, int m)
{
//如果该物品足够塞满背包 说明对于该物品来说是个完全背包问题
if( m*c >= v)
{
CompletePack(a, c, w, v);
return;
}
//不能塞满 用二进制拆分 转化为01背包问题
int k = 1;
while(k < m)
{
ZeroOnePack(a, k*c, k*w, v);
m = m - k;
k = k * 2;
}
ZeroOnePack(a, m*c, m*w, v);
}
int main ()
{
int m[7],c[7],w[7];
int cases = 0, sum;
while(scanf("%d%d%d%d%d%d", &m[1], &m[2], &m[3], &m[4], &m[5], &m[6]))
{
cases++;
if(m[1]==0 && m[2]==0 && m[3]==0 && m[4]==0 && m[5]==0 && m[6]==0)
break;
sum = 0;
for(int i = 1; i <= 6; i++)
{
c[i] = w[i] = i;
sum += c[i] * m[i];
}
printf("Collection #%d:\n", cases);
if(sum % 2 == 1)
printf("Can't be divided.\n\n");
else
{
vol = sum / 2;
f[0] = 0; //恰好装满的条件
for(int j = 1; j <= vol; j++)
f[j] = -inf;
for(int k = 1; k <= 6; k++)
MultiPack(f, c[k], w[k], vol, m[k]);
if(f[vol] < 0)
printf("Can't be divided.\n\n");
else
printf("Can be divided.\n\n");
}
}
return 0;
}