POJ2976 Dropping tests

最新推荐文章于 2022-06-06 12:21:21 发布

forezxl

最新推荐文章于 2022-06-06 12:21:21 发布

阅读量132

点赞数 1

分类专栏： POJ 其他---分数规划蒟蒻zxl的Blog专栏文章标签： POJ 01分数规划

本文链接：https://blog.csdn.net/a1799342217/article/details/80385820

版权

蒟蒻zxl的Blog专栏同时被 3 个专栏收录

460 篇文章 6 订阅

订阅专栏

POJ

22 篇文章 0 订阅

订阅专栏

其他---分数规划

6 篇文章 0 订阅

订阅专栏

这篇博客介绍了POJ2976题目的解题思路，主要涉及01分数规划的运用。通过二分法寻找最大值，将a[i]-r*b[i]进行排序，并选取n-k个满足不等式的元素，从而达到最大化目标。

摘要由CSDN通过智能技术生成

01分数规划

题目传送门

题目大意：给你两串数 $a[i],b[i]$ ，选 $n-k$ 个数使得 $\frac {\sum a[i]}{\sum b[i]}$ 最大。

这道题就是最裸的01分数规划了。01分数规划就是求 $\max/\min(\frac{\sum a[i]*x_i}{\sum b[i]*x_i})$ ，其中 $x_i\in\{0,1\}$

二分一个答案 $r$ ，即 $\frac {\sum a[i]}{\sum b[i]}\geq r$ ，化简后得 $\sum a[i]-r\sum b[i]\geq 0$

那么我们把 $a[i]-r*b[i]$ 排序，然后取前 $n-k$ 个看看是不是 $\geq0$ 就行了。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
#define eps 1e-8
using namespace std;
typedef double DB;
int n,k,a[N],b[N];
DB s[N],sum;
inline bool pd(DB x){
    for (int i=1;i<=n;i++) s[i]=a[i]-x*b[i];
    sort(s+1,s+n+1),sum=0;
    for (int i=k+1;i<=n;i++) sum+=s[i];
    return sum>-eps;
}
int main(){
    while (~scanf("%d%d",&n,&k)&&(n||k)){
        for (int i=1;i<=n;i++) scanf("%d",&a[i]);
        for (int i=1;i<=n;i++) scanf("%d",&b[i]);
        DB l=0,r=1,mid;
        while (r-l>-eps)
            if (pd(mid=(l+r)/2)) l=mid+eps*100;
            else r=mid-eps*100;
        printf("%.0f\n",100*l);
    }
}