/*======================================================================
求一元二次方程的根
总时间限制: 1000ms 内存限制: 65536kB
描述
利用公式x1 = (-b + sqrt(b*b-4*a*c))/(2*a), x2 = (-b - sqrt(b*b-4*a*c))/(2*a)求一元二次方程ax^2 + bx + c =0的根,其中a不等于0。
输入
第一行是待解方程的数目n。
其余n行每行含三个浮点数a, b, c(它们之间用空格隔开),分别表示方程ax2 + bx + c =0的系数。
输出
输出共有n行,每行是一个方程的根:
若是两个实根,则输出:x1=...;x2 = ...
若两个实根相等,则输出:x1=x2=...
若是两个虚根,则输出:x1=实部+虚部i; x2=实部-虚部i 所有实数部分要求精确到小数点后5位,数字、符号之间没有空格。
x1和x2的顺序:x1的实部>x2的实部||(x1的实部==x2的实部&&x1的虚部>=x2的虚部)
样例输入
3
1.0 3.0 1.0
2.0 -4.0 2.0
1.0 2.0 8.0
样例输出
x1=-0.38197;x2=-2.61803
x1=x2=1.00000
x1=-1.00000+2.64575i;x2=-1.00000-2.64575i
提示
1、需要严格按照题目描述的顺序求解x1、x2。
2、方程的根以及其它中间变量用double类型变量表示。
3、函数sqrt()在头文件math.h中。
4、要输出浮点数、双精度数小数点后5位数字,可以用下面这种形式: printf("%.5f", num); 注意,在使用Java做此题时,可能会出现x1或x2等于-0的情形,此时,需要把负号去掉
========================================================================*/
#include<stdio.h>
#include<math.h>
int comparToZero(double x);//当x与0之差的绝对值小于0.00001(即:1e-5)时 认为x等于0
int main()
{
int n,i;
double a,b,c,delt;
double x1,x2,real,image;
freopen("222.in","r",stdin);
scanf("%d",&n);
for(i=;i<n;i++)
{
scanf("%lf%lf%lf",&a,&b,&c);
delt=b*b-*a*c;
if(comparToZero(delt)==) //if(delt>0)
{
x1=(-b+sqrt(delt))//a;
x2=(-b-sqrt(delt))//a;
if(comparToZero(x1)==)
x1=;
if(comparToZero(x2)==)
x2=;
printf("x1=%.5lf;x2=%.5lf\n",x1,x2);
}
else if(comparToZero(delt)==) //if(delt==0)
{
x1=-b//a;
if(comparToZero(x1)==)
x1=;
printf("x1=x2=%.5lf\n",x1);
}
else //delt<0
{
delt=sqrt(-delt);
real=-b//a;
image=delt//a;
if(comparToZero(real)==) real=;
if(comparToZero(image)==) image==;
printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi\n",real,image,real,image);
}
}
return ;
}
int comparToZero(double x)//当x与0之差的绝对值小于0.00001(即:1e-5)时 认为x等于0
{
/*if(x>0) return 1;
else if(x<0) return -1;
else return 0;*/
//上面的写法会出现浮点数计算精度的误差 。需要改用下面的方法来比较浮点数x和0之间的大小关系。
if((x-)>1e-) return ;
else if((x-)<-(1e-)) return -;
else return ;
}
/*==============================================================================
实数运算中,经常需要判断实数x和y是否相等。编程者往往把判断的条件简单设成x==y或者y-x==0.
实际上,这种写法是有失偏颇的,可能会产生精度误差。
避免精度误差的办法就是设置一个精度常量delta。
若y-x的实数值与0之间的区间长度小于delta,则认为x和y相等,这样就可以把误差控制在delta范围内。
显然,判断实数x和y是否相等的条件应设成|x-y|<=delta。下面是一个实例。
===============================================================================*/
这个题目需要注意的两个地方:
/*==========================================
这个题需要注意的地方:
1. 当一个数如果在(-0.000005,0)之间,输出精确到小数点后5位,
就是0而不是期望的某个小数0.00……,所以需要判断处理一下输出数据
2. 浮点数无法和0比较,浮点数本身不是一个精确值,
可以用fabs(x-0)<1e-5来比较(这里精确度为5)
============================================*/