【三分】Error Curves
来源:互联网 发布:英文表达祝福知乎 编辑:程序博客网 时间:2024/05/17 23:19
In order to test the algorithm's efficiency, she collects many datasets. What's more, each data is divided into two parts: training data and test data. She gets the parameters of the model on training data and test the model on test data.
To her surprise, she finds each dataset's test error curve is just a parabolic curve. A parabolic curve corresponds to a quadratic function. In mathematics, a quadratic function is a polynomial function of the formf(x) = ax2 + bx + c. The quadratic will degrade to linear function ifa = 0.
It's very easy to calculate the minimal error if there is only one test error curve. However, there are several datasets, which means Josephina will obtain many parabolic curves. Josephina wants to get the tuned parameters that make the best performance on all datasets. So she should take all error curves into account, i.e., she has to deal with many quadric functions and make a new error definition to represent the total error. Now, she focuses on the following new function's minimal which related to multiple quadric functions.
The new function F(x) is defined as follow:
F(x) = max(Si(x)), i = 1...n. The domain ofx is [0, 1000]. Si(x) is a quadric function.
Josephina wonders the minimum of F(x). Unfortunately, it's too hard for her to solve this problem. As a super programmer, can you help her?
Input
The input contains multiple test cases. The first line is the number of cases T (T < 100). Each case begins with a number n(n ≤ 10000). Followingn lines, each line contains three integers a (0 ≤ a ≤ 100),b (|b| ≤ 5000), c (|c| ≤ 5000), which mean the corresponding coefficients of a quadratic function.
Output
For each test case, output the answer in a line. Round to 4 digits after the decimal point.
Sample Input
212 0 022 0 02 -4 2
Sample Output
0.00000.5000
容易证明,F(x)是下凸函数(也即是一种单峰函数),可以用三分解决。
小心最后是输出函数值而不是x值,还有二分次数可以多一点,二分区间为10^6,才二分二十多次,可以二分100次。
#include <cstdio>double f(int a,int b,int c,double x){return x*x*a+x*b+c;}int a[10010];int b[10010];int c[10010];int n;double max(double a,double b){return a>b?a:b;}double F(double x){double ff = -1e21;for (int i=1;i<=n;i++){ff = max(ff,f(a[i],b[i],c[i],x));}return ff;}int main(){freopen("5009.in","r",stdin);freopen("5009.out","w",stdout);int T;scanf("%d",&T);while (T--){scanf("%d",&n);for (int i=1;i<=n;i++){scanf("%d%d%d",a+i,b+i,c+i);}double l = 0;double r = 1000;for (int i=0;i<100;i++){double lm = l+(r-l)/3;double rm = r-(r-l)/3;double f1 = F(lm);double f2 = F(rm);if (f1 < f2)r = rm;elsel = lm;}printf("%.4lf\n",F(l));}return 0;}
- 【三分】Error Curves
- UVA5009 Error Curves三分
- hdu3714 Error Curves---三分复习
- [三分]HDOJ 3714Error Curves
- la 5009 - Error Curves(三分)
- uva 1476 - Error Curves(三分)
- zoj 3421 Error Curves(三分)
- LA5009 Error Curves(三分)
- HDU3714:Error Curves(三分)
- hdu 3714 Error Curves (三分)
- 三分搜索--hdu3714 Error Curves
- uva-1476 Error Curves 三分
- HDOJ 3714 Error Curves 三分搜索
- hdu 3714 Error Curves(三分)
- Uva - 1476 - Error Curves(三分)
- 数学-三分->HDU 3714 Error Curves
- HDU 3714 Error Curves(三分)
- nyoj 1029/hdu 3714 Error Curves 三分
- 关于在ie下无法读写cookie的解决
- 基于jquery ajax 无刷新 文件批量上传插件 GooUploader整合struts2实现上传
- 北漂
- Oracle中的&号问题
- 聊聊并发(二)——Java SE1.6中的Synchronized
- 【三分】Error Curves
- 编程面试的10大算法概念汇总
- 数字图像处理(四) 彩色图像拼接
- 下载索引
- 测试朋友的东西!
- Fix协议封装挑战-消息处理
- 功率谱
- 聊聊并发(四)——深入分析ConcurrentHashMap
- java的Date和Calendar的一些方法