[分数规划] [USACO18OPEN]Talent Show G

最新推荐文章于 2023-09-02 20:33:59 发布

锑元素使者

最新推荐文章于 2023-09-02 20:33:59 发布

阅读量138

点赞数

本文链接：https://blog.csdn.net/Antimonysbguy/article/details/107821719

版权

Portkey

分数考虑分数规划
$\begin{aligned} &\frac{\sum t_i}{\sum w_i}\ge mid\\ \Rightarrow&\sum (t_i-w_i\times mid)\ge0 \end{aligned}$

转化成背包：
花费 $w_i$ ，价值 $t_i-w_i\times mid$ 和要 $> 0$

$O(n^2\log n)$

#include<bits/stdc++.h>
using namespace std;
#define in Read()
int in{
	int i=0,f=1;char ch=0;
	while(!isdigit(ch)&&ch!='-') ch=getchar();
	if(ch=='-') ch=getchar(),f=-1;
	while(isdigit(ch)) i=(i<<1)+(i<<3)+ch-48,ch=getchar();
	return i*f;
}

const int N=1e3+5,INF=1e9;
const double eps=1e-9;
int n,W,w[N],t[N];
double v[N],f[N];

bool check(double mid){
	for(int i=1;i<=n;++i)
		v[i]=(double)t[i]-w[i]*mid;
	for(int i=1;i<=W;++i)
		f[i]=-INF;
	for(int j=1;j<=n;++j)
		for(int i=W;i>=0;--i){
			if(i+w[j]>=W) f[W]=max(f[W],f[i]+v[j]);
			else f[i+w[j]]=max(f[i+w[j]],f[i]+v[j]);
		}
	// cout<<mid<<endl;
	// for(int i=1;i<=n;++i) cout<<v[i]<<" ";cout<<endl;
	// for(int i=1;i<=n;++i) cout<<f[i]<<" ";cout<<endl;
	return f[W]>=0;
}

int main(){
	n=in,W=in;
	for(int i=1;i<=n;++i)
		w[i]=in,t[i]=in;
	double l=0,r=INF;
	while(r-l>eps){
		double mid=(l+r)/2;
		if(check(mid)) l=mid;
		else r=mid;
	}
	printf("%d\n",(int)floor(l*1000));
	return 0;
}