欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

洛谷题解——P1024:一元三次方程求解

程序员文章站 2022-05-09 10:27:53
...

视频讲解可以直接点击这个 B 站链接,https://www.bilibili.com/video/BV1qT4y13717/

题目相关

题目链接

洛谷,https://www.luogu.com.cn/problem/P1024

我的 OJ,http://47.110.135.197/problem.php?id=4242

题目描述

有形如:洛谷题解——P1024:一元三次方程求解 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,d 为实数),并约定该方程存在三个不同实根(根的范围在 -100 − 100 之间),且根与根之差的绝对值 ≥ 1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后 2 位。

提示:记方程 f(x) = 0,若存在 2 个数 x1 ​和 x2​,且 x1 ​< x2​,f(x1​)×f(x2​) < 0,则在 (x1​, x2​) 之间一定有一个根。

输入格式

一行,4 个实数 A,B,C,D。

输出格式

一行,3 个实根,并精确到小数点后 2 位。

输入样例

1 -5 -4 20

输出样例

-2.00 2.00 5.00

题目分析

题意分析

给了一个特殊的一元三次方程,该方程在 [-100, 100] 之间存在 3 个不同的实根,而且实根之间差的绝对值大于等于 1。要求我们求出这 3 个实根。

样例数据分析

本题是一个方程求解问题,样例数据没什么需要分析的,跳过。

数据规模分析

[-100, 100] 之间求解,而且要求的精度为小数点后 2 位,也就是 1e-2。

这句话意味着什么?

洛谷题解——P1024:一元三次方程求解,看懂没有?

算法思路

从上面的数据规模分析可以得出,最大的数据个数就是 20,000 个中筛选出 3 个数据。数据量小得可怜。因此本题的解法有很多种。

枚举算法

所谓暴力出奇迹。我们可以看到,要求精度是 0.01,数据范围为 [-100, 100],在这个范围一共有 20,000 个数据,完全可以考虑用暴力枚举。

也就是从 -100 开始,到 100 结束,每间隔 0.01,对数据进行枚举。

AC 参考代码

#include <bits/stdc++.h>

const double EPS = 1e-6;

int main() {
    double a,b,c,d;
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);

    //枚举算法
    for (double i=-100; i<=100; i+=0.01) {
        double ans = a*i*i*i+b*i*i+c*i+d;
        if (fabs(ans)<EPS) {
            printf("%.2lf ", (double)i);
        }
    }
    printf("\n");

    return 0;
}

二分查找

本题和标准二分查找有些不一样。我们知道标准二分查找都是在某个范围内,查找一个数据。而本题是在 -100 到 100 之间,要求查找出所有答案。

因此为了套用标准浮点数二分查找,我们需要将 [-100, 100] 这个区间进行细分,然后在这个分区间里进行二分查找。因此本题的核心就是如何找到有解的区间,然后在这个区间进行二分。

查找有解区间

其实题目已经告诉我们如何查找。关注以下几个部分:

1、根与根之差的绝对值 ≥ 1。

2、记方程 f(x) = 0,若存在 2 个数 x1 ​和 x2​,且 x1 ​< x2​,f(x1​)×f(x2​) < 0,则在 (x1​, x2​) 之间一定有一个根。

浮点二分的难点

就是 check() 函数如何写。

本题 check() 函数还是需要参考“记方程 f(x) = 0,若存在 2 个数 x1 ​和 x2​,且 x1 ​< x2​,f(x1​)×f(x2​) < 0,则在 (x1​, x2​) 之间一定有一个根”这句话。

AC 参考代码

直接套用浮点数二分查找模板代码即可。

#include <bits/stdc++.h>

double a,b,c,d;
const double EPS = 1e-4;

double calc(double mid) {
    return a*mid*mid*mid+b*mid*mid+c*mid+d;
}

bool check(double left, double mid) {
    double x1 = calc(left);
    double x2 = calc(mid);
    return x1*x2<0 ? true : false;
}

double bsearch(double left, double right) {
    double mid;
    while (fabs(right-left)>EPS) {
        mid = (left+right)/2;

        if (check(left, mid)) {
            right = mid;
        } else {
            left = mid;
        }
    }

    return left;
}

int main() {
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);

    for (int i=-100; i<=100; i++) {
        double x1 = calc(i);
        double x2 = calc(i+1-EPS);

        if (fabs(x1)<EPS) {
            printf("%.2lf ", (double)i);
        } else if (x1*x2<0) {
            printf("%.2lf ", bsearch(i, i+1));
        }
    }
    printf("\n");

    return 0;
}

代码细节

1、EPS 定义。根据经验公式,只需要比题目要求精度大两个数量级即可。因此 1e-4 对于本题而言,足够了。

2、check() 函数。我们知道该函数的作用是缩小查找的范围,本来我们是在 [left, right] 这个范围找答案。因此根据题目提供的有根条件,我们可以尝试在 [left, mid] 这个范围查找是否存在根,然后根据反馈进行范围调整。当然也可以在 [mid, right] 这个范围查找。

3、由于根与根之差的绝对值 ≥ 1,这样查找区间总长不能为 2。如果查找的区间总长为 2,依据上面的代码,输入的数据为 1 -4.65 2.25 1.4。当 i=0 的时候,x1!=0 && x1*x2<0,我们用二分查找区间 [0, 1]之间可以查找到一个根为 1.00。当 i=1 的时候,x1=0,又会输出一个根 1.00。导致出现重复答案。