zju 3421 Error Curves(凸函数最值,三分法)

来源:互联网 发布:淘宝宝贝详情怎么设置 编辑:程序博客网 时间:2024/04/30 17:15

题意:
在[0, 1000]区间内,给n个二次曲线。
定义:f(x)=max(si(x)), si为第i条二次曲线。。
求f(x)最小值。
思路:
经过观察(其实是一个结论),我们发现f(x)是凸函数!
所以用三分法求解。
tip 三分法:
mid = (l+r)/2, mmid = (mid+r)/2
1)如果mmid更靠近最值点,有两种情况,mmid和mid在最值点同侧和异侧,这两种情况下,都应该取l=mid。
2)mid更靠近最值点,则r=mmid。
PS:精度1e-8wa,1e-10ac。。。

double get(double x, int id) {    double a = p[id][0], b = p[id][1], c = p[id][2];    return a*x*x+b*x+c;}double get(double x) {    double ret = get(x, 1);    rep(i, 2, n) ret = max (ret, get(x, i));    return ret;}double solve() {    double l = 0, r = 1000;    while (r-l>1e-10) {        double mid = (l+r)/2;        double mmid = (mid+r)/2;        if (get(mmid) < get(mid)) l = mid;        else r = mmid;    }    return get(l);}
0 0
原创粉丝点击