h0053. 求一元二次方程的根 (10 分)

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

输入格式:

第一行是待解方程的数目n。 其余n行每行含三个浮点数a, b, c（它们之间用空格隔开），分别表示方程ax^2 + bx + c =0的系数。

输出格式:

输出共有n行，每行是一个方程的根： 若是两个实根，则输出：x1=...;x2 = ... 若两个实根相等，则输出：x1=x2=... 若是两个虚根，则输出：x1=实部+虚部i; x2=实部-虚部i

所有实数部分要求精确到小数点后5位，数字、符号之间没有空格。 x1和x2的顺序：x1的实部>Re的实部||(x1的实部==x2的实部&&x1的虚部>=x2的虚部)

输入样例:

3
1.0 3.0 1.0
2.0 -4.0 2.0
1.0 2.0 8.0

输出样例:

在这里给出相应的输出。例如：

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

#include<bits/stdc++.h>

using namespace std;

int main()
{
    double a,b,c;
    int n;
    cin>>n;
    while(n--)
    {
        cin>>a>>b>>c;
        double p = b * b - 4 * a * c;
        double x1 = -b / (2 * a);
        if(!x1) x1 = 0;//防止出现 -0 的情况
        if(p > 0) printf("x1=%.5f;x2=%.5f\n",(-b + sqrt(p)) / (2 * a),(-b - sqrt(p)) / (2 * a));
        else if(p == 0) printf("x1=x2=%.5f\n",x1);
        else
        {
            p = -p;
            double y1 = sqrt(p) / (2 * a);
            double y2 = -sqrt(p) / (2 * a);
            printf("x1=%.5f%+.5fi;x2=%.5f%+.5fi\n",x1,y1,x1,y2);
        }
    }
    return 0;
}