#include <stdio.h>
#include <math.h>
typedef float R;
R sq(R x) {return x*x;}
R p1(R x, R y, R z){R r2 = x*x+y*y+z*z; return log((sqrt(r2-2*x+1)-x+1)/(sqrt(r2)-x));}
R p2(R x, R y, R z){return p1((x+1)/2, y/2, z/2);}
R p3(R L, R x, R y, R z){R r = 1./L; return p2(x*r, y*r, z*r);}
#define lr 100
#define lm (2*lr+1)
const float r = 1., dx=1/(R)lr;
int main(){R L[lm][lm]; R pma=0, pmi=1.e10; int ima, jma, imi, jmi;
for(int i=0; i<lm; ++i) for(int j=0; j<lm; ++j)
   L[i][j] =({R z = 1.- sq((i-lr)*dx) - sq((j-lr)*dx); z<0?0:sqrt(z);});
for(int i=0; i<lm; ++i) for(int j=0; j<lm; ++j) if(L[i][j]>0) {R p = 0;
  for(int k=0; k<lm; ++k) for(int l=0; l<lm; ++l) if(L[k][l]>0 && (i!=k || j!=l)) 
    p += p3(L[k][l], L[i][j], (i-k)*dx, (j-l)*dx);
  p *= dx*dx;
  if(p>pma) {pma=p; ima=i; jma=j;}
  if(p<pmi) {pmi=p; imi=i; jmi=j;}}
printf("max %d %d %17.14f\n", ima, jma, pma);
printf("min %d %d %17.14f\n", imi, jmi, pmi);
return 0;}
/*
max 20 40     4.1884
min 105 1     4.1878

real	0m35.038s
user	0m34.889s
sys	0m0.044s
Norms-MBP:Viscous norm$ 
*/