解方程

【链接】nbu1042点击打开链接

【题目】

解方程

Time Limit:1000MS  Memory Limit:65536K

Description

给定一个方程的系数，求出它的解

Input

每行三个整数
代表方程
a*x^2+b*x+c=0，保证有解
0 0 0代表结束

Output

输出所有可能的解，每个一行
如果解数值唯一，则输出两次，值为实数按绝对值从大到小排序，相同时先输正数。
复数的输出格式如SAMPLE，保留三位小数，
虚部为正的先输出

Sample Input

1 -2 1

1 2 2

0 0 0

Sample Output

1.000

1.000

-1.000+1.000i

-1.000-1.000i

【代码】

#include<stdio.h>
#include<math.h>
int main()
{
	double a,b,c,x1,x2,t;
	while(scanf("%lf%lf%lf",&a,&b,&c)!=EOF)
	{
		if(a==0&&b==0&&c==0) break;
		if(b*b-4*a*c>=0)
		{
			if(a==0)
			{x1=-c/b;x2=-c/b;}
			else
			{
				x1=(-b+sqrt(b*b-4*a*c))/(2*a);
				x2=(-b-sqrt(b*b-4*a*c))/(2*a);
				if(fabs(x1)<fabs(x2))
				{t=x1;x1=x2;x2=t;}
				else if(fabs(x1)==fabs(x2)&&x1<x2) 
				{t=x1;x1=x2;x2=t;}
			}
			printf("%.3lf\n%.3lf\n",x1,x2);
		}
		else
		{
			if(b==0)
				printf("+%.3lfi\n-%.3lfi",sqrt(c/a),sqrt(c/a));//+没有也能AC
			else
			{
				x1=-b/(2*a);
				x2=fabs(sqrt(4*a*c-b*b)/(2*a));
				printf("%.3lf+%.3lfi\n%.3lf-%.3lfi\n",x1,x2,x1,x2);
			}
		}
	}
	return 0;
}