HDU Error Curves

//这题有坑, 用C++交题才能A, G++TLE

//把图画在纸上, 直觉会告诉你, max(si(x))构成一个有凸多边形, 极值在最低处~,
//然后可以用三分, 精度至少在1e-9

#include <iostream>
#include <cstdio>
#include <queue>
#include <cstring>
#include <cmath>
#include <vector>
#include <set>
#include <stack>
#include <algorithm>
//#include "myAlgorithm.h"

#define MAX 10005
#define INF (1e8 + 5)
#define eps 1e-9
#define Rep(s, e) for( int i = s; i <= e; i++)
#define Cep(e, s) for( int i = e; i >= s; i --)
using namespace std;
int T, n;
double a[MAX], b[MAX], c[MAX];
double func(int i,  double val){
    return a[i] * val * val + b[i] * val + c[i];
}
double get_min(){
    double d = 0, u = 1000, mid, mmid, res1, res2;
    while(u - d > eps){
        mid = (u + d)/2; mmid = (mid + u)/2;
        res1 = func(0, mid);
        res2 = func(0, mmid);
        Rep(1, n - 1){
                res1 = max(res1, func(i, mid));
                res2 = max(res2, func(i, mmid));
        }
        if(res1 < res2) {
            u = mmid;
        } else {
            d = mid;
        }
    }
    return res1;
}

int main() {
    cin>>T;
    //freopen("in.tx0t", "w", stdout);
    while(T--){
        cin>>n;
        Rep(0, n - 1){
            cin>>a[i]>>b[i]>>c[i];
        }
        printf("%.4f\n", get_min());
    }
    return 0;
}
/*
*/




 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值