POJ1014 多重背包问题

这道题大意是有六种石子，分别价值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;
}