HDU 1171 Big Event in HDU【多重背包||母函数】

题目：http://acm.hdu.edu.cn/showproblem.php?pid=1171

题意：将财产尽可能对半分

两种解法：一是背包，二是母函数
            可以以财产的总和的一半为背包上限，求背包最大值。
            因为每件物品可以有很多件，所以本题为多重背包问题.

背包解法（比较快）：

#define max(a,b) ( (a) > (b) ? (a) : (b) )
#include<stdio.h>
#include<string.h>
int dp[250000],v[53],m[53];
int main()
{
    int n,i,j,sum,he,k;
	while( scanf("%d",&n)!=EOF)	
	{
		if( n <= 0 ) break;	
		he = 0;
		memset(dp,0,sizeof(dp));
		for(i=0;i<n;i++)
		{
			scanf("%d%d",&v[i],&m[i]);	
			he += v[i]*m[i];
		}
		sum = he/2;	
		for(i=0;i<n;i++)
		{
			if( v[i] * m[i] > sum)	//完全背包
			{
				for(j=v[i];j<=sum;j++)
					dp[j] = max ( dp[j],dp[ j-v[i] ] + v[i] );
				continue;
			}
			k = 1;
			while(k < m[i])	
			{
				for(j = sum;j>=k*v[i];j--)		//01
					dp[j] = max ( dp[j],dp[ j - k*v[i] ] + k*v[i] );
				m[i] -= k;
				k <<= 1;
			}
			for(j = sum;j >= m[i]*v[i];j--)		//继续01
				dp[j] = max ( dp[j],dp[ j - m[i] * v[i] ] + m[i] * v[i] );
		}
		printf("%d %d\n",he-dp[sum],dp[sum]);	
	}
	return 0;
}

母函数解法（够慢）列出情况后，从总值的一半开始往下找，第1个系数非0就是答案

母函数：

(1+x^v[0]+x^2v[0]+...+x^num[0]*v[0])(1+x^v[1]+x^2v[1]+...+x^num[1]*v[1])...(1+x^v[n]+ x ^2v[n] +x ^3v[n] +...+x^v[n]*num[n])

#define MAX 255556	//杭电题目上在数组大小上总是个坑
int c[MAX],tc[MAX];
#include<stdio.h>
#include<string.h>
int main(){
	int time,max,i,j,k;
	int v[55],num[55];
//	freopen("in.txt","r",stdin);
	while(scanf("%d",&time)!=EOF && time >= 0){
		max = 0;
		memset(v,0,sizeof(v));
		memset(num,0,sizeof(num));

		for(i = 0; i < time; i++){
			scanf("%d%d",&v[i],&num[i]);
			max += v[i] * num[i];
		}
		for(i = 0; i <= max; i++){
			c[i] = 0;
			tc[i] = 0;
		}
		for(i = 0; i <= num[0] * v[0]; i+=v[0]){
			c[i] = 1;
		}
		for(i = 1; i < time; i++){
			for(j = 0; j <= max; j++){
				for(k = 0; k + j <= max && k <= num[i] * v[i]; k+=v[i]){
					tc[j+k] += c[j];
				}
			}
			for(j = 0; j <= max; j++){
				c[j] = tc[j];
				tc[j] = 0;
			}
		}
		for(i = max / 2; i >= 0; i--){
			if(c[i]){
				printf("%d %d\n",max-i,i);
				break;
			}
		}
	}
	return 0;
}