描述
题解
这个题官方题解是用反演搞的,但是我不是特别会,看到我 佐学姐
用 平面几何 + 代数方法
解的这个题感觉十分巧妙,也更加容易理解一些。
佐理慧’s blog,虽然学姐写的十分详细,但是我依然因为高中的数学知识都不知道丢哪儿了而看了好一会儿,数学这东西,说忘就忘啊。利用椭圆逼近获取最近距离即可。
代码
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;
typedef pair<double, double> point;
const double eps = 1e-8;
double r;
point p1, p2;
double get_dis(const point &P1, const point &P2)
{
double dx = P1.first - P2.first;
double dy = P1.second - P2.second;
return sqrt(dx * dx + dy * dy);
}
int main()
{
int T;
scanf("%d", &T);
while (T--)
{
scanf("%lf", &r);
scanf("%lf%lf", &p1.first, &p1.second);
scanf("%lf%lf", &p2.first, &p2.second);
double C = get_dis(p1, p2) / 2.0;
double D = get_dis(point(0, 0), p1);
if (D < eps)
{
printf("%.7f\n", r * 2.0);
continue;
}
double d = sqrt(D * D - C * C);
double a = C * r / D;
double y = d * (a * a - C * C) / (C * C);
if (y < r - d)
{
printf("%.7f\n", a * 2.0);
continue;
}
printf("%.7f\n", sqrt(C * C + (r - d) * (r - d)) * 2.0);
}
return 0;
}