分数规划

#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
*/