[LUOGU1024] 一元三次方程求解

题目描述

有形如：ax3+bx2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a，b，c，d 均为实数)，并约定该方程存在三个不同实根(根的范围在-100至100之间)，且根与根之差的绝对值>=1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后2位。

提示：记方程f(x)=0，若存在2个数x1和x2，且x1

#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
double a,b,c,d,ax[105];
int t; 
double f(double x)
{
    if(abs(a*x*x*x+b*x*x+c*x+d)<0.00000000001)
        return 0.0;
    else 
        return a*x*x*x+b*x*x+c*x+d;
}
double rst;
void search(double l,double r)
{
    double mid=(l+r)/2;
    double fmid=f(mid);
    double fl=f(l);
    double fr=f(r);
    if(fmid==0) {rst=mid;return;}
    if((fl>0 && fmid>0)||(fl<0 && fmid<0))
        search(mid,r);
    else
        search(l,mid);
}
int main()
{
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
    for(double i=-100;i<=100;i=i+1.0)
    {
        double j=i+1.0;
        double fi=f(i);
        double fj=f(j);
        if(fi==0.0) {ax[++t]=i;continue;}//且根与根之差的绝对值>=1
        if(fi*fj<0)
        {
            search(i,j);
            ax[++t]=rst;
        }
        if(t==3) break;
    }
    for(int i=1;i<=t;i++)   
        printf("%0.2lf ",ax[i]);
}