P1024 [NOIP2001 提高组] 一元三次方程求解（数学，二分，枚举，分治）

题目描述

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

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

输入格式

一行，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<bits/stdc++.h>
using namespace std;
int a, b, c, d;
double f(double x){
    return a*x*x*x + b *x *x + c * x + d;//一元三次方程
}
int main(){
    cin >> a >> b >> c >> d;
    for(int i = -100; i <= 100;i++){
        double x1= i;
        double x2 = x1+1;
        if(f(x1) == 0) printf("%.2lf ",x1);
        else if(f(x1) * f(x2) < 0){
            double l = x1, r = x2,mid;
            while(r-l>0.001){
                mid=(l+r)/2;
                if(f(l)*f(mid)<=0)
                    r = mid;
                else 
                    l = mid;
                }
            printf("%0.2lf ",l);
        }
    }
    return 0;
}