Pagini recente » Cod sursa (job #517418) | Cod sursa (job #410740) | Cod sursa (job #2632146) | Cod sursa (job #2401386) | Cod sursa (job #2032095)
# 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;
}