hdu 3714

来源:互联网 发布:内蒙古网络春晚 编辑:程序博客网 时间:2024/06/16 13:14

参考厚白书

题意:求出F(x)的最小值。x的取值范围为[0, 1000]。F(x) = max(Si(x))。Si = ai * x *x + bi * x + ci。其中ai恒定>= 0。

也就是求0 -> n - 1 的每个Si的最小值,再在每个最小值之中取最大值。

因为ai恒大于零,所以Si必定能用三分法求出。因为F(x)=max(S0, S1, S2, S3.....) i = 0 -> n  - 1, x 属于[0, 1000]。因为每个S都是下凹函数,组成的函数也必定是下凹的函数。所以F(x)也就是下凹的函数。(宝宝心里苦,希望有大神看到可以告诉我为什么)所以,三分法。


#include<cstdio>#include<algorithm>#include<cmath>using namespace std;const int maxn = 1e5 + 5;const double eps= 1e-7;int a[maxn], b[maxn], c[maxn];int n;double F(double x){    double ans = a[0] * x * x + b[0] * x + c[0];    for(int i = 1; i< n; i++)    {        ans = max(ans, a[i] * x * x + b[i] * x + c[i]);    }    return ans;}int main(){    int T;    scanf("%d", &T);    for(int kase = 1; kase <= T; kase++)    {        scanf("%d", &n);        for(int i = 0; i < n; i++)            scanf("%d%d%d", &a[i], &b[i], &c[i]);        double L = 0, R = 1000;        while(fabs(F(R) - F(L)) >= eps)        {            double mid = R - L;            //printf("%lf   %lf\n", R, L);            if(F(L + mid / 3) > F(R - mid / 3)) L = L + mid / 3;            else R = R - mid / 3;        }        //printf("R = %lf\n", R);        printf("%.4lf\n", F(R));    }    return 0;}


原创粉丝点击