POJ2976-01分数规划&二分-Dropping tests

https://vjudge.net/problem/POJ-2976
给定n个人，有两种分数,a,b
要求你选取一部分,要求 sigma(a)/sigma(b) *100得值最大。
化简式子

枚举的mid 不断的逼近最优解。 设定一个精度。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;
/* 01分数规划 可以用二分，或者用一种递推。
 思路是 不断的维护一个最优解。
 因为暴力计算最优解，是不可行的，因为他要找那么那么多那么多组合
  而01分数规划是通过对 式子的转化。
  如果 sigma(x1-y1*mid) 大于0，如果当前选择的 分数构成的比例已经大于mid
  也就是说，必然存在一个更优的解（因为我们当前计算就已经比他优了）
  那为什么不用我们 现在计算的那个点呢，
  因为可能有些选择比他更优，
  （如果没有比他更优的话，最后不断的二分，mid会不断的回退到这里）
*/
const int maxn=3000;
const double eps=1e-5;
double  a[maxn];
double  b[maxn];
double aa[maxn];
bool cmp2(double a,double b){
    return a<b;
}
int main()
{   int m,k;
    while(~scanf("%d%d",&m,&k)){
         if(!m&&!k) break;
         for(int i=0;i<m;i++){
              scanf("%lf",&a[i]);
         }
         for(int i=0;i<m;i++){
             scanf("%lf",&b[i]);
         }
         double l=0.0;
         double r=1.0;
         double ans=0;
         double mid;
         while(r-l>=eps){
               mid=(r+l)/2;
               //cout<<mid<<endl;
              double sum=0;
              for(int i=0;i<m;i++){
                 aa[i]=a[i]*1.0-b[i]*mid;
              }
              sort(aa,aa+m,cmp2);
              for(int i=k;i<m;i++){
                  sum+=aa[i];
              }
              if(sum>=0){
                 l=mid;
                 ans=mid;
              }
              else
                r=mid;
         }
         printf("%.0f\n",ans*100);
    }
    return 0;
}