POJ1015 Jury Compromise

类似于背包问题，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]);

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

#include
using namespace std;
#define max(a,b) (a>b?a:b)
const int N=810;
const int oo=-(1<<30);
int dp[210][21]
,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;
else
dp[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
[m][k+400]>=0&&dp
[m][k+400]>=dp
[m][400-k])
{
ansk=400+k;
break;
}
elseif(dp
[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;
}