poj 1015 Jury Compromise(神一般的DP)

先说题意

一共n个人 每人有两个值 p 和 d。现在要从中选出m个人

使得这m个人的p 的和 与 d的和的差的绝对值最小

如果有相同情况输出p的总和最大的那个

要求输出这种方法选的人的编号，生序输出

这道题。半年前我就看过，我想了不会，我看代码看不懂。

这段时间来一直想搞这道题 每次想想 想不通 都放弃了

刚才，就在刚才。我终于过了。

我决定把这道神一般的dp再敲一遍。

这种题肯定是看了以后完全没有头绪的 先看下数据范围，n是 1-200 m是1-20 p和d 是0-20；

接着想 dp数组肯定不能是一维的。

那么想 如果是二维的 dp[i][j] 怎么表达能算出最终的结果？

因为要选m个人 而且选m在20内 理所应当的想 i应该是21 i表示选了i个人的情况

至于j 要考虑有递推关系 也就是说我们如果知道了dp[i][j]，能知道下一个值 

最简单的情况就是dp[i+1] 在这种情况下再加入一个人，而加入一个人后 j的改变。。。。

怎么处理？

我们选j为当前的pd差！！！

那么dp[i+1][j+p-d] 就可以得到更新！

总结：

dp[i][j]就表示 选i个人pd差是j的情况下pd的最大的和

还要解决一个问题 就是pd差 可能为负数

我们整体右移m*20个单位

而j的范围就从-400到400 变成0-800

至于输出选的人的编号

用path[i][j]表示i j情况下选的最后一个人

那么上一个人就是 path[i-1][j-p+d];

一路追溯到第一个人就行了。

ok

代码

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
using namespace std;

int p[205],d[205];
int dp[21][805],path[21][805];
int ans[21];
int n,m,N;

void init()
{
    memset(dp,-1,sizeof(dp));
    memset(path,0,sizeof(path));
    N=m*20;
    dp[0][N]=0;
}
void dp_work()
{
    int i,j,k,t1,t2,t3;
    for(i=0;i<m;++i)
    {
        for(j=0;j<=N*2;++j)
        {
            if(dp[i][j]>=0)
            {
                for(k=1;k<=n;++k)
                {
                    if(dp[i][j]+p[k]+d[k]>dp[i+1][j+p[k]-d[k]])
                    {
                        t1=i;t2=j;
                        while(t1>0 && path[t1][t2]!=k)
                        {
                            t3=path[t1][t2];
                            t2-=p[t3]-d[t3];
                            t1--;
                        }
                        if(t1==0)
                        {
                            dp[i+1][j+p[k]-d[k]]=dp[i][j]+p[k]+d[k];
                            path[i+1][j+p[k]-d[k]]=k;
                        }
                    }
                }
            }
        }
    }
}

int main()
{
    int i,j,k,pans,dans,ca=1;
    while(scanf("%d %d",&n,&m),(n||m))
    {
        init();
        for(i=1;i<=n;++i)
        scanf("%d %d",&p[i],&d[i]);
        dp_work();
        i=0;
        while(dp[m][N-i]<0 && dp[m][N+i]<0)
        i++;

        if(dp[m][N+i]>dp[m][N-i])
        k=N+i;
        else
        k=N-i;

        pans=(dp[m][k]+k-N)/2;
        dans=(dp[m][k]-k+N)/2;
        printf("Jury #%d\n",ca++);
        printf("Best jury has value %d for prosecution and value %d for defence:\n",pans,dans);

        for(i=0;i<m;++i)
        {
            ans[i]=path[m-i][k];
            k-=p[ans[i]]-d[ans[i]];
        }

        sort(ans,ans+m);

        for(i=0;i<m;++i)
        printf(" %d",ans[i]);
        printf("\n\n");
    }
    return 0;
}