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

最新推荐文章于 2024-01-18 17:08:46 发布

乄卝⑦

最新推荐文章于 2024-01-18 17:08:46 发布

阅读量237

点赞数 1

文章标签：算法 c++

本文链接：https://blog.csdn.net/weixin_66162442/article/details/128397640

版权

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

题目描述

有形如： $a x^3 + b x^2 + c x + d = 0$ 这样的一个一元三次方程。给出该方程中各项的系数（a,b,c,d均为实数），并约定该方程存在三个不同实根（根的范围在 -100 至 100 之间），且根与根之差的绝对值 ≥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

题目解析

因为根与根之差的绝对值 ≥1，所以在长度为1的区间内只可能有一个解。因此，我们可以把区间[-100,100]分成长度为1的区间（[-100,-99),[-99,-98),……,[99,100)），依次枚举来找方程的解。
如何找解？首先题目给了提示就是我们熟知的零点定理（利用二分来完成）显然这里用到的是小数二分，每一个区间[i,i+1) 对应的最初的l=i,r=i+1。如果f(l) * f(mid) < 0，即方程的解在l和mid中间。反之在r和mid间。切记，l和r不能+1。到这里看似已经完成，可我们忽略了区间的端点，零点定理的范围是在开区间 $(x_1, x_2)$ 上，我们还需再对每个区间的端点进行判断即可。由于100这一点最后没有枚举，在最后进行判断即可。
为了方便可写一个函数计算函数的值。此外，我们枚举时i应该枚举到99，因为我们枚举的是区间[99,100)。

AC代码

#include<iostream>
using namespace std;

double 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;
    //从[-100,-99),[-99,-98)···[99,100)依次枚举
    for (double i = -100; i <= 99; i++)
    {
        //判断左端点是否为方程的解
        if (f(i) == 0) printf("%.2lf ",i);
        //如果在此区间有界，利用二分来找答案
        else if (f(i) * f(i + 1) < 0)
        {
            double l = i, r = i + 1;
            while (r-l>=1e-4)
            {
                double mid = (l + r) / 2;
                if (f(l) * f(mid) < 0) r = mid;
                else l = mid;
            }
            printf("%.2lf ", l);
        }
    }
    //最后判断100是否为解
    if (f(100) == 0) cout << "100.00 ";

    return 0;
}