[NOIP2001 提高组] 一元三次方程求解_一元三次方程求解 noip-CSDN博客

本文链接：https://blog.csdn.net/naexting/article/details/120060621

该博客介绍了一种解决一元三次方程的方法，通过二分法寻找方程在特定区间内的实根。代码示例展示了如何在C++中实现这个算法，确保找到的根满足条件且精确到小数点后两位。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目描述

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

输入输出格式

输入格式

一行， $4$ 个实数 $a, b, c, d$ 。

输出格式

一行， $3$ 个实根，从小到大输出，并精确到小数点后 $2$ 位。

输入输出样例

输入样例

1 -5 -4 20

输出样例

-2.00 2.00 5.00

#include<iostream>
#include<cstdio>
using namespace std;
double a,b,c,d;
double fun(double x)
{
    double y=((a*x+b)*x+c)*x+d; 
    return y;
}
int main()
{
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);  //输入
    for(double i=-100;i<=100;i++)
    {
        double f1=fun(i);
        double f2=fun(i+1);
        // [i,i+1]之间存在根
        if(f1==0)
            printf("%.2lf ",i); 
        if(f1*f2<0)
        {
            double l=i,r=i+1;
            double mid;
                //在整数区间内求近似值，精度在0.001即可
            while(r-l>1e-3)
            {
                mid=l+(r-l)/2.0;
 
                //二分法的技巧
                if(fun(l)*fun(mid)<=0)
                    r=mid;
                else
                    l=mid;
            }
            printf("%.2lf ",mid);
        }
    }
    return 0;
}
/*
1 -5 -4 20

-2.00 2.00 5.00
*/