在一个2维平面上有两条传送带,每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。lxhgww在AB上的移动速度为P,在CD上的移动速度为Q,在平面上的移动速度R。现在lxhgww想从A点走到D点,他想知道最少需要走多长时间
这道题是一道三分,大家可以yy一下,本蒟蒻就不给出证明了(其实是因为不会)。
那接下来就是普通的三分套三分,注意一下精度问题,要开到1e-4,然后就没了吧。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define eps 1e-4
using namespace std;
double ax,ay,bx,by,cx,cy,dx,dy,Ak,Ab,Bk,Bb,A,B,C;
inline double getf(double x){return x*x;}
inline double getval(double xx,double yy,int id,int cnt)
{
double xy=xx*Ak+Ab,wy=yy*Bk+Bb;
if(id==1)xy=xx,xx=ax;
if(cnt==1)wy=yy,yy=cx;
double s1=sqrt(getf(xx-ax)+getf(xy-ay))/A;
double s2=sqrt(getf(yy-dx)+getf(wy-dy))/B;
double s3=sqrt(getf(xx-yy)+getf(xy-wy))/C;
return s1+s2+s3;
}
inline double sanfen2(double x,int id)
{
int cnt=0;
double l=min(cx,dx),r=max(cx,dx);
if(cx==dx)l=min(cy,dy),r=max(cy,dy),cnt=1;
while(r-l>=eps)
{
double lm=l+(r-l)/3.0,rm=r-(r-l)/3.0;
double wy1=getval(x,lm,id,cnt),wy2=getval(x,rm,id,cnt);
if(wy1>=wy2)l=lm;
else r=rm;
}
return getval(x,l,id,cnt);
}
double sanfen()
{
int id=0;
double l=min(ax,bx),r=max(ax,bx);
if(ax==bx)l=min(ay,by),r=max(ay,by),id=1;
while(r-l>=eps)
{
double lm=l+(r-l)/3.0,rm=r-(r-l)/3.0;
double wy1=sanfen2(lm,id),wy2=sanfen2(rm,id);
if(wy1>=wy2)l=lm;
else r=rm;
}
return sanfen2(l,id);
}
int main()
{
scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&ax,&ay,&bx,&by,&cx,&cy,&dx,&dy);
scanf("%lf%lf%lf",&A,&B,&C);
Ak=(ay-by)/(ax-bx);
if(ax==bx)Ak=0.0;Ab=ay-Ak*ax;
Bk=(cy-dy)/(cx-dx);
if(cx==dx)Bk=0.0;Bb=cy-Bk*cx;
printf("%.2lf\n",sanfen());
return 0;
}