Cod sursa(job #930415)

Utilizator DEYDEY2Tudorica Andrei DEYDEY2 Data 27 martie 2013 17:12:25
Problema Adapost 2 Scor 90
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.89 kb
#include<iostream>
#include<fstream>
#include<cmath>
using namespace std;
ifstream f("adapost2.in");
ofstream g("adapost2.out");
const int dx[]={0,1,0,-1},dy[]={1,0,-1,0};
int n,nr=38;
double x[50010],y[50010],px,py,dmin=1000000000,t,pas=512;

inline double dist(double pozx, double pozy)
{
    double rez = 0;
    for(int i=1;i<=n;++i) rez+=sqrt((pozx-x[i])*(pozx-x[i])+(pozy-y[i])*(pozy-y[i]));
    return rez;
}

int main()
{
    f>>n;
    for(int i=1;i<=n;++i){f>>x[i]>>y[i]; px+=x[i]; py+=y[i];}
    px/=n; py/=n;
    int nr;
    if(n<10000) nr=50;
    for(int i=1;i<=nr;pas/=2,++i)
        for(int j=0;j!=4;++j)
        {
            t=dist(px+pas*dx[j],py+pas*dy[j]);
            if(t<dmin)
            {
                dmin=t;
                px+=pas*dx[j];
                py+=pas*dy[j];
                pas*=2;
                break;
            }
        }
    g<<px<<" "<<py<<"\n";
    return 0;
}