Решите квадратное уравнение в C++

Что-то вроде этого будет работать:

struct complex { double r,i; }
struct pair<T> { T p1, p2; }

pair<complex> GetResults(double a, double b, double c)
{
  pair<complex> result={0};

  if(a<0.000001)    // ==0
  {
    if(b>0.000001)  // !=0
      result.p1.r=result.p2.r=-c/b;
    else
      if(c>0.00001) throw exception("no solutions");
    return result;
  }

  double delta=b*b-4*a*c;
  if(delta>=0)
  {
    result.p1.r=(-b-sqrt(delta))/2/a;
    result.p2.r=(-b+sqrt(delta))/2/a;
  }
  else
  {
    result.p1.r=result.p2.r=-b/2/a;
    result.p1.i=sqrt(-delta)/2/a;
    result.p2.i=-sqrt(-delta)/2/a;
  }

  return result;
}

Таким образом, вы получаете одинаковые результаты как для реальных, так и для комплексных результатов (у реальных результатов просто мнимая часть установлена ​​на 0). С наддувом смотрелось бы еще красивее!

edit: исправлено для дельты и добавлена ​​проверка вырожденных случаев, таких как a=0. Бессонная ночь ftl!

Важное замечание ко всему этому. Решения, показанные в этих ответах и ​​в исходном вопросе, не являются надежными.

Хорошо известное решение (-b +- sqrt(b^2 - 4ac)) / 2a известно как ненадежное в вычислениях, когда ac очень мало по сравнению с b^2, потому что нужно вычесть два очень похожих значения. Для другого корня лучше использовать менее известное решение 2c / (-b --+ sqrt(b^2 -4ac)).

Устойчивое решение может быть рассчитано как:

temp = -0.5 * (b + sign(b) * sqrt(b*b - 4*a*c);
x1 = temp / a;
x2 = c / temp;

Использование знака (b) гарантирует, что мы не вычитаем два одинаковых значения.

Для OP измените это для комплексных чисел, как показано другими плакатами.

У вас более или менее это есть, просто проверьте, является ли часть, которая находится внутри квадратного корня, отрицательной, а затем отслеживайте это отдельно в своих сокращениях.

Вы могли бы просто использовать std::complex<float> вместо float, чтобы получить поддержку комплексных чисел.

Никая идея от Blindy:

typedef std::complex<double> complex;
using std::pair;
pair<complex> GetResults(double a, double b, double c)
{
  double delta=(b*b-4*a*c);
  double inv_2a = 1/2/a;
  if(delta >= 0) {
    double root = sqrt(delta);
    return std::make_pair(
        complex((-b-root)*inv_2a),
        complex((-b+root)*inv_2a);
  } else {
    double root = sqrt(-delta);
    return std::make_pair(
        complex(-b*inv_2a, -root*inv_2a)),
        complex(-b*inv_2a, +root*inv_2a)));
  }
}

Я пробовал программу без использования заголовка «math.h», а также пробовал другую логику ... но моя программа может отвечать только на те квадратные уравнения, которые имеют коэффициент «x квадрат» как единицу ..... и где коэффициент «x ' можно выразить как сложение двух чисел, которые являются множителями постоянного члена. например. х квадрат +8х+16; х квадрат +7х+12; и т.д. здесь 8=4+4 и 16=4*4; здесь коэффициент при х может быть выражен как сложение двух чисел, которые являются множителями постоянного члена 16... Сам я им не полностью доволен, но попробовал что-то другое, не используя формулу решения квадратного уравнения. код есть;

        #include<iostream.h>
        #include<conio.h>
         class quadratic
              {
                int b,c ;
                float l,k;
                public:
               void solution();
              };
        void quadratic::solution()
             {
                 cout<<"Enter coefficient of x and the constant term of the quadratic eqn where coefficient of x square is one";
                 cin>>b>>c;

                 for(l=1;l<b;l++)
                  {
                   for(k=1;k<b;k++)
                    {
                     if(l+k==b&&l*k==c)
                        {
                          cout<<"x="<<-l<<"\t"<<"or"<<"\t"<<"x="<<-k;
                          cout<<"\n";
                         }
                    }
                }
            }
              void main()
                 {
                  quadratic a;
                   clrscr();
                  a.solution();
                  getch();
                 }