[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]);
}