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

最新推荐文章于 2024-07-21 13:38:21 发布

Ryan1229

最新推荐文章于 2024-07-21 13:38:21 发布

阅读量477

点赞数 16

分类专栏：洛谷【普及-】文章标签：算法 c++ c语言开发语言

本文链接：https://blog.csdn.net/Ryan1229/article/details/137418016

版权

洛谷【普及-】专栏收录该内容

68 篇文章 0 订阅

订阅专栏

本文介绍了一种求解一元三次方程ax^3+bx^2+cx+d=0的方法，利用二分法确定三个不同实根，保证根与根之差大于1，结果精确到小数点后两位。

摘要由CSDN通过智能技术生成

时间限制1.00s 内存限制125.00MB 难易度：普及−

【题目描述】

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

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

【输入格式】

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

【输出格式】

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

【输入输出样例】

输入 #1

1 -5 -4 20

输出 #1

-2.00 2.00 5.00

【说明/提示】

【题目来源】

NOIP 2001 提高组第一题

【算法分析】

这道题用二分法，若f(x1)=0，则确定x1为f(x)的根；若f(x1)*f(x2)>0，则确定根x不在[x1,x2]内，设[x2,x2+1]为下一个搜索区间；若f(x1)*f(x2)<0，则确定根x在[x1,x2]内。

【参考代码】

#include <bits/stdc++.h>
using namespace std;
double a,b,c,d,x,x1,x2,xx;
double f(double x)
{
    return (x*x*x*a+x*x*b+x*c+d);
}
int main()
{
    cin>>a>>b>>c>>d;
    for(x=-100;x<=100;x++)
    {
        x1=x;x2=x+1;
        if(f(x1)==0) printf("%.2lf ",x1);
        else if(f(x1)*f(x2)<0)
        {
            while(x2-x1>=0.001)
            {
                xx=(x1+x2)/2;
                if((f(x1)*f(xx))<=0) x2=xx;
                else x1=xx;
            }
            printf("%.2lf ",x1);
        }
    }
    return 0;
}