Pagini recente » Cod sursa (job #297244) | Cod sursa (job #2047312) | Cod sursa (job #2694898) | Cod sursa (job #2582875) | Cod sursa (job #762296)
Cod sursa(job #762296)
#include<fstream>
#include<iomanip>
#include<cmath>
using namespace std;
const double eps=0.0001;
int N;
double x[50005],y[50005],daux;
double verific(double midx,double midy){
daux=0;
for(int i=1;i<=N;++i)
daux+=sqrt((x[i]-midx)*(x[i]-midx)+(y[i]-midy)*(y[i]-midy));
return daux;
}
int main(void){
ifstream fin("adapost2.in");
ofstream fout("adapost2.out");
int i;
double midx,midy,pasx,W,xx=0,yy=0;
fin>>N;
for(i=1;i<=N;++i){
fin>>x[i]>>y[i];
midx+=x[i]; midy+=y[i];
}
midx/=N; midy/=N; pasx=512;
while(pasx>=eps){
W=verific(midx,midy);
if(verific(midx+pasx,midy)<W){ midx+=pasx; continue; }
if(verific(midx,midy+pasx)<W){ midy+=pasx; continue; }
if(verific(midx-pasx,midy)<W){ midx-=pasx; continue; }
if(verific(midx,midy-pasx)<W){ midy-=pasx; continue; }
pasx/=2;
}
fout<<setprecision(4)<<fixed;
fout<<midx<<' '<<midy<<'\n';
return 0;
}