002:求一元二次方程的根

总时间限制: 

1000ms

 

内存限制: 

65536kB

描述

利用公式x1 = (-b + sqrt(b*b-4*a*c))/(2*a), x2 = (-b - sqrt(b*b-4*a*c))/(2*a)求一元二次方程ax²+ bx + c =0的根，其中a不等于0。

输入

输入一行，包含三个浮点数a, b, c（它们之间以一个空格分开），分别表示方程ax ² + bx + c =0的系数。

输出

输出一行，表示方程的解。
若b ² = 4 * a * c,则两个实根相等，则输出形式为：x1=x2=...。
若b ² > 4 * a * c,则两个实根不等，则输出形式为：x1=...;x2 = ...，其中x1>x2。
若b ² < 4 * a * c，则有两个虚根，则输出：x1=实部+虚部i; x2=实部-虚部i，即x1的虚部系数大于等于x2的虚部系数，实部为0时不可省略。实部 = -b / (2*a), 虚部 = sqrt(4*a*c-b*b) / (2*a)

所有实数部分要求精确到小数点后5位，数字、符号之间没有空格。

样例输入

样例输入1
1.0 2.0 8.0

样例输入2
1 0 1

样例输出

样例输出1
x1=-1.00000+2.64575i;x2=-1.00000-2.64575i

样例输出2
x1=0.00000+1.00000i;x2=0.00000-1.00000i

 1 #include<iostream>
 2 #include <math.h>
 3 #include <cmath>
 4 using namespace std;
 5 
 6 int main(){    
 7     double a=1.0,b=2.0,c=8.0;
 8     double x1=0.0,x2=0.0;
 9     double eps=0.000001;
10     //cin>>a>>b>>c;
11     scanf("%lf%lf%lf",&a,&b,&c);
12     double tmp = b*b - 4*a*c;
13     if(tmp<eps && tmp > -eps){
14         printf("x1=x2=%.5f",(-b)/(2*a)+eps);
15     }else if(tmp >eps){
16         x1=(-b+sqrt(tmp))/(2*a);
17         x2=(-b-sqrt(tmp))/(2*a)+eps;
18         if((x1-x2)>eps){
19             printf("x1=%.5f;x2=%.5f",x1+eps,x2+eps);
20         }else{
21             printf("x1=%.5f;x2=%.5f",x2+eps,x1+eps);
22         }
23     }else{
24         printf("x1=%.5f+%.5fi;x2=%.5f-%.5fi",
25         (-b)/(2*a)+eps,sqrt(-tmp)/(2*a)+eps,(-b)/(2*a)+eps,sqrt(-tmp)/(2*a)+eps);
26 
27     }
28 
29 
30 
31 
32 return 0;
33 }