poj1015【DP.......无奈了】
首先,读题,真是一窍不通。后来看完程序的意思,才明白吧。。
题意:
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;
}