ordinal COM_QQ_D(const polynom& left, const polynom& right) { if (left.degree() < right.degree()) return ordinal::LT; if (left.degree() > right.degree()) return ordinal::GT; for (int i = left.degree(); i >= 0; i--) { auto currComp = COM_QQ_D(left[i], right[i]); if (currComp != ordinal::EQ) return currComp; } return ordinal::EQ; }
polynom DER_P_P(const polynom& polynom_1) { // вспомогательные переменные + степень многочлена polynom polynom_2 = polynom_1; unsigned size_polynom = polynom_1.degree(); fraction Number_1(0), Number_2(1); polynom_2[size_polynom] = Number_1; // высшей степени мн-на присваиваем 0 //дифференцирование многочлена for (unsigned p = 0 ; p < size_polynom; ++p) { Number_1 = ADD_QQ_Q(Number_1, Number_2); polynom_2[p] = MUL_QQ_Q(polynom_1[p+1], Number_1); } polynom_2.reduce(); // уменьшение степени на 1 return polynom_2; }
natural DEG_P_N(const polynom &pn) { natural result(pn.degree()); return result; }