【POJ 3621】[01分数规划]Dropping tests

最新推荐文章于 2021-02-01 19:59:44 发布

ORZSunIsMe

最新推荐文章于 2021-02-01 19:59:44 发布

阅读量355

点赞数

分类专栏： 01分数规划文章标签： 01分数规划

本文链接：https://blog.csdn.net/qq_35479641/article/details/59548594

版权

01分数规划专栏收录该内容

2 篇文章 0 订阅

订阅专栏

题意

给出n，k和n对二元组（a，b），要求从中选出n-k个二元组使得 $100\times\frac{\Sigma{a_i}}{\Sigma{b_i}}$ 最大。
原题链接–POJ

题目分析

我们假设答案为ans，那么一定对于任意方案有 $ans \geq 100\times\frac{\Sigma{a_i}}{\Sigma{b_i}}$ ，
而 $ans = max(100\times \frac{\Sigma{a_i}}{\Sigma{b_i}})$ 时 $ans$ 为最优解，
$ans > max(100\times \frac{\Sigma{a_i}}{\Sigma{b_i}})$ 时 $ans$ 一定不是解，
故解在 $ans \leq max(100\times\frac{\Sigma{a_i}}{\Sigma{b_i}})$ 中
变一下形就可以得到 $max(\Sigma{(100\times a_i-ans\cdot b_i)}) \geq 0$ 。
若令 $g(ans)=max(\Sigma{(100\times a_i-ans\cdot b_i)})$ ,我们很容易就可以发现其为单调的： $g(ans)$ 随着 $ans$ 的增大而减小。
于是我们不妨二分 $ans$ 再预处理出 $d_i=a_i-ans\cdot b_i$ ，从大到小排序后贪心算出 $g(ans)$ 判定一下就行了。

代码

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

#define MAXN 1000
#define MAXM
#define INF 0x3f3f3f3f
typedef long long int LL;

template<class T>
void Read(T &x){
    x=0;char c=getchar();bool flag=0;
    while(c<'0'||'9'<c){if(c=='-')flag=1;c=getchar();}
    while('0'<=c&&c<='9'){x=x*10+c-'0';c=getchar();}
    if(flag)x=-x;
}

const double EPS = 1e-3;
int Sign(const double x){
    if(x>EPS)return 1;
    else if(x<-EPS)return -1;
    else return 0;
}

bool dcmp(const double &x,const double &y){
    return Sign(x-y)>0;
}

double a[MAXN+10],b[MAXN+10];
double d[MAXN+10];
int n,k;

bool check(double l){
    for(int i=1;i<=n;++i)
        d[i]=a[i]*100.0-l*b[i];
    sort(d+1,d+n+1,dcmp);
    double sum=0;
    for(int i=1;i<=k;++i)
        sum+=d[i];
    return Sign(sum)>=0;
}

int main(){
    while(~scanf("%d%d",&n,&k)){
        if(!n&&!k)break;
        k=n-k;

        for(int i=1;i<=n;++i)scanf("%lf",&a[i]);
        for(int i=1;i<=n;++i)scanf("%lf",&b[i]);

        double l=0,r=100,mid;
        double ans=0;
        while(r-l>EPS){
            mid=(l+r)/2.0;
            if(check(mid)){
                ans=max(ans,mid);
                l=mid+EPS;
            }
            else r=mid-EPS;
        }

        printf("%0.0lf\n",ans);
    }
}