18、第四课：求一元二次方程的根（作业）

                                            求一元二次方程的根（作业）

                                                                （时间限制：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。
若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)。

输入

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

输出

输出一行，表示方程的解。所有实数部分要求精确到小数点后5位，数字、符号之间没有空格。

难度

一般

输入示例

1.0 2.0 8.0

输出示例

x1=-1.00000+2.64575i; x2=-1.00000-2.64575i

---

源代码：

#include <cstdio>    
    #include <iostream>    
    #include <cmath>    
    int main() {    
        double a,b,c,d;    
        double x1,x2;    
        double m,n;     
        double esp=pow(10,-7);    
        scanf("%lf%lf%lf",&a,&b,&c);    
        d=b*b-4*a*c;    
        if(d==0||d>0)    
        {    
            if(d==0)    
            {    
                x1=x2=(-b+sqrt(d))/(2*a);    
                printf("x1=x2=%.5lf\n",x1,x2);    
            }    
            else    
            {    
                x1=(-b+sqrt(d))/(2*a);    
                x2=(-b-sqrt(d))/(2*a);    
                printf("x1=%.5lf;x2=%.5lf\n",x1,x2);        
            }    
        }    
        else    
        {    
            m=-b/(2*a)+esp;    
            n=sqrt(-d)/(2*a);    
            printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi\n",m,n,m,n);            
        }       
        return 0;    
    }