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.
The Euclidean algorithm for computing the polynomial greatest common divisor is a special case of Buchberger's algorithm restricted to polynomials of a single variable. Gaussian elimination of a system of linear equations is another special case where the degree of all polynomials equals one.
For other Gröbner basis algorithms, see.

Algorithm

A crude version of this algorithm to find a basis for an ideal of a polynomial ring R proceeds as follows:
The polynomial S_ij is commonly referred to as the S-polynomial, where S refers to subtraction or syzygy. The pair of polynomials with which it is associated is commonly referred to as critical pair.
There are many ways to improve this algorithm beyond what has been stated above. For example, all the new elements of F can be reduced relative to each other before adding them. If the leading terms of f_i and f_j share no variables in common, then S_ij will always reduce to 0, so it needn't be calculated 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 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 that the degrees of the elements of a reduced Gröbner basis are always bounded by
where is the number of variables, and 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 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.

Implementations

At least one implementation of Buchberger’s algorithm has been proved correct within the proof assistant Coq.
In the SymPy library for Python, the Buchberger algorithm is implemented as sympy.polys.polytools.groebner.