2018CCPC吉林E The Tower

题意

给你一个锥心位于(0,0,0)的圆锥体 给出r h

再给你一个人的坐标(x0,y0,z0)和他三个方向的速度(vx,vy,vz) 问你人和锥体相撞的最短时间是？

思路

列出各种不等式 然后各种代换解二元一次方程 这玩意儿是真难算啊 过程如图（字丑 写给自己随意看的）：

代码

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>

using namespace std;

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

int main()
{
	int n;
	scanf("%d",&n);
	int cnt = 0;
	while(n--){
		cnt++;
		scanf("%lf %lf",&r,&h);
		scanf("%lf %lf %lf",&x,&y,&z);
		scanf("%lf %lf %lf",&vx,&vy,&vz);
		double K = h - z;
		double C = r * r / h / h;
		double a = vx * vx + vy * vy - C * vz * vz;
		double b = 2.0*(vx*x + vy*y + C * K * vz);
		double c = x * x + y * y - C * K * K;
		double t = (-1.0*b - sqrt(b*b - 4.0*a*c)) / (2.0*a);
		printf("Case %d: %.10f\n",cnt,t);
	}
	return 0;
}