Pagini recente » Cod sursa (job #2928737) | Cod sursa (job #2187015) | Cod sursa (job #1379803) | Cod sursa (job #670205) | Cod sursa (job #1054125)
#include<fstream>
#include<iomanip>
#include<cmath>
using namespace std;
typedef struct { double x,y; } punct;
int n,i;
punct a[50005];
double px,py;
double dist(double x, double y) {
double rez=0;
for (int i=1; i<=n; ++i) rez+=sqrt( (x-a[i].x)*(x-a[i].x) + (y-a[i].y)*(y-a[i].y) );
return(rez);
}
int main(void) {
ifstream fin("adapost2.in");
ofstream fout("adapost2.out");
fin>>n;
for (i=1; i<=n; ++i) {
fin>>a[i].x>>a[i].y;
px+=a[i].x; py=a[i].y;
}
px/=n; py/=n; double prec=1.0, s=dist(px,py);
while (prec>=0.0001) {
bool ok=1;
while (ok) {
ok=0;
double dn=dist(px,py+prec), ds=dist(px,py-prec),dv=dist(px-prec,py),de=dist(px+prec,py);
if (dn<s) { s=dn; py+=prec; ok=1;}
else if (ds<s) { s=ds; py-=prec; ok=1; }
else if (de<s) { s=de; px+=prec; ok=1; }
else if (dv<s) { s=dv; px-=prec; ok=1; }
}
prec/=10;
}
fout<<setprecision(4)<<fixed<<px<<" "<<py;
return(0);
}