c++中的cardano方法

Question

我试图用cardano的方法制作一个计算三次函数根的程序

这是我的密码：

#include<iostream>
#include<cmath>
using namespace std;

int main()
{
    double a, b, c, d, value;
    double f, g, h;
    double i, j, k, l, m, n, p, po;
    double r, s, t, u;
    double x1, x2, x2re, x2im, x3re, x3im, x3;

    cin >> value;
    for(int w=1; w <= value; w++){
       cin >> a >> b >> c >> d;
       cout << "CUBIC EQUATION : " << a << " x^3 + " << b << " x^2 + " << c <<" x + " << d << " = 0" << endl;

       f = ((3*c/a)-((b*b)/(a*a)))/3;
       g = ((2*(b*b*b)/(a*a*a))-(9*b*c/(a*a))+(27*d/a))/27;   
       h = ((g*g)/4)+((f*f*f)/27);


    if(f==0 && g==0 && h==0){     // all roots are real and equal
       x1 = pow((d/a),0.33333333333333333333333333333333);
       x2 = pow((d/a),0.33333333333333333333333333333333);
       x3 = pow((d/a),0.33333333333333333333333333333333);
       cout << "x = " << x1 << endl;
       cout << "x = " << x2 << endl;
       cout << "x = " << x3 << endl;
       }
    else if(h<=0){         // all 3 roots are real
       i = pow((((g*g)/4)-h),0.5);
       j = pow(i,0.33333333333333333333333333333333);
       k = acos((g/(2*i))*-1);
       l = j * -1;
       m = cos(k/3);
       n = sqrt(3) * sin(k/3);
       p = (b/(3*a))*-1;
       x1 = (2*j)*m-(b/(3*a));
       cout << "x = " << x1 << endl;
       x2 = l * (m+n) + p;
       cout << "x = " << x2 << endl;
       x3 = l * (m-n) + p;
       cout << "x = " << x3 << endl;
       }
    else if(h>0){
       r = ((g/2)*-1)+pow(h,0.5);
       s = pow(r,0.33333333333333333333333333333333);
       t = ((g/2)*-1)-pow(h,0.5);
       u = pow((t),0.33333333333333333333333333333333);
       x1 = (s+u) - (b/(3*a));
       cout << "x = " << x1 << endl;
       x2re = (((s+u)*-1)/2) - (b/(3*a));
       x2im = -(s-u)*pow(3,0.5)/2;
       cout << "x = (" << x2re << "," << x2im << ")" << endl;
       x3re = (((s+u)*-1)/2) - (b/(3*a));
       x3im = (s-u)*pow(3,0.5)/2;
       cout << "x = (" << x3re << "," << x3im << ")" << endl;
       }
       }

    return 0;
}

有人能在这段代码中帮助实现用户定义的ComplexNumber吗？我想使用这个链接function，但我无法理解这一点。

Answer 1

有人能在这段代码中帮助实现用户定义的ComplexNumber吗？

不要实现用户定义的ComplexNumber类型.使用语言提供的语言。只有#include <complex>。这样，您就可以通过complex<double> variable_name拥有一个复杂的变量。

代码注释：

1. 最好使用std::sqrt(x)而不是std::pow(x, 0.5)。
2. 如果您在符合POSIX标准的机器上，则数学库中有一个立方体根函数，头中有cbrt(double)。(它可能不会导出到C++头。)这也比std::pow(x,0.33333333333333333333333333333333)更可取。