c语言解一元二次方程虚根oj,OJ 1040 Problem J 求一元二次方程的根

学校OJ作业题，挺吓人的。

一元二次方程的标准形式为ax2+bx+c=0(a≠0)，其中a、b、c为常数。求解一元二次方程的根x时有三种情况，分别为(记Δ=b2-4ac)：

1. Δ>0，有两个不等的实根；

2. Δ=0，有两个相同的实根；

3. Δ<0，有两个共轭的虚根。

Input

输入为多行，每行为一元二次方程的三个常数a,b,c，在double类型范围之内。当输入的a为0时，表示输入结束。

Output

每行输入的样例对应三行输出。

第一行输出为样例的编号。

第二行输出为所输入常数a,b,c对应的一元二次方程的标准形式，要求输出满足a>0。

第三行输出为所输入方程的根，分为三种情况：

1. 若方程满足Δ>0，即有两不等实根x1、x2，则按顺序(先小后大)输出这两个实根。

2. 若方程满足Δ=0，即有两相同实根x，则输出一个实根。

3. 若方程满足Δ<0，即有两共轭的虚根x1、x2，则输出两个虚根，虚部符号为正的(即u+vi形式)先输出，虚部符号为负的(x-yi形式)后输出。

以上输出均不输出数学上无意义或可省略的的符号，所有数值最多保留6位有效数字。每个样例之后都有一个空行分隔。

Sample Input

1 2 1

-1 2 -1

-5 2 -0.2

-3 2 0

3 0 12

2 4 4

0

Sample Output

Case 1 :

x^2 + 2x + 1 = 0

only one real root : -1

Case 2 :

x^2 - 2x + 1 = 0

only one real root : 1

Case 3 :

5x^2 - 2x + 0.2 = 0

only one real root : 0.2

Case 4 :

3x^2 - 2x = 0

two real roots : 0, 0.666667

Case 5 :

3x^2 + 12 = 0

two imaginary roots : 2i, -2i

Case 6 :

2x^2 + 4x + 4 = 0

two imaginary roots : -1+i, -1-i

HINT

输出方程格式的各种情况要想清楚，这一部分测试数据给的很全面。另一个就是浮点数的精度控制，这一部分sample给出了例子。

值得注意的是，linux下gcc编译的浮点数运算结果有-0，这是OJ系统Judge端使用的系统；而windows XP下的minGW编译器和VC6不会产生-0，只会输出0；但windows 7下的minGW编译器是能够产生-0的(确实很诡异)。因此使用windows XP的同学忽略了对结果为0的检测，程序需要对结果为0的情况进行全面考虑，确保正确的输出0。这个问题卡了好些同学好几天。

关于是否会产生-0，输出表达式0.0/-1的结果就能测试出来。浮点数从负数方向运算出结果为0，则浮点值为-0是符合C语言浮点数运算规则的，目前尚不清楚windows XP系统不能产生-0的原因。

#include

#include

int main()

{

double a, b, c, dt;

int case_ = 0;

while(scanf("%lf", &a) == 1&& fabs(a) > 0.0000001)

{

scanf("%lf%lf", &b, &c);

if(a < -0.0000001)

{

a = -a;

b = -b;

c = -c;

}

if(case_)

printf("\n");

++case_;

printf("Case %d :\n", case_);

(a == 1)?printf("x^2"):printf("%lgx^2", a);

if(fabs(b) < 0.0000001);

else if(b < -0.0000001&& fabs(fabs(b) - 1)> 0.0000001)

printf(" - %lgx", fabs(b));

else if(b > 0.0000001&& fabs(fabs(b) - 1)> 0.0000001)

printf(" + %lgx", fabs(b));

else if(b < 0.0000001&& fabs(fabs(b) - 1)< 0.0000001)

printf(" - x");

else if(b > 0.0000001&& fabs(fabs(b) - 1)< 0.0000001)

printf(" + x");

if(fabs(c) < 0.0000001);

else if(c > 0.0000001)

printf(" + %lg", c);

else if(c < -0.0000001)

printf(" - %lg", fabs(c));

printf(" = 0\n");

double x1, x2, temp;

dt = b * b - 4 * a * c;

if(dt > 0.0000001)

{

x1 = (- b - sqrt(dt)) / (2 * a);

x2 = (- b + sqrt(dt)) / (2 * a);

if(x1 - x2 > 0.0000001)

{

temp = x1;

x1 = x2;

x2 = temp;

}

if(fabs(x1) < 0.0000001)

x1 = 0;

if(fabs(x2) < 0.0000001)

x2 = 0;

printf("two real roots : %lg, %lg\n", x1, x2);

}

if(fabs(dt) < 0.0000001)

{

x1 = (-b) / (2 * a);

if(fabs(x1) < 0.0000001)

x1 = 0;

printf("only one real root : %lg\n", x1);

}

if(dt < -0.0000001)

{

if((-b / (2 * a)) == 0&& (fabs(sqrt(4 * a * c-b * b) / (2 * a) - 1 ) > 0.0000001))

{

printf("two imaginary roots : %lgi, ", sqrt(4 * a * c-b * b) / (2 * a));

printf("-%lgi\n", sqrt(4 * a * c-b * b) / (2 * a));

}

else if((-b / (2 * a)) == 0&& (fabs(sqrt(4 * a * c-b * b) / (2 * a) - 1 ) < 0.0000001))

{

printf("two imaginary roots : i, ");

printf("-i\n");

}

else if(fabs(sqrt(4 * a * c-b * b) / (2 * a) - 1 ) < 0.0000001)

{

printf("two imaginary roots : %lg+i, ", (-b / (2 * a)));

printf("%lg-i\n", (-b / (2 * a)));

}

else

{

printf("two imaginary roots : %lg+%lgi, ", (-b / (2 * a)), sqrt(4 * a * c-b * b) / (2 * a));

printf("%lg-%lgi\n", (-b / (2 * a)), sqrt(4 * a * c-b * b) / (2 * a));

}

}

}

return 0;

}