一元三次方程求解

题目描述
有形如：a x^3 + b x^2 + c x + d = 0ax
3
+bx
2
+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数（a,b,c,da,b,c,d 均为实数），并约定该方程存在三个不同实根（根的范围在 -100−100 至 100100 之间），且根与根之差的绝对值 \ge 1≥1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后 22 位。

提示：记方程 f(x) = 0f(x)=0，若存在 22 个数 x_1x
1
​
和 x_2x
2
​
，且 x_1 < x_2x
1
​
<x
2
​
，f(x_1) \times f(x_2) < 0f(x
1
​
)×f(x
2
​
)<0，则在 (x_1, x_2)(x
1
​
,x
2
​
) 之间一定有一个根。

输入格式
一行，44 个实数 a, b, c, da,b,c,d。

输出格式
一行，33 个实根，从小到大输出，并精确到小数点后 22 位。

输入输出样例
输入 #1复制
1 -5 -4 20
输出 #1复制
-2.00 2.00 5.00
说明/提示
【题目来源】

NOIP 2001 提高组第一题

#include <algorithm>
#include <iostream>
#include <queue>
#include<map>
#include <set>
#include <stack>
#include <string>
#include <vector>
#include <cstring>
#include <cstdio>
#include<cmath>
#define ms(a,b) memset(a, b, sizeof a)
#define _for(i,a,b) for(int i=a;i<b;i++)
#define sc1(a) scanf("%d", &a)
#define sc2(a,b) scanf("%d%d", &a,&b)
#define sc3(a,b,c) scanf("%d%d%d", &a,&b,&c)
#define PII pair<int,int>
using namespace std;
double a, b, c, d;
int k;
double l[5], r[5];
set<double> st;
double con(double x){
    return a * pow(x, 3) + b * pow(x, 2) + c * x + d;
}


int main(){
    //freopen("in.txt", "r", stdin);
    cin >>a >> b >> c >> d;
    //查区间
    double pre = -101;
    for (double i = -100; i <= 101;i+=0.5){
        if(con(pre)*con(i)<=0){
            l[k] = pre;
            r[k++] = i;
        }
        pre = i;
    }

    //查
    _for(i,0,k){
        double ll = l[i]-0.1, rr = r[i]+0.1;
        while(rr-ll>=1e-5){
           
            double mid = (ll + rr) / 2;
            if((con(r[i]+0.1)>=0 ? (con(mid)>= 0)  : (con(mid)<=0) )){
                rr = mid;
            }
            else
                ll = mid;
        
        }
        
        //近似，判重
        ll = round(ll * 100) / 100;
        st.insert(ll);
    }

    for(auto i:st){
        printf("%.2lf ", i);
    }
}