CF489E Hiking 01分数规划

一.题目

在坐标轴上有$n(1\le n\le 1000)$个景点，第i个景点坐标为$x_i(1\le x_i\le 10^6)$，美丽度为$b_i(1\le b_i\le 10^6)$。数据保证xi严格递增。一个人从坐标0出发，每天行进一段大于0的距离到达坐标更大的景点，他希望每天正好行进$L(1\le L\le 10^5)$距离，但他只能在景点停下。如果他某天的行进距离为r，那么当天的不满足度为$\sqrt{|L-r|}$，同时他会获得当天旅程结束所停在的景点的美丽度。最终他要停在n号节点。求使得$\frac{\sum \sqrt{|L-r|}}{\sum b}$取最小值的任意方案，其中$\sqrt{|L-r|}$表示每天的不满足度之和，$\sum b$表示每天停留的景点的美丽度之和。
传送门

二.Solution

这道题目是01分数规划的板题。
首先考虑二分，我们二分一个$mid$，使
$$\frac{\sum \sqrt{|L-r|}}{\sum b}\le mid$$
$$\to \sum \sqrt{|L-r|}-mid\ast \sum b\le 0$$
$$\to \sum \sqrt{|L-(x[i]-x[j])|}-mid\ast \sum b\le 0$$
于是我们写出一个方程式：
$$f[i]=f[j]+\sqrt{|L-(x[i]-x[j])|}-mid\ast b[i]$$
当$f[n]\le 0$时$mid$符合条件。
对于输出方案，我们用存前驱的方法，递归输出答案就行了。
初始值$f[0]=0$，其他为最小值。

三.Code

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
using namespace std;

#define M 1005

int n, pre[M], ans[M];
double L, x[M], b[M], f[M];

void print (int a){
	if (! a)
		return ;
	print (pre[a]);
	printf ("%d ", a);
}
bool check (double lim){
	memset (pre, 0, sizeof pre);
	f[0] = 0;
	for (int i = 1; i <= n; i ++){
		f[i] = -0x7f7f7f7f;
		for (int j = 0; j < i; j ++){
			if (f[j] + lim * b[i] - sqrt (fabs (L - (x[i] - x[j]))) > f[i]){
				f[i] = f[j] + lim * b[i] - sqrt (fabs (L - (x[i] - x[j])));
				pre[i] = j;
			}
		}
	}
	if (f[n] >= 0){
		for (int i = 0; i <= n; i ++)
			ans[i] = pre[i];
		return 1;
	}
	return 0;
}
int main (){
	scanf ("%d %lf", &n, &L);
	for (int i = 1; i <= n; i ++)
		scanf ("%lf %lf", &x[i], &b[i]);
	double l = 0, r = 1e9, mid;
	while (r - l > 1e-9){
		mid = (l + r) / 2.0;
		if (! check (mid))
			l = mid + 1e-9;
		else
			r = mid - 1e-9;
	}
	print (n);
	return 0;
}

Thanks！