pdf
A Gröbner basis is a set of multivariate polynomials that has desirable algorithmic properties. Every set of polynomials can be transformed into a Gröbner basis. This process generalizes three familiar techniques: Gaussian elimination for solving linear systems of equations, the Euclidean algorithm for computing the greatest common divisor of two univariate polynomials, and the Simplex Algorithm for linear programming; see [3]. For example, the input for Gaussian elimination is a collection of linear forms such as and the algorithm transforms into the Gröbner basis .
Let be any field, such as the real numbers , the complex numbers , the rational numbers , or a finite field . We write for the ring of polynomials in variables with coefficients in the field . If is any set of polynomials, then the ideal generated by is the set consisting of all polynomial linear combinations:In our example the set and its Gröbner basis generate the same ideal: . By Hilbert’s Basis Theorem, every ideal in has the form ; i.e., it is generated by some finite set of polynomials.
A term order on is a total order ≺ on the set of all monomials which has the following two properties:
(1) It is multiplicative; i.e., implies for all .
(2) The constant monomial is the smallest; i.e., for all .
An example of a term order (for ) is the degree lexicographic order
If we fix a term order ≺, then every polynomial has a unique initial term . This is the ≺-largest monomial which occurs with nonzero coefficient in the expansion of . We write the terms of in ≺-decreasing order, and we often underline the initial term. For instance, a quadratic polynomial is written . Suppose now that is an ideal in . Then its initial ideal is the ideal generated by the initial terms of all the polynomials in :
A finite subset of is a Gröbner basis with respect to the term order ≺ if the initial terms of the elements in suffice to generate the initial ideal:
There is no minimality requirement for being a Gröbner basis. If is a Gröbner basis for , then any finite subset of that contains is also a Gröbner basis. To remedy this nonminimality, we say that is a reduced Gröbner basis if
(1) for each , the coefficient of in is 1,
(2) the set minimally generates , and
(3) no trailing term of any lies in .
With this definition, we have the following theorem:
The reduced Gröbner basis can be computed from any generating set of by a method that was introduced in Bruno Buchberger’s 1965 dissertation.
Buchberger named his method after his advisor, Wolfgang Gröbner. With hindsight, the idea of Gröbner bases can be traced back to earlier sources, including a paper written in 1900 by the invariant theorist Paul Gordan. But Buchberger was the first to give an algorithm for computing Gröbner bases.
Gröbner bases are very useful for solving systems of polynomial equations. Suppose , and let be a finite set of polynomials in .
The variety of is the set of all common complex zeros:The variety does not change if we replace by another set of polynomials that generates the same ideal in . In particular, the reduced Gröbner basis for the ideal specifies the same variety:The advantage of is that it reveals geometric properties of the variety that are not visible from . The first question that one might ask about a variety is whether it is empty. Hilbert’s Nullstellensatz implies that
The variety is finite if and only if the set of standard monomials is finite, and the number of standard monomials equals the cardinality of , when zeros are counted with multiplicity.
For this is the Fundamental Theorem of Algebra, which states that the variety of a univariate polynomial of degree consists of complex numbers. Here the singleton is a Gröbner basis, and the standard monomials are .
Our criterion for deciding whether a variety is finite generalizes to the following formula for the dimension of a variety. Consider a subset of the variables such that no monomial in the variables in appears in , and suppose that has maximal cardinality among all subsets with this property. That maximal cardinality equals the dimension of .
The set of standard monomials is a -vectorspace basis for the residue ring .
The image of a polynomial modulo can be expressed uniquely as a -linear combination of standard monomials. This expression is the normal form of . The process of computing the normal form is the division algorithm. In the familar case of only one variable , where and has degree , the division algorithm writes any polynomial as a -linear combination of . But the division algorithm works relative to any Gröbner basis in any number of variables.
How can we test whether a given set of polynomials is a Gröbner basis or not? Consider any two polynomials and in , and form their S-polynomial . Here and are monomials of smallest possible degree such that . The S-polynomial lies in the ideal . We apply the division algorithm with respect to the tentative Gröbner basis to .
The resulting normal form is a -linear combination of monomials none of which is divisible by an initial monomial from . A necessary condition for to be a Gröbner basis is
Buchberger’s Criterion states that this necessary condition is sufficient: a set of polynomials is a Gröbner basis if and only if all its S-polynomials have normal form zero. From this criterion, one derives Buchberger’s Algorithm [1] for computing the reduced Gröbner basis from any given input set .
In summary, Gröbner bases and the Buchberger Algorithm for finding them are fundamental notions in algebra. They furnish the engine for more advanced computations in algebraic geometry, such as elimination theory, computing cohomology, resolving singularities, etc. Given that polynomial models are ubiquitous across the sciences and engineering, Gröbner bases have been used by researchers in optimization, coding, robotics, control theory, statistics, molecular biology, and many other fields. We invite the reader to experiment with one of the many implementations of Buchberger’s algorithm (e.g., in CoCoA, Macaulay2, Magma, Maple, Mathematica, or Singular).
References
[1] David Cox, John Little, and Donal O’ Shea, Ideals, Varieties and Algorithms. An Introduction to Computational Algebraic Geometry and Commutative Algebra, second ed., Undergraduate Texts in Mathematics, Springer-Verlag, New York, 1997.
[2] Niels Lauritzen, Concrete Abstract Algebra: From Numbers to Gröbner Bases, Cambridge University Press, 2003.
[3] Bernd Sturmfels, Two Lectures on Gröbner Bases, New Horizons in Undergraduate Mathematics, VMath Lecture Series, Mathematical Sciences Research Institute, Berkeley, California, 2005, http://www.msri.org/communications/vmath/special_productions/.
A Gröbner basis is a set of multivariate polynomials that has desirable algorithmic properties. Every set of polynomials can be transformed into a Gröbner basis. This process generalizes three familiar techniques: Gaussian elimination for solving linear systems of equations, the Euclidean algorithm for computing the greatest common divisor of two univariate polynomials, and the Simplex Algorithm for linear programming; see [3]. For example, the input for Gaussian elimination is a collection of linear forms such as and the algorithm transforms into the Gröbner basis .
Let be any field, such as the real numbers , the complex numbers , the rational numbers , or a finite field . We write for the ring of polynomials in variables with coefficients in the field . If is any set of polynomials, then the ideal generated by is the set consisting of all polynomial linear combinations:In our example the set and its Gröbner basis generate the same ideal: . By Hilbert’s Basis Theorem, every ideal in has the form ; i.e., it is generated by some finite set of polynomials.
A term order on is a total order ≺ on the set of all monomials which has the following two properties:
(1) It is multiplicative; i.e., implies for all .
(2) The constant monomial is the smallest; i.e., for all .
An example of a term order (for ) is the degree lexicographic order
If we fix a term order ≺, then every polynomial has a unique initial term . This is the ≺-largest monomial which occurs with nonzero coefficient in the expansion of . We write the terms of in ≺-decreasing order, and we often underline the initial term. For instance, a quadratic polynomial is written . Suppose now that is an ideal in . Then its initial ideal is the ideal generated by the initial terms of all the polynomials in :
A finite subset of is a Gröbner basis with respect to the term order ≺ if the initial terms of the elements in suffice to generate the initial ideal:
There is no minimality requirement for being a Gröbner basis. If is a Gröbner basis for , then any finite subset of that contains is also a Gröbner basis. To remedy this nonminimality, we say that is a reduced Gröbner basis if
(1) for each , the coefficient of in is 1,
(2) the set minimally generates , and
(3) no trailing term of any lies in .
With this definition, we have the following theorem:
If the term order ≺ is fixed, then every ideal in has a unique reduced Gröbner basis.
The reduced Gröbner basis can be computed from any generating set of by a method that was introduced in Bruno Buchberger’s 1965 dissertation.
Buchberger named his method after his advisor, Wolfgang Gröbner. With hindsight, the idea of Gröbner bases can be traced back to earlier sources, including a paper written in 1900 by the invariant theorist Paul Gordan. But Buchberger was the first to give an algorithm for computing Gröbner bases.
Gröbner bases are very useful for solving systems of polynomial equations. Suppose , and let be a finite set of polynomials in .
The variety of is the set of all common complex zeros:The variety does not change if we replace by another set of polynomials that generates the same ideal in . In particular, the reduced Gröbner basis for the ideal specifies the same variety:The advantage of is that it reveals geometric properties of the variety that are not visible from . The first question that one might ask about a variety is whether it is empty. Hilbert’s Nullstellensatz implies that
the variety is empty if and only if equals {1}.
How can one count the number of zeros of a given system of equations? To answer this, we need one more definition. Given a fixed ideal in and a term order ≺, a monomial is called standard if it is not in the initial ideal . The number of standard monomials is finite if and only if every variable appears to some power in the initial ideal. For example, if , then there are sixty standard monomials, but if , then the set of standard monomials is infinite.The variety is finite if and only if the set of standard monomials is finite, and the number of standard monomials equals the cardinality of , when zeros are counted with multiplicity.
For this is the Fundamental Theorem of Algebra, which states that the variety of a univariate polynomial of degree consists of complex numbers. Here the singleton is a Gröbner basis, and the standard monomials are .
Our criterion for deciding whether a variety is finite generalizes to the following formula for the dimension of a variety. Consider a subset of the variables such that no monomial in the variables in appears in , and suppose that has maximal cardinality among all subsets with this property. That maximal cardinality equals the dimension of .
The set of standard monomials is a -vectorspace basis for the residue ring .
The image of a polynomial modulo can be expressed uniquely as a -linear combination of standard monomials. This expression is the normal form of . The process of computing the normal form is the division algorithm. In the familar case of only one variable , where and has degree , the division algorithm writes any polynomial as a -linear combination of . But the division algorithm works relative to any Gröbner basis in any number of variables.
How can we test whether a given set of polynomials is a Gröbner basis or not? Consider any two polynomials and in , and form their S-polynomial . Here and are monomials of smallest possible degree such that . The S-polynomial lies in the ideal . We apply the division algorithm with respect to the tentative Gröbner basis to .
The resulting normal form is a -linear combination of monomials none of which is divisible by an initial monomial from . A necessary condition for to be a Gröbner basis is
for all .
Buchberger’s Criterion states that this necessary condition is sufficient: a set of polynomials is a Gröbner basis if and only if all its S-polynomials have normal form zero. From this criterion, one derives Buchberger’s Algorithm [1] for computing the reduced Gröbner basis from any given input set .
In summary, Gröbner bases and the Buchberger Algorithm for finding them are fundamental notions in algebra. They furnish the engine for more advanced computations in algebraic geometry, such as elimination theory, computing cohomology, resolving singularities, etc. Given that polynomial models are ubiquitous across the sciences and engineering, Gröbner bases have been used by researchers in optimization, coding, robotics, control theory, statistics, molecular biology, and many other fields. We invite the reader to experiment with one of the many implementations of Buchberger’s algorithm (e.g., in CoCoA, Macaulay2, Magma, Maple, Mathematica, or Singular).
References
[1] David Cox, John Little, and Donal O’ Shea, Ideals, Varieties and Algorithms. An Introduction to Computational Algebraic Geometry and Commutative Algebra, second ed., Undergraduate Texts in Mathematics, Springer-Verlag, New York, 1997.
[2] Niels Lauritzen, Concrete Abstract Algebra: From Numbers to Gröbner Bases, Cambridge University Press, 2003.
[3] Bernd Sturmfels, Two Lectures on Gröbner Bases, New Horizons in Undergraduate Mathematics, VMath Lecture Series, Mathematical Sciences Research Institute, Berkeley, California, 2005, http://www.msri.org/communications/vmath/special_productions/.
PREVIOUSHomogeneous Coordinates
NEXTAbel's Test