总时间限制: 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的系数。输出输出一行,表示方程的解。
若两个实根相等,则输出形式为:x1=x2=...。
若两个实根不等,则输出形式为:x1=...;x2 = ...,其中x1若是两个虚根,则输出:x1=实部+虚部i; x2=实部-虚部i,其中x1,x2满足以下两个条件中的一个:
样例输入1.0 2.0 8.0
样例输出x1=-1.00000+2.64575i;x2=-1.00000-2.64575i
我是在OPENJUDGE上面做的这题,拼死了还是只有7分不能ACCEPTED;
我的代码如下:
#include<cmath>
#include<cstdio>
#include<iostream>
using namespace std;
int main(){
double a , b , c , x1 , x2;
double L;
double derta;
scanf("%lf%lf%lf",&a,&b,&c);
L = b*b - 4*a*c;
if(L >= 0)
derta = sqrt(L);
else
derta = sqrt(-L);
x1 = (-b + derta)/2*a;
x2 = (-b - derta)/2*a;
if(L > 0){
double temp;
if(x1 < x2){
temp = x1;
x1 = x2;
x2 = temp;
}
printf("x1=%.5lf;x2=%.5lf", x1 , x2);
}
if(L == 0)
printf("x1=x2=%.5lf",x1);
if(L < 0){
if(b != 0){
double x , y , temp;
x = derta/(2*a);
y = (derta)/(2*a);
if(x < y){
temp = x;
x = y;
y = temp;
}
printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi",(-b)/(2*a),x,(-b)/(2*a),y);
}
else
printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi",b/(2*a),(derta)/(2*a),b/(2*a),(derta)/(2*a));
}
return 0;
}
#include <iostream>
#include<iomanip>
#include <cmath>
#include <complex>
using namespace std;
static const double e = 1e-12;
bool operator == (complex<double> c1, complex<double> c2) { return abs(c1-c2) < e;}
int main()
{
complex<double> a,b,c;
complex<double> x1,x2;
cin >> a >> b >> c;
x1 = (-b + sqrt(b*b-a*c*4.0))/(a*2.0);
x2 = (-b - sqrt(b*b-a*c*4.0))/(a*2.0);
cout << setiosflags(ios::fixed);
cout.precision(6);
if ( abs(x1.imag()) < e ) // real
{
if (x1 == x2) {
cout << "x1=x2=" << x1.real();
} else {
cout << "x1=" << x1.real() <<";x2=" << x1.real();
}
}
else {
cout << "x1=" << x1.real()<<"+"<<x1.imag()<<"i;"
<<"x2=" << x2.real()<<"+"<<x2.imag()<<"i";
}
return 0;
}
版权声明:本文内容由阿里云实名注册用户自发贡献,版权归原作者所有,阿里云开发者社区不拥有其著作权,亦不承担相应法律责任。具体规则请查看《阿里云开发者社区用户服务协议》和《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容,填写侵权投诉表单进行举报,一经查实,本社区将立刻删除涉嫌侵权内容。