Pagini recente » Cod sursa (job #1767997) | Cod sursa (job #2568217) | Cod sursa (job #388039) | Cod sursa (job #1505255) | Cod sursa (job #2082369)
#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;
}