HDU 3968 Angry Birds Again（计算几何）

Description
给出小鸟的位置以及猪的位置，小鸟做斜抛运动经过ts后砸到猪，问小鸟经过的路程
Input
第一行为一整数T表示你用例组数，每组用例包括四个浮点数x0,y0,x1,y1,t分别表示小鸟和猪的横纵坐标以及小鸟砸到猪的用时，保证x0和x1不相等
Output
输出小鸟经过的路程，如果小鸟砸不到猪则输出Impossible.（重力加速度g=9.80m/s^2）
Sample Input
1
0 0 2 2 1
Sample Output
3.687
Solution
简单计算几何，首先计算出小鸟的轨迹方程，然后求出积分值即可
Code

#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        double x0,y0,x1,y1,t,g=9.8;
        scanf("%lf%lf%lf%lf%lf",&x0,&y0,&x1,&y1,&t);
        if(t==0)
        {
            printf("Impossible.\n");
            continue;
        }
        double vx=(x1-x0)/t,vy=(y1-y0+0.5*g*t*t)/t;
        double a=-0.5*g/(vx*vx),b=vy/vx;
        x1-=x0,x0=0;
        double t0=b,t1=2*a*x1+b;
        double l0=(t0*sqrt(1+t0*t0)+log(t0+sqrt(1+t0*t0)))/(4*a);
        double l1=(t1*sqrt(1+t1*t1)+log(t1+sqrt(1+t1*t1)))/(4*a);
        double ans=fabs(l1-l0);
        printf("%.3lf\n",ans);
    }
    return 0;
}