大清早的我就来刷题了:D

http://noi.openjudge.cn/ch0104/20/

20:求一元二次方程的根

总时间限制:
1000ms
内存限制:
65536kB
描述
利用公式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。

输入
输入一行，包含三个浮点数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 / (2*a), 虚部 = sqrt(4*a*c-b*b) / (2*a)所有实数部分要求精确到小数点后5位，数字、符号之间没有空格。
样例输入
样例输入1
1.0 2.0 8.0

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

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

乍看之下很水么。但对我这种弱鸡来说陷阱还是蛮多的。

QQ图片20160220103530

四个月前..4WA1CE的记录

没心思再去另外写一个代码了，就看看历史代码有没有什么可以查漏补缺吧..

先把样例数据带进去..

QQ图片20160220103831

卧槽样例都过不了？！那么看来我还真是个傻缺..

也没有想到别的办法了，只好用笨办法写了个函数：

double chk(double x){
	if (x<0.00001&&x>-0.00001) return 0.00000;
	return x;
}

嗯，轻喷...

把所有的输出都套上了这个函数，兴冲冲地跑去提交：

结果还是WA，⑨分！

不能忍啊。继续排查。到网上找了几道初中生的一元二次方程题目(⊙﹏⊙)b代入，果然出错了，把两根相同的结果ABS掉了...

改之，AC。

QQ图片20160220104507

附上AC代码：

#include<cstdio>
#include <cmath>
double chk(double x){
if (x<0.00001&&x>-0.00001) return 0.00000;
return x;
}
double abs(double x){
return x<0?-x:x;
}
int main(){
double a,b,c,d;
scanf("%lf %lf %lf",&a,&b,&c);
d=b*b-4*a*c;
if (d<0.00001&&d>-0.00001) d=0;
if (d<0.0){ printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi",chk(-b/(2*a)),chk(sqrt(-(b*b-4*a*c))/(2*a)),chk(-b/(2*a)),chk(sqrt(-(b*b-4*a*c))/(2*a))); }else if(d==0.0){ printf("x1=x2=%.5lf",chk((-b)/(2*a)) ); }else if(d>0.0){
printf("x1=%.5lf;x2=%.5lf",chk((-b + sqrt(b*b-4*a*c))/(2*a)),chk((-b - sqrt(b*b-4*a*c))/(2*a)));
}
return 0;
}

To be continued...