GYM 100971 F.Two Points（三分）

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-8
struct 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;
}