hdu3264

给 n 个圆，作一个新的圆（以 n 个圆里面的一个圆心为圆心），要求将所有的圆覆盖至少一半的面积，求出最小半径。

const  int  maxn = 28  ;

const  double  eps = 1e-8 ;

const  double  pi =  acos(-1.0) ;

struct point{
       double x , y ;
       point(){}
       point(double x , double y){
            this->x = x ;
            this->y = y ;
       }
};

double  dist(const point &a , const point &b){
        return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y)) ;
}

double circle_area_inst(point o1 , double r1 ,  point o2 , double r2){
       double d = dist(o1 , o2) ;
       if(d > r1 + r2 - eps)  return 0 ;
       if(d < r2 - r1 + eps)  return pi * r1 * r1 ;
       if(d < r1 - r2 + eps)  return pi * r2 * r2 ;
       double a1 = acos( (r1*r1 + d*d - r2*r2) / (2.0*r1*d) ) ;
       double a2 = acos( (r2*r2 + d*d - r1*r1) / (2.0*r2*d) ) ;
       return r1*r1*(a1 - 0.5*sin(2.0*a1) ) + r2*r2*(a2 - 0.5*sin(2.0*a2) ) ;
}

point o[maxn] ;
double  r[maxn] ;
int   n ;

int     judge(double R){
        int i ,  j ;
        for(i = 1 ; i <= n ; i++){
             for(j = 1 ; j <= n ; j++){
                  if(circle_area_inst(o[i] , R , o[j] , r[j]) < pi*r[j]*r[j]*0.5) break ;
             }
             if(j > n) return 1 ;
        }
        return 0 ;
}

double  getr(){
        double l = 0 , r = 1e6  , m  , s ;
        while(l + eps < r){
              m = (l + r) * 0.5 ;
              if(judge(m)){
                    s = m ; r = m - eps ; }
              else   l = m + eps ;
        }
        return s ;
}

int  main(){
     int t , i ;
     cin>>t ;
     while(t--){
          cin>>n ;
          for(i = 1 ; i <= n ; i++) scanf("%lf%lf%lf" , &o[i].x , &o[i].y , &r[i]) ;
          printf("%.4lf\n" , getr()) ;
     }
     return 0 ;
}