poj1015

题目有难度。大致题意为要选定n个群众，从中选取m个作为陪审团。每个群众都有两个得分，分别为defence and prosecution给出的得分（0——20），表示对该群众的满意度（与得分成正比）。先要使得选出的m个群众的所有|D得分-P得分|最小，若存在相等（同时最小），则选择所有D+P之和最大的一组，并输出所有D之和，和所有P之和，以及陪审团的各群众编号。

题目意思很明确了。这里我最先的思路是：令dp[i][j]表示前i个人中选择了j个人时的最小D-P之和（或者|D-P|之和）。状态转移方程为：

dp[i][j]=Max{dp[i-1][j-1]+v[i],dp[i-1][j]},求解dp,然后按p+d为关键字从小到大排序，再逆序判定，确定满足题意的m个群众。若该题改为所有群众D-P之和最小而不是|D-P|的绝对值最小则上述思想就是正确的dp,然而题目条件是后者，所以这个状态转移方程是错误（很明显）。

下面讲述正确的dp思想：(v=d-p,s=d+p)

令dp[i][j]表示选了i个群众，所有i个群众D-P得分为j时最大的D+P之和。则分析如下：令k不在dp[i-1][j-v[k]]所选群众中，则有dp[i][j]=max{s[k]+dp[i-1][j-v[k]]}。同时令path[i][j]=k表示选了i个群众，所有群众D-P得分为j且D+P最大时，选了第k个群众。则可以依据path[i][j]可以递推出所有的选择群众。例如下一个群众编号为path[i-1][j-v[path[i][j]]]，依次类推。由于D-P之和可能为负数，故必须加上一个偏移量x，m的范围为1——20，而评分范围在0——20，故所有m个人D-P范围为-20*20——20*20,得-400<=j<=400，故可以加上偏移量400(20*m)使之在0——800之间（0——40*m之间）。然后就是dp求解了。求得dp[m][0——800] (dp[m][0——40*m])，然后从中心400（20*m)向两侧枚举，知道有有意义的dp出现，然后比较相对的两个值，取大值，其j即为D-P之和，而dp为D+P之和。那么D=j+dp-400(j+dp-20*m)，P=dp-j+400(dp+j+20*m)。最后利用path数组将m个群众编号保存，排序，然后从小到大输出即可。

下面是程序：320K+735MS

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <algorithm>
#define Max 210 //n限制
#define Maxx 30 //m限制
#define Inf 100000010 // 无穷大
using namespace std;
int d[Max]; 
int p[Max];
int v[Max]; // d-p
int s[Max]; // d+p
int record[Max]; //记录最后m个群众编号
int dp[Maxx][Maxx*30]; //dp[i][j]表示i个人，D-P之和为j时最大的D+P之和
int path[Maxx][Maxx*30]; //path[i][j]表示i个人，D-P之和为j时最大的D+P之和时选择了第k个人
int n,m; 
bool Is_exist(int i,int j,int k){ // 判断第k个人是否在dp[i][j]所选的群众中
	while(i>0 && path[i][j]!=k){ //i若为0，则直接返回false，即不再。
		if(path[i][j]==0) // 注意当i>0时，若path[i][j]=0，则说明不是合法的情况，直接返回false
			return false;
		j-=v[path[i][j]]; // 正常情况，下一个编号，直到i为0
		i--;
	}
	return (i>0)?true:false; // 若退出while时，是由于path[i][j]!=k，则说明在其中，否则不再其中
}
int main(){
	int Case=1;
	while(scanf("%d%d",&n,&m),n){ 
		int i,j,k;
		for(i=1;i<=n;i++){
			scanf("%d%d",&p[i],&d[i]);
			v[i]=d[i]-p[i]; //求d-p
			s[i]=d[i]+p[i]; //求d+p
		}
		for(i=0;i<=m;i++)
			for(j=0;j<30*Maxx;j++)
				dp[i][j]=-Inf; //初始化为负无穷
		memset(path,0,sizeof(path)); // 初始化为0
		int line=20*m; //最大D-P
		dp[0][line]=0; //0的偏移为line,初始化为0，其余为负无穷
		for(i=1;i<=m;i++){ // 求解dp[m][0——2*line],注意枚举范围为0——2*line
			for(j=0;j<=2*line;j++){
				for(k=1;k<=n;k++){  //枚举k(1——n)
					int temp=j-v[k];
					if(temp>=0 && temp<=2*line){ //若D-P合法
						if(!Is_exist(i-1,temp,k)){ //且不再dp[i-1][temp]群众中
							if(dp[i-1][temp]+s[k]>dp[i][j]){ // 状态转移方程求最大值和此时的群众编号
								dp[i][j]=dp[i-1][temp]+s[k];
								path[i][j]=k;
							}
						}
					}
				}
			}
		}
		for(i=0;i<=line;i++) // 从中心向两边试探，只到有>=0情况出现
			if(dp[m][line+i]>=0 || dp[m][line-i]>=0)
				break;
	    int dd,pp;
		if(dp[m][line+i]>dp[m][line-i]){ // 取大值
			dd=line+i;
		    pp=dp[m][line+i];
		}
		else{
			dd=line-i;
			pp=dp[m][line-i];
		}
		printf("Jury #%d\nBest jury has value ",Case++);
		int dsum=(dd+pp-line)/2; // 计算所有m个人的D和，注意要-line
		int psum=(pp-dd+line)/2; // 计算所有m个人的P和,注意要+line
		printf("%d for prosecution and value %d for defence:\n",psum,dsum);
		int pivot=0;
		for(i=0;i<m;i++){ // 利用path计算所有选择的群众编号，保存在record数组中
			int index=path[m-i][dd];
		    record[pivot++]=index;
			dd-=v[index];
		}
		sort(record,record+pivot); //从小到大排序
		for(i=0;i<pivot;i++) // 顺序输出
		    printf(" %d",record[i]);
		printf("\n"); //空格
	}
	return 0;
}