Fifa and Fafa are sharing a flat. Fifa loves video games and wants to download a new soccer game. Unfortunately, Fafa heavily uses the internet which consumes the quota. Fifa can access the internet through his Wi-Fi access point. This access point can be accessed within a range of r meters (this range can be chosen by Fifa) from its position. Fifa must put the access point inside the flat which has a circular shape of radius R. Fifa wants to minimize the area that is not covered by the access point inside the flat without letting Fafa or anyone outside the flat to get access to the internet.
The world is represented as an infinite 2D plane. The flat is centered at (x1, y1) and has radius R and Fafa's laptop is located at (x2, y2), not necessarily inside the flat. Find the position and the radius chosen by Fifa for his access point which minimizes the uncovered area.
The single line of the input contains 5 space-separated integers R, x1, y1, x2, y2 (1 ≤ R ≤ 105, |x1|, |y1|, |x2|, |y2| ≤ 105).
Print three space-separated numbers xap, yap, r where (xap, yap) is the position which Fifa chose for the access point and r is the radius of its range.
Your answer will be considered correct if the radius does not differ from optimal more than 10 - 6 absolutely or relatively, and also the radius you printed can be changed by no more than 10 - 6 (absolutely or relatively) in such a way that all points outside the flat and Fafa's laptop position are outside circle of the access point range.
5 3 3 1 1
3.7677669529663684 3.7677669529663684 3.914213562373095
10 5 5 5 15
5.0 5.0 10.0
题意:已知一个圆,点P可能在圆内圆外,求最大的圆T,使得圆T不包含P.
思路:
分类讨论 :若P在圆内,明显使得直径最大的圆即可, 在圆外就是该圆面积.半径很容易算.坐标用向量法[高中知识]
补充一下过程:
(x2-x1,y2-y1)/(x-x1,y-y1)=dis/r =k
化简一下得到 x=(x2-x1)/k+x1 y=(y2-y1)/k+y1
#include <bits/stdc++.h>
#pragma comment(linker, "/STACK:102400000,102400000")
#define INF 0x3f3f3f3f
#define MOD 1000000007
#define bug1 cout <<"bug1"<<endl
#define bug2 cout <<"bug2"<<endl
#define bug3 cout <<"bug3"<<endl
using namespace std;
typedef long long ll;
const int MAX_N=1e5+5;
int main(void){
double R,x1,y1,x2,y2;// x1,y1是P点坐标 x2,y2是家的圆心坐标
cin >>R>>x2>>y2>>x1>>y1;
if(x1==x2&&y1==y2){// 重合的特殊情况
double r=(double)R/2;
double x=x1;
double y=y1+r;
printf("%.12lf %.10lf %.10lf\n",x,y,r);
return 0;
}
double dis=sqrt((y2-y1)*(y2-y1)+(x2-x1)*(x2-x1));
if(dis-R>0){
printf("%.12lf %.12lf %.12lf\n",x2,y2,R);
return 0;
}
double r=(dis+R)/2;
double k=dis/r;
double x=(double)(x2-x1)/k+(double)x1;
double y=(double)(y2-y1)/k+(double)y1;
printf("%.12lf %.12lf %.12lf\n",x,y,r);
}