第三周作业——求一元二次方程的根

Description

利用公式x1 = (-b + sqrt(b*b-4*a*c))/(2*a), x2 = (-b - sqrt(b*b-4*a*c))/(2*a)求一元二次方程ax2+ bx + c =0的根，其中a不等于0。

Input

输入一行，包含三个浮点数a, b, c（它们之间以一个空格分开），分别表示方程ax2 + bx + c =0的系数。

Output

输出一行，表示方程的解。

若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 / (2*a), 虚部 = sqrt(4*a*c-b*b) / (2*a)
 

所有实数部分要求精确到小数点后5位，数字、符号之间没有空格。

Sample Input

样例输入1
1.0 2.0 8.0

样例输入2
1 0 1

 

Sample Output

样例输出1
x1=-1.00000+2.64575i;x2=-1.00000-2.64575i

样例输出2
x1=0.00000+1.00000i;x2=0.00000-1.00000i

 

Solution

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
#define EPS 1e-7
int main()
{
	double a,b,c; //不要用 float,精度不够
	scanf("%lf%lf%lf",&a,&b,&c);
	double tmp = b*b - 4*a*c;
	if( tmp < EPS && tmp > - EPS)
		printf("x1=x2=%.5f",(-b)/(2*a)+EPS);  // + EPS是为了防止出现 -0.00000 
	else if( tmp > EPS) {
		double x1 = (-b+sqrt(tmp))/(2*a);
		double x2 = (-b-sqrt(tmp))/(2*a)+EPS;
		if( x1 - x2 > EPS)
			printf("x1=%.5f;x2=%.5f",x1+EPS,x2+EPS);
		else
			printf("x1=%.5f;x2=%.5f",x2+EPS,x1+EPS);
	}
	else {
		printf("x1=%.5f+%.5fi;x2=%.5f-%.5fi",(-b)/(2*a)+EPS,sqrt(-tmp)/(2*a)+EPS,(-b)/(2*a)+EPS,sqrt(-tmp)/(2*a)+EPS);
	}
	return 0;
}

 

做浮点数计算相关的题目要注意以下几点：

 

1) 用double ，不要用 float，float精度不够

2) double 类型变量是用 %lf 读入的，不是 %f。输出double值可以用 %f

3)不要用 a == b的办法判断两个浮点数是否相等，包括不要用  a== 0的办法判断浮点数 a是否等于0，因为浮点数是有误差的。 应该用  a-b > -eps && a-b < eps  ，即a和b的差的绝对值小于某个很小值 eps的办法来判断a和b是否相等。如果结果要保留小数点后面n位，那么 eps可以取 10的-(n+1)次方

4) 为避免出现 -0.000这样的输出，对要输出的值要先判断是不是足够接近0，如果是，就直接输出 0.0。 例如如果要保留小数点后面 n位，那么 eps 就取 10的-(n+1)次方，然后 认为 数的绝对值小于 eps时，就是0

5) 最好不要用 a!=b 判断a,b是否相等, a > b, a < b 比较a,b大小。和4)类似，用两者的差的绝对值和 eps比较