hdu 5105——Math Problem

题意：计算a*x^3+b*x*x+cx+d的绝对值的最大值

思路：当a为0的时候，这是一个二次方程，取了绝对值以后最大值要么在L或者R的地方，要么在二次函数的对称轴上（当a和b都是0的时候在L或者R上）。当a不是0的时候，这是一个三次函数。首先求导以后得到 3ax^2+2bx+c，当判别式小于0的时候，导数与x轴没有交点，原函数是一个单调函数，最大值就在L或者R上，当判别式大于0的时候，求原函数的导数与坐标轴的两个焦点，然后在L，R，这两个焦点的位置求最大值（判断焦点是否在范围内）。

错误：比赛的时候求导求错了。。。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
double a,b,c,d,L,R;

double f(double x)
{
        return fabs(a*x*x*x+b*x*x+c*x+d);
}

const double eps=1e-9;
int dcmp(double x)
{
        if(fabs(x)<=eps)return 0;
        else return x<0?-1:1;
}

bool judge(double x)
{
        return x<=R&&x>=L;
}

int main()
{
//freopen("data.txt","r",stdin);

        while(scanf("%lf%lf%lf%lf%lf%lf",&a,&b,&c,&d,&L,&R)!=EOF)
        {
                double delta=4*b*b-12*a*c;
                double ans;
                if(dcmp(a)==0)
                {
                        double mid=-c/2/b;
                        if(judge(mid))
                                ans=max(f(mid),f(L));
                        else ans=f(L);
                        ans=max(ans,f(R));
                }
                else if(delta<=0)
                {
                         ans=max(f(L),f(R));
                }
                else
                {
                        delta=sqrt(delta);
                        double x1=-2*b+delta;
                        x1=x1/6/a;
                        double x2=-2*b-delta;
                        x2=x2/6/a;
                        ans=f(L);
                        ans=max(ans,f(R));
                        if(judge(x1))
                                ans=max(ans,f(x1));
                        if(judge(x2))
                                ans=max(ans,f(x2));

                }
                printf("%.2lf\n",ans);
        }
        return 0;
}