TOPICS
Search

Subresultant


A subresultant is a polynomial generalizing the resultant of two univariate polynomials. Subresultants provide a way to compute their polynomial greatest common divisor using determinants formed from their coefficients (Collins 1967, Brown and Traub 1971).

Let f,g in K[x] have polynomial degrees m>=n>0, where K is a field. For 0<=j<n, form the matrix M_j whose rows are the coefficients of

 x^(n-j-1)f,...,f,x^(m-j-1)g,...,g,
(1)

with columns in descending powers from x^(m+n-j-1) to 1. This matrix has r=m+n-2j rows and r+j columns. Let M_(j,k) consist of its first r-1 columns followed by its column for x^k. The jth subresultant is

 S_j(f,g)=sum_(k=0)^jdet(M_(j,k))x^k.
(2)

It has degree at most j, and each coefficient is a polynomial in the coefficients of f and g (El Kahoui 2003, pp. 283-285). In particular, M_0 is the Sylvester matrix, so S_0(f,g) is the resultant.

With b_n the leading coefficient of g, the terminal convention

 S_n(f,g)=b_n^(m-n-1)g
(3)

gives a sequence with n+1 members. When m=n, this last member is g/b_n. Other normalizations of the subresultant sequence are also used.

The principal subresultant coefficient s_j is the coefficient of x^j in S_j(f,g). Writing C(f,g)={s_0,...,s_n}, some simple examples are

C(x-a,x-b)={a-b,1}
(4)
C((x-a)(x-b),x-c)={(a-c)(b-c),1}
(5)
C((x-a)(x-b),(x-c)(x-d))={(a-c)(b-c)(a-d)(b-d),a+b-c-d,1}.
(6)

These are lists of principal subresultant coefficients, rather than lists of the subresultant polynomials themselves. For example, the subresultant polynomial at index 1 for the two quadratic polynomials is

 S_1((x-a)(x-b),(x-c)(x-d))=(a+b-c-d)x+cd-ab.
(7)

If d is the degree of the polynomial greatest common divisor, then all subresultants of index less than d vanish, and S_d is a nonzero scalar multiple of that divisor (El Kahoui 2003, pp. 290-291). Equivalently, the first nonzero principal subresultant coefficient has index d. Over an algebraic closure of K, this counts common roots with each multiplicity equal to the smaller of its multiplicities in f and g.

The subresultant polynomials are computed in the Wolfram Language using SubresultantPolynomials[f, g, x], while their principal subresultant coefficients are computed using Subresultants[f, g, x].


See also

Polynomial Discriminant, Polynomial Greatest Common Divisor, Principal Subresultant Coefficient, Resultant, Sylvester Matrix

Explore with Wolfram|Alpha

References

Akritas, A. G. Elements of Computer Algebra with Applications. New York: Wiley, 1989.Brown, W. S. and Traub, J. F. "On Euclid's Algorithm and the Theory of Subresultants." J. ACM 18, 505-514, 1971. https://doi.org/10.1145/321662.321665.Cohen, H. Ch. 3 in A Course in Computational Algebraic Number Theory. Berlin, Germany: Springer-Verlag, 1993.Collins, G. E. "Subresultants and Reduced Polynomial Remainder Sequences." J. ACM 14, 128-142, 1967. https://doi.org/10.1145/321371.321381.D'Andrea, C.; Krick, T.; and Szanto, A. "Multivariate Subresultants in Roots." J. Algebra 302, 16-36, 2006. https://doi.org/10.1016/j.jalgebra.2005.08.016.Ducos, L. "Optimizations of the Subresultant Algorithm." J. Pure Appl. Algebra 145, 149-163, 2000.El Kahoui, M. "An Elementary Approach to Subresultants Theory." J. Symb. Comput. 35, 281-292, 2003. https://doi.org/10.1016/S0747-7171(02)00135-9.Geddes, K. O.; Czapor, S. R.; and Labahn, G. Algorithms for Computer Algebra. Amsterdam, Netherlands: Kluwer, 1992.Hong, H. "Subresultants Under Composition." J. Symb. Comput. 23, 355-365, 1997.Hong, H. "Subresultants in Roots." Submitted 1999. http://www4.ncsu.edu/~hong/papers/Hong99a.html.Trott, M. "Polynomials in Equations." §1.2.2 in The Mathematica GuideBook for Symbolics. New York: Springer-Verlag, 2006. https://www.mathematicaguidebooks.org/.

Referenced on Wolfram|Alpha

Subresultant

Cite this as:

Weisstein, Eric W. "Subresultant." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/Subresultant.html

Subject classifications