This function computes the monic greatest common divisor of two polynomials with coefficients in an algebraic function field or an algebraic number field.

Since the ordering of objects may vary from a session to another, the leading coefficient may change accordingly.

Algebraic functions and algebraic numbers may be represented by radicals or with the RootOf notation (see type,algnum, type,algfun, type,radnum, type,radfun).

The optional arguments p and q are respectively assigned the co-factors of P and Q. Namely, $p=\frac{P}{G}$, $q=\frac{Q}{G}$ where G is the Gcd.

The computation is performed in the domain $K_x$, where:

x is the set of names in P and Q which do not appear inside a RootOf or a radical,

K is a field generated over the rational numbers by the coefficients of P and Q.

The arguments P and Q must be polynomials in x.

Partial factorizations are preserved. The result may be a product of powers of expanded polynomials.

Algebraic numbers and functions occurring in the results are reduced modulo their minimal polynomial (see Normal).

If the RootOf and the radicals defining the algebraic numbers are not independent, a case discussion T of the gcd of P and Q over the direct product of rings F defined by the towers of algebraic extensions is computed first. Then, the result of T is combined. If the gcd values over all the branches of T can be combined into one, namely G, then G is the output. Otherwise, an error is reported.

If the expression contains only algebraic numbers in radical notation (2^(1/2), 3^(1/2), 6^(1/2)), the algebraic quantities are not necessarily independent. A basis over Q for the radicals can be computed by Maple in this case.

To skip the independence checking, use the option 'independent'. However, if a zero-divisor is hit when normalizing P and Q during the computation, an error is reported. Furthermore, the result may be incorrect if the option 'independent' is used but the RootOfs are not actually independent.

If a or b contains functions, their arguments are normalized recursively and the functions are frozen before the computation proceeds.

Other objects are frozen and considered as variables.

$\mathrm{alias}\left(\mathrm{\alpha }=\mathrm{RootOf}\left(x^3+x+1\right)\right)$

$\mathrm{\alpha }$

$\mathrm{P1}≔\left(x-\mathrm{\alpha }\right){\left(x^2-{\mathrm{\alpha }}^2\right)}^2\left(x-1\right)$

$\mathrm{P1}≔\left(x-\mathrm{\alpha }\right){\left(-{\mathrm{\alpha }}^2+x^2\right)}^2\left(x-1\right)$

$\mathrm{Q1}≔\mathrm{evala}\left(\mathrm{Expand}\left({\left(x+\mathrm{\alpha }\right)}^3\right)\right)$

$\mathrm{Q1}≔3{\mathrm{\alpha }}^{2}x+3\mathrm{\alpha }{x}^2+x^3-\mathrm{\alpha }-1$

$\mathrm{G1}≔\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P1}\,\mathrm{Q1}\,'\mathrm{p1}'\,'\mathrm{q1}'\right)\right)$

$\mathrm{G1}≔{\left(x+\mathrm{\alpha }\right)}^2$

$\mathrm{p1}$

${\left(x-\mathrm{\alpha }\right)}^3\left(x-1\right)$

$\mathrm{q1}$

$x+\mathrm{\alpha }$

$\mathrm{P2}≔\left(x-\mathrm{sqrt}\left(t\right)y\right){\left(\mathrm{sqrt}\left(t\right)x+1\right)}^2$

$\mathrm{P2}≔\left(x-\sqrt{t}y\right){\left(\sqrt{t}x+1\right)}^2$

$\mathrm{Q2}≔\mathrm{expand}\left(\mathrm{diff}\left(\mathrm{P2}\,x\right)\right)$

$\mathrm{Q2}≔3t{x}^2+4\sqrt{t}x+1-2t^{\raisebox{1ex}{$3$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}yx-2ty$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P2}\,\mathrm{Q2}\right)\right)$

$x+\frac{1}{\sqrt{t}}$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P2}\,\mathrm{Q2}+\frac{1}{x}\right)\right)$

Error, (in `evala/Gcd/preproc0`) arguments should be polynomials

$\mathrm{infolevel}\left[\mathrm{Gcd}\right]≔1$

${\mathrm{infolevel}}_{\mathrm{Gcd}}≔1$

$\mathrm{P3}≔x^2-\mathrm{RootOf}\left(x^2-x\right)$

$\mathrm{P3}≔x^2-\mathrm{RootOf}\left({\mathrm{\_Z}}^2-\mathrm{\_Z}\right)$

$\mathrm{Q3}≔x^2-2x+1$

$\mathrm{Q3}≔x^2-2x+1$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P3}\,\mathrm{Q3}\right)\right)$

evala/Gcd/preproc0:   "a case discussion result: Array(1..2, {1 = [[1], [RootOf(_Z^2-_Z) = 0]], 2 = [[x-1], [RootOf(_Z^2-_Z) = 1]]}) "
MergeGcds:   "Merge result: #cases: 2; gcd(s): [[[1], Record('property' = isPrime, 'polynomials' = [[_z[1], 1, 1, _z[1], isPrimeIrreducible]], 'type' = regular_chain, 'ModulePrint' = module_print_inner)], [[x-1], Record('property' = isPrime, 'polynomials' = [[_z[1], 1, 1, _z[1]-1, isPrimeIrreducible]], 'type' = regular_chain, 'ModulePrint' = module_print_inner)]]"

Error, (in `evala/Gcd/preproc`) reducible RootOf detected. Substitutions are {RootOf(_Z^2-_Z) = 0, RootOf(_Z^2-_Z) = 1}

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P3}\,\mathrm{Q3}\right)\,'\mathrm{independent}'\right)$

$1$

$\mathrm{P31}≔x^2-\mathrm{RootOf}\left(x^2-x\,\mathrm{index}=1\right)$

$\mathrm{P31}≔x^2$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P31}\,\mathrm{Q3}\right)\right)$

$1$

$\mathrm{P32}≔x^2-\mathrm{RootOf}\left(x^2-x\,\mathrm{index}=2\right)$

$\mathrm{P32}≔x^2-1$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P32}\,\mathrm{Q3}\right)\right)$

$x-1$

$\mathrm{Q4}≔x^2+x$

$\mathrm{Q4}≔x^2+x$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P3}\,\mathrm{Q4}\right)\right)$

evala/Gcd/preproc0:   "a case discussion result: Array(1..2, {1 = [[x], [RootOf(_Z^2-_Z) = 0]], 2 = [[x+1], [RootOf(_Z^2-_Z) = 1]]}) "
MergeGcds:   "Merge result: #cases: 1; gcd(s): [[x+RootOf(_Z^2-_Z), [RootOf(_Z^2-_Z)]]]"

$x+\mathrm{RootOf}\left({\mathrm{\_Z}}^2-\mathrm{\_Z}\right)$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P3}\,\mathrm{Q4}\right)\,'\mathrm{independent}'\right)$

$1$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P31}\,\mathrm{Q4}\right)\right)$

$x$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P32}\,\mathrm{Q4}\right)\right)$

$x+1$

$\mathrm{alias}\left(s=\mathrm{RootOf}\left({\mathrm{\_Z}}^2-2\right)\,t=\mathrm{RootOf}\left({\mathrm{\_Z}}^2-\mathrm{\_Z}+\mathrm{RootOf}\left({\mathrm{\_Z}}^2-2\right)-2\right)\right)\:$

$\mathrm{P5}≔\left(t-s\right)\left(y+1\right)x\cdot 3+\left(s+t-1\right)y^{2}x\cdot 4$

$\mathrm{P5}≔3\left(t-s\right)\left(y+1\right)x+4\left(s+t-1\right)y^{2}x$

$\mathrm{Q5}≔\left(t-s\right)\left(y+1\right)x\cdot 6+\left(s+t-1\right)y{\left(y-1\right)}^{2}x\cdot 8$

$\mathrm{Q5}≔6\left(t-s\right)\left(y+1\right)x+8\left(s+t-1\right)y{\left(y-1\right)}^{2}x$

$\mathrm{evala}\left(\mathrm{Gcd}\left(\mathrm{P5}\,\mathrm{Q5}\right)\right)$

evala/Gcd/preproc0:   "a case discussion result: Array(1..2, {1 = [[(y+1)*x], [RootOf(_Z^2-_Z+RootOf(_Z^2-2)-2) = -RootOf(_Z^2-2)+1]], 2 = [[y*x], [RootOf(_Z^2-_Z+RootOf(_Z^2-2)-2) = RootOf(_Z^2-2)]]}) "
MergeGcds:   "Merge result: #cases: 1; gcd(s): [[y*x+(1/7*(-2*s-1)*t+1/7*s+4/7)*x, [s, t]]]"

$yx+\left(\frac{\left(-2s-1\right)t}{7}+\frac{s}{7}+\frac{4}{7}\right)x$