Умножить два полинома заданных степеней

#include <iostream>
 
int main()
{
    double   a[] = {-7,4,2}, // 2*x^2 + 4x - 7
             b[] = {1,1,-3,4}; // 4*x^3 - 3*x^2 + x + 1
    unsigned A = sizeof(a)/sizeof(*a),
             B = sizeof(b)/sizeof(*b),
             R = (A-1)*(B-1);
    double *r = new double[R];
    
    for (unsigned i=0; i < R; i++) {
        r[i] = 0;
    }
    
    for (unsigned i=0; i < A; i++) {
        for (unsigned j=0; j < B; j++) {
            if (i+j < R) {
                r[i+j] += a[i]*b[j];
            }
        }
    }
    
    for (int i = R-1; i >= 0; i--) {
        std::cout << r[i] << ' ';
    }
}

#include <iostream>
 
void product(double *a, unsigned A, double *b, unsigned B, double *&r, unsigned &R)
{
    R = (A-1)*(B-1);
    r = new double[R];
    for (unsigned i=0; i < R; i++)
        r[i] = 0;
    for (unsigned i=0; i < A; i++)
        for (unsigned j=0; j < B; j++)
            if (i+j < R) r[i+j] += a[i]*b[j];
}
 
int main()
{
    double   a[] = {-7,4,2},
             b[] = {1,1,-3,4},
             *r = NULL;
    unsigned R = 0;
    product(a,3,b,4,r,R);
    for (int i = R-1; i >= 0; i--) {
        std::cout << r[i] << ' ';
    }
}