poj3111 K Best (最大化平均值,二分/牛顿迭代)

题意:

n个物品,每个物品有v(i)和w(i),要求选出k个物品,问sigma(v)/sigma(w)的最大值是多少

思路:

二分解法:

设答案为x,则x=sigma(v)/sigma(w)

如果要满足式子,那么x<=sigma(v)/sigma(w)

把右边的w乘到左边,则式子变为:
x*sigma(w)<=sigma(v)

移项得:
sigma(v)-x*sigma(w)>=0

二分x,check方法为令d(i)=v(i)-x*w(i),排序一下贪心第选择前k大得,判断和是否大于0。

牛顿迭代:

待更

code1:

二分
ps:g++提交ac,c++提交超时

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxm=1e5+5;
const double eps=1e-6;
struct Node{
    int id;
    double v;
}e[maxm];
int v[maxm],w[maxm];
int n,k;
bool cmp(Node a,Node b){
    return a.v>b.v;
}
bool check(double mid){
    for(int i=1;i<=n;i++){
        e[i].id=i;
        e[i].v=v[i]-mid*w[i];
    }
    sort(e+1,e+1+n,cmp);
    double sum=0;
    for(int i=1;i<=k;i++){
        sum+=e[i].v;
    }
    return sum>eps;
}
signed main(){
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++){
        scanf("%d%d",&v[i],&w[i]);
    }
    double l=0,r=1e6;
    while(r-l>eps){
        double mid=(l+r)/2;
        if(check(mid)){
            l=mid;
        }else{
            r=mid;
        }
    }
    for(int i=1;i<=k;i++){
        printf("%d",e[i].id);
        if(i!=k)printf(" ");
        else printf("\n");
    }
    return 0;
}

发布了439 篇原创文章 · 获赞 36 · 访问量 3万+
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 技术黑板 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览