алгоритм деления в клен - PullRequest
Новая
       10

алгоритм деления в клен

0 голосов
/ 25 октября 2018

Я пытаюсь написать функцию, которая определяет остаток от алгоритма длинного деления. В частности, найдите остаток от многочлена f на делении на F, где F - это семейство многочленов.Я застрял, и вот мой код. 

    with(Groebner); LT := proc (f, order) 
    return LeadingCoefficient(f, order)*LeadingMonomial(f, order) 
    end proc

    with(Groebner);  
    CHIA:=proc(f,Lf,order)  local Lq,p,i,j,r,divisionoccurred;    
    for j from 1 to nops(Lf) do  
    Lq[j]:=0;  
    od;  
    r:=0;p:=f;  
    while p<>0 do  
    i:=1  
    divisionoccurred:=false  
    while i<=nops(Lf) and divisionoccurred=false do  
    if divide(LT(p,order),LT(Lf[i],order)) then  
    Lq[i]:=Lq[i]+(LT(p,order))/(LT(Lf[i],order));  
    p:=p-((LT(p,order))/(LT(Lf[i],order)))*Lf[i]; 
    divisionoccurred:=true;  
    else  
    i:=i+1;
    fi;
    end do;
    if divisionoccurred=false then   
    r:=r+LT(p,order);  
    p:=p-LT(p,order);
    fi;
    end do;
    return Lq,r   
    end proc;
...