Codeforces Round #465 (Div. 2) Fifa and Fafa

C. Fifa and Fafa

time limit per test

1 second

memory limit per test

256 megabytes

input

standard input

output

standard output

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.

Input

The single line of the input contains 5 space-separated integers R, x1, y1, x2, y2 (1 ≤ R ≤ 105, |x1|, |y1|, |x2|, |y2| ≤ 105).

Output

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.

Examples

input

Copy

5 3 3 1 1

output

3.7677669529663684 3.7677669529663684 3.914213562373095

input

Copy

10 5 5 5 15

output

5.0 5.0 10.0

 

 一道数学题

但是想练练模拟退火，万幸的A了，如果wrong answer了，就把m写大点儿~ 逃~~~~

#include<iostream>
using namespace std;
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<ctime>
#include<cmath>
typedef long long LL;
double T;
double R,x1,x2,yy1,y2;
const int M = 20;
double ax[M],ay[M],dr[M];
/*inline long long randt(){
    static long long seed = 2333;
    return seed = (long long)((((seed ^ 998244353) + 19260817ll) * 19890604ll) % 10000000000000007);
}*/
double Getrand(){
    double ans = rand();
//    cout<<ans<<endl;
    ans/=RAND_MAX;
    if(rand()%2)
        ans*=-1.0;
    ans*=3.145926535;
    return ans;
}
double get_dis(double x1,double yy1,double x2,double y2){
    double d1 = x1-x2;
    double d2 = yy1-y2;
    return sqrt(d1*d1+d2*d2);
}
double get(double X,double Y){
    double r1 = R-get_dis(x1,yy1,X,Y);
    double r2 = get_dis(x2,y2,X,Y);
    double r = r1;
    if(r>r2)
        r=r2;
    return r;
}
int main(){
    srand(time(0));
    scanf("%lf%lf%lf%lf%lf",&R,&x1,&yy1,&x2,&y2);
    for(int i=0;i<M;i++){
        ax[i]=Getrand()*R+x1;
        ay[i]=Getrand()*R+yy1;
        dr[i]=get(ax[i],ay[i]);
    }
    T=1.2*R;
    while(T>0.000009){
        for(int i=0;i<M;i++){
            for(int j=0;j<500;j++){
            double nx,ny,nd;
            double pp =Getrand() ,qq=Getrand();
            nx = cos(Getrand())*T+ax[i];
            ny = sin(Getrand())*T+ay[i];
            if(get_dis(x1,yy1,nx,ny)>R)
                continue;
            nd = get(nx,ny);
        //    cout<<nx<<" "<<ny<<" "<<nd<<" "<<pp<<" "<<qq<<" "<<T<<endl;
        //    getchar();
            if(nd>dr[i]){
                dr[i]=nd;
                ax[i]=nx;
                ay[i]=ny;
                continue;
            }
            else{
                 int p=rand()%10000;
                 if(p<=50){
                    dr[i]=nd;
                    ax[i]=nx;
                    ay[i]=ny;
                    continue;    
                 }
            }
            }
        }
        T*=0.9;
    }
    double ansx=ax[0],ansy=ay[0],ansr=dr[0];
    for (int i=1;i<M;i++){
        if(ansr<dr[i]){
            ansx=ax[i],ansy=ay[i],ansr=dr[i];
        }
    }
    printf("%.10f %.10f %.10f",ansx,ansy,ansr);
    return 0;
}