poj1015【DP.......无奈了】

来源:互联网 发布:华为网络机顶盒怎么样使用 编辑:程序博客网 时间:2024/05/22 13:49

首先,读题,真是一窍不通。后来看完程序的意思,才明白吧。。
题意:
n个人中选m个,条件是取sum|D-P|最小,当有|D-P|相同的时候取|D+P|最大的。然后输出那些m个人的sumD,sumP,最后升序输出m个人的编号。看完题意,想想好像也不难?。。。瞎想想可能会想出来。
然后看了题解:http://blog.csdn.net/lyy289065406/article/details/6671105
说说感觉吧:
可能一开始去想真的会死命地盯着目标“取最小”不放了。先说说题解DP( j , k )的好处吧,用k代表|D-P|的最小,如果我说用k代表|D+P|的最大情况,似乎可行,而且总共范围是[0 , 20],k也就是[0 , 800]这反倒是取消了题解的转换|D-P|的负数,那么状态转化的公式,我们可以看到题解中说了:
“显然,方案dp(j, k)是由某个可行的方案dp(j-1, x)( -20×m ≤ x ≤ 20×m)演化而来的。”
而且这样真的非常完美“存在某个候选人i,i 在方案dp(j-1, x)中没有被选上,且x+V(i) = k。在所有满足该必要条件的dp(j-1, x)中,选出 dp(j-1, x) + S(i) 的值最大的那个,那么方案dp(j-1, x)再加上候选人i,就演变成了方案 dp(j, k)。”
中间的处理也是很清楚的,加一个path[ j ][ k ]去记录编号。如果要考察编号,就是类似这样嘛path[j-1][k-V[path[j][k]]]
那么现在如果说
我把k看成是最大|D+P|,依样画葫芦一下:方案dp(j, k)是由某个可行的方案dp(j-1, x)( 0≤ x ≤ 40×m)演化而来的,“存在某个候选人i,i 在方案dp(j-1, x)中没有被选上,且x+S(i) = k。在所有满足该必要条件的dp(j-1, x)中,选出 dp(j-1, x) + v(i) 的值最小的那个,那么方案dp(j-1, x)再加上候选人i,就演变成了方案 dp(j, k)。”
中间的处理也是很清楚的,加一个path[ j ][ k ]去记录编号。如果要考察编号,就是类似这样嘛path[j-1][k-S[path[j][k]]]
让我去验证一下。
2333333333,后来证明这种想法是不可取的,透析一下j,k,dp[ j ] [ k ],j是人,但其实k是sum(d-p),以及DP[ i ] [ j ] 是sum(d+p),他的比较就是越大越好;那么DP[i][j]就是第i个人辩控差最小的辩控和最大,最后怎么找出来呢?
从两边遍历当辩控和一旦大于等于0的时候出来,取个辩控差最小的时候。
然后就是那样了。这题主要还是01背包的转化,以及DP[ i ][ j ] 的确定,这种还是多练,多搞吧。

//#include <bits/stdc++.h>#include <iostream>#include <cstdio>#include <math.h>#include <string.h>#include <algorithm>using namespace std;typedef long long LL;typedef unsigned long long ULL;const double eps=1e-5;const double pi=acos(-1.0);const int mod=1e9+7;const int INF=0x3f3f3f3f;const int N=1e2+7;int n;int m;int dp[25][807];int path[24][807];int p[N*2];int d[N*2];int s[N*2];int v[N*2];bool select(int j,int k,int i,int *v){    while(j>0&&path[j][k]!=i){        k-=v[path[j][k]];        j--;    }    return j?0:1;}int main(){    int time=1;    while(cin>>n>>m&&n)    {        int j,k,i;        int *p=new int[n+1];        int *d=new int[n+1];        int *s=new int[n+1];        int *v=new int[n+1];        memset(dp,-1,sizeof(dp));        memset(path,0,sizeof(path));        for(i=1;i<=n;i++)        {            cin>>p[i]>>d[i];            s[i]=p[i]+d[i];            v[i]=p[i]-d[i];        }        int fix=m*20;        dp[0][fix]=0;        for(j=1;j<=m;j++)            for(k=0;k<=2*fix;k++)            {                if(dp[j-1][k]>=0)                {                    for(i=1;i<=n;i++)                        if(dp[j][k+v[i]]<dp[j-1][k]+s[i])                        {                            if(select(j-1,k,i,v))                            {                                dp[j][k+v[i]]=dp[j-1][k]+s[i];                                path[j][k+v[i]]=i;                            }                        }                }            }        for(k=0;k<=fix;k++)            if(dp[m][fix-k]>=0||dp[m][fix+k]>=0)                break;        int div=dp[m][fix-k]>dp[m][fix+k]?(fix-k):(fix+k);        cout<<"Jury #"<<time++<<endl;        printf("Best jury has value %d for prosecution and value %d for defence\n",(dp[m][div]+div-fix)/2,(dp[m][div]-div+fix)/2);    //  cout<<"Best jury has value ";    //  cout<<(dp[m][div]+div-fix)/2<<" for prosecution and value ";    //  cout<<(dp[m][div]-div+fix)/2<<" for defence:"<<endl;        /*int id[N<<1];        for(i=0,j=m,k=div;i<m;i++){            id[i]=path[j][k];            k-=v[id[i]];            j--;        }        sort(id,id+m);        for(i=0;i<m;i++)            cout<<" "<<id[i];*/        cout<<endl;    }    return 0;}
0 0