Error Curves（三分习题）（精度）

之前对double类型的数据的比较一直一无所知，A了这道题后感觉很神奇

突然很想写个博客记录一下，以便后面查看

 

题目：http://acm.bnu.edu.cn/v3/external/50/5009.pdf

下面代码中的宏定义“#define eps 1e-12”

刚开始，我是“#define eps 1e-8”，提交是：

然后网上找了找，看到有位大佬写的是“#define eps 1e-15”，我也改成了1e-15，提交结果是：

 

然后我看了一下我和他的代码，思路一样，测试了几条数据也是一样的结果

我就改成了“#define eps 1e-12”，然后就

 

 

#include<iostream>
#include<math.h>
#include<stdio.h>
using namespace std;
int Data[10005][3];
#define eps 1e-12
double Cal(double x,int a,int b, int c)
{
    return a*x*x+b*x+c;
}
double Max(double *x,int n)
{
    int maxid=0;
    for(int i=1;i<n;++i)
        if(x[i]>x[maxid])
            maxid=i;
    return x[maxid];
}
int main()
{
    int T;
    int n;
    cin>>T;
    while(T>0)
    {
        cin>>n;
        for(int i=0;i<n;++i)
            cin>>Data[i][0]>>Data[i][1]>>Data[i][2];
        double mid1,mid2;
        double res1[n]={0};
        double res2[n]={0};
        double f1,f2;
        double left=0,right=1000;
        while(left<right-eps){
            mid1=left+(right-left)/3;
            mid2=mid1+(right-left)/3;
            for(int i=0;i<n;++i)
            {
                res1[i]=Cal(mid1,Data[i][0],Data[i][1],Data[i][2]);
            }
            f1=Max(res1,n);
            for(int i=0;i<n;++i)
            {
                res2[i]=Cal(mid2,Data[i][0],Data[i][1],Data[i][2]);
            }
            f2=Max(res2,n);
            if(f1<f2) right=mid2;
            else left=mid1;
        }
        printf("%.4lf\n",f2);
        T--;
    }
    return 0;
}