P1570 KC喝咖啡

题目背景

KC喝咖啡

题目描述

话说KC和SH在福州的时候常常跑去85°C喝咖啡或者其他的一些什么东西。

这天，KC想要喝一杯咖啡，服务员告诉他，现在有n种调料，这杯咖啡只可以加入其中的m种(当然KC一定会加入m种，不会加入少于m种的调料)，根据加入的调料不同，制成这杯咖啡要用的时间也不同，得到的咖啡的美味度也不同。

KC在得知所有的n种调料后，作为曾经的化竞之神的他，马上就知道了所有调料消耗的时间ci以及调料的美味度vi。由于KC急着回去刷(play)题(games)，所以他想尽快喝到这杯咖啡，但他又想喝到美味的咖啡，所以他想出了一个办法，他要喝到(Σvi)/(Σci)最大的咖啡，也就是单位时间的美味度最大的咖啡。

现在，KC把调料信息告诉了SH，要SH帮他算出喝到的咖啡的(Σvi)/(Σci)，但SH不想帮KC算这东西，于是KC就只能拜托你来算了。

注释：Σ表示求和，所以(Σvi)/(Σci)表示 美味度的总和 除以 消耗时间的总和。

输入输出格式

输入格式：

输入数据第一行为一个整数n，m，表示调料种数和能加入的调料数。接下来2行，每行为n个数，第一行第i个整数表示调料i的美味度vi，第二行第i个整数表示调料i消耗的时间ci。

输出格式：

一个整数T，表示KC喝的咖啡的最大(Σvi)/(Σci)，保留三位小数。

输入输出样例

输入样例#1：

3 2
1 2 3
3 2 1

输出样例#1：

1.667
（选2号和3号调料,（2+3）/(2+1)=1.667）

说明

数据范围：

20%：1<=n<=5。

50%：1<=n<=10。

80%：1<=n<=50。

100%：1<=n<=200，1<=m<=n，1<=c[i],v[i]<=10000。

保证答案不超过1000。

思路：二分答案，题中条件sigma(vi)/sigma(ci)>=x可以转换为sigma(vi-x*ci)>=0，这样我们对 yi= vi - x*ci {1<=i<=n}从大到小排序，计算sum(yi) {1<=i<=k} 如果sum(yi){1<=i<=k}>=0 ，则说明 sigma(vi) / sigma(ci) >= x, 成立。

#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
#define mem(ar,num) memset(ar,num,sizeof(ar))
#define me(ar) memset(ar,0,sizeof(ar))
#define lowbit(x) (x&(-x))
#define IOS ios::sync_with_stdio(false)
#define DEBuG cout<<endl<<"DEBuG"<<endl;
using namespace std;
int n,m;//调料种数和能加入的调料数
struct node{
    int v,c;double w;
}no[210];
bool cmp(node a,node b){return a.w>b.w;}
bool check(double x){
    for(int i=1;i<=n;i++)no[i].w=no[i].v-x*no[i].c;
    sort(no+1,no+n+1,cmp);
    double tot=0;
    for(int i=1;i<=m;i++)tot+=no[i].w;
    return tot>=0;
}
int main(){IOS;
    cin>>n>>m;
    double l=0,r,mid;
    for(int i=1;i<=n;i++)
        cin>>no[i].v;
    for(int i=1;i<=n;i++){
        cin>>no[i].c;
        r=r<no[i].v*1.0/no[i].c?no[i].v*1.0/no[i].c:r;
    }
    while(r-l>1e-6){
        mid=(l+r)/2;
        if(check(mid))l=mid;
        else r=mid;
    }
    printf("%.3lf",l);
    return 0;
}