数轴line

输入：

10 4974
3636 3
6679 70
7182 93
10618 98
14768 23
15242 99
16077 35
23368 46
27723 15
32404 81

输出：

0.465308

类似于求最小平均值问题

首先二分答案mid

然后根据一个式子：∑w/∑a[i]<=mid

∑w-mid*a[i]<=0说明答案可以更小

这个过程用一个n^2的dp实现

f[i]=min(f[i],f[j]+dis(j,i)-mid*a[i])

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<cmath>
 6 using namespace std;
 7 double eps=0.0000001;
 8 double f[1001],x[1001],a[1001],L,inf;
 9 int n;
10 double abss(double p)
11 {
12   if (p<0) return -p;
13   else return p;
14 }
15 double check(double k)
16 {int i,j;
17   f[0]=0;
18   for (i=1;i<=n;i++)
19     {
20       f[i]=inf;
21       for (j=0;j<=i-1;j++)
22     if (f[i]>f[j]+sqrt(abss(x[i]-x[j]-L))-k*a[i])
23       f[i]=f[j]+sqrt(abss(x[i]-x[j]-L))-k*a[i];
24     }
25   return f[n];
26 }
27 int main()
28 {int i;
29   cin>>n>>L;
30   inf=1000000000.0;
31   for (i=1;i<=n;i++)
32     {
33       scanf("%lf%lf",&x[i],&a[i]);
34     }
35   double l=0,r=1000000.0,mid;
36   while (r-l>eps)
37     {
38       mid=(l+r)/2.0;
39       if (check(mid)<eps) r=mid;
40       else l=mid;
41     }
42   printf("%.6lf\n",r);
43 }