码迷,mamicode.com
首页 > 其他好文 > 详细

zoj 3386 Trick or Treat 三分 求最大值的 最小值

时间:2014-05-15 21:39:24      阅读:314      评论:0      收藏:0      [点我收藏+]

标签:style   blog   class   code   java   c   

题目来源:

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=3963

题意:  给定 N 个不同的点, 求在x轴上的 一点,  使 这点到N个点的 距离 最大 的 最小值。

f(x) =  max(i){ (xi - x) ^2 + yi ^2 }

求 x  使  min(f(x)) , f(x)为凹函数   ,  采用三分的形式

代码如下:

bubuko.com,布布扣
const double EPS = 1e-10 ;
const int Max_N =  50005 ;
int n;
double add(double a, double b){
    return (fabs(a + b) < EPS * (fabs(a) + fabs(b))) ? 0 : (a + b) ;
}
struct Point {
    double x, y;
    double dist(double a){
        return sqrt(add((x - a)*(x -a) ,(y)*(y) )) ;
    }
};
Point pt[Max_N] ;
double f(double x){
    int i ;
    double Max = 0 ;
    for( i = 0  ; i <  n;  i++){
        Max = pt[i].dist(x) > Max ? pt[i].dist(x) : Max ;
    }
    return Max ;
}
double tri_search(){
    double Mid, Midmid , L,  R ;
    L = -400000.0  , R = 400000.0 ;
    while(L + EPS < R){
        Mid = (L + R) * 0.5 ;
        Midmid = (Mid + R) *0.5 ;
        if(f(Mid) <= f(Midmid ) )
            R = Midmid ;
        else L = Mid ;
    }
    return L ;
}
int main(){
    while(scanf("%d", &n) && n){
        for(int i =0 ; i < n ; i++)
            scanf("%lf%lf" , &pt[i].x ,&pt[i].y ) ;
        double xx = tri_search() ;
        double Max = f(xx) ;
        printf("%.9lf %.9lf\n" , xx + EPS , Max) ;
    }
}
bubuko.com,布布扣

 

zoj 3386 Trick or Treat 三分 求最大值的 最小值,布布扣,bubuko.com

zoj 3386 Trick or Treat 三分 求最大值的 最小值

标签:style   blog   class   code   java   c   

原文地址:http://www.cnblogs.com/zn505119020/p/3729693.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!