【zzulioj 2127】 tmk射气球

比较简单的题，直接求空间中一个点到直线的距离而已，这道题说了直线和水平的平面 平行，我们可以先求投影到直线的距离，然后再算当前点到直线的距离。

Description

有一天TMK在做一个飞艇环游世界，突然他发现有一个气球匀速沿直线飘过，tmk想起了他飞艇上有一把弓，他打算拿弓去射气球，为了提高射击的准确性，他首先在飞艇上找到一个离气球最近的一个点，然后射击（即使气球在飞船的正上方），现在求某些时刻飞艇上的点和气球的距离最小是多少（这个最小距离我们简称为飞艇到气球的距离）。

Input

第一行一个整数T(T<=20)，表示有T组测试数据

每组测试数据，有两行。

第一行有5个整数，h,x1,y1,x2,y2，其中h表示飞船的高度，飞船可抽象为一个线段，(x1,y1)(x2,y2)分别是这个线段的端点（有可能会有(x1,y1)(x2,y2)重合的情况）

第二行有6个整数，x,y,z,X,Y,Z分别表示气球的在第0秒的时候的横坐标，纵坐标，高度，一秒时间气球横坐标的变化量，一秒时间气球纵坐标的变化量，一秒时间气球高度的变化量(如果现在气球在(x0,y0,z0)下一秒坐标就为(x0+X,y0+Y,z0+Z))

第三行1个整数n，表示询问组数

接下来的n行，每行一个整数，表示询问的秒数t

题目涉及的整数除了T以外，范围均为[0,1000]

Output

每组询问输出n行，每行输出一个数，表示在t秒的时候飞艇与气球的距离最小是多少，保留两位小数

Sample Input

1 1 1 1 2 2 0 0 0 4 4 4 2 0 3

Sample Output

1.73 17.92

 

 #include <iostream>
 #include<string>
 #include<cstring>
 #include<vector>
 #include<set>
 #include<map>
 #include<stack>
 #include<queue>
 #include<list>
 #include<cmath>
 #include<algorithm>
 #include<cstdio>
 #include <bitset>
 #include <climits>
 #include <time.h>
 #include<iomanip>
 #define INF 0x3f3f3f3f
 using namespace std;
 #define ll long long
 #define eps 1e-8
 #define zero(x) (((x>0)?(x):-(x))<eps)
 struct point{double x,y;};

 double distancen(point p1,point p2)
 {
     return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
 }
 ll xmult(point p1,point p2,point p0)
 {
     return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);

 }
 point intersection(point u1,point u2,point v1,point v2)
 {
     point ret = u1;
     double t = ((u1.x-v1.x)*(v1.y-v2.y)-(u1.y-v1.y)*(v1.x-v2.x))
                 /((u1.x-u2.x)*(v1.y-v2.y)-(u1.y-u2.y)*(v1.x-v2.x));
     ret.x+=(u2.x-u1.x)*t;
     ret.y+=(u2.y-u1.y)*t;
     return ret;

 }

 point ptoseg(point p,point l1,point l2)
 {
     point t = p;
     t.x+=l1.y-l2.y;
     t.y+=l2.x-l1.x;

     if(xmult(l1,t,p)*xmult(l2,t,p)>eps)
     {
         return (distancen(p,l1)<distancen(p,l2))?l1:l2;
     }
     return intersection(p,t,l1,l2);
 }

 int dot_online(point p,point l1,point l2)
 {
     return (zero(xmult(p,l1,l2))&&(l1.x-p.x)*(l2.x-p.x)==&&(l1.y-p.y)*(l2.y-p.y)==);
 }

 int main()
 {
     //freopen("D://in.txt","r",stdin);
     //freopen("D://out.txt","w",stdout);
     int T;
     scanf("%d",&T);
     while(T--)
     {
         bool mark = false;
         point ll1,ll2;
         double h,x1,y1,x2,y2;
         scanf("%lf%lf%lf%lf%lf",&h,&x1,&y1,&x2,&y2);
         ll1.x = x1;ll1.y = y1;
         ll2.x = x2;ll2.y = y2;
         if(x1==x2&&y1==y2)
             mark = true;

         double x,y,z,X,Y,Z;
         scanf("%lf%lf%lf%lf%lf%lf",&x,&y,&z,&X,&Y,&Z);
         int n;
         scanf("%d",&n);
         while(n--)
         {
             int t;
             scanf("%d",&t);
             point newp;
             newp.x = x+t*X;
             newp.y = y+t*Y;
             double zz = z+t*Z;
             if(mark)
             {
                 double res = distancen(newp,ll1);
                 double ans = sqrt((h-zz)*(h-zz)+res*res);
                 printf("%.2lf\n",ans);
             }
             else if(dot_online(newp,ll1,ll2))
             {
                 double ans = fabs(h-zz);
                 printf("%.2lf\n",ans);
             }
             else
             {
                 point pseg = ptoseg(newp,ll1,ll2);
                 double res = distancen(pseg,newp);
                 double ans = sqrt(res*res+(zz-h)*(zz-h));
                 printf("%.2lf\n",ans);
             }
         }
     }
     return ;
 }