Cod sursa(job #2032097)

Utilizator RaduMirceaAndreiRadu Mircea Andrei RaduMirceaAndrei Data 4 octombrie 2017 15:52:22
Problema Adapost 2 Scor 90
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.18 kb
# include <fstream>
# include <cmath>
# include <iomanip>
# define DIM 50010
# define f first
# define s second
using namespace std;
ifstream fin("adapost2.in");
ofstream fout("adapost2.out");
int di[]={1,-1,0,0};
int dj[]={0,0,1,-1};
pair<double,double> v[DIM];
int n,i;
double x,y,xp,yp,val,a;
inline double distanta(double a,double b,double c,double d){
    return sqrt((a-c)*(a-c)+(b-d)*(b-d));
}
inline double dist(double x,double y){
    double sum=0;
    for(int i=1;i<=n;i++)
        sum+=distanta(x,y,v[i].f,v[i].s);
    return sum;
}
int main () {
    fin>>n;
    for(int i=1;i<=n;i++){
        fin>>v[i].f>>v[i].s;
        x+=v[i].f;
        y+=v[i].s;
    }
    x/=n;
    y/=n;
    val=dist(x,y);
    for(double step=1024;step>0.0001;){
        xp=yp=-1000;
        for(int i=0;i<4;i++){
            a=dist(x+step*di[i],y+step*dj[i]);
            if(val>a){
                xp=x+step*di[i];
                yp=y+step*dj[i];
                val=a;
            }
        }
        if(xp==-1000)
            step/=2;
        else{
            x=xp;
            y=yp;
        }
    }
    fout<<setprecision(4)<<fixed<<x<<" "<<y<<"\n";
    return 0;
}