HDU-2017 多校训练赛6-1002-Mindis

ACM模版

描述

题解

这个题官方题解是用反演搞的，但是我不是特别会，看到我 佐学姐 用 平面几何 + 代数方法 解的这个题感觉十分巧妙，也更加容易理解一些。

佐理慧’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;}