Buchberger's algorithm: Difference between revisions

Revision as of 18:37, 13 December 2022

Latest revision as of 14:35, 1 June 2025

(8 intermediate revisions by 7 users not shown)

Line 1:

{{Short description|Algorithm for computing Gröbner bases}}

In the theory of [[multivariate polynomial]]s, '''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]] computationis anda special case of Buchberger's algorithm restricted to polynomials of a single variable. [[Gaussian elimination]] of a [[system of linear equations|linear system]]s areis another special casescase of Buchberger's algorithm whenwhere the numberdegree of variables or the degrees of theall polynomials are respectively equal toequals one.

For other Gröbner basis algorithms, see {{slink|Gröbner basis#Algorithms and implementations}}.

Line 15:

:# Choose two polynomials in ''G'' and let {{math|1=''S''_''ij'' = {{sfrac|''a''_''ij'' | ''g''_{{var|i}}}} ''f''_{{var|i}} − {{sfrac|''a''_''ij'' | ''g''_''j''}} ''f''_''j''}} ''(Note that the leading terms here will cancel by construction)''.

:# Reduce ''S''_''ij'', 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-42–4 until all possible pairs are considered, including those involving the new polynomials added in step 4.

:# Output ''G''

Line 30:

<math>d^{2^{n+o(1)}}</math>.

On the other hand, there are examples<ref>{{cite journal|doi=10.1016/0001-8708(82)90048-2|doi-access=free|title=The complexity of the word problems for commutative semigroups and polynomial ideals|journal=[[Advances in Mathematics]]|volume=46|issue=3|pages=305305–329|year=1982|last1=Mayr|first1=Ernst W|last2=Meyer|first2=Albert R|hdl=1721.1/149010|hdl-access=free}}</ref> where the Gröbner basis contains elements of degree

:<math>d^{2^{\Omega(n)}}</math>,

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 ==

In the [[SymPy|SymPy library]] for [[Python (programming language)|Python]], the (improved) Buchberger algorithm is implemented as <code>sympy.polys.polytools.groebner()</code>.<ref>{{cite web |title=Polynomials Manipulation Module Reference - SymPy 1.14.0 documentation |url=https://docs.sympy.org/latest/modules/polys/reference.html#sympy.polys.polytools.groebner |website=docs.sympy.org}}</ref>

There is an implementation of Buchberger’s algorithm that has been proved correct

within the proof assistant [[Coq (proof assistant)|Coq]].<ref>{{cite journal |last1=Théry |first1=Laurent |title=A Machine-Checked Implementation of Buchberger's Algorithm |journal=Journal of Automated Reasoning |date=2001 |volume=26 |issue=2 |pages=107–137 |doi=10.1023/A:1026518331905}}</ref>

== See also ==