poj 2184 Cow Exhibition

参考:

http://www.cnblogs.com/rainydays/archive/2012/07/04/2576077.html

题意：给定一些奶牛，每个牛有s和f两个属性值，有正有负，要求选出一些牛，使得这些牛的两种属性的和的加和最大，且这些牛的两种属性分别求加和不能为负。

分析：dp，开始想到dp[i][s][f]，表示前i头牛能否实现属性和分别为s,f。空间和时间都不允许，要将f从状态中拿出來，让f的属性和作为所求的值。即变为d[i][s]=f的形式。表示用前i头牛构成s属性和为s的情况下f属性和最大为多少。状态转移从两种情况来，即用或者不用当前的牛。dp[i][j] = max(dp[i - 1][j - s[i]] + f[i], dp[i - 1][j])。在实际操作的时候可以将第一维去掉，进行空间上的优化。但是由于s[i]的值有正有负，所以在填写数组的顺序要根据s[i]的值来决定，若为正则从右到左（类似01背包的空间优化），若为负则从左到右。

注意：动态规划中状态维和值是可以相互转化的。状态维过多，效率低的时候，可以把将其转化为数组值；同理，数组值不唯一无法规划时，可以增加状态维使状态更详细。

仔细看了背包九讲 1.4里提及的初始化的细节问题:

有的题目要求 1)恰好装满背包时的最优解  是==关系，而另外一些则要求 2) 没有“恰"字，即不要求装满背包，是<=关系

解此2184就利用 == 关系。

<=关系，经典背包初始化，DP[0..v] = 0

==关系，则按照D[1..v] = -inf D[0] =0的初始化方式

可以这样理解：初始化的F数组事实上就是在没有任何物品可以放入背包时的合法状态。如果要求背包恰好装满，那么此时只有容量为0的背包可以在什么也不装且价值为0的情况下被“恰好装满”，其它容量的背包均没有合法的解，属于未定义的状态，应该被赋值为-∞ 了。如果背包并非必须被装满，那么任何容量的背包
都有一个合法解“什么都不装”，这个解的价值为0，所以初始时状态的值也就全部为0了。

#include <iostream>
#include <stdio.h>
#include <cstring>
#include <algorithm>


using namespace std;

#define MAX_N (100+1)
int S[MAX_N],F[MAX_N];

int n;

#define inf 200000000

#define MAX_F (2000 * 100)
#define SHIFT (1000 * 100)
//F[i][s] = f //表示放入第i个cow 和为s的情况下,所有达到的最大值f
int DP[MAX_F+100];

int l=0,r=0;
#define printf //
void ZeroOnePack(int s, int f)
{
	l=min(l,l+s);
	r=max(r,r+s);

	int begin =l;
	int end = r;
	int step = 1;

	

	if(s > 0)
	{
		swap(begin,end);
		step = -1;
	}
	//    0 <= v-s <<MAX_F ->    s<= v <= s+MAX_F
	//    0 <= v << MAX_F
	printf("s:%d f:%d begin=%d,end=%d\n",s,f,begin,end);
	for(int v =  begin; v != end+step ; v += step)
	{
		printf("v[%d]-s[%d] = %d\n",v,s,v-s);
		int ldp = DP[ v -s +  SHIFT ] +f; 

		if( DP[v  + SHIFT] == -inf || ldp > DP[v + SHIFT])
		{
			//printf("DP[%d] = %d\n",v,ldp);
			DP[v + SHIFT] = ldp;
		}
	}
}

int main()
{
	while(cin>>n)
	{
		int maxp = 0;

		for(int i=1;i<=n;i++)
		{
			cin>>S[i]>>F[i];
		}

		for(int v=0;v<= MAX_F;v++)
			DP[v] = -inf;
		
		DP[0+SHIFT] = 0;

		l=0;
		r=0;
		printf("n=%d\n",n);
		for(int _i=1;_i <= n;_i++)
		{
			printf("i=%d n=%d\n",_i,n);
			ZeroOnePack(S[_i],F[_i]);
			printf("i=%d n=%d\n",_i,n);
		}


		for(int v=0;v<=r;v++)
		{
			if(DP[v+SHIFT] >= 0 )
				maxp = max(maxp, v + DP[v+SHIFT]);
		}

		cout<<maxp<<endl;
	}

}