【poj2976 Dropping tests】0/1分数规划裸题

01分数规划的裸题,式子可以变形为 ∑ ( a i − x ∗ b i ) \sum(a_i-x*b_i) (aixbi)

之后对x进行二分即可

#include<cstdio>
#include<climits>
#include<iostream>
#include<algorithm>
using namespace std;
const long double eps = 1e-7;
typedef struct ass {
   
	long long a;
	long long b;
}ass;
bool cmp(long double a, long double b) {
   
	return a > b;
}
long long n, k;
ass temp[10000];
long double tt[10000];
bool check(long double x) {
   
	for (int i = 0; i < n; i++) {
   
		tt[i] = (long double)temp[i].a - x*(temp[i].b);
	}
	sort(tt, tt + n, cmp);
	long double sum = 0;
	for (int i = 0; i < n - k; i++) {
   
		sum += tt[i];
	}
	return sum+eps >= 0;
}
int main(void) {
   
	while (scanf("%lld%lld", &n, &k) && !(n == 0&& k == 0)) {
   
		for (int i = 0; i < n; i++) {
   
			scanf("%lld", &temp[i].a);
		}
		for (int i = 0; i < n; i++) {
   
			scanf("%lld", &temp[i].b);
		}
		long double left = 0.0;
		long double right = 0x3f3f3f3f3f3f;
		while (left + eps < right)
		{
   
			long double mid = (left + right) / 2.0;
			if (check(mid)) {
   
				left = mid;
			}
			else {
   
				right = mid;
			}
		}
		printf("%.0lf\n", left * 100);
	}
	return 0;
}
全部评论

相关推荐

不愿透露姓名的神秘牛友
06-29 17:30
找实习找着找着就要进入7月了,马上秋招也要开始了,找实习还有意义吗?
绝迹的星:有面就面, 没面上就当日薪4位数大佬免费培训, 面上了再考虑要不要实习
点赞 评论 收藏
分享
在瑞幸干两年,奥特曼都得闪灯
不知名的牛友:奥特曼每天只上3分钟班
点赞 评论 收藏
分享
06-23 11:28
门头沟学院 Java
牛客91966197...:也有可能是点拒绝的时候自动弹的话术
点赞 评论 收藏
分享
每晚夜里独自颤抖:你cet6就cet6,cet4就cet4,你写个cet证书等是什么意思。专业技能快赶上项目行数,你做的这2个项目哪里能提现你有这么多技能呢
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

更多
牛客网
牛客网在线编程
牛客网题解
牛客企业服务