POJ 1015 Jury Compromise【DP】

罗大神说这题很简单,,,,然而我着实写的很难过。。。


题目链接:

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=110495#problem/K

题意:

给定n个人对罪犯的d值和p值,从中选m个人使得abs(sum(d)−sum(p))最小,如果有多种情况则输出sum(p)+sum(d)最大的情况,并升序输出选择的人的编号。

分析:

这题是既要考虑差又要考虑和。还是绝对值最小。。刚看见有点懵逼。

其实不方,这题d和p最大只有20,m最大又只有20,就是说差的绝对值最大只有20 * 20,

那么我们设dp[i][j]为已经选出i个人,sum(p)−sum(d)为j的最大sum(p)+sum(d)。dp[i][j]=−1 表示这个状态无法达到。

那么问题来了,绝对值怎么处理,如果sum(p)−sum(d)小于0怎么办。。。

我们知道差的绝对值最大为20 * 20,那么我们给每个差都加上20 * 20 ,保证j大于等于0,这样就相当于以20 * 20 为原点,最终只要找到距离20 * 20 最近的dp值大于等于0的点就好了。。

想到这里就差不多了。。

很容易得到状态转移方程。

对于每个i和j的状态枚举给定的n个人,注意判断之前当前这个人之前是否出现过。

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define sa(a) scanf("%d", &a)
#define sal(a) scanf("%I64d", &a)
const int maxn = 1000 + 5, maxm = 200 + 5, INF = 0x3f3f3f3f;
int p[maxm], d[maxm];
struct DP{int a; int now;};
DP dp[20 + 5][maxn];
int tt[maxn];
int cnt;
void output(int he, int ans)
{
if(he == 0) return;
int index = dp[he][ans].now;
tt[he--] = index;
output(he, ans - p[index] + d[index]);
}
bool vis(int i, int j, int k)
{
while(i >= 0){
int index = dp[i][j].now;
if(index == k) return true;
i--;
j -= p[index] - d[index];
}
return false;
}
int main (void)
{
int n, m;
int cnt = 1;
while(scanf("%d%d", &n, &m) && n + m){
for(int i = 0; i < n; i++){
sa(p[i]); sa(d[i]);
}
int t = m * 20;
for(int i = 0; i <= m; i++){
for(int j = 0; j <= 2 * t; j++){
dp[i][j].a = -1;
dp[i][j].now = -1;
}
}
dp[0][t].a = 0;
for(int i = 0; i < m; i++){
for(int j = 0; j <= t * 2; j++){
if(dp[i][j].a < 0) continue;//状态无法到达
for(int k = 0; k < n; k++){
if(dp[i][j].a + p[k] + d[k] > dp[i + 1][j + p[k] - d[k]].a){
if(vis(i, j, k)) continue;//之前出现过
dp[i + 1][j + p[k] - d[k]].a = dp[i][j].a + p[k] + d[k];
dp[i + 1][j + p[k] - d[k]].now = k;
}
}
}
}
int i = 0;
while(dp[m][t + i].a < 0 && dp[m][t - i].a < 0) i++;
int ans;
if(dp[m][t + i].a > dp[m][t - i].a) ans = t + i;
else ans = t - i;
printf("Jury #%d\n", cnt++);
printf("Best jury has value %d for prosecution and value %d for defence:\n",
(ans - t + dp[m][ans].a)/2, (dp[m][ans].a - ans + t)/2);
output(m, ans);
sort(tt + 1, tt + m + 1);
for(int i = 1; i <= m; i++)
printf(" %d%c", tt[i] + 1, i == m?'\n':' ');
printf("\n");
}
return 0;
}
上一篇:POJ 1015 Jury Compromise(dp坑)


下一篇:1.5.6、CDH 搭建Hadoop在安装之前(定制安装解决方案---使用Cloudera Manager模板创建CDH群集)