Алгоритм умножения произвольных многочленов длины?

Question

Допустим, мы представим полиномы как массив поплавки, где степень каждого элемента многочлена соответствует индексу в массиве (например, 4.2x^5+x^2-1.4 будет представлена ​​как {-1.4, 0, 1, 0, 0, 4.2}.

вопрос Моя задача заключается в написании метод в C, который умножает два многочлена произвольной длины и выводит результат (вместо того, чтобы возвращать его).

Обычно, когда я задаю вопросы по этому вопросу, я включаю то, что я пытался сделать до сих пор, но я действительно совершенно незнакомый с этим. Это все, что у меня есть:

void multpoly(float *a, int len_a, float *b, int len_b) 
{ 
    for(i = 0; i < len_result; i++) 
    { 
     printf(" %.5f, ", product[i]); 
    } 
} 

Любая помощь будет принята с благодарностью!

Answer 1

Я считаю, что это то, что вы хотите:

// constraints: 'result' must have space for at least len1 + len2 - 1 elements. 
void multpoly(const float *poly1, int len1, const float *poly2, int len2, float *result) 
{ 
    int i, p1i, p2i; 
    int len_result = len1 + len2 - 1; 

    for (i = 0; i < len_result; i++) result[i] = 0.0; 

    for (p1i = 0; p1i < len1; ++p1i) 
     for (p2i = 0; p2i < len2; ++p2i) 
      result[p1i + p2i] += poly1[p1i] * poly2[p2i]; 
} 

Ideone example of this function

Answer 2

#include <stdio.h> 
#include <string.h> 

main() 
{ 
    float a[5]={1, 2, 3}; 
    float b[5]={2, 0, 1}; 

    printPol(a, b, 3, 3); 

} 

void printPol(float*a, float*b, int len1, int len2) 
{ 
    int i, j; 

    // order of resulting poly is o1+o2 
    // o1 = len1 -1 
    // o2 = len2 -1 
    // length is order + 1 (+1 is the constant number) 
    int len = (len1-1)+(len2-1)+1; 

    float res[len]; 
    //initialize 
    for(i=0;i<len; i++) res[i] = 0; 


    for(i=0; i<len1; i++)  
     for(j=0; j<len2; j++) 
     { 
      // mutually multiply all elements 

      res[i+j] += a[i]*b[j]; 
     } 
    printf("%f ", res[0]); 
    for(i=1;i<len; i++) printf("+%f*x^%d ", res[i], i); 


}