#10017 「一本通 1.2 练习 4」传送带+三分套三分

题目描述

原题来自：SCOI 2010

在一个 2 维平面上有两条传送带，每一条传送带可以看成是一条线段。两条传送带分别为线段  AB和线段CD 。lxhgww 在  AB上的移动速度为 P ，在  CD上的移动速度为 Q ，在平面上的移动速度R 。

现在 lxhgww 想从A  点走到 D 点，他想知道最少需要走多长时间。

输入格式

输入数据第一行是4  个整数，表示  A和  B的坐标，分别为 ；;

第二行是 4 个整数，表示  和  的坐标，分别为 ，；;

第三行是  3个整数，分别是 。P,Q,R.

输出格式

输出数据为一行，表示 lxhgww 从A  点走到D  点的最短时间，保留到小数点后 2 位。

样例

样例输入

0 0 0 100
100 0 100 100
2 2 1

样例输出

136.60

描述：emmmm，看完之后，第一反应是暴力做，wa了几发之后，才晓得用三分套三分，果然不能轻视每一个算法啊。

分析：

给到四个点，我们可以分到最多三个线段。三分ab上的点到得到线段，将改点代入再三分cd上的点得到线段；用返回值（P,Q,R，分别为线段AB，线段CD，平地上的速度）,来判断两个点是不是最优点。

上代码：

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<algorithm>
using namespace std;
const double eps=1e-4;
int m,n,k,a1,a2,b1,b2,c1,c2,d1,d2;
double dfs(double x1,double y1,double x2,double y2)
{
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
double bfs(double x,double y)
{
    double cx=c1, cy=c2,dx=d1, dy=d2,ab=dfs(a1,a2,x,y)/m;
    while(fabs(dx-cx)>eps||fabs(dy-cy)>eps)
    {
        double x1=cx+(dx-cx)/3,x2=dx-(dx-cx)/3;
        double y1=cy+(dy-cy)/3,y2=dy-(dy-cy)/3;
        double s1=ab+dfs(x,y,x1,y1)/k+dfs(x1,y1,d1,d2)/n;
        double s2=ab+dfs(x,y,x2,y2)/k+dfs(x2,y2,d1,d2)/n;
        if(s1<=s2)
            dx=x2,dy=y2;
        else
            cx=x1,cy=y1;
    }
    return ab+dfs(x,y,cx,cy)/k+dfs(cx,cy,d1,d2)/n;
}
int main()
{
    scanf("%d%d%d%d",&a1,&a2,&b1,&b2);
    scanf("%d%d%d%d",&c1,&c2,&d1,&d2);
    scanf("%d%d%d",&m,&n,&k);
    double  ax=a1,ay=a2,bx=b1,by=b2;
    while(fabs(bx-ax)>eps||fabs(by-ay)>eps)
    {
        double x1=ax+(bx-ax)/3,x2=bx-(bx-ax)/3;
        double y1=ay+(by-ay)/3,y2=by-(by-ay)/3;
        if(bfs(x1,y1)<=bfs(x2,y2))
            bx=x2,by=y2;
        else
            ax=x1,ay=y1;
    }
    printf("%.2lf\n",bfs(ax,ay));
    return 0;
}