GYM 100971 F.Two Points（三分）

来源：互联网发布：淘宝晚班客服兼职编辑：程序博客网时间：2024/05/05 05:31

Description
给出两个点的初始二维坐标和速度向量，问这两个点运动过程中的最短距离
Input
第一行四整数ax,ay,bx,by表示两个点的初始二维坐标，第二行四个整数vax,vay,vbx,vby表示两个点的速度向量
(-1e4<=ax,ay,bx,by,vax,vay,vbx,vby<=1e4))
Output
输出两个点运动过程中的最短距离
Sample Input
1 1 2 2
0 0 -1 0
Sample Output
1.000000000000000
Solution
三分时间即可
Code

#include<cstdio>#include<iostream>#include<cstring>#include<algorithm>#include<cmath>#include<vector>#include<queue>#include<map>#include<set>#include<ctime>using namespace std;typedef long long ll;#define INF 0x3f3f3f3f#define maxn 1111#define eps 1e-8struct node{    double x,y;    node operator*(const double &n)const    {        node a;        a.x=x*n,a.y=y*n;        return a;    }    node operator+(const node &b)const    {        node a;        a.x=x+b.x,a.y=y+b.y;        return a;    }}a,b,va,vb;int sign(double x){    if(fabs(x)<eps)return 0;    if(x>eps)return 1;    return -1;}double dis(node a,node b){    double x=b.x-a.x,y=b.y-a.y;    return sqrt(x*x+y*y);}int main(){    while(~scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y,&va.x,&va.y,&vb.x,&vb.y))    {        double l=0,r=1e12,mid1,mid2,ans=dis(a,b);        int T=100000;        while(T--)        {            mid1=0.5*(l+r),mid2=0.5*(mid1+r);            double d1=dis(a+va*mid1,b+vb*mid1),d2=dis(a+va*mid2,b+vb*mid2);            ans=min(ans,d1),ans=min(ans,d2);            if(sign(d1-d2)<=0)r=mid2;            else l=mid1;        }        printf("%.10f\n",ans);    }    return 0;}

0 0