HDU 6559The Tower(计算几何)

http://acm.hdu.edu.cn/showproblem.php?pid=6559

题意：一个圆锥的地面圆的圆心在原点，在给定一个质点和一个初速度，问过多久能与圆锥相交。

做法：设交点为（x,y,z)根据参数方程

在根据几何关系，比如相似三角形可以得到一些关系，关系不同，可能有一些特殊情况，我觉得下面这个方程最简单

然后得到一个关于t的二次方程，暴力求解就可以了。

#include "bits/stdc++.h"

using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const ull hash1 = 201326611;
const ull hash2 = 50331653;
const int N = 100000 + 10;
const int M = 20;
const double eps = 1e-8;

int dcmp(double x) {
    if (fabs(x) < eps) return 0;
    return (x > 0) ? 1 : -1;
}

int main() {
    int T;
    double r, h, x, y, z, vx, vy, vz;
    cin >> T;
    for (int Case = 1; Case <= T; Case++) {
        cin >> r >> h >> x >> y >> z >> vx >> vy >> vz;
        double a = (vx * vx + vy * vy) / (r * r) - vz * vz / (h * h);
        double b = (2 * vx * x + 2 * vy * y) / (r * r) - 2 * vz * z / (h * h) + 2 * vz / h;
        double c = (x * x + y * y) / (r * r) - z * z / (h * h) + 2 * z / h - 1;
        double d = b * b - 4 * a * c;
        double t1 = (-b - sqrt(d)) / (2 * a);
        double t2 = (-b + sqrt(d)) / (2 * a);
        if (dcmp(t1) <= 0) printf("Case %d: %.8f\n", Case, t2);
        else printf("Case %d: %.8f\n", Case, t1);
    }
    return 0;
}