hdu 6097 Mindis（多校联赛）

#include <bits/stdc++.h>
#define eps 1e-8 //确定一个最小值 用于后续的比较
int main()
{
int T;
scanf("%d", &T);
while (T--)
{
double r, x1, y1, x2, y2;
scanf("%lf%lf%lf%lf%lf", &r, &x1, &y1, &x2, &y2);//输入半径和点的坐标等信息
double d0 = sqrt(pow(x1, 2) + pow(y1, 2));
if (fabs(d0) < eps)//p和q两个点无限接近于原点
{
printf("%.7f\n", 2 * r);//此时的pd+qd就相当与两条半径
continue;
}
double k = r * r / (d0 * d0);//用这个比值 确定p和q的反演点
double x3 = x1 * k, x4 = x2 * k, y3 = y1 * k, y4 = y2 * k;//x3 y3 x4 y4分别是p和q对应的反演点
double mx = (x3+x4)/2, my = (y3+y4)/2, ans;//mx和my是 p和q反演点连线的中点坐标（以下用连线中点代替）
double d = sqrt(pow(mx,2)+pow(my,2));//原点（圆心）到pq反演点的中点的距离
if (d <= r)//判断反演点和半径的关系 如果两个反演点的中点到圆心的距离小于半径
{//即反演点之间的连线与圆相交或者相切 此时dp+qd的最小时的d点即是线与圆的交点
double dist = sqrt(pow(x3 - x4, 2) + pow(y3 - y4, 2));//长度符合反演点与原来点的比例关系
ans = dist * d0 / r;
}
else//其他的即是连线与圆相离时的状态 这时候的d点是p和q的反演点的连线的中垂线与圆的交点
{
double kk = r / d;//找出半径和连线中点到圆心距离的比例 即是现在的d点坐标与连线中点坐标的比例
double smx = mx * kk, smy = my * kk;//求出连线中点的坐标
ans = 2 *sqrt(pow(smx - x1, 2) + pow(smy - y1, 2));//求出最后连线的长度 乘2是因为此时pd=qd
}
printf("%.7f\n", ans);//注意保留小数的位数
}
return 0;
}