01分数规划

给定 $n$ 个二元组 $(valu{e}_i，cos{t}_i)$，分别表示价值和代价
设 x i ∈ { 0 , 1 } x_i∈\{0, 1\} xi​∈{0,1}，表示第 $i$ 个二元组是否选取
要求最大(小)化下式：
$$r=\frac{\sum valu{e}_i\times x_i}{\sum cos{t}_i\times x_i}$$

---

二分 $mid$ ，若最优值 $\ge mid:$
$$\frac{\sum valu{e}_i\times x_i}{\sum cos{t}_i\times x_i}\ge mid$$
$$\sum (valu{e}_i-mid\ast cos{t}_i)\times x_i\ge 0$$
也就是找到一组 { x i } \{x_i\} {xi​} 使得 $\sum (valu{e}_i-mid\ast cos{t}_i)$ 最大
若找到的最大值 $\ge 0$，也就是 最优值大于 $mid$，反之就是小于 $mid$

详细证明请见01分数规划问题相关算法与题目讲解（二分法与Dinkelbach算法）

---

例题一：POJ 2976 Dropping tests

选取 $n-k$ 个二元组使得 $\frac{\sum a_i}{\sum b_i}$ 最大

#include<bits/stdc++.h>
#define rint register int
#define deb(x) cerr<<#x<<" = "<<(x)<<'\n';
using namespace std;
typedef long long ll;
typedef pair <int,int> pii;
const double eps = 1e-5;
const int maxn = 1e3 + 5;
int n, k, a[maxn], b[maxn];
double c[maxn];

bool ck(double x){
	for(int i=1; i<=n; i++) c[i] = a[i] - x * b[i];
	sort(c+1, c+1+n);
	double res = 0;
	for(int i=k+1; i<=n; i++) res += c[i];
	return res >= eps;
}

signed main() {
	while(scanf("%d%d", &n, &k) && n){
		for(int i=1; i<=n; i++) scanf("%d", a+i);
		for(int i=1; i<=n; i++) scanf("%d", b+i);
		double l = 0, r = 1e5, mid;
		while((r - l) >= eps){
			mid = (l + r) / 2;
			if(ck(mid)) l = mid + eps;
			else r = mid - eps;
		}
		printf("%.0f\n", r * 100);
	}
}

---

例题二：P4322 [JSOI2016]最佳团体

选取 $k$ 个二元组使得 $\frac{\sum p_i}{\sum s_i}$ 最大
同时选取一个二元组有可能有一个必须先选的依赖二元组

因此二分后，用 $n^2$ 的树形背包，每个点的点权为 $\sum (p_i-mid\ast s_i)$
我看题解里有在 $dfs$ 里 $dp$，然后说是 $n^2$ 的，但是代码却写成 $n^3$ 。。。
树形背包 $n^2$ 的写法，合并子树时，必须先 $dp$，然后再加上这个子树的 $size$
或者在 $dfs$ 序上 $dp$ 也可，常数更小

#include<bits/stdc++.h>
#define rint register int
#define deb(x) cerr<<#x<<" = "<<(x)<<'\n';
using namespace std;
typedef long long ll;
typedef pair <int,int> pii;
const int maxn = 2505;
const double eps = 1e-5;
int n, m, s[maxn], p[maxn], r[maxn];
int dfn[maxn], tot, nxt[maxn];
double dp[maxn][maxn], v[maxn];
vector <int> g[maxn];

void dfs(int u){
	dfn[u] = tot++;
	for(auto v : g[u]) dfs(v);
	nxt[dfn[u]] = tot;
}

inline double ck(double x){
	for(rint i=1; i<=n+1; i++)
		for(rint j=0; j<=m; j++)
			dp[i][j] = -1e9;
	for(rint i=1; i<=n; i++) v[dfn[i]] = p[i] - s[i] * x;
	for(rint i=0; i<=n; i++)
		for(rint j=0; j<=min(m, i); j++){
			dp[i+1][j+1] = max(dp[i+1][j+1], dp[i][j] + v[i]);
			dp[nxt[i]][j] = max(dp[nxt[i]][j], dp[i][j]);
		}
	return dp[n+1][m];
}

signed main() {
	scanf("%d%d", &m, &n); m++;
	for(rint i=1; i<=n; i++) {
		scanf("%d%d%d", s+i, p+i, r+i);
		g[r[i]].push_back(i);
	}
	dfs(0);
	double l = 0, r = 1e4, mid;
	while((r - l) >= eps){
		mid = (l + r) * 0.5;
		if(ck(mid) >= eps) l = mid + eps;
		else r = mid - eps;
	}
	printf("%.3f\n", r);
}