Pagini recente » Cod sursa (job #250088) | Cod sursa (job #1416850) | Cod sursa (job #2297767) | Cod sursa (job #982184) | Cod sursa (job #1846010)
#include <iostream>
#include <cmath>
#include <map>
using namespace std;
const int Q=50007;
int n;
struct point
{
double x,y;
}v[Q];
std::map<double, map<double,double> > prec;
double dist(double a, double b)
{
//if(prec[a][b]!=0)
// return prec[a][b];
double rez=0;
for(int i=1; i<=n; i++)
{
rez+=sqrt((a-v[i].x)*(a-v[i].x) + (b-v[i].y)*(b-v[i].y));
}
//prec[a][b]=rez;
return rez;
}
double aux,dist_act,gx,gy,pas;
void bra1()
{
aux=dist(gx,gy+pas);
if(aux<dist_act)
{
gy+=pas;
dist_act=aux;
}
else
{
aux=dist(gx,gy-pas);
if(aux<dist_act)
{
gy-=pas;
dist_act=aux;
}
}
}
void bra2()
{
aux=dist(gx+pas,gy);
if(aux<dist_act)
{
gx+=pas;
dist_act=aux;
}
else
{
aux=dist(gx-pas,gy);
if(aux<dist_act)
{
gx-=pas;
dist_act=aux;
}
}
}
int main()
{
freopen("adapost2.in","r",stdin);
freopen("adapost2.out","w",stdout);
cin>>n;
gx=0,gy=0;
for(int i=1; i<=n; i++)
{
cin>>v[i].x>>v[i].y;
gx+=v[i].x;
gy+=v[i].y;
}
gx/=n;
gy/=n;
pas=1<<9;
aux,dist_act=dist(gx,gy);
for(int tmp=1; tmp<=40; tmp++)
{
if(rand()&1)
{
bra1();
bra2();
}
else
{
bra2();
bra1();
}
if(tmp%2==0)
pas/=2;
}
printf("%.6lf %.6lf\n",gx,gy);
return 0;
}