Buchberger's algorithm
In the theory of multivariate polynomials, Buchberger's algorithm is a method for transforming a given set of polynomials into a Gröbner basis, which is another set of polynomials that have the same common zeros and are more convenient for extracting information on these common zeros. It was introduced by Bruno Buchberger simultaneously with the definition of Gröbner bases.
Euclidean algorithm for polynomial Greatest common divisor computation and Gaussian elimination of linear systems are special cases of Buchberger's algorithm when the number of variables or the degrees of the polynomials are respectively equal to one.
For other Gröbner basis algorithms, see Gröbner basis § Algorithms and implementations.
Algorithm
A crude version of this algorithm to find a basis for an ideal I of a polynomial ring R proceeds as follows:
- Input A set of polynomials F that generates I
- Output A Gröbner basis G for I
- G := F
- For every fi, fj in G, denote by gi the leading term of fi with respect to the given ordering, and by aij the least common multiple of gi and gj.
- Choose two polynomials in G and let Sij = aij/ gi fi − aij/ gj fj (Note that the leading terms here will cancel by construction).
- Reduce Sij, with the multivariate division algorithm relative to the set G until the result is not further reducible. If the result is non-zero, add it to G.
- Repeat steps 2-4 until all possible pairs are considered, including those involving the new polynomials added in step 4.
- Output G
The polynomial Sij is commonly referred to as the S-polynomial, where S refers to subtraction (Buchberger) or Syzygy (others). The pair of polynomials with which it is associated is commonly referred to as critical pair.
There are numerous ways to improve this algorithm beyond what has been stated above. For example, one could reduce all the new elements of F relative to each other before adding them. If the leading terms of fi and fj share no variables in common, then Sij will always reduce to 0 (if we use only fi and fj for reduction), so we needn't calculate it at all.
The algorithm terminates because it is consistently increasing the size of the monomial ideal generated by the leading terms of our set F, and Dickson's lemma (or the Hilbert basis theorem) guarantees that any such ascending chain must eventually become constant.
Complexity
The computational complexity of Buchberger's algorithm is very difficult to estimate, because of the number of choices that may dramatically change the computation time. Nevertheless, T. W. Dubé has proved[1] that the degrees of the elements of a reduced Gröbner basis are always bounded by
- ,
where n is the number of variables, and d the maximal total degree of the input polynomials. This allows, in theory, to use linear algebra over the vector space of the polynomials of degree bounded by this value, for getting an algorithm of complexity .
On the other hand, there are examples[2] where the Gröbner basis contains elements of degree
- ,
and the above upper bound of complexity is optimal. Nevertheless, such examples are extremely rare.
Since its discovery, many variants of Buchberger's have been introduced to improve its efficiency. Faugère's F4 and F5 algorithms are presently the most efficient algorithms for computing Gröbner bases, and allow to compute routinely Gröbner bases consisting of several hundreds of polynomials, having each several hundreds of terms and coefficients of several hundreds of digits.
See also
- Knuth–Bendix completion algorithm
- Quine–McCluskey algorithm – analogous algorithm for Boolean algebra
References
- ^ Dubé, Thomas W. (1990). "The Structure of Polynomial Ideals and Gröbner Bases". SIAM Journal on Computing. 19 (4): 750–773. doi:10.1137/0219053.
- ^ Mayr, Ernst W; Meyer, Albert R (1982). "The complexity of the word problems for commutative semigroups and polynomial ideals". Advances in Mathematics. 46 (3): 305. doi:10.1016/0001-8708(82)90048-2.
Further reading
- Buchberger, B. (August 1976). "Theoretical Basis for the Reduction of Polynomials to Canonical Forms". ACM SIGSAM Bulletin. ACM. 10 (3): 19–29. doi:10.1145/1088216.1088219. MR 0463136. S2CID 15179417.
- David Cox, John Little, and Donald O'Shea (1997). Ideals, Varieties, and Algorithms: An Introduction to Computational Algebraic Geometry and Commutative Algebra, Springer. ISBN 0-387-94680-2.
- Vladimir P. Gerdt, Yuri A. Blinkov (1998). Involutive Bases of Polynomial Ideals, Mathematics and Computers in Simulation, 45:519ff
External links
- "Buchberger algorithm", Encyclopedia of Mathematics, EMS Press, 2001 [1994]
- Buchberger's algorithm on Scholarpedia
- Weisstein, Eric W. "Buchberger's Algorithm". MathWorld.