急！一元二方程根求解，OpenJudge测试不通过，有没有Baby施以援手~

Open Judge NOI组：1.4编程基础之逻辑表达式与条件分支：
测试总是不通过，有没有Baby看看问题在哪里！

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

C语言：
#include <stdio.h> 
#include <stdlib.h>
#include <math.h> 
int main(){
    double a,b,c;
    scanf("%lf %lf %lf",&a,&b,&c);
    
    if(b*b==4*a*c) {
        printf("x1=x2=%lf\n",-b/(2*a));
    }
    else if(b*b>4*a*c){
        printf("x1=%.5lf;x2=%.5lf\n",(-b+sqrt(b*b-4*a*c))/(2*a),(-b-sqrt(b*b-4*a*c))/(2*a));
    }
    else if(b*b<4*a*c){
        if(b==0){
            printf("x1=0.00000+%.5lfi;x2=0.00000-%.5lfi\n",sqrt(4*a*c-b*b)/(2*a),sqrt(4*a*c-b*b)/(2*a));
                } else{
                    printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi\n",-b/(2*a),sqrt(4*a*c-b*b)/(2*a),-b/(2*a),sqrt(4*a*c-b*b)/(2*a));
                }
    }
    return 0;
}

你看看这篇博客：


里面的代码是这样的，是博主zoe_ed的代码，引用：

#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
#define ll long long 
int main()
{
    freopen("in.txt","r",stdin);
    double y;char ch;
    double a,b,c;
    scanf("%lf %lf %lf",&a,&b,&c);
    double dt=b*b-4*a*c;
    if(dt==0) 
    {
        double x1=-b/(2*a);
        if (x1==0) x1=fabs(x1);
        printf("x1=x2=%.5f\n",x1);
        return 0;
    }
    if(dt>0)
    {   double x1=(-b + sqrt(dt))/(2*a),x2=(-b - sqrt(dt))/(2*a);
        if(x1==0) x1=fabs(x1);
        if(x2==0) x2=fabs(x2);
        printf("x1=%.5f;x2=%.5f\n", x1,x2);
    }
    else 
    {
        double x=sqrt(fabs(dt))/(2*a);
        double y=-b/(2*a);
        if(y==0) y=fabs(y);

        printf("x1=%.5f+%.5fi;x2=%.5f-%.5fi\n",y,x,y,x);
    }
    return 0;
}

你稍微改一下输入

这样写，供参考：

#include <stdio.h>
#include <math.h>
int main()
{
    double a, b, c, disc, x1, x2, realpart, imagpart;
    scanf("%lf %lf %lf", &a, &b, &c);
    if (fabs(a) <= 1e-6)
        printf("x1=x2=%.5lf", -(c / b));
    else{
        disc = b * b - 4 * a * c;
        if (fabs(disc) <= 1e-6)
        {
            printf("x1=x2=%.5lf\n", -b / (2 * a));
        }
        else
        {
            if (disc > 1e-6)
            {
                x1 = (-b + sqrt(disc)) / (2 * a);
                x2 = (-b - sqrt(disc)) / (2 * a);
                printf("x1=%.5lf;x2=%.5lf\n", x1, x2);
            }
            else
            {
                realpart = -b / (2 * a);
                if (fabs(realpart) <= 1e-6)
                    realpart = fabs(realpart);
                imagpart = sqrt(-disc) / (2 * a);
                printf("x1=%.5lf+%.5lfi;",  realpart, imagpart);
                printf("x2=%.5lf-%.5lfi\n", realpart, imagpart);
            }
        }
    }
    return 0;
}