每日一题 4.6 P1024 一元三次方程求解【二分】

题目描述

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

提示：记方程f(x)=0f(x)=0，若存在22个数x_1x1和x_2x2，且x_1<x_2x1<x2，f(x_1) \times f(x_2)<0f(x1)×f(x2)<0，则在(x_1,x_2)(x1,x2)之间一定有一个根。

输入格式

一行，44个实数A,B,C,DA,B,C,D。

输出格式

一行，33个实根，并精确到小数点后22位。

输入输出样例

输入 #1

1 -5 -4 20

输出 #1

-2.00 2.00 5.00

题解

限定条件比较多，大大降低了题目难度，加上数据水，混过了。。。

自己运了几次都是奇奇怪怪的东西。。。

后续再加强代码。。。

#include <bits/stdc++.h>
using namespace std;
double a,b,c,d;
double f(double x)//将函数提出来
{
    return a * x * x * x + b * x * x + c * x + d;
}

double dfs(double l,double r)//二分
{
    if(r - l <= 0.001)
        return l;
    double mid = (l + r)/2;
    if(f(l) * f(mid) <= 0)
        return dfs(l,mid);
    else
        return dfs(mid,r);
}
int main()
{
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
    double ans[105];
    int num = 0;
    for(int i = -100;i < 100;++i){//遍历
        double l = i,r = i + 1;
        if(f(l) == 0)
            ans[++num] = l;
        else if(f(l)*f(r) < 0)
            ans[++num] = dfs(l,r);
        if(num >= 3)
            break;
    }
    printf("%0.2lf %0.2lf %0.2lf\n",ans[1],ans[2],ans[3]);
    return 0;
}