总时间限制: 

1000ms

内存限制: 

65536kB

描述

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

输入

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

输出

输出一行，表示方程的解。
若b2 = 4 * a * c,则两个实根相等，则输出形式为：x1=x2=…。
若b2 > 4 * a * c,则两个实根不等，则输出形式为：x1=…;x2 = …，其中x1>x2。
若b2 < 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
2 #include
3 #include
4 #include
5 using namespace std;
6 int main()
7 {
8 double a,b,c;
9 cin>>a>>b>>c;
10 double d=b*b-4*a*c;
11 double der=sqrt(d);
12 double x1=(-b+der)/(2*a);
13 double x2=(-b-der)/(2*a);
14 if(d==0)
15 {
16 cout<0)
19 {
20 if (a>0)
21 {
22 cout << setprecision(5) << fixed << "x1=" << x1 << ";" << "x2=" << x2; 23 } 24 else 25 { 26 cout << setprecision(5) << fixed << "x2=" << x2 << ";" << "x1=" << x1; 27 } 28 } 29 else 30 { 31 double m, n; 32 double esp = pow(10, -7); 33 m = -b / (2 * a) + esp; 34 n = sqrt(-d) / (2 * a); 35 printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi\n", m, n, m, n); 36 /*der = sqrt(-d); 37 double x3 = -b / 2.0 / a; 38 double x4 = der / 2.0 / a; 39 if (b==0) 40 { 41 x3 = 0; 42 } 43 if (a>0)
44 {
45 cout << setprecision(5) << fixed << "x1=" << x3 << "+" << x4 << "i" << ";" << "x2=" << x3 << "-" << x4 << "i" << endl;
46 }
47 else
48 {
49 cout<<setprecision(5)<< fixed<<"x1="<<x3<<"+"<<-x4<<"i"<<";"<<"x2="<<x3<<x4<<"i"<<endl;
50 }*/
51
52
53 }
54
55 return 0;
56 }