【西南交大峨眉校区INT杯决赛 --- C. INT擅长捉弄的高木同学2】
Description
这天放学,西片和高木同学一起回家,在路上,他们捡到一张藏宝图,上面记录了寻宝的路线和宝藏的位置,于是他们就去寻宝了。经过不懈的努力,终于他们找到了宝藏位置的最后一个标志地点,一个圆形的花坛,宝藏就在花坛的周围。此时,西片心生一计,他刁难高木同学,让她计算一下他们到宝藏的最近的距离,请你帮助高木同学计算距离。
Input
输入第一行包含三个整数x,y,r(-105≤x,y≤105 ,0<r≤105)代表花坛圆心(x,y)和半径r。
第二行包含四个整数x1,y1,x2,y2(-105≤x1,y1,x2,y2≤105 )代表高木和西片的坐标(x1,y1)和宝藏坐标(x2,y2),且(x1,y1)和(x2,y2)不会在花坛内。
Output
输出他们的位置到宝藏的最短距离(不能穿过花坛,小数点后保留两位)。
Sample Input
1 1 1
3 0 0 0
Sample Output
3.00
解题思路
思路:判断线段与圆是否相交,三种情况。 设起点为 A,终点为 B,圆心为 O,则
- AB 与圆相交, AB 到圆心距离小于半径,最短距离为切线 AA’+BB’+弧 A’B’。
证明:若将 AA‘上移,则从 A 到达图上 A’的距离会变长,若将 AA’点右移,则会被圆挡住,
无法直线到达, 故 AA’最短,同理 BB’最短,得最短距离即为切线 AA’+BB’+弧 A’B’。
- 坑点: O 到 AB 距离小于半径但 AB 未与圆相交, 利用余弦定理判断角 A 或角 B 是否是钝角
即可排除此种情况。 - AB 与圆不相交且 O 到 AB 距离大于半径,直接计算 AB 距离即可
AC代码:
#include <iostream>
#include <algorithm>
#include <cmath>
#include <cstdio>
using namespace std;
double dis(double x1,double y1,double x2,double y2)
{
return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
int main()
{
double x,y,r,x1,y1,x2,y2;
cin >> x >> y >> r >> x1 >> y1 >> x2 >> y2;
double ao=dis(x1,y1,x,y);
double bo=dis(x2,y2,x,y);
double ab=dis(x1,y1,x2,y2);
double cosb=(ab*ab+bo*bo-ao*ao)/(2*ab*bo);
double cosa=(ab*ab+ao*ao-bo*bo)/(2*ab*ao);
double p=(ao+bo+ab)/2;
double s=sqrt(p*(p-ao)*(p-bo)*(p-ab));
double h=2*s/ab;
if(h>=r || h<r&&(cosa<0 || cosb<0))
printf("%.2f\n",ab);
else
{
double len1=sqrt(ao*ao-r*r);
double len3=sqrt(bo*bo-r*r);
double d1=acos(r/ao),d2=acos(r/bo),d3=acos((ao*ao+bo*bo-ab*ab)/(2*ao*bo));
double len2=(d3-d1-d2)*r;
printf("%.2f\n",len1+len2+len3);
}
return 0;
}