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;
}