POJ1015 Jury Compromise

来源：互联网发布：淘宝俄罗斯教父代购编辑：程序博客网时间：2024/05/13 03:15

类似于背包问题，n个人每一个都有取或者不取。在0-1背包中，定义一个递归式dp[i][j]，表示前i个物品，在背包容量大小为j的情况下，最大装载量。动态规划方程如下：dp[i][j]= max(dp[i-1][j], dp[i-1][j-c[i]+v[i])

在这个问题中，d[i]-p[i]看成背包中每个物品占用的体积。设dp[i][j][k]前i个人中选择j个人差值为k的总分最大和。状态转移方程：dp[i][j][k]=max(dp[i-1][j][k],dp[i-1][j-1][k-(d[i]-p[i])]+d[i]+p[i]);

代码一直改不对，贴一段好了。。。

#includeusing namespace std;#define max(a,b) (a>b?a:b)const int N=810;const int oo=-(1<<30);int dp[210][21][N],d[210],p[210],q[210],end;void dfs(int i,int j,int k){if(i==0||j==0) return;if(dp[i-1][j][k]{ q[end++]=i; dfs(i-1,j-1,k-d[i]+p[i]);}else dfs(i-1,j,k);}int main(){int n,m,num=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]); for(i=0;i<=n;i++)for(j=0;j<=m;j++)for(k=0;k<=800;k++)if(j==0&&k==400)dp[i][j][k]=0;elsedp[i][j][k]=oo; for(i=1;i<=n;i++) for(j=1;j<=m;j++)  {for(k=0;k<=800;k++)if(k>=d[i]-p[i]&&k-d[i]+p[i]<=800) dp[i][j][k]=max(dp[i-1][j][k],dp[i-1][j-1][k-d[i]+p[i]]+d[i]+p[i]);  } int ansk; for(k=0;;k++) { if(dp[n][m][k+400]>=0&&dp[n][m][k+400]>=dp[n][m][400-k])  {ansk=400+k;break;  }  elseif(dp[n][m][400-k]>=0)  {ansk=400-k;break;  } } end=0; dfs(n,m,ansk); int px=0,dx=0; for(i=0;i  px+=p[q[i]],dx+=d[q[i]]; printf("Jury #%d\nBest jury has value %d forprosecution and value %d for defence:\n",num++,px,dx); for(i=end-1;i>=0;i--)  printf(" %d",q[i]); printf("\n\n");}return 0;}