传送门：HDU - 6559

题意

在一个三维空间，给定一个点和他的三维速度，给定一个圆锥，问这个点最早什么时候能撞上圆锥。

题解

本来一直想着怎么求圆锥的方程，然后....队友：这不是二分吗！然后问题就转换成了要怎么求当前时间是不是已经穿过了圆锥了....然后就gg了。

正解就是联立:

① r ' = (h-z)/h*r

② x=x0+vxt; y=y0+vyt; z=z0+vz*t;

③ x2+y2=r ' 2

然后就可以求解出 t 了（这个一元二次方程及其麻烦.....）不知道为啥直接输出t1就可以了，而不是输出非负的最小值....（逃

代码

#include<bits/stdc++.h>
#define eps 1e-12
using namespace std;

double r,h;
double x,y,z;
double vx,vy,vz;

int main()
{
    int t;
    scanf("%d",&t);
    int tt=0;
    while(t--){
        scanf("%lf%lf",&r,&h);
        scanf("%lf%lf%lf",&x,&y,&z);
        scanf("%lf%lf%lf",&vx,&vy,&vz);
        printf("Case %d: ",++tt);
        double a=(r*r/h/h*vz*vz-vx*vx-vy*vy);
        double b=-r*r/h/h*(2*h*vz-2*z*vz)-2*x*vx-2*y*vy;
        double c=-(x*x+y*y-r*r-r*r*z*z/h/h+2*h*z*r*r/h/h);
        double derta=b*b-4*a*c;
        double t1=(-b+sqrt(derta))/(2*a);
        double t2=(-b-sqrt(derta))/(2*a);
        printf("%.10f\n",t1);
    }
    return 0;
}