Computing greatest common divisor of several parametric univariate polynomials via generalized subresultant polynomials
Hoon Hong, Jing Yang
TL;DR
The paper tackles the problem of computing the gcd of several univariate polynomials with parametric coefficients by partitioning the parameter space into cells where the gcd expression is uniform. It introduces a natural extension of subresultants to multiple polynomials and develops generalized Sylvester's and Habicht's theories, yielding determinant-based formulas that remain non-nested and easier to analyze. The authors propose two core algorithms, PGCD and EPGCD, which compute parametric gcds by leveraging generalized Sylvester and Habicht relations to produce explicit gcd expressions $G_i$ on cells $C_i$ with conditions $C_i$. Performance assessments show substantial improvements in output size and running time compared with recursive and matrix-based approaches, making the method scalable to larger families of polynomials. The work provides a concrete, efficient framework for parametric gcds with multiple polynomials, enabling structural analysis and robust symbolic computation in parameterized settings.
Abstract
In this paper, we tackle the following problem: compute the gcd for several univariate polynomials with parametric coefficients. It amounts to partitioning the parameter space into ``cells'' so that the gcd has a uniform expression over each cell and constructing a uniform expression of gcd in each cell. We tackle the problem as follows. We begin by making a natural and obvious extension of subresultant polynomials of two polynomials to several polynomials. Then we develop the following structural theories about them. 1. We generalize Sylvester's theory to several polynomials, in order to obtain an elegant relationship between generalized subresultant polynomials and the gcd of several polynomials, yielding an elegant algorithm. 2. We generalize Habicht's theory to several polynomials, in order to obtain a systematic relationship between generalized subresultant polynomials and pseudo-remainders, yielding an efficient algorithm. Using the generalized theories, we present a simple (structurally elegant) algorithm which is significantly more efficient (both in the output size and computing time) than algorithms based on previous approaches.
