20:求一元二次方程的根

​​原题链接​​

总时间限制: 1000ms 内存限制: 65536kB

描述

利用公式x1 = (-b + sqrt(bb-4ac))/(2a), x2 = (-b - sqrt(bb-4ac))/(2a)求一元二次方程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 / (2a), 虚部 = sqrt(4ac-bb) / (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<cstdio>
#include<algorithm>
using namespace std;
int main()
{
  double a,b,c,x1,x2,s;
  cin>>a>>b>>c;
  s=sqrt(b*b-4*a*c);
  x1=(-1*b+s)*1.0/(2*a);
  x2=(-1*b-s)*1.0/(2*a);
  if(b*b==4*a*c) 
  {
    cout<<"x1=x2=";
    printf("%.5lf",x1);
    cout<<endl;
  }
  if(b*b>4*a*c) 
  {
    cout<<"x1=";
    printf("%.5lf",x1);
    cout<<";";
    cout<<"x2=";
    printf("%.5lf",x2);
    cout<<endl;
  }
  if(b*b<4*a*c) 
  {
    cout<<"x1=";
  if(-1.0*b/(2*a)<0)
  {
    cout<<"-";
    printf("%.5lf",abs(-1.0*b/(2*a)));
  }
  else
    printf("%.5lf",abs(-1.0*b/(2*a)));
    cout<<"+";
    printf("%.5lf",abs(sqrt(4*a*c-b*b)*1.0/(2*a)));
    cout<<"i;";
    cout<<"x2=";
  if(-1.0*b/(2*a)<0)
  {
    cout<<"-";
    printf("%.5lf",abs(-1.0*b/(2*a)));
  }
  else
    printf("%.5lf",abs(-1.0*b/(2*a)));
    cout<<"-";
    printf("%.5lf",abs(sqrt(4*a*c-b*b)*1.0/(2*a)));
    cout<<"i"<<endl;
  }
  return 0;
}

以上代码仅供参考