The command SquarefreeFactorization(p, v, rc, R) returns a list of pairs [sqf_i, rc_i].  For each pair, the list sqf_i is a square-free decomposition of p modulo the saturated ideal of rc_i; each element in the list sqf_i is a pair as [s_j,e_j], where s_j is a square-free polynomial modulo rc_i and e_j is the exponent of s_j in p.

All the regular chains from the output pairs form a triangular decomposition of rc in the sense of Kalkbrener.

The option 'method' specifies which gcd algorithm to use. The default option, 'method'='evala', uses a modular algorithm. The other option is 'method'='src', which uses a subresultant-based approach. This method is generally slower, but can be faster in some cases, for instance, if the dimension of the saturated ideal of rc is high, say greater than 4.

Assumptions: the polynomial ring is assumed to have characteristic zero; the initial of p is regular w.r.t. rc; v is greater than the main variables of the regular chain rc in R.

$\mathrm{with}\left(\mathrm{RegularChains}\right)\:$$\mathrm{with}\left(\mathrm{ChainTools}\right)\:$

$R≔\mathrm{PolynomialRing}\left(\left[y\,x\right]\right)$

$R≔\mathrm{polynomial\_ring}$

$q≔\mathrm{expand}\left({\left(x-2\right)}^2{\left(x-1\right)}^2\right)$

$q≔x^4-6x^3+13x^2-12x+4$

$p≔\mathrm{prem}\left(\left(y+6\right)\left(y+2-x\right)\left(x{y}^2-1\right)\left(y-x\right)\,q\,x\right)$

$p≔x^3{y}^3-2x^2{y}^4+x{y}^5+6x^3{y}^2-14x^2{y}^3+8x{y}^4-12x^2{y}^2+12x{y}^3-x^{2}y+2x{y}^2-y^3-6x^2+14yx-8y^2+12x-12y$

$C≔\mathrm{Chain}\left(\left[q\right]\,\mathrm{Empty}\left(R\right)\,R\right)$

$C≔\mathrm{regular\_chain}$

$\mathrm{SquarefreeFactorization}\left(p\,y\,C\,R\right)$

$\left[\left[\left[\left[2y^5+8y^4-25y^3-4y^2+12y\,1\right]\right]\,\mathrm{regular\_chain}\right]\,\left[\left[\left[y+6\,1\right]\,\left[y^2-1\,2\right]\right]\,\mathrm{regular\_chain}\right]\right]$

$\mathrm{SquarefreeFactorization}\left(p\,y\,C\,R\,\mathrm{method}='\mathrm{src}'\right)$

$\left[\left[\left[\left[2y^5+8y^4-25y^3-4y^2+12y\,1\right]\right]\,\mathrm{regular\_chain}\right]\,\left[\left[\left[y+6\,1\right]\,\left[y^2-1\,2\right]\right]\,\mathrm{regular\_chain}\right]\right]$

$R≔\mathrm{PolynomialRing}\left(\left[x\,y\right]\right)$

$R≔\mathrm{polynomial\_ring}$

$p≔\left(x+y\right){\left(x+2y\right)}^2{\left(x+3y\right)}^3$

$p≔\left(x+y\right){\left(x+2y\right)}^2{\left(x+3y\right)}^3$

$\mathrm{rc}≔\mathrm{Chain}\left(\left[y^2+1\right]\,\mathrm{Empty}\left(R\right)\,R\right)$

$\mathrm{rc}≔\mathrm{regular\_chain}$

$p≔\mathrm{SparsePseudoRemainder}\left(p\,\mathrm{rc}\,R\right)$

$p≔x^6+14x^{5}y-80x^4-238x^{3}y+387x^2+324xy-108$

$\mathrm{SquarefreeFactorization}\left(p\,x\,\mathrm{rc}\,R\right)$

$\left[\left[\left[\left[x+y\,1\right]\,\left[x+2y\,2\right]\,\left[x+3y\,3\right]\right]\,\mathrm{rc}\right]\right]$

$\mathrm{rc}≔\mathrm{Chain}\left(\left[y^2+y\right]\,\mathrm{Empty}\left(R\right)\,R\right)$

$\mathrm{rc}≔\mathrm{regular\_chain}$

$p≔\left(x+y\right){\left(x+2y\right)}^2{\left(x+3y\right)}^3$

$p≔\left(x+y\right){\left(x+2y\right)}^2{\left(x+3y\right)}^3$

$p≔\mathrm{SparsePseudoRemainder}\left(p\,\mathrm{rc}\,R\right)$

$p≔x^6+14x^{5}y-80x^{4}y+238x^{3}y-387x^{2}y+324xy-108y$

$\mathrm{SquarefreeFactorization}\left(p\,x\,\mathrm{rc}\,R\right)$

$\left[\left[\left[\left[x-1\,1\right]\,\left[x-2\,2\right]\,\left[x-3\,3\right]\right]\,\mathrm{regular\_chain}\right]\,\left[\left[\left[x\,6\right]\right]\,\mathrm{regular\_chain}\right]\right]$

$\mathrm{SquarefreeFactorization}\left(p\,x\,\mathrm{rc}\,R\,\mathrm{method}='\mathrm{src}'\right)$

$\left[\left[\left[\left[x-1\,1\right]\,\left[x-2\,2\right]\,\left[x-3\,3\right]\right]\,\mathrm{regular\_chain}\right]\,\left[\left[\left[x\,6\right]\right]\,\mathrm{regular\_chain}\right]\right]$

The RegularChains[ChainTools][SquarefreeFactorization] command was introduced in Maple 16.

For more information on Maple 16 changes, see Updates in Maple 16.