洛谷-1024 一元三次方程求解

题目描述
有形如:ax3+bx2+cx1+dx0=0这样的一个一元三次方程。给出该方程中各项的系数(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个实根,并精确到小数点后222位。

输入输出样例
输入样例#1:
1 -5 -4 20

输出样例#1:
-2.00 2.00 5.00

解释:线性搜索每次从-100开始每次步长0.5,一直到100,如果出现符号不同就开始二分得到答案就好了,注意精度

#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
double A,B,C,D;
double f(double x){
    return A*x*x*x+B*x*x+C*x+D;
}
double ok(double l,double r){
    double mid=0;
    while(fabs(r-l)>=1e-7){
        mid=(l+r)/2;
        if(f(mid)*f(l)<=0){
            r=mid;
        }else{
            l=mid;
        }
    }
    return mid;
}
int main(){
    cin>>A>>B>>C>>D;
    for(double i=-100.0;i<=100.0;i+=0.5){
        if(f(i)*f(i+0.5)<=0&&fabs(f(i+0.5)!=0)) printf("%.2lf ",ok(i,i+0.5));
    }
    return 0;
}

猜你喜欢

转载自blog.csdn.net/mkopvec/article/details/90029654
今日推荐