BZOJ 5281--[Usaco2018 Open]Talent Show(分数规划&单调队列&DP)

5281: [Usaco2018 Open]Talent Show

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 79  Solved: 58
[Submit][Status][Discuss]

Description

FarmerJohn要带着他的N头奶牛,方便起见编号为1…N,到农业展览会上去,参加每年的达牛秀!他的第i头奶牛重量为wi,才艺水平为ti,两者都是整数。在到达时,FarmerJohn就被今年达牛秀的新规则吓到了:
 
(一)参加比赛的一组奶牛必须总重量至少为W
 
(这是为了确保是强大的队伍在比赛,而不仅是强大的某头奶牛),并且
 
(二)总才艺值与总重量的比值最大的一组获得胜利。
 
FJ注意到他的所有奶牛的总重量不小于W,所以他能够派出符合规则(一)的队伍。帮助他确定这样的队伍中能够达到的最佳的才艺与重量的比值。
 

Input

输入的第一行包含N和W。下面N行,每行用两个整数wi和ti描述了一头奶牛。
1≤N≤250
1≤W≤1000
1≤wi≤10^6
1≤ti≤10^3
 

Output

请求出Farmer用一组总重量最少为W的奶牛最大可能达到的总才艺值与总重量的比值。
如果你的答案是A,输出1000A向下取整的值,以使得输出是整数
(当问题中的数不是一个整数的时候,向下取整操作在向下舍入到整数的时候去除所有小数部分)。
 

Sample Input

3 15
20 21
10 11
30 31

Sample Output

1066
在这个例子中,总体来看最佳的才艺与重量的比值应该是仅用一头才艺值为11、重量为10的奶牛,但是由于我们需
要至少15单位的重量,最优解最终为使用这头奶牛加上才艺值为21、重量为20的奶牛。这样的话才艺与重量的比值
为(11+21)/(10+20)=32/30=1.0666666...,乘以1000向下取整之后得到1066。
 

题目链接:

    http://www.lydsy.com/JudgeOnline/problem.php?id=5281 

Solution

  一看题面就知道要01分数规划。。。

  假设答案为c,式子比较显然 ti的总和 / wi的总和 >= c

  ti的总和 >=  wi的总和 * c

  然后就是二分c验证正确性。。

  将每只牛的价值赋值为 ti - wi * c

  显然如果有一种方案使得总的牛的价值不小于0,就说明c存在合法方案。。

  由于W<=1000,我们可以直接DP。。

  f [ i ] 表示wi总和不小于i时的最大价值总和。。

  f [ i ] = max(f[j]+val)其中val表示当前牛的价值,并且i-j<=w

  由于i和j的关系,要用单调队列维护。。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<map>
#include<cmath>
#include<set>
#define pa pair<LL,LL>
#define LL long long
#define ept 1e-5
using namespace std;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
inline void Out(int a){
    if(a>9) Out(a/10);
    putchar(a%10+'0');
}
const double inf=1e9;
const LL mod=1e9+7;
const int N=300;
int n,m,cnt=1,W;
struct cow{
	int w,s;
}a[N];
double f[2000],t[2000];
int q[2000];

int main(){
    n=read();W=read();
	for(int i=1;i<=n;++i){
		a[i].w=read();a[i].s=read();
	}
	double l=0,r=10000,mid,y;
	int x,L=1,R=0;
	while(fabs(r-l)>ept){
		//cout<<l<<" "<<r<<endl;
		mid=(l+r)/2;
		for(int i=0;i<=W;++i) f[i]=-inf;
		f[0]=0;
		for(int i=1;i<=n;++i){
			x=a[i].w<=W?a[i].w:W;
			y=(double)a[i].s-(double)mid*(double)a[i].w;
			L=1;R=1;q[1]=0;
			for(int j=1;j<=W;++j){
				while(L<=R&&f[j]>f[q[R]]) --R;
				q[++R]=j;
				while(j-q[L]>x) ++L;
				t[j]=f[q[L]]+y;
			}
			for(int j=1;j<=W;++j)
				if(t[j]>f[j]) f[j]=t[j];
		}
		if(f[W]>=0) l=mid;
		else r=mid; 
	}
	l=l*1000;
	int ans=l;
	printf("%d\n",ans);
	return 0;
}

  

  

This passage is made by Iscream-2001.

 

posted @ 2018-09-20 13:45 Iscream-2001 阅读(...) 评论(...) 编辑 收藏