Dropping tests (01分数规划)

Dropping tests

题目

给一个二元组，可以扔掉这个二元组的k个元素，求取扔掉后的二元组的比值最大值。

思路

新学的01分数规划 写写题解理理思路。
想学习的参考这篇博客即可

• 首先把这个分式化成便于数学语言表达的形式，所以引入一个辅助符号函数$sgn$,取值只有0，1.
• 引入辅助参数r，显然$r_{m}ax$为答案，那就是说找到$r$的值使得这个方程成立，这个过程可以二分。
• 以上述的核心思想，对等式进行变形，由于要删去k个元素，在引入$sgn$函数后，可以得到最后要满足的等式为$d[i]\ast x[i]$，其中$d[i]=a[i]-r\ast b[i]$，$d[i]$的大小和r有关，并且是一个和r线性相关的单调函数，而它们整体则表示从$d[i]$中找到$n-k$个（它们显然是满足某种约束的）如果我们每次贪心的选出的前$n-k$个$d[i]$的求和大于0，这就说明r可以继续增大使这个函数变小，反之说明太大，即使用最大的$sum(d[i]...)$也无法使函数值大于0，所以要减小r
• 由上可以以$sum(d[i])$大于还是小于0为标志二分答案。

---

• 理解障碍还是在$sgn$这个辅助数组上，时刻紧扣定义的话，如果在不删数的时候，就是求$sum(d[i])=0$这玩意的解，普通的浮点二分即可，但是当要删数后，就要引入一个辅助数组来表达这个等式，也就是$d[i]\ast x[i]$，这时候我们选择前n-k大的$d[i]$ 来源于$d[i]$这个函数的单调性，和所求答案的性质，如果r大，d必然小，所以为了让r尽可能大，我们在一开始选择d的组件上，就尽可能往大了选。

AC代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<fstream>
#include<functional>
#define yes puts("yes");
#define inf 0x3f3f3f3f
#define linf 0x3f3f3f3f3f3f3f3f
#define ll long long
#define ull unsigned long long
#define debug(x) cout<<"> "<< x<<endl;
#define endl '\n'
#define lowbit(x) x&-x
//#define int long long
using namespace std;
typedef pair<int,int> PII;
const int N =10 + 1e5, mod = 1e9 + 7;
int n,k;
bool check(double mid,vector<int>&a,vector<int>&b){
	vector<double>d;
	for(int i=0;i<n;i++)
		d.push_back(a[i]-mid*b[i]);
	sort(d.rbegin(),d.rend());
	double sum = 0;
	for(int i=0;i<n-k;i++) sum += d[i];
	return sum>=0;
}
void solve()
{
	vector<int>a(n),b(n);
	for(int i=0;i<n;i++)cin>>a[i];
	for(int i=0;i<n;i++)cin>>b[i];
	double l=0.0,r=1.0;
	double eps = 1e-7;
	while(r-l>eps){
		double mid = (l+r)/2;
		if(check(mid,a,b)) l = mid;
		else r = mid;
	}    
	int res = 100*r +0.5;
	cout<<res<<endl;
}
signed main()
{
	ios::sync_with_stdio();cin.tie();cout.tie();
	while(cin>>n>>k&&(n||k))
	solve();

	return 0;
}