C++分支结构,求一元二次方程的根

总时间限制: 

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

#include<iostream>
#include<cmath>
#include<iomanip>
#include <cstdio>
using namespace std;
int main()
{
    double a,b,c;
    cin>>a>>b>>c;
    double d=b*b-4*a*c;
    double der=sqrt(d);
    double x1=(-b+der)/(2*a);
    double x2=(-b-der)/(2*a);
    if(d==0)
    {
        cout<<setprecision(5)<<fixed<<"x1=x2="<<x1;
    }
    else if (d>0)
    {
        if (a>0)
        {
            cout << setprecision(5) << fixed << "x1=" << x1 << ";" << "x2=" << x2;
        }
        else
        {
            cout << setprecision(5) << fixed << "x2=" << x2 << ";" << "x1=" << x1;
        }
    }
    else
    {
        double m, n;
        double esp = pow(10, -7);
        m = -b / (2 * a) + esp;
        n = sqrt(-d) / (2 * a);
        printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi\n", m, n, m, n);
        /*der = sqrt(-d);
        double x3 = -b / 2.0 / a;
        double x4 = der / 2.0 / a;
        if (b==0)
        {
            x3 = 0;
        }
        if (a>0)
        {
            cout << setprecision(5) << fixed << "x1=" << x3 << "+" << x4 << "i" << ";" << "x2=" << x3 << "-" << x4 << "i" << endl;
        } 
        else
        {
            cout<<setprecision(5)<< fixed<<"x1="<<x3<<"+"<<-x4<<"i"<<";"<<"x2="<<x3<<x4<<"i"<<endl;
        }*/

        
    }

    return 0;
}