分数规划

#include<algorithm> 
#include<iostream>
#include<cstdio>
using namespace std;
int a[1000+10],b[1000+10];
int n,m;
double c[1000+10];
bool cmp(double a,double b){
    return a>b;
}
bool judge(double k){//判断当前k是否满足要求，即是否<=ans 
    for(int i=1;i<=n;i++)
        c[i]=a[i]-k*b[i];
    sort(c+1,c+n+1,cmp);
    double sum=0.00;
    for(int i=1;i<=m;i++)
        sum+=c[i];
    if(sum>=0)
        return true;
    return false;
} 
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%d%d",&a[i],&b[i]);//数对（ai，bi） 
    double k,l=0,r=100000;
    for(int i=1;i<=80;i++){//i枚举二分次数
        double mid=(l+r)/2;
        if(judge(mid))
            l=mid;//如果k满足小于ans则枚举比k大的数 
        else
            r=mid;  
    } 
    printf("%lf\n",l);
    return 0;
}
/*2 7 10 3 5 6 100 101 2----1.666667*/ 
/* ans>=k ∑ai/∑bi>=k ∑ai>=∑k*bi ∑(ai-k*bi)>=0 */