20:求一元二次方程的根

描述
利用公式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
来源
1709

#include <iostream>
#include <cstdio>
#include <cmath>


//浮点类型比较只能通过之间是否小于某个精度判断
const float precision = 1e-6;

using namespace std;
int main() {
    float a, b, c;
    float x1,x2;
    cin >> a >> b >> c;
    if (a) {
        float q = pow(b,2) - 4 * a * c;
        if (abs(q) < precision) {
            x1 = (-b)/(2 * a);
            x1 = (abs(x1) < precision)?0:x1;
            printf("x1=x2=%.5f", x1);
        } else if (q > 0) {
            x1 = (-b + sqrt(q))/(2 * a);
            x2 = (-b - sqrt(q))/(2 * a);
            x1 = (abs(x1) < precision)?0:x1;
            x2 = (abs(x2) < precision)?0:x2;
            printf("x1=%.5f;x2=%.5f", x1, x2);
        } else {
            float r = -b / (2 * a);
            float i = sqrt(-q)/(2 * a);
            r = (abs(r) < precision)?0:r;
            i = (abs(i) < precision)?0:i;
            printf("x1=%.5f+%.5fi;x2=%.5f-%.5fi",r,i,r,i);
        }
    }
}