一元三次方程求解

Problem Description

有形如:ax^3+bx^2+cx+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)之间一定有一个根。

Input

输入有多行测试数据,每行为四个系数a,b,c,d,输入以0 0 0 0结束。

Output

对于每组测试数据,输出一个三个实根(根与根之间留有空格),并精确到小数点后2位。

Sample Input

1 -5 -4 20
0 0 0 0

Sample Output

-2.00 2.00 5.00
/*思路:若存在2个数x1和x2,且x1<x2,f(x1)*f(x2)<0,则在(x1,x2)之间一定有一个根。
可以用二分法来写,查找。范围在-100到100,三个根。
先对一元三次方程求导,得到极值(驻点)利用高数方法就分成三部分
,分别在三部分里面求根。缩小范围*/ 
#include<cstdio>
#include<cmath>
#include<iostream>
using namespace std;
double a,b,c,d;
double fun(double x)
{
    return a*x*x*x+b*x*x+c*x+d;
}
double f(double x1,double x2)//二分法查找
{
     double mid;
     while(1)
    {
      mid=(x1+x2)/2;      
      if( fabs(x1-x2 )<1E-8) return mid;              
           else
                {
                     if(fun(mid) ==0)  return mid;
                       else  if(fun(x1)*fun(mid )<0) x2=mid;
                       else  x1=mid;
               }
    }
 }
int main()
{
   double x1,x2;
    while( scanf("%lf%lf%lf%lf",&a,&b,&c,&d)!=EOF)
   {
        if(!(a&&b&&c&&d)) break;
                 //一元三次变为一元二次,求导,降幂,得f'(x)=3*a*x*x+2*b*x+c;
               //解方程f'(x) 一元二次求根公式
               x1=(-2*b-sqrt(4*b*b-12*a*c))/(6*a);
         x2=(-2*b+sqrt(4*b*b-12*a*c))/(6*a);
                  //得到两个实根x1、x2,也就是f(x)的极值点(驻点)。           
             //f(x)=0的3个根应在区间[-100,x1]、[x1,x2]、[x2,100]中。  
         printf("%.2lf %.2lf %.2lf\n",f(-100,x1),f(x1,x2),f(x2,100));
    }
    return 0;
}

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值