analysis
题目有两个条件:
- 总重量至少为W
- 总才艺值与总重量的比值最大
由于出现了比值,这个题一定是01分数规划了
那么这个比值可以表示为
∑
i
=
1
n
C
i
×
x
i
W
i
×
x
i
\sum_{i=1}^{n}\frac{C_i\times x_i}{W_i\times x_i}
i=1∑nWi×xiCi×xi
且
∑
i
=
1
n
W
i
>
=
W
且\sum_{i=1}^{n}W_i>=W
且i=1∑nWi>=W
开始套用模板分析:
首先假设一个L,然后假设存在一组解x使得
∑ i = 1 n C i × x i W i × x i > L ∑ i = 1 n C i × x i > L × ∑ i = 1 n ( W i × x i ) ∑ i = 1 n C i × x i > ∑ i = 1 n L × ( W i × x i ) ∑ i = 1 n C i × x i − L × ( W i × x i ) > 0 ∑ i = 1 n ( C i − L × W i ) x i > 0 \sum_{i=1}^{n}\frac{C_i\times x_i}{W_i\times x_i}>L\\ \sum_{i=1}^{n}{C_i\times x_i}>L\times\sum_{i=1}^{n}({W_i\times x_i})\\ \sum_{i=1}^{n}{C_i\times x_i}>\sum_{i=1}^{n}L\times ({W_i\times x_i})\\ \sum_{i=1}^{n}{C_i\times x_i-L\times ({W_i\times x_i})}>0\\ \sum_{i=1}^{n}{(C_i-L\times W_i)x_i}>0\\ i=1∑nWi×xiCi×xi>Li=1∑nCi×xi>L×i=1∑n(Wi×xi)i=1∑nCi×xi>i=1∑nL×(Wi×xi)i=1∑nCi×xi−L×(Wi×xi)>0i=1∑n(Ci−L×Wi)xi>0
同理,如果任意解都有:
∑
i
=
1
n
C
i
×
x
i
W
i
×
x
i
<
=
L
\sum_{i=1}^{n}\frac{C_i\times x_i}{W_i\times x_i}<=L\\
i=1∑nWi×xiCi×xi<=L
则
∑
i
=
1
n
C
i
×
x
i
−
L
×
(
W
i
×
x
i
)
<
=
0
\sum_{i=1}^{n}{C_i\times x_i-L\times ({W_i\times x_i})}<=0\\
i=1∑nCi×xi−L×(Wi×xi)<=0
也就是要找出最大的 ( C i − L × W i ) (C_i-L\times W_i) (Ci−L×Wi)并且 ∑ W i > = W \sum W_i>=W ∑Wi>=W
为了满足总重量大于W(这是必须满足也一定可以满足的前提条件),我们可能需要选一些权值为负的牛来加大所有牛的总重量,但是如果仅仅考虑重量的话也不行,因为可能用几头重量更小的牛来加大重量对最后的答案的损失要比一头很重的牛少
这就必须用到背包了啊
相当于是有n个物品,每个物品i的贡献为 C i − L × W i C_i-L\times W_i Ci−L×Wi,体积为 W i W_i Wi,然后现在要装满一个体积为 W W W的背包,这个背包可以装满,甚至可以溢出,但是就只不能留有空间,求最大的贡献
上面那种背包用刷表法可以解决
由于背包体积可以超限使用,我们就把超过w的方案归在体积为 f [ w ] f[w] f[w]里面,这样子的话 f [ w ] f[w] f[w]就相当于是存了大于等于w的最后方案,于是最后只需要看f[w]就可以得到二分的答案了
code
#include<bits/stdc++.h>
using namespace std;
#define loop(i,start,end) for(register int i=start;i<=end;++i)
#define clean(arry,num) memset(arry,num,sizeof(arry))
#define anti_loop(i,start,end) for(register int i=start;i>=end;--i)
#define ll long long
template<typename T>void read(T &x){
x=0;char r=getchar();T neg=1;
while(r>'9'||r<'0'){if(r=='-')neg=-1;r=getchar();}
while(r>='0'&&r<='9'){x=(x<<1)+(x<<3)+r-'0';r=getchar();}
x*=neg;
}
const int maxn=300;
const int maxw=1000+10;
const int maxwi=1e6+10;
const int maxti=1e3;
int n,w;
int W[maxn],T[maxn];
double f[maxwi];
double beibao(double L){
loop(i,1,w)f[i]=-1e9;
loop(i,1,n){
anti_loop(j,w,0){
int k=min(j+W[i],w);
f[k]=max(f[k],f[j]+(double)T[i]-L*(double)W[i]);
}
}
return f[w];
}
void bin(){
double L=0,R=1000000,eps=1e-6;
while(R-L>eps){
double mid=(L+R)/2;
if(beibao(mid)>0)L=mid;
else R=mid;
}
printf("%d",(int)floor(L*1000));
}
int main(){
#ifndef ONLINE_JUDGE
freopen("datain.txt","r",stdin);
#endif
read(n);
read(w);
double maxx=-1;
loop(i,1,n)
read(W[i]),
read(T[i]),
maxx=max(maxx,(double)T[i]);
bin();
return 0;
}