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][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;
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[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;
}