开发者社区> 问答> 正文

c++编写求解一元二次方程组根问题

总时间限制: 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. x1的实部大于x2的实部
  2. x1的实部等于x2的实部且x1的虚部大于等于x2的虚部
    所有实数部分要求精确到小数点后5位,数字、符号之间没有空格。

样例输入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;
}

展开
收起
a123456678 2016-03-09 09:20:13 2939 0
1 条回答
写回答
取消 提交回答
  •  #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;
    }
    2019-07-17 18:55:20
    赞同 展开评论 打赏
问答分类:
C++
问答地址:
问答排行榜
最热
最新

相关电子书

更多
使用C++11开发PHP7扩展 立即下载
GPON Class C++ SFP O;T Transce 立即下载
GPON Class C++ SFP OLT Transce 立即下载