题意:给出一个轮子,上面有一个随着它转动的传感器在圆周上,给出一个指定距离m,和轮子向前行进的速度v以及轮子的半径r。问让传感器通过该距离最少需要多少时间。
分析:首先我们列出传感器行进距离与时间的轮子行进距离的关系:f(c) = c+r*sin(c/r)。其中c是距离,r是半径。该公式表示传感器从轮子正上方开始运行索性进的距离。
我们可以先去除一种情况,就是距离m是直径2r的整数倍,那么传感器在开始点和结束点所处的角度是相同的,无论从哪个角度开始所需的时间都是一样的。
我们只考虑非整数倍的情况。
要看传感器行进速度的变化,我们需要对这个函数f关于c进行求导。
求得导数是:f'(c) = 1+cos(c/r)。这就是传感器速度的函数。而这个函数图像与x轴所夹的面积就是行进距离。
那么,我们现在的任务就变成了在这个图像上截取一段,使得与x轴所夹的面积是m,且要保证我们截取的横向长度最短。(因为速度是一定的,距离和时间成正比)
假设我们截取了x=a到x=b。那么一定有f'(a)=f'(b)。如果两者不等,例如f'(a)<f'(b),我只要将右边缘向右平移一点点,保证面积不变的话,就需要将左边缘向右平移更多。
因为我们如果把这个增量看成小长方体,如果左右平移量相等,f'(a)<f'(b)导致右边进来的就比左边进来的长方体高,那面积就变了。
既然f'(a)=f'(b),就一定a在图像的递增区域,b在递减区域,否则还会出现上面的情况,或者就是m是2r的整数倍。
也就是说两边缘必须对称,也就是说选区的中间是最高点或者最低点。我们只需要对两种情况分别进行二分查找即可。
二分查找的时候查找轮子的行进长度,每次观察f(c)>=m/2是否成立。
二分查找有个需要特别注意的地方。如果每次只是限定L和R的差值是不行的。差值设大了就会出现WA,设小了会TLE。
需要限定二分迭代的次数为50次左右。不知道是不是其他的比赛或者OJ也有这种情况。
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std; #define d(x) #define zero(x) (((x)>0?(x):-(x))<eps)
#define eps 1.0E-8 int double_cmp(double a)
{
if (zero(a))
return ;
return a > ? : -;
} int n, radius, v;
int dist; bool top_ok(double bike_dist)
{
double ret = bike_dist + radius * sin(bike_dist / radius);
return double_cmp(ret * 2.0 - dist) >= ;
} bool bottom_ok(double bike_dist)
{
double ret = bike_dist - radius * sin(bike_dist / radius);
return double_cmp(ret * 2.0 - dist) >= ;
} double binary_search(bool (*ok)(double))
{
double l = max(0.0, dist / 2.0 - * radius);
double r = dist / 2.0 + * radius;
int cnt = ;
while (double_cmp(r - l) > )
{
double mid = (l + r) / ;
if (ok(mid))
r = mid;
else
l = mid;
cnt++;
if (cnt > )
break;
}
d(printf("%.2f\n", l / radius / / 3.14));
return l;
} double work()
{
double time_top = binary_search(&top_ok) / v;
double time_bottom = binary_search(&bottom_ok) / v;
d(printf("%.2f %.2f\n", time_top, time_bottom));
return min(time_top, time_bottom) * ;
} int main()
{
scanf("%d", &n);
scanf("%d%d", &radius, &v);
for (int i = ; i < n; i++)
{
int a, b;
scanf("%d%d", &a, &b);
dist = b - a;
printf("%.12f\n", work());
}
return ;
}