一元二次方程求解

[問題描述]

       解一元二次方程a x²+ b x + c = 0

[輸入格式]

       輸入的第一行為三個(gè)實(shí)數(shù)a, b, c  (|a|,|b|,|c|<10³)，相鄰兩個(gè)數(shù)之間用一個(gè)空格分隔。

[輸出格式]

       若方程有兩個(gè)相同的解，則只輸出一個(gè)數(shù)，為所求得的解。

       若方程有兩個(gè)實(shí)數(shù)解，則在同一行輸出兩個(gè)實(shí)數(shù)，先輸出較小的，后輸出較大的，中間用Tab(“\t”)分隔。

       若方程有兩個(gè)虛數(shù)解，則在同一行輸出兩個(gè)解，前坐標(biāo)軸上方的解先輸出，然后再輸出坐標(biāo)軸下方的解，中間用Tab分隔。

[樣例輸入1]

1 2 3

[樣例輸出1]

-1+i1.41421     -1-i1.41421

[樣例輸入2]

1 2 1

[樣例輸出2]

-1

#include<iostream>
#include<cmath>
using namespace std;
int main()
{
    float a,b,c;
    double    d;
    char i;
    cin>>a>>b>>c;
    d=b*b-4*a*c;
    if (d>0)
    {
        if (a>0)
        {
            cout<<(-b-sqrt(d))/(2*a)<<'\t'<<(-b+sqrt(d))/(2*a)<<endl;}
        if (a<0)
        {    cout<<(-b+sqrt(d))/(2*a)<<'\t'<<(-b-sqrt(d))/(2*a)<<endl;}
    }
    if (d==0)
    {
        cout<<-b/(2*a)<<endl;
    }
    if (d<0)
    {
    if (a>0)
    {
        d=-d;
        double e=sqrt(d)/(2*a);
        cout<<-b/(2*a)<<"+i"<<e<<'\t'<<-b/(2*a)<<"-i"<<e<<endl;
    }
    if (a<0)
    { d=-d;
      double e=-sqrt(d)/(2*a);
      cout<<-b/(2*a)<<"+i"<<e<<'\t'<<-b/(2*a)<<"-i"<<e<<endl;
    }
    }
    return 0;
}