求一元二次方程的根——慕课程序设计与算法（一）练习题

题目来源

程序设计实习MOOC / 程序设计与算法（一）测验汇总(2022春季)的第十三题。

---

题目：求一元二次方程的解

`
总时间限制: 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位，数字、符号之间没有空格。

最初按照题目的写法如下：

#include<iostream>
#include<cmath>
using namespace std;
int main()
{	
	double a,b,c;
	cin>>a>>b>>c;
	if(b*b==4*a*c)
	{
		double k=(-b)/2/a;
		printf("x1=x2=%.5lf",k);
	}
	if(b*b>4*a*c)
	{
		double k1=(-b+sqrt(b*b-4*a*c))/2/a;
		double k2=(-b-sqrt(b*b-4*a*c))/2/a;
		printf("x1=%.5lf;x2=%.5lf",k1,k2);
	}
	if(b*b<4*a*c)
	{
		double m=(-b)/2/a;
		double k1=sqrt(4*a*c-b*b)/2/a;
		printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi",m,k1,m,k1);
	}
	return 0;
} 

但是提交后发现有样例错误，其中输入为 1 0 7.3
注意到原来程序跑出的结果重有0 并且 0带了前缀负号。
即把0.00000输出成了-0.00000，从而导致出错.
下面给出更正解法：

#include<iostream>
#include<cmath>
using namespace std;
int main()
{	
	double a,b,c;
	cin>>a>>b>>c;
	if(b==-b)
	b=-b;
	if(b*b==4*a*c)
	{
		double k=(-b)/2/a;
		printf("x1=x2=%.5lf",k);
	}
	if(b*b>4*a*c)
	{
		double k1=(-b+sqrt(b*b-4*a*c))/2/a;
		double k2=(-b-sqrt(b*b-4*a*c))/2/a;
		printf("x1=%.5lf;x2=%.5lf",k1,k2);
	}
	if(b*b<4*a*c)
	{
		double m=(-b)/2/a;
		double k1=sqrt(4*a*c-b*b)/2/a;
		printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi",m,k1,m,k1);
	}
	return 0;
}

比起第一次提交多了下列步骤：

if(b==-b)
	b=-b;

判断b是否为零，为零则取消之后的-b。

---

总结
在输出为 0 的情况下 若前置有更改数据 0 的负号 则其有可能输出-0；