- 总时间限制:
- 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
1 #include<iostream> 2 #include<cmath> 3 #include<iomanip> 4 #include <cstdio> 5 using namespace std; 6 int main() 7 { 8 double a,b,c; 9 cin>>a>>b>>c; 10 double d=b*b-4*a*c; 11 double der=sqrt(d); 12 double x1=(-b+der)/(2*a); 13 double x2=(-b-der)/(2*a); 14 if(d==0) 15 { 16 cout<<setprecision(5)<<fixed<<"x1=x2="<<x1; 17 } 18 else if (d>0) 19 { 20 if (a>0) 21 { 22 cout << setprecision(5) << fixed << "x1=" << x1 << ";" << "x2=" << x2; 23 } 24 else 25 { 26 cout << setprecision(5) << fixed << "x2=" << x2 << ";" << "x1=" << x1; 27 } 28 } 29 else 30 { 31 double m, n; 32 double esp = pow(10, -7); 33 m = -b / (2 * a) + esp; 34 n = sqrt(-d) / (2 * a); 35 printf("x1=%.5lf+%.5lfi;x2=%.5lf-%.5lfi\n", m, n, m, n); 36 /*der = sqrt(-d); 37 double x3 = -b / 2.0 / a; 38 double x4 = der / 2.0 / a; 39 if (b==0) 40 { 41 x3 = 0; 42 } 43 if (a>0) 44 { 45 cout << setprecision(5) << fixed << "x1=" << x3 << "+" << x4 << "i" << ";" << "x2=" << x3 << "-" << x4 << "i" << endl; 46 } 47 else 48 { 49 cout<<setprecision(5)<< fixed<<"x1="<<x3<<"+"<<-x4<<"i"<<";"<<"x2="<<x3<<x4<<"i"<<endl; 50 }*/ 51 52 53 } 54 55 return 0; 56 }
原文链接: https://www.cnblogs.com/dss-99/p/14021341.html
欢迎关注
微信关注下方公众号,第一时间获取干货硬货;公众号内回复【pdf】免费获取数百本计算机经典书籍
原创文章受到原创版权保护。转载请注明出处:https://www.ccppcoding.com/archives/205340
非原创文章文中已经注明原地址,如有侵权,联系删除
关注公众号【高性能架构探索】,第一时间获取最新文章
转载文章受原作者版权保护。转载请注明原作者出处!