二分法 最大化平均值

摘自 挑战程序设计竞赛

和最大化最小值类似，最大化平均值也可以通过二分法求得。

比如下面这个经典的问题：
有n个物品的重量和价值分别是wi和vi，从中选出k个物品使得单位重量价值最大。

样例输入：

3 2
2 2
5 3
2 1
1
2
3
4
样例输出：

0.75
1
分析：

一般先想到的是将每个物品的单位重量价值算出来，然后排个序，从大到小贪心进行选择，可惜这样是不对的，这样不能保证最后一定是最大平均值，直接用贪心对于这类要涉及多个因素比如求最大平均值的问题就显得不那么正确了。

那么我们这样分析这个问题：
令C(x)为可以选择使得单位重量的价值不小于x。
这样就可以用二分法来解决，不断二分x进行判断，取最大。

我们继续分析：

∑vi / ∑wi 这个式子是我们需要求的单位重量的价值。
i∈s i∈s
1
2
那么就求是否满足：

∑vi / ∑wi >= x
i∈s i∈s
1
2
转换一下得到：

∑(vi - x*wi) >= 0
i∈s
1
2
判断这个式子是否成立即可。这下就可以用一个数组来保存vi - x * wi 的值，并进行排序，从大到小贪心地进行选择求和，如果求和的值大于0，那么此时 的x就是成立的。

代码如下：

#include<cstdio>
#include<algorithm> 
using namespace std;
const int maxn = 10000+10;
const int INF = 1000000;
int w[maxn],v[maxn];
double y[maxn];
int n,k;

bool C(double d)
{
    for(int i=0; i<n; i++)
    {
        y[i] = v[i] - d * w[i];
    }
    sort(y, y+n);

    double sum = 0;
    for(int i=0; i<k; i++)
    {
        sum += y[n-i-1];
    }
    return sum >= 0;
}
void solve()
{
    double lb = 0, ub = INF;
    for(int i=0; i<100; i++)
    {
        double mid = (lb + ub) / 2;
        if(C(mid)) lb = mid;
        else ub = mid;
    }
    printf("%.2f
",lb);
}
int main()
{
    scanf("%d%d",&n,&k);
    for(int i=0; i<n; i++)
    {
        scanf("%d%d",&w[i],&v[i]);
    }
    solve();
    return 0;
}