01分数规划

题套：
有一堆物品，每一个物品有一个收益ai，一个代价bi，我们要求一个方案使选择的∑ai/∑bi为最值。

套路解题：
我们可用二分。设我们求的是最大，则我们要使∑ai/∑bi>=x ,等价于∑a[i]-x*∑b[i]>=0

例一：
Dropping Tests poj2976

#include<cstdio>
#include<algorithm>
using namespace std;
int n,k,a[2000],b[2000];
double p[2000];
bool ok(double x)
{
    for(int i=1;i<=n;i++)p[i]=a[i]-x*b[i];
    sort(p+1,p+n+1);
    double ans=0;
    for(int i=n;i>=k+1;i--) ans+=p[i];
    return ans>=0;
}
void solve()
{
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(int i=1;i<=n;i++) scanf("%d",&b[i]);
    double l=0,r=0x3f;
    while(r-l>1e-8)
    {
        double mid=(l+r)/2;
        if(ok(mid)) l=mid;else r=mid;
    }
    printf("%.0lf\n",l*100);
}
int main()
{
    while(scanf("%d%d",&n,&k),n|k)
        solve();
}

例二：
Desert King poj2728
给出一个n个点的完全图，每条边有两个权值cost和len。
是二分答案x，把每条边边权设为cost-x*len，跑最小生成树，判答案的正负。建议用prim跑，复杂度比较科学。

#include<cstdio>
#include<cmath>
#include<iostream>
#include<algorithm> 
#define ld long double
using namespace std;
#define SZ 1066
int n,x[SZ],y[SZ],h[SZ],fm[SZ];
bool del[SZ];
ld mind[SZ],d[SZ][SZ],d1[SZ][SZ],d2[SZ][SZ];
ld Abs(ld x) 
{
    return (x>=0)? x:-x;
}
ld sqr(ld x) {return x*x;}
bool ok(ld x)
{
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++) 
     d[i][j]=d1[i][j]-x*d2[i][j];
    for(int i=0;i<=n;i++) mind[i]=1000000000,del[i]=0;
    mind[1]=0;fm[1]=0;
    ld ans=0;
    for(int i=1;i<=n;i++)
    {
        int mp=0;
        for(int j=1;j<=n;j++)
        {
            if(mind[j]<mind[mp]&&!del[j]) mp=j;
        }
        ans+=d[fm[mp]][mp];del[mp]=1;
        for(int j=1;j<=n;j++)
        {
            if(d[mp][j]>=mind[j]||del[j]) continue;
            mind[j]=d[mp][j];fm[j]=mp;
        }
    }
    return ans<=0;
}
void solve()
{
    for(int i=1;i<=n;i++) scanf("%d%d%d",&x[i],&y[i],&h[i]);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
     {
        d1[i][j]=Abs(h[i]-h[j]);//
        d2[i][j]=sqrt(sqr(x[i]-x[j])+sqr(y[i]-y[j]));
     }
    ld l=0,r=1000000000;
    while(r-l>1e-6)
    {
        double mid=(l+r)/2;////
        if(ok(mid)) r=mid;else l=mid;//qiu min
    }
    printf("%.3lf\n",(double)l);
}
int main()
{
    while(scanf("%d",&n),n) solve();
}