[洛谷P1024] [NOIP2001 提高组] 一元三次方程求解 题解

原题链接

题目描述

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

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

输入格式

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

输出格式

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

样例输入

1 -5 -4 20

样例输出

-2.00 2.00 5.00

思路：函数求导 + 二分

• 首先将一元三次方程进行求导 $f(x)$，得到一个一元二次方程 $f^{\prime }(x)$，题目给出函数 $f(x)$必定有3个不同的实根，那么 $f^{\prime }(x)$求导后得到的一元二次方程求得 $\Delta$ > 0，设 $f^{\prime }(x)$的两个不同解分别为$x_1,x_2$，那么这两个解分别是 $f(x)$图像的最高点和最低点，则 $f(x)$的解分别处于 $-100\sim x_1,x_1\sim x_2,x_2\sim 100$，要求出解，直接对三个区间分别二分即可
  

static double a, b, c, d;
public static void solve() throws IOException{
    a = readDouble(); b = readDouble(); c = readDouble(); d = readDouble();
    double p = 3.0 * a, q = 2.0 * b;//求导后方程 f'的系数
    double x1 = (-q - Math.sqrt(q * q - 4 * p * c)) / (2 * p),
           x2 = (-q + Math.sqrt(q * q - 4 * p * c)) / (2 * p);//方程 f'的两个根
    List<Double> list = new ArrayList<>();
    // 区间 -100 ~ x1
    double l = -100, r = x1;
    while (r - l > 1e-4) {
        double mid = (l + r) / 2;
        if (fnc(mid) <= 0) {
            l = mid;
        } else {
            r = mid;
        }
    }
    list.add(l);
    // 区间 x1 ~ x2
    l = x1; r = x2;
    while (r - l > 1e-4) {
        double mid = (l + r) / 2;
        if (fnc(mid) <= 0) {
            r = mid;
        } else {
            l = mid;
        }
    }
    list.add(r);
    // 区间 x2 ~ 100
    l = x2; r = 101;
    while (r - l > 1e-4) {
        double mid = (l + r) / 2;
        if (fnc(mid) <= 0) {
            l = mid;
        } else {
            r = mid;
        }
    }
    list.add(l);
    for (double res : list) {
        printWriter.printf("%.2f ", res);
    }
}
//计算函数值
public static double fnc(double x) {
    return 1.0 * a * x * x * x + 1.0 * b * x * x + c * x + d;
}