1038 一元三次方程求解

2018-06-17 23:00:50来源:未知 阅读 ()

新老客户大回馈,云服务器低至5折

1038 一元三次方程求解

2001年NOIP全国联赛提高组

时间限制: 1 s
空间限制: 128000 KB
题目等级 : 白银 Silver
 
 
 
题目描述 Description

有形如:ax3+bx2+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 Description

一个三次方程的各项系数

输出描述 Output Description

三个解

样例输入 Sample Input

1   -5   -4   20

样例输出 Sample Output

-2.00   2.00   5.00

数据范围及提示 Data Size & Hint
 1 #include<iostream>
 2 #include<cstdio>
 3 using namespace std;
 4 double a,b,c,d;
 5 double mid;
 6 int flag=0;
 7 double qz(double p)
 8 {
 9     return a*p*p*p+b*p*p+c*p+d;
10 }
11 
12 void find(double x,double y)
13 {
14 
15     if(y-x>=0.001)
16     {
17         mid=(x+y)/2;
18         if(qz(mid)*qz(x)<0)
19         {
20             find(x,mid);
21         }// 
22         else
23         {
24             find(mid,y);
25         }
26     }
27     if(flag==1)return;
28     printf("%.2lf ",mid);
29     flag=1;
30     return;
31 }
32 int main()
33 {
34     
35     cin>>a>>b>>c>>d;
36     for(int i=-100;i<=100;i++)
37     {
38         double x1=i;
39         double x2=i+1;
40         if(qz(x1)==0)
41         {
42             printf("%.2lf ",x1);
43         }
44         if(qz(x1)*qz(x2)<0)
45         {
46             flag=0;
47             find(x1,x2);
48         }//两个数之间有一个根
49         else
50         continue;    
51     }
52     return 0;
53 }

 

标签:

版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有

上一篇:C++拷贝构造函数(深拷贝与浅拷贝)

下一篇:循环比赛日程表