JZOJ 4230. 【五校联考4day1】淬炼神体

最新推荐文章于 2022-10-08 18:48:17 发布

滑稽大佬

最新推荐文章于 2022-10-08 18:48:17 发布

阅读量136

点赞数

分类专栏： 01分数规划文章标签： 01分数规划

本文链接：https://blog.csdn.net/qq_35786326/article/details/86564757

版权

01分数规划专栏收录该内容

1 篇文章 0 订阅

订阅专栏

...

题目：
分析：
代码：

题目：

传送门

分析：

$01$ 分数规划裸题( $y i d o n g$ 说得好：)逃

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<set>
#include<queue>
#include<vector>
#include<map>
#include<list>
#include<ctime>
#include<iomanip>
#include<string>
#include<bitset>
#include<deque>
#include<set>
#define LL long long
using namespace std;
inline LL read(){
	LL d=0,f=1;char s=getchar();
	while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
	while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
	return d*f;
}
double d[100005],maxx;
int x[100005],y[100005]; 
bool cmp(double a,double b)
{
	return a>b;
}
int main()
{
	int n=read(),m=read();
	for(int i=1;i<=n;i++) x[i]=read();
	maxx=0;
	for(int i=1;i<=n;i++) 
	{
		y[i]=read();
		if(maxx<x[i]/y[i]) maxx=x[i]/y[i];
	}
	double l=0,r=maxx;
	while(abs(r-l)>=0.0001)
	{
		double mid=(l+r+0.0001)/2;
		for(int i=1;i<=n;i++) d[i]=x[i]-mid*y[i];
		sort(d+1,d+1+n,cmp);
		double sum=0;
		for(int i=1;i<=m;i++) sum+=d[i];
		if(sum>=0) l=mid; else r=mid-0.0001;
	}
	printf("%.3lf",l);
	return 0;
}