HDU 6559 The Tower（一元二次方程）

大致题意

给一个圆锥，底部圆心在(0,0,0)，高为 H , 底部圆半径为 R ，再给一个质点 (x0,y0,z0)，给定其单位时间的速度速度 (vx,vy,vz)，保证质点一定能撞上圆锥，求质点撞上圆锥表面的最小单位时间。

思路

感觉列个式子就是解一个一元二次方程，首先当高度为 h，圆锥的横截圆半径通过相似三角形易得 r = R x (H - h)/H 。 圆上的点有性质 x ^ 2 + y ^ 2 = r^2
设 t 时刻质点的位置为 (x0 + vx t , y0 + vy t , z0 + vz t)，带入上述两式，解一元二次方程，求较小的解，不必担心无解。

代码

代码不能，主要就是推公式算出 A , B , C 即可。

#include<bits/stdc++.h>
using namespace std;
#define maxn 1000005
#define maxm 1000006
#define ll long long int
#define INF 0x3f3f3f3f
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
#define mem(a) memset(a,0,sizeof(a))
#define sqr(x) (x*x)
#define inf (ll)2e18+1
#define mod 998244353
#define db double
int read(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
     while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
     return f*x;
}
int T;
int main()
{
    db R,H,x0,y0,z0,vx,vy,vz;
    T=read();
    int kase=0;
    while(T--){
        scanf("%lf%lf",&R,&H);
        scanf("%lf%lf%lf%lf%lf%lf",&x0,&y0,&z0,&vx,&vy,&vz);
        db a,b,c;
        a=vx*vx+vy*vy-(R*R*vz*vz/(H*H));
        b=2.0*x0*vx+2.0*y0*vy-2.0*(R*R*z0*vz/(H*H))+2.0*(R*R*vz/H);
        c=x0*x0+y0*y0-(R*R*z0*z0/(H*H))+2.0*(R*R*z0/H)-R*R;
        db delt=b*b-4.0*a*c;
        delt=sqrt(delt);
        db x1,x2;
        x1=(-b+delt)/(2.0*a);
        x2=(-b-delt)/(2.0*a);
        printf("Case %d: %lf\n",++kase,x2);
    }
    return 0;
}