题目
题目大意
一个数轴上有n个点,每个点的坐标均为正整数,每个点有一个美丽值。
给定常数p,从原点出发,只能向正方向移动,每次必须停在一个点上,假设经过的距离为d,那么这次hike的疲劳值为 a b s ( d − p ) \sqrt{abs(d-p)} abs(d−p)。现在要从原点前往最左端的点,求总疲劳值/总美丽值的最小值。
总美丽值为停留的点的美丽值之和。
题目分析
分数规划
注意到求两个和式的比值,并且以一个点为终点就要对疲劳值、美丽值都产生贡献。
所以可以考虑01分数规划
1.二分一个答案,然后令新点权为 上个点到这里的距离-二分答案*该点美丽值
2.dp求解判断是否合法
关于01分数规划
动态规划
现在问题就转化为了在跳跃的过程中获得最小点权。
考虑最简单的动态规划: d p [ i ] dp[i] dp[i]表示到i号点的时候的最小点权。
那么 d p [ i ] = m i n ( d p [ k ] + c [ i ] ) dp[i]=min(dp[k]+c[i]) dp[i]=min(dp[k]+c[i])
其中c表示点权, k ∈ ( 1 , i − 1 ) k∈(1,i-1) k∈(1,i−1)
代码
下附AC代码:
#include<bits/stdc++.h>
using namespace std;
int read(){
char s;
int x=0,f=1;
s=getchar();
while(s<'0'||s>'9'){
if(s=='-')f=-1;
s=getchar();
}
while(s>='0'&&s<='9'){
x*=10;
x+=s-'0';
s=getchar();
}
return x*f;
}
const int N=1010;
int n,len;
int x[N],c[N];
double dp[N];
int pre[N];//记录答案路径
bool check(double limit){
for(int i=1;i<=n;i++){
dp[i]=1e9;
for(int j=0;j<i;j++){
double tmp=dp[j]+sqrt(abs(len-(x[i]-x[j])))-limit*c[i];
if(tmp<dp[i]){
pre[i]=j;
dp[i]=tmp;
}
}
}
return dp[n]>0;
}
int ans[N];
int main(){
n=read(),len=read();
for(int i=1;i<=n;i++){
x[i]=read(),c[i]=read();
}
double l=0,r=1e10;
while(r-l>1e-8){
double mid=(l+r)/2;
if(check(mid))l=mid;
else r=mid;
}
check(l);
int cnt=1;
int id=n;
while(id>0){
ans[cnt++]=id;
id=pre[id];
}
for(int i=cnt-1;i>0;i--){
printf("%d ",ans[i]);
}
}