题目来源
程序设计实习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;