Cod sursa(job #2082369)

Utilizator anca.sotirAnca Sotir anca.sotir Data 6 decembrie 2017 01:53:11
Problema Adapost 2 Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.43 kb
#include <fstream>
#include <math.h>
#define nmax 50002
#define err 0.001
using namespace std;
ifstream f("adapost2.in");
ofstream g("adapost2.out");

struct punct{
    float x,y;
} p[nmax];
float pas=100;
int n;

float distanta(punct a, punct b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
float suma_distante(punct x)
{
    float s=0;
    for(int i=1;i<=n;++i)
        s+=distanta(x,p[i]);
    return s;
}
punct CDG()
{
    punct cdg;
    float x,y;
    for(int i=1;i<=n;++i)
    {
        x+=p[i].x;
        y+=p[i].y;
    }
    cdg.x=x/n;
    cdg.y=y/n;
    return cdg;
}
punct gaseste_loc(punct x)
{
    if(pas<err) return x;

    punct vecin;

    vecin.x=x.x;
    vecin.y=x.y+pas;
    if(suma_distante(vecin)<suma_distante(x))
        return gaseste_loc(vecin);

    vecin.x=x.x;
    vecin.y=x.y-pas;
    if(suma_distante(vecin)<suma_distante(x))
        return gaseste_loc(vecin);

    vecin.x=x.x+pas;
    vecin.y=x.y;
    if(suma_distante(vecin)<suma_distante(x))
        return gaseste_loc(vecin);

    vecin.x=x.x-pas;
    vecin.y=x.y;
    if(suma_distante(vecin)<suma_distante(x))
        return gaseste_loc(vecin);

    pas/=2;
    return gaseste_loc(x);
}
int main()
{
    f>>n;
    for(int i=1;i<=n;++i)
    {
        f>>p[i].x;
        f>>p[i].y;
    }
    punct cdg=CDG();
    punct loc=gaseste_loc(cdg);
    g<<loc.x<<' '<<loc.y;
    return 0;
}