因为这里l和R 给我们的范围不大,然后考虑暴力的方法。。真的可以
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <vector>
#include <set>
#include <map>
#include <queue>
using namespace std;
double a,b,c,d,l,r;
double f(double x)
{
return abs(a*x*x*x+b*x*x+c*x+d);
}
int main()
{
while(scanf("%lf %lf %lf %lf %lf %lf",&a,&b,&c,&d,&l,&r)!=EOF)
{
double ans=0;
while(r-l>=0.01)
{
ans=max(ans,f(l));
l+=0.01;
}
printf("%.2f\n",max(ans,f(r)));
}
return 0;
}
或者也可以把可能的点列出来。
如果a=0时,看对称轴和端点。如果不为0 ,则求导看delta,如果delta小于等于0,则就是端点,如果大于0,就是极点。
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <vector>
#include <set>
#include <map>
#include <queue>
using namespace std;
double a,b,c,d,l,r;
#define eps 1e-5
double f(double x)
{
return abs(a*x*x*x+b*x*x+c*x+d);
}
double x[10];
bool ok(double x1)
{
if(x1>0) return 1;
if(abs(x1)<eps) return 1;
else return 0;
}
int main()
{
while(scanf("%lf %lf %lf %lf %lf %lf",&a,&b,&c,&d,&l,&r)!=EOF)
{
double ans=0;
x[1]=l,x[2]=r;
int h=3;
if(b!=0) x[h++]=-c/(2*b);
double delta=b*b-3*a*c;
if(delta>0){
x[h++]=(-b+sqrt(delta))/(3*a);
x[h++]=(-b-sqrt(delta))/(3*a);
}
for(int i=1;i<h;i++)
{
if(ok(x[i]-l)>0&&ok(r-x[i])>0) ans=max(ans,f(x[i]));
}
printf("%.2f\n",max(ans,f(r)));
}
return 0;
}