【贪心】【01规划】luogu 1570 KC喝咖啡

题目链接：http://oj.luogu.org:8888/problemshow.php?pid=1570

分析：
1、裸搞。

按比值降序排列，然后去前m个

60分。

2、0 1 分数规划

 关于 0 1 规划：http://blog.csdn.net/hhaile/article/details/8883652

设最终答案为R，x[i]为0或1，则 Σ(v[i]*x[i]) / Σ(c[i]*x[i]) = R

即 Σ(v[i]*x[i]) - Σ(c[i]*x[i]*R) =0
即 Σ[(v[i]*x[i]) - (c[i]*x[i]*R)] =0

设 f(R)= Σ[(v[i]*x[i]) - (c[i]*x[i]*R)]

二分R，求f(R)=0

当f(R)<0时，R偏大

当f(R)>0时，R偏小

计算每一对 v[i] - (c[i]*R)，排序求前m大的和即为f(R) 。

代码：

60：

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
typedef struct
{
int a,b;
double c;
} nc;
nc t[201];
int n,m,i;
bool cmp(nc a1,nc a2)
{
if (fabs(a1.c-a2.c)<0.00000001) return a1.a<a2.a;
else  return a1.c>a2.c;
}
int main()
{
scanf("%d%d",&n,&m);
for (i=1;i<=n;i++) scanf("%d",&t[i].a);
for (i=1;i<=n;i++)
{
scanf("%d",&t[i].b);
t[i].c=t[i].a/(t[i].b*1.0);
}
sort(t+1,t+1+n,cmp);
int m1=0,m2=0;
for (i=1;i<=m;i++)
{
m1+=t[i].a;
m2+=t[i].b;
}
printf("%.3lf",m1/(m2*1.0));
return 0;

}

100：

#include<cstdio>
using namespace std;
int n,m,i,j,k;double ma,ans1,ans2;bool b[205];
struct jgtt
{
double c,v;
}a[205];
int main()
{
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)
scanf("%lf",&a[i].v);
for(i=1;i<=n;i++)
scanf("%lf",&a[i].c);
ans1=0;ans2=0;
for(i=1;i<=m;i++)
{
ma=0;int ji;
for(j=1;j<=n;j++)
if(!b[j]&&((a[j].v+ans1)/(a[j].c+ans2))>ma)
{
ma=(a[j].v+ans1)/(a[j].c+ans2);
ji=j;
}
b[ji]=1;
ans1+=a[ji].v;
ans2+=a[ji].c;
}
printf("%0.3f",ans1/ans2);
}