poj3111

来源：互联网发布：森林优化编辑：程序博客网时间：2024/06/13 20:45

#include <iostream>#include<cstdio>#include<algorithm>using namespace std;const double eps=1e-6;int n,k;struct node2{    double res;    int num;}value[100005];struct node1{    int v, w, num;}j[100005];bool qrt(node2 f,node2 s){    return (f.res)>(s.res);}int main(){   while(scanf("%d%d",&n,&k)!=EOF)   {       for(int i=0;i<n;i++)       {            scanf("%d%d",&j[i].v,&j[i].w);            j[i].num=i;       }           double ub=1<<29,lb=0;           while(ub-lb>eps)           {               double mid=(ub+lb)/2;               for(int i=0;i<n;i++)               {                   value[i].res=j[i].v-mid*j[i].w;        value[i].num=i;               }               sort(value,value+n,qrt);               double sum=0;               for(int i=0;i<k;i++)                sum+=value[i].res;                if(sum>=0)                lb=mid;                else ub=mid;           }           for(int i=0;i<k;i++)            printf("%d%c",value[i].num+1,(i==k-1)?'\n':' ');   }    return 0;}

二分法

据说迭代可以更快的实现，暂时水平有限就这样吧

0 0