2021-01-27

洛谷 一元三次方程求解

数论,数学
二分答案
枚举,暴力
分治

原题

题目描述

有形如：a x^3 + b x^2 + c x + d = 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 个实根，从小到大输出，并精确到小数点后 2 位。
输入输出样例
输入 #1

1 -5 -4 20

输出 #1

-2.00 2.00 5.00

我的题解

有一个麻烦的地方就是始终要保持两位小数输出。

#include <iostream>
#include <iomanip>
using namespace std;
double a,b,c,d;
double zhi(double y)
{
    double y1;
    y1=a*y*y*y+b*y*y+c*y+d;
    return y1;
}
int main()
{
    double n=0,x,y,y1,y2,y3,z;
    cin>>a>>b>>c>>d;
    for(y=-100;y<100;y++)
    {
        z=y+1;
            y1=zhi(y);
            y2=zhi(z);
            if(y1*y2==0)
            {
                if(y1==0)
                {
                    x=y;
                    cout<<setiosflags(ios::fixed)<<setprecision(2)<<x<<" ";//设置保留两位小数输出
                    n++;
                }
          /*      if(y2==0)//只取区间左边的零点防止重复
                {
                    x=z;
                    cout<<x<<" ";
                    n++;
                }*/
            }
            if(y1*y2<0)
            {        while(z-y>=0.001)//二分控制精度
        {
                y3=zhi((y+z)/2);
                if(y1*y3<0) z=(y+z)/2;
                if(y2*y3<0) y=(y+z)/2;
            /*    if(y3==0)
                {
                    x=(y+z)/2;
                    cout<<x<<" ";
                    n++;
                    break;
                }*/
            }
            x=y;
            cout<<setiosflags(ios::fixed)<<setprecision(2)<<x<<" ";//设置保留两位小数输出
            n++;}
        if(n>2) break;

    }
    return 0;
}