poj 1015 动态规划 难

这道题我是看别人的思路写的，因为有个绝对值，无法转移状态，自己傻得非要把绝对值表示出来，

其实可以在计算完成后筛选，我这个菜鸟估计是受背包的的影响，想状态的时候还想着前几个中选出多少个。

这道题的状态是f[i][j],i选了的人数，边界为零，目标为M，j为d-p(a[])的差值和，f值为d+p(b[])值。记如何状态转移呢，

f[i][j]是由f[i-1][x]而来，而x+a[k]=j；要求k没有被用到，因此用path[i][j]表示该状态时，最后一个选着的编号，前面的用

到的编号可以通过path[i-1][j-a[path[i][j]]得到，由此判断k是否可用。而且，由于j值可能为负，因此加上一个可能取值的最大值，

这儿为400，边界为f[0][0+400]=0，在递推的时候，如果j为负或大于2*400，都是不合法的。ps:这道题我调试了一个下午

，居然是数组开的太小，结果以为开的很大，结果出现完全没解释的奇怪现象却不报错。

/*
poj 1015
*/
#include <stdio.h>
#include <algorithm>
using namespace std;
#define M 25
#define N 10000
#define move 400
#define min -9999999
int f[M+10][N],path[M+10][N],a[N],b[N];
int n,m;
int main(int argc, char const *argv[])
{
	int h=0;
	while(scanf("%d%d",&n,&m),h++,n && m)
	{
		for(int i=0;i<n;i++)
		{
			int c,d;
			scanf("%d%d",&c,&d);
			a[i]=c-d;
			b[i]=c+d;
		}
		for(int i=0;i<=M;i++)
			for(int j=0;j<=2*move;j++)
				f[i][j]=min;

		memset(path,-1,sizeof(path));
		f[0][0+move]=0;
		for(int i=1;i<=m;i++)
			for(int j=0;j<=2*move;j++)
			{
				
				//对1到n个依次扫描，看是否可以找到f[i-1][x]不是负无穷，x为j-a[k],k没有使用
				for(int k=0;k<n;k++)
				{
					if(j-a[k]>=0 && j-a[k]<=2*move && f[i-1][j-a[k]]!=min )
					{

						//看k是否用过
						int sig=1;
						int v=i-1;
						int y=j-a[k];
						while(v)
						{
							if(path[v][y]==k)
							{

								sig=0;
								break;
							}
							y=y-a[path[v][y]];
							v--;
						}
						if(sig==1)
						{
							if(f[i][j]<f[i-1][j-a[k]]+b[k])
							{
								f[i][j]=f[i-1][j-a[k]]+b[k];
									path[i][j]=k;
							}
						}
						
					}
				}
			}

		int mm=99999999;
		int z;
		for(int j=0;j<=2*move;j++)
			if(f[m][j]!=min)
			{
				if(abs(j-move)<mm)
				{
					mm=abs(j-move);
					 z=j-move;
				}
				if(abs(j-move)==mm)
				{
					if(f[m][j]>f[m][z+move])
						z=j-move;
				}
			}

		printf("Jury #%d\n",h );
		printf("Best jury has value %d for prosecution and value %d for defence: \n",
			(z+f[m][z+move])/2,(z+f[m][z+move])/2-z );

		int g[M];
		int l=0;
		fill(g,g+M,999999999);

		int j=z+move;
		int i=m;
		while(i)
		{
			g[l++]=path[i][j];
			j=j-a[path[i][j]];
			i--;
		}

		sort(g,g+M);
		for(int i=0;i<m;i++)
			printf("%d ",g[i]+1 );
		printf("\n");
		
	}
	return 0;
}