hdu_1059

/*
 * 使用二进制优化,普通多重背包超时
 * The use of binary optimization, ordinary multiple backpack timeout
*/
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

#define MAXTYPE		7
#define MAXN		120001

int dp[MAXN], cnt[MAXTYPE];

void complete_knapsack(int max_v, int v)
{
	for(int i = 0; i <= max_v-v; i ++) {
		if( dp[i] ) {
			dp[i+v] = 1;
		}
	}
}

void zero_one_knapsack(int max_v, int v)
{
	for(int i = max_v; i >= v; i --) {
		if( dp[i-v] ) {
			dp[i] = 1;
		}
	}
}

int can_divided(int max_v)
{
	int binary, tot;
	memset(dp, 0, sizeof(dp));	dp[0] = 1;
	for(int i = 1; i < MAXTYPE; i ++) {
		if( cnt[i]*i >= max_v ) {
			complete_knapsack(max_v>>1, i);
		}
		else {
			binary = 1; tot = cnt[i];
			while( binary < tot ) {
				zero_one_knapsack(max_v>>1, binary*i);
				tot -= binary;
				binary <<= 1;
			}
			zero_one_knapsack(max_v, tot*i);
		}
	}
	for(int v = max_v>>1; v >= 0; v --) {
		if( dp[v] && 2*v == max_v ) {
			return 1;
		}
	}
	return 0;
}

int main(int argc, char const *argv[])
{
#ifndef ONLINE_JUDGE
	freopen("test.in", "r", stdin);
#endif
	int cas(1), max_v;
	while( true ) {
		max_v = 0;
		for(int i = 1; i < MAXTYPE; i ++) {
			scanf("%d", &cnt[i]);
			max_v += i*cnt[i];
		}
		if( !max_v ) {
			break;
		}
		printf("Collection #%d:\n", cas ++);
		if( can_divided(max_v) ) {
			printf("Can be divided.\n\n");
		}
		else {
			printf("Can't be divided.\n\n");
		}
	}
	return 0;
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值