Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Gröbner basis

From Wikipedia, the free encyclopedia
Mathematical construct in computer algebra

Inmathematics, and more specifically incomputer algebra, computationalalgebraic geometry, and computationalcommutative algebra, aGröbner basis is a particular kind ofgenerating set of an ideal in apolynomial ringK[x1,,xn]{\displaystyle K[x_{1},\ldots ,x_{n}]} over afieldK{\displaystyle K}. A Gröbner basis allows many important properties of the ideal and the associatedalgebraic variety to be deduced easily, such as thedimension and the number of zeros when it is finite. Gröbner basis computation is one of the main practical tools for solvingsystems of polynomial equations and computing the images ofalgebraic varieties underprojections orrational maps.

Gröbner basis computation can be seen as a multivariate, non-linear generalization of bothEuclid's algorithm for computingpolynomial greatest common divisors, andGaussian elimination for linear systems.[1]

Gröbner bases were introduced byBruno Buchberger in his 1965 Ph.D. thesis, which also included an algorithm to compute them (Buchberger's algorithm). He named them after his advisorWolfgang Gröbner. In 2007, Buchberger received theAssociation for Computing Machinery'sParis Kanellakis Theory and Practice Award for this work.However, the Russian mathematicianNikolai Günther had introduced a similar notion in 1913, published in various Russian mathematical journals. These papers were largely ignored by the mathematical community until their rediscovery in 1987 by Bodo Renschuchet al.[2] An analogous concept for multivariatepower series was developed independently byHeisuke Hironaka in 1964, who named themstandard bases. This term has been used by some authors to also denote Gröbner bases.

The theory of Gröbner bases has been extended by many authors in various directions. It has been generalized to other structures such as polynomials overprincipal ideal rings orpolynomial rings, and also some classes of non-commutative rings and algebras, likeOre algebras.

Tools

[edit]

Polynomial ring

[edit]
Main article:Polynomial ring

Gröbner bases are primarily defined forideals in apolynomial ringR=K[x1,,xn]{\displaystyle R=K[x_{1},\ldots ,x_{n}]} over afieldK. Although the theory works for any field, most Gröbner basis computations are done either whenK is thefield of rationals or the integers modulo a prime number.

In the context of Gröbner bases, a nonzeropolynomial inR=K[x1,,xn]{\displaystyle R=K[x_{1},\ldots ,x_{n}]} is commonly represented as a sumc1M1++cmMm,{\displaystyle c_{1}M_{1}+\cdots +c_{m}M_{m},} where theci{\displaystyle c_{i}} are nonzero elements ofK, calledcoefficients, and theMi{\displaystyle M_{i}} aremonomials (calledpower products by Buchberger and some of his followers) of the formx1a1xnan,{\displaystyle x_{1}^{a_{1}}\cdots x_{n}^{a_{n}},} where theai{\displaystyle a_{i}} are nonnegative integers. The vectorA=[a1,,an]{\displaystyle A=[a_{1},\ldots ,a_{n}]} is called theexponent vector of the monomial. When the listX=[x1,,xn]{\displaystyle X=[x_{1},\ldots ,x_{n}]} of the variables is fixed, the notation of monomials is often abbreviated asx1a1xnan=XA.{\displaystyle x_{1}^{a_{1}}\cdots x_{n}^{a_{n}}=X^{A}.}

Monomials are uniquely defined by their exponent vectors, and, when amonomial ordering (see below) is fixed, a polynomial is uniquely represented by theordered list of theordered pairs formed by an exponent vector and the corresponding coefficient. This representation of polynomials is especially efficient for Gröbner basis computation in computers, although it is less convenient for other computations such aspolynomial factorization andpolynomial greatest common divisor.

IfF={f1,,fk}{\displaystyle F=\{f_{1},\ldots ,f_{k}\}} is a finite set of polynomials in the polynomial ringR, theideal generated byF is the set of linear combinations of elements ofF with coefficients inR; that is the set of polynomials that can be writteni=1kgifi{\textstyle \sum _{i=1}^{k}g_{i}f_{i}} withg1,,gkR.{\displaystyle g_{1},\ldots ,g_{k}\in R.}

Monomial ordering

[edit]
Main article:Monomial order

All operations related to Gröbner bases require the choice of atotal order on the monomials, with the following properties of compatibility with multiplication. For all monomialsM,N,P,

  1. MNMPNP{\displaystyle M\leq N\Longleftrightarrow MP\leq NP}
  2. MMP{\displaystyle M\leq MP}.

A total order satisfying these condition is sometimes called anadmissible ordering.

These conditions imply that the order is awell-order, that is, every strictly decreasing sequence of monomials is finite.

Although Gröbner basis theory does not depend on a particular choice of an admissible monomial ordering, three monomial orderings are especially important for the applications:

  • Lexicographical ordering, commonly calledlex orplex (for pure lexical ordering).
  • Total degree reverse lexicographical ordering, commonly calleddegrevlex.
  • Elimination ordering,lexdeg.

Gröbner basis theory was initially introduced for the lexicographical ordering. It was soon realised that the Gröbner basis for degrevlex is almost always much easier to compute, and that it is almost always easier to compute a lex Gröbner basis by first computing the degrevlex basis and then using a "change of ordering algorithm". Whenelimination is needed, degrevlex is not convenient; both lex and lexdeg may be used but, again, many computations are relatively easy with lexdeg and almost impossible with lex.

Basic operations

[edit]

Leading term, coefficient and monomial

[edit]

Once a monomial ordering is fixed, the terms of a polynomial (product of a monomial with its nonzero coefficient) are naturally ordered by decreasing monomials (for this order). This makes the representation of a polynomial as a sorted list of pairs coefficient–exponent vector acanonical representation of the polynomials (that is, two polynomials are equal if and only if they have the same representation).

The first (greatest) term of a polynomialp for this ordering and the corresponding monomial and coefficient are respectively called theleading term,leading monomial andleading coefficient and denoted, in this article,lt(p), lm(p) andlc(p).

Most polynomial operations related to Gröbner bases involve the leading terms. So, the representation of polynomials as sorted lists make these operations particularly efficient (reading the first element of a list takes a constant time, independently of the length of the list).

Polynomial operations

[edit]

The other polynomial operations involved in Gröbner basis computations are also compatible with the monomial ordering; that is, they can be performed without reordering the result:

  • The addition of two polynomials consists in amerge of the two corresponding lists of terms, with a special treatment in the case of a conflict (that is, when the same monomial appears in the two polynomials).
  • The multiplication of a polynomial by a scalar consists of multiplying each coefficient by this scalar, without any other change in the representation.
  • The multiplication of a polynomial by a monomialm consists of multiplying each monomial of the polynomial bym. This does not change the term ordering by definition of a monomial ordering.

Divisibility of monomials

[edit]

LetM=x1a1xnan{\displaystyle M=x_{1}^{a_{1}}\cdots x_{n}^{a_{n}}} andN=x1b1xnbn{\displaystyle N=x_{1}^{b_{1}}\cdots x_{n}^{b_{n}}} be two monomials, with exponent vectorsA=[a1,,an]{\displaystyle A=[a_{1},\ldots ,a_{n}]} andB=[b1,,bn].{\displaystyle B=[b_{1},\ldots ,b_{n}].}

One says thatMdividesN, or thatN is amultiple ofM, ifaibi{\displaystyle a_{i}\leq b_{i}} for everyi; that is, ifA is componentwise not greater thanB. In this case, the quotientNM{\textstyle {\frac {N}{M}}} is defined asNM=x1b1a1xnbnan.{\textstyle {\frac {N}{M}}=x_{1}^{b_{1}-a_{1}}\cdots x_{n}^{b_{n}-a_{n}}.} In other words, the exponent vector ofNM{\textstyle {\frac {N}{M}}} is the componentwise subtraction of the exponent vectors ofN andM.

Thegreatest common divisorgcd(M,N) ofM andN is the monomialx1min(a1,b1)xnmin(an,bn){\textstyle x_{1}^{\min(a_{1},b_{1})}\cdots x_{n}^{\min(a_{n},b_{n})}} whose exponent vector is the componentwise minimum ofA andB. Theleast common multiplelcm(M,N) is defined similarly withmax instead ofmin.

One has

lcm(M,N)=MNgcd(M,N).{\displaystyle \operatorname {lcm} (M,N)={\frac {MN}{\gcd(M,N)}}.}

Reduction

[edit]

Thereduction of a polynomial by other polynomials with respect to a monomial ordering is central to Gröbner basis theory. It is a generalization of bothrow reduction occurring inGaussian elimination and division steps of theEuclidean division of univariate polynomials.[1] When completed as much as possible, it is sometimes calledmultivariate division although its result is not uniquely defined.

Lead-reduction is a special case of reduction that is easier to compute. It is fundamental for Gröbner basis computation, since general reduction is needed only at the end of a Gröbner basis computation, for getting a reduced Gröbner basis from a non-reduced one.

Let an admissible monomial ordering be fixed, to which refers every monomial comparison that will occur in this section.

A polynomialf islead-reducible by another polynomialg if the leading monomiallm(f) is a multiple oflm(g). The polynomialf isreducible byg if some monomial off is a multiple oflm(g). (So, iff is lead-reducible byg, it is also reducible, butf may be reducible without being lead-reducible.)

Suppose thatf isreducible byg, and letcm be a term off such that the monomialm is a multiple oflm(g). Aone-step reduction off byg consists of replacingf by

red1(f,g)=fclc(g)mlm(g)g.{\displaystyle \operatorname {red} _{1}(f,g)=f-{\frac {c}{\operatorname {lc} (g)}}\,{\frac {m}{\operatorname {lm} (g)}}\,g.}

This operation removes the monomialm fromf without changing the terms with a monomial greater thanm (for the monomial ordering). In particular, aone step lead-reduction off produces a polynomial all of whose monomials are smaller thanlm(f).

Given a finite setG of polynomials, one says thatf isreducible orlead-reducible byG if it is reducible or lead-reducible, respectively, by at least one elementg ofG. In this case, a one-step reduction (resp. one-step lead-reduction) off byG is any one-step reduction (resp. one-step lead-reduction) off by an element ofG.

The (complete) reduction (resp. lead-reduction) off byG consists of iterating one-step reductions (respect. one-step lead reductions) until getting a polynomial that is irreducible (resp. lead-irreducible) byG. It is sometimes called anormal form off byG. In general this form is not uniquely defined because there are, in general, several elements ofG that can be used for reducingf; this non-uniqueness is the starting point of Gröbner basis theory.

The definition of the reduction shows immediately that, ifh is a normal form off byG, one has

f=h+gGqgg,{\displaystyle f=h+\sum _{g\in G}q_{g}\,g,}

whereh is irreducible byG and theqg{\displaystyle q_{g}} are polynomials such thatlm(qgg)lm(f).{\displaystyle \operatorname {lm} (q_{g}\,g)\leq \operatorname {lm} (f).} In the case of univariate polynomials, ifG consists of a single elementg, thenh is the remainder of theEuclidean division off byg, andqg is the quotient. Moreover, the division algorithm is exactly the process of lead-reduction. For this reason, some authors use the termmultivariate division instead of reduction.

Non uniqueness of reduction

[edit]

In the example that follows, there are exactly two complete lead-reductions that produce two very different results. The fact that the results are irreducible (not only lead-irreducible) is specific to the example, although this is rather common with such small examples.

In this two variable example, the monomial ordering that is used is thelexicographic order withx>y,{\displaystyle x>y,} and we consider the reduction off=2x3x2y+y3+3y{\displaystyle f=2x^{3}-x^{2}y+y^{3}+3y}, byG={g1,g2},{\displaystyle G=\{g_{1},g_{2}\},} withg1=x2+y21,g2=xy2.{\displaystyle {\begin{aligned}g_{1}&=x^{2}+y^{2}-1,\\g_{2}&=xy-2.\end{aligned}}}

For the first reduction step, either the first or the second term off may be reduced. However, the reduction of a term amounts to removing this term at the cost of adding new lower terms; if it is not the first reducible term that is reduced, it may occur that a further reduction adds a similar term, which must be reduced again. It is therefore always better to reduce first the largest (for the monomial order) reducible term; that is, in particular, to lead-reduce first until getting a lead-irreducible polynomial.

The leading term2x3{\displaystyle 2x^{3}} off is reducible byg1{\displaystyle g_{1}} and not byg2.{\displaystyle g_{2}.} So the first reduction step consists of multiplyingg1{\displaystyle g_{1}} by−2x and adding the result tof:f2xg1f1=f2xg1=x2y2xy2+2x+y3+3y.{\displaystyle f\;\xrightarrow {\overset {}{-2xg_{1}}} \;f_{1}=f-2xg_{1}=-x^{2}y-2xy^{2}+2x+y^{3}+3y.}

The leading termx2y{\displaystyle -x^{2}y} off1{\displaystyle f_{1}} is a multiple of the leading monomials of bothg1{\displaystyle g_{1}} andg2,{\displaystyle g_{2},} So, one has two choices for the second reduction step. If one choosesg2,{\displaystyle g_{2},} one gets a polynomial that can be reduced again byg2:{\displaystyle g_{2}\colon }f2xg1f1xg22xy2+y3+3y2yg2f2=y3y.{\displaystyle f\;\xrightarrow {\overset {}{-2xg_{1}}} \;f_{1}\;\xrightarrow {xg_{2}} \;-2xy^{2}+y^{3}+3y\;\xrightarrow {2yg_{2}} \;f_{2}=y^{3}-y.}No further reduction is possible, sof2{\displaystyle f_{2}} is a complete reduction off.

One gets a different result with the other choice for the second step:f2xg1f1yg12xy2+2x+2y3+2y2yg2f3=2x+2y32y.{\displaystyle f\;\xrightarrow {\overset {}{-2xg_{1}}} \;f_{1}\;\xrightarrow {yg_{1}} \;-2xy^{2}+2x+2y^{3}+2y\;\xrightarrow {2yg_{2}} \;f_{3}=2x+2y^{3}-2y.}Again, the resultf3{\displaystyle f_{3}} is irreducible, although only lead reductions were done.

In summary, the complete reduction off can result in eitherf2=y3y{\displaystyle f_{2}=y^{3}-y} orf3=2x+2y32y.{\displaystyle f_{3}=2x+2y^{3}-2y.}

It is for dealing with the problems set by this non-uniqueness thatBuchberger introduced Gröbner bases andS-polynomials. Intuitively,0=ff{\displaystyle 0=f-f} may be reduced tof2f3.{\displaystyle f_{2}-f_{3}.} This implies thatf2f3{\displaystyle f_{2}-f_{3}} belongs to the ideal generated byG. So, this ideal is not changed by addingf3f2{\displaystyle f_{3}-f_{2}} toG, and this allows more reductions. In particular,f3{\displaystyle f_{3}} can be reduced tof2{\displaystyle f_{2}} byf3f2{\displaystyle f_{3}-f_{2}} and this restores the uniqueness of the reduced form.

Here Buchberger's algorithm for Gröbner bases would begin by adding toG the polynomial

g3=yg1xg2=2x+y3y.{\displaystyle g_{3}=yg_{1}-xg_{2}=2x+y^{3}-y.}

This polynomial, calledS-polynomial by Buchberger, is the difference of the one-step reductions of the least common multiplex2y{\displaystyle x^{2}y} of the leading monomials ofg1{\displaystyle g_{1}} andg2{\displaystyle g_{2}}, byg2{\displaystyle g_{2}} andg1{\displaystyle g_{1}} respectively:

g3=(x2yx2ylt(g2)g2)(x2yx2ylt(g1)g1)=x2ylt(g1)g1x2ylt(g2)g2{\displaystyle g_{3}=\left(x^{2}y-{\frac {x^{2}y}{lt(g_{2})}}g_{2}\right)-\left(x^{2}y-{\frac {x^{2}y}{lt(g_{1})}}g_{1}\right)={\frac {x^{2}y}{lt(g_{1})}}g_{1}-{\frac {x^{2}y}{lt(g_{2})}}g_{2}}.

In this example, one hasg3=f3f2.{\displaystyle g_{3}=f_{3}-f_{2}.} This does not complete Buchberger's algorithm, asxy gives different results, when reduced byg2{\displaystyle g_{2}} org3.{\displaystyle g_{3}.}

S-polynomial

[edit]

Given monomial ordering, theS-polynomial orcritical pair of two polynomialsf andg is the polynomial

S(f,g)=red1(lcm,g)red1(lcm,f){\displaystyle S(f,g)=\operatorname {red} _{1}(\mathrm {lcm} ,g)-\operatorname {red} _{1}(\mathrm {lcm} ,f)};

wherelcm denotes the least common multiple of the leading monomials off andg.Using the definition ofred1{\displaystyle \operatorname {red} _{1}}, this translates to:

S(f,g)=(lcm1lc(g)lcmlm(g)g)(lcm1lc(f)lcmlm(f)f)=1lc(f)lcmlm(f)f1lc(g)lcmlm(g)g.{\displaystyle {\begin{aligned}S(f,g)&=\left(\mathrm {lcm} -{\frac {1}{\operatorname {lc} (g)}}\,{\frac {\mathrm {lcm} }{\operatorname {lm} (g)}}\,g\right)-\left(\mathrm {lcm} -{\frac {1}{\operatorname {lc} (f)}}\,{\frac {\mathrm {lcm} }{\operatorname {lm} (f)}}\,f\right)\\&={\frac {1}{\operatorname {lc} (f)}}\,{\frac {\mathrm {lcm} }{\operatorname {lm} (f)}}\,f-{\frac {1}{\operatorname {lc} (g)}}\,{\frac {\mathrm {lcm} }{\operatorname {lm} (g)}}\,g\\\end{aligned}}.}

Using the property that relates thelcm and thegcd, theS-polynomial can also be written as:

S(f,g)=1lc(f)lm(g)gcdf1lc(g)lm(f)gcdg;{\displaystyle S(f,g)={\frac {1}{\operatorname {lc} (f)}}\,{\frac {\operatorname {lm} (g)}{\mathrm {gcd} }}\,f-{\frac {1}{\operatorname {lc} (g)}}\,{\frac {\operatorname {lm} (f)}{\mathrm {gcd} }}\,g;}

wheregcd denotes the greatest common divisor of the leading monomials off andg.

As the monomials that are reducible by bothf andg are exactly the multiples oflcm, one can deal with all cases of non-uniqueness of the reduction by considering only theS-polynomials. This is a fundamental fact for Gröbner basis theory and all algorithms for computing them.

For avoiding fractions when dealing with polynomials with integer coefficients, theS polynomial is often defined as

S(f,g)=lc(g)lm(g)gcdflc(f)lm(f)gcdg;{\displaystyle S(f,g)=\operatorname {lc} (g)\,{\frac {\operatorname {lm} (g)}{\mathrm {gcd} }}\,f-\operatorname {lc} (f)\,{\frac {\operatorname {lm} (f)}{\mathrm {gcd} }}\,g;}

This does not change anything to the theory since the two polynomials areassociates.

Definition

[edit]

LetR=F[x1,,xn]{\displaystyle R=F[x_{1},\ldots ,x_{n}]} be apolynomial ring over a fieldF. In this section, we suppose that an admissible monomial ordering has been fixed.

LetG be a finite set of polynomials inR thatgenerates anidealI. The setG is a Gröbner basis (with respect to the monomial ordering), or, more precisely, a Gröbner basis ofI if

  1. the ideal generated by the leading monomials of the polynomials inI equals the ideal generated by the leading monomials ofG,

or, equivalently,

  1. the leading monomial of every polynomial inI is a multiple of the leading monomial of some polynomial inG.

There are many characterizing properties, which can each be taken as an equivalent definition of Gröbner bases. For conciseness, in the following list, the notation "one-word/another word" means that one can take either "one-word" or "another word" for having two different characterizations of Gröbner bases. All the following assertions are characterizations of Gröbner bases:

  1. a polynomialf is inI, if and only if some/every complete lead-reduction/reduction off byG produces the zero polynomial;
  2. for everyS-polynomials of elements ofG, some/every complete lead-reduction/reduction ofs byG produces zero;
  3. all complete reductions of an element ofR produce the same result;
  4. the monomials that are irreducible byG form a basis of theF-vector spaceR/I.{\displaystyle R/I.}

Counting the above definition, this provides 12 characterizations of Gröbner bases. The fact that so many characterizations are possible makes Gröbner bases very useful. For example, condition 3 provides an algorithm for testingideal membership; condition 4 provides an algorithm for testing whether a set of polynomials is a Gröbner basis and forms the basis ofBuchberger's algorithm for computing Gröbner bases; conditions 5 and 6 allow computing inR/I{\displaystyle R/I} in a way that is very similar tomodular arithmetic.

Existence

[edit]

For every admissible monomial ordering and every finite setG of polynomials, there is a Gröbner basis that containsG and generates the same ideal. Moreover, such a Gröbner basis may be computed withBuchberger's algorithm.

This algorithm uses condition 4, and proceeds roughly as follows: for any two elements ofG, compute the complete reduction byG of theirS-polynomial, and add the result toG if it is not zero; repeat this operation with the new elements ofG included until, eventually, all reductions produce zero.

The algorithm terminates always because ofDickson's lemma or because polynomial rings areNoetherian (Hilbert's basis theorem). Condition 4 ensures that the result is a Gröbner basis, and the definitions ofS-polynomials and reduction ensure that the generated ideal is not changed.

The above method is an algorithm for computing Gröbner bases; however, it is very inefficient. Many improvements of the original Buchberger's algorithm, and several other algorithms have been proposed and implemented, which dramatically improve the efficiency. See§ Algorithms and implementations, below.

Reduced Gröbner bases

[edit]

A Gröbner basis isminimal if all leading monomials of its elements are irreducible by the other elements of the basis. Given a Gröbner basis of an idealI, one gets a minimal Gröbner basis ofI by removing the polynomials whose leading monomials are multiple of the leading monomial of another element of the Gröbner basis. However, if two polynomials of the basis have the same leading monomial, only one must be removed. So, every Gröbner basis contains a minimal Gröbner basis as a subset.

All minimal Gröbner bases of a given ideal (for a fixed monomial ordering) have the same number of elements, and the same leading monomials, and the non-minimal Gröbner bases have more elements than the minimal ones.

A Gröbner basis isreduced if every polynomial in it is irreducible by the other elements of the basis, and has1 as leading coefficient. So, every reduced Gröbner basis is minimal, but a minimal Gröbner basis need not be reduced.

Given a Gröbner basis of an idealI, one gets a reduced Gröbner basis ofI by first removing the polynomials that are lead-reducible by other elements of the basis (for getting a minimal basis); then replacing each element of the basis by the result of the complete reduction by the other elements of the basis; and, finally, by dividing each element of the basis by its leading coefficient.

All reduced Gröbner bases of an ideal (for a fixed monomial ordering) are equal. It follows that two ideals are equal if and only if they have the same reduced Gröbner basis.

Sometimes, reduced Gröbner bases are defined without the condition on the leading coefficients. In this case, the uniqueness of reduced Gröbner bases is true onlyup to the multiplication of polynomials by a nonzero constant.

When working with polynomials over the fieldQ{\displaystyle \mathbb {Q} } of therational numbers, it is useful to work only with polynomials with integer coefficients. In this case, the condition on the leading coefficients in the definition of a reduced basis may be replaced by the condition that all elements of the basis areprimitive polynomials with integer coefficients, with positive leading coefficients. This restores the uniqueness of reduced bases.

Special cases

[edit]

For every monomial ordering, theempty set of polynomials is the unique Gröbner basis of thezero ideal.

For every monomial ordering, a set of polynomials that contains a nonzero constant is a Gröbner basis of theunit ideal (the whole polynomial ring). Conversely, every Gröbner basis of the unit ideal contains a nonzero constant. The reduced Gröbner basis of the unit is formed by the single polynomial1.

In the case of polynomials in a single variable, there is a unique admissible monomial ordering, the ordering by the degree. The minimal Gröbner bases are thesingletons consisting of a single polynomial. The reduced Gröbner bases are themonic polynomials.

Example and counterexample

[edit]
The zeroes off{\displaystyle f} form the red parabola; the zeroes ofg{\displaystyle g} form the three blue vertical lines. Their intersection consists of three points.

LetR=Q[x,y]{\displaystyle R=\mathbb {Q} [x,y]} be the ring of bivariate polynomials with rational coefficients and consider the idealI=f,g{\displaystyle I=\langle f,g\rangle } generated by the polynomials

f=x2y{\displaystyle f=x^{2}-y},
g=x3x{\displaystyle g=x^{3}-x}.

By reducingg byf, one obtains a new polynomialk such thatI=f,k:{\displaystyle I=\langle f,k\rangle :}

k=gxf=xyx.{\displaystyle k=g-xf=xy-x.}

None off andk is reducible by the other, butxk is reducible byf, which gives another polynomial inI:

h=xk(y1)f=y2y.{\displaystyle h=xk-(y-1)f=y^{2}-y.}

Under lexicographic ordering withx>y{\displaystyle x>y} we have

lt(f) =x2
lt(k) =xy
lt(h) =y2

Asf,k andh belong toI, and none of them is reducible by the others, none of{f,k},{\displaystyle \{f,k\},}{f,h},{\displaystyle \{f,h\},} and{h,k}{\displaystyle \{h,k\}} is a Gröbner basis ofI.

On the other hand,{f,k,h} is a Gröbner basis ofI, since the S-polynomials

yfxk=y(x2y)x(xyx)=fhykxh=y(xyx)x(y2y)=0y2fx2h=y(yfxk)+x(ykxh){\displaystyle {\begin{aligned}yf-xk&=y(x^{2}-y)-x(xy-x)=f-h\\yk-xh&=y(xy-x)-x(y^{2}-y)=0\\y^{2}f-x^{2}h&=y(yf-xk)+x(yk-xh)\end{aligned}}}

can be reduced to zero byf,k andh.

The method that has been used here for findingh andk, and proving that{f,k,h} is a Gröbner basis is a direct application ofBuchberger's algorithm. So, it can be applied mechanically to any similar example, although, in general, there are many polynomials and S-polynomials to consider, and the computation is generally too large for being done without a computer.

Properties and applications of Gröbner bases

[edit]

Unless explicitly stated, all the results that follow[3] are true for anymonomial ordering (see that article for the definitions of the different orders that are mentioned below).

It is a common misconception that the lexicographical order is needed for some of these results. On the contrary, the lexicographical order is, almost always, the most difficult to compute, and using it makes impractical many computations that are relatively easy with graded reverse lexicographic order (grevlex), or, when elimination is needed, the elimination order (lexdeg) which restricts to grevlex on each block of variables.

Equality of ideals

[edit]

Reduced Gröbner bases areunique for any given ideal and any monomial ordering. Thus two ideals are equal if and only if they have the same (reduced) Gröbner basis (usually a Gröbner basis software always produces reduced Gröbner bases).

Membership and inclusion of ideals

[edit]

Thereduction of a polynomialf by the Gröbner basisG of an idealI yields 0if and only iff is inI. This allows to test the membership of an element in an ideal. Another method consists in verifying that the Gröbner basis ofG∪{f} is equal toG.

To test if the idealI generated byf1, ...,fk is contained in the idealJ, it suffices to test that everyfI is inJ. One may also test the equality of the reduced Gröbner bases ofJ andJ ∪ {f1, ...,fk}.

Solutions of a system of algebraic equations

[edit]
Main article:System of polynomial equations

Any set of polynomials may be viewed as asystem of polynomial equations by equating the polynomials to zero. The set of the solutions of such a system depends only on the generated ideal, and, therefore does not change when the given generating set is replaced by the Gröbner basis, for any ordering, of the generated ideal. Such a solution, with coordinates in analgebraically closed field containing the coefficients of the polynomials, is called azero of the ideal. In the usual case ofrational coefficients, this algebraically closed field is chosen as thecomplex field.

An ideal does not have any zero (the system of equations isinconsistent) if and only if 1 belongs to the ideal (this isHilbert's Nullstellensatz), or, equivalently, if its Gröbner basis (for any monomial ordering) contains 1, or, also, if the corresponding reduced Gröbner basis is [1].

Given the Gröbner basisG of an idealI, it has only a finite number of zeros, if and only if, for each variablex,G contains a polynomial with a leading monomial that is a power ofx (without any other variable appearing in the leading term). If this is the case, then the number of zeros, counted with multiplicity, is equal to the number of monomials that are not multiples of any leading monomial ofG. This number is called thedegree of the ideal.

When the number of zeros is finite, the Gröbner basis for a lexicographical monomial ordering provides, theoretically, a solution: the first coordinate of a solution is a root of thegreatest common divisor of polynomials of the basis that depend only on the first variable. After substituting this root in the basis, the second coordinate of this solution is a root of the greatest common divisor of the resulting polynomials that depend only on the second variable, and so on. This solving process is only theoretical, because it implies GCD computation and root-finding of polynomials with approximate coefficients, which are not practicable because of numeric instability. Therefore, other methods have been developed to solve polynomial systems through Gröbner bases (seeSystem of polynomial equations for more details).

Dimension, degree and Hilbert series

[edit]

Thedimension of an idealI in a polynomial ringR is theKrull dimension of the ringR/I and is equal to thedimension of the algebraic set of the zeros ofI. It is also equal to number ofhyperplanes ingeneral position which are needed to have an intersection with the algebraic set, which is a finite number of points. Thedegree of the ideal and of its associated algebraic set is the number of points of this finite intersection, counted with multiplicity. In particular, the degree of ahypersurface is equal to the degree of its definition polynomial.

The dimension depend only on the set of the leading monomials of the Gröbner basis of the ideal for any monomial ordering. The same is true for the degree and degree-compatible monomial orderings; a monomial ordering isdegree compatible is smaller for the degree implies smaller for the monomial ordering.

The dimension is the maximal size of a subsetS of the variables such that there is no leading monomial depending only on the variables inS. Thus, if the ideal has dimension 0, then for each variablex there is a leading monomial in the Gröbner basis that is a power ofx.

Both dimension anddegree may be deduced from theHilbert series of the ideal, which is the seriesi=0diti{\textstyle \sum _{i=0}^{\infty }d_{i}t^{i}}, wheredi{\displaystyle d_{i}} is the number of monomials of degreei that are not multiple of any leading monomial in the Gröbner basis.[4] The Hilbert series may be summed into a rational fraction

i=0diti=P(t)(1t)d,{\displaystyle \sum _{i=0}^{\infty }d_{i}t^{i}={\frac {P(t)}{(1-t)^{d}}},}

whered is the dimension of the ideal andP(t){\displaystyle P(t)} is a polynomial. The numberP(1){\displaystyle P(1)} is thedegree of the algebraic set defined by the ideal, in the case of ahomogeneous ideal or a monomial ordering compatible with the degree; that is, to compare two monomials, one compares their total degrees first.

The dimension does not depend on the choice of a monomial ordering, although the Hilbert series and the polynomialP(t){\displaystyle P(t)} may change with changes of the monomial ordering. However, for homogeneous ideals or monomial orderings compatible with the degree, the Hilbert series and the polynomialP(t){\displaystyle P(t)} do not depend on the choice of monomial ordering.[5]

Mostcomputer algebra systems that provide functions to compute Gröbner bases provide also functions for computing the Hilbert series, and thus also the dimension and the degree.

Elimination

[edit]

The computation of Gröbner bases for anelimination monomial ordering allows computationalelimination theory. This is based on the following theorem.

Consider a polynomial ringK[x1,,xn,y1,,ym]=K[X,Y],{\displaystyle K[x_{1},\ldots ,x_{n},y_{1},\ldots ,y_{m}]=K[X,Y],} in which the variables are split into two subsetsX andY. Let us also choose an elimination monomial ordering "eliminating"X, that is a monomial ordering for which two monomials are compared by comparing first theX-parts, and, in case of equality only, considering theY-parts. This implies that a monomial containing anX-variable is greater than every monomial independent ofX.IfG is a Gröbner basis of an idealI for this monomial ordering, thenGK[Y]{\displaystyle G\cap K[Y]} is a Gröbner basis ofIK[Y]{\displaystyle I\cap K[Y]} (this ideal is often called theelimination ideal). Moreover,GK[Y]{\displaystyle G\cap K[Y]} consists exactly of the polynomials ofG whose leading terms belong toK[Y] (this makes the computation ofGK[Y]{\displaystyle G\cap K[Y]} very easy, as only the leading monomials need to be checked).

Thiselimination property has many applications, some described in the next sections.

Another application, inalgebraic geometry, is thatelimination realizes the geometric operation ofprojection of anaffine algebraic set into a subspace of the ambient space: with above notation, the (Zariski closure of) the projection of the algebraic set defined by the idealI into theY-subspace is defined by the idealIK[Y].{\displaystyle I\cap K[Y].}

The lexicographical ordering such thatx1>>xn{\displaystyle x_{1}>\cdots >x_{n}} is an elimination ordering for every partition{x1,,xk},{xk+1,,xn}.{\displaystyle \{x_{1},\ldots ,x_{k}\},\{x_{k+1},\ldots ,x_{n}\}.} Thus a Gröbner basis for this ordering carries much more information than usually necessary. This may explain why Gröbner bases for the lexicographical ordering are usually the most difficult to compute.

Intersecting ideals

[edit]

IfI andJ are two ideals generated respectively by {f1, ...,fm} and {g1, ...,gk}, then a single Gröbner basis computation produces a Gröbner basis of their intersectionIJ. For this, one introduces a new indeterminatet, and one uses an elimination ordering such that the first block contains onlyt and the other block contains all the other variables (this means that a monomial containingt is greater than every monomial that does not containt). With this monomial ordering, a Gröbner basis ofIJ consists in the polynomials that do not containt, in the Gröbner basis of the ideal

K=tf1,,tfm,(1t)g1,,(1t)gk.{\displaystyle K=\langle tf_{1},\ldots ,tf_{m},(1-t)g_{1},\ldots ,(1-t)g_{k}\rangle .}

In other words,IJ is obtained byeliminatingt inK.This may be proven by observing that the idealK consists of the polynomials(ab)t+b{\displaystyle (a-b)t+b} such thataI{\displaystyle a\in I} andbJ{\displaystyle b\in J}. Such a polynomial is independent oft if and only ifa =b, which means thatbIJ.{\displaystyle b\in I\cap J.}

Implicitization of a rational curve

[edit]

Arational curve is analgebraic curve that has a set ofparametric equations of the form

x1=f1(t)g1(t)xn=fn(t)gn(t),{\displaystyle {\begin{aligned}x_{1}&={\frac {f_{1}(t)}{g_{1}(t)}}\\&\;\;\vdots \\x_{n}&={\frac {f_{n}(t)}{g_{n}(t)}},\end{aligned}}}

wherefi(t){\displaystyle f_{i}(t)} andgi(t){\displaystyle g_{i}(t)} are univariate polynomials for 1 ≤in. One may (and will) suppose thatfi(t){\displaystyle f_{i}(t)} andgi(t){\displaystyle g_{i}(t)} arecoprime (they have no non-constant common factors).

Implicitization consists in computing theimplicit equations of such a curve. In case ofn = 2, that is for plane curves, this may be computed with theresultant. The implicit equation is the following resultant:

Rest(g1x1f1,g2x2f2).{\displaystyle {\text{Res}}_{t}(g_{1}x_{1}-f_{1},g_{2}x_{2}-f_{2}).}

Elimination with Gröbner bases allows to implicitize for any value ofn, simply by eliminatingt in the idealg1x1f1,,gnxnfn.{\displaystyle \langle g_{1}x_{1}-f_{1},\ldots ,g_{n}x_{n}-f_{n}\rangle .}Ifn = 2, the result is the same as with the resultant, if the mapt(x1,x2){\displaystyle t\mapsto (x_{1},x_{2})} is injective for almost everyt. In the other case, the resultant is a power of the result of the elimination.

Saturation

[edit]

When modeling a problem by polynomial equations, it is often assumed that some quantities are non-zero, so as to avoid degenerate cases. For example, when dealing withtriangles, many properties become false if the triangle degenerates to a line segment, i.e. the length of one side is equal to the sum of the lengths of the other sides. In such situations, one cannot deduce relevant information from the polynomial system unless the degenerate solutions are ignored. More precisely, the system of equations defines analgebraic set which may have severalirreducible components, and one must remove the components on which the degeneracy conditions are everywhere zero.

This is done bysaturating the equations by the degeneracy conditions, which may be done via the elimination property of Gröbner bases.

Definition of the saturation

[edit]

Thelocalization of a ring consists in adjoining to it the formal inverses of some elements. This section concerns only the case of a single element, or equivalently a finite number of elements (adjoining the inverses of several elements is equivalent to adjoining the inverse of their product). Thelocalization of a ringR by an elementf is the ringRf=R[t]/(1ft),{\displaystyle R_{f}=R[t]/(1-ft),} wheret is a new indeterminate representing the inverse off. Thelocalization of an idealI ofR is the idealIf=RfI{\displaystyle I_{f}=R_{f}I} ofRf.{\displaystyle R_{f}.} WhenR is a polynomial ring, computing inRf{\displaystyle R_{f}} is not efficient because of the need to manage the denominators. Therefore, localization is usually replaced by the operation ofsaturation.

Thesaturation with respect tof of an idealI inR is the inverse image ofRfI{\displaystyle R_{f}I} under the canonical map fromR toRf.{\displaystyle R_{f}.} It is the idealI:f={gR(kN)fkgI}{\displaystyle I:f^{\infty }=\{g\in R\mid (\exists k\in \mathbb {N} )f^{k}g\in I\}} consisting in all elements ofR whose product with some power off belongs toI.

IfJ is the ideal generated byI and 1−ft inR[t], thenI:f=JR.{\displaystyle I:f^{\infty }=J\cap R.} It follows that, ifR is a polynomial ring, a Gröbner basis computation eliminatingt produces a Gröbner basis of the saturation of an ideal by a polynomial.

The important property of the saturation, which ensures that it removes from the algebraic set defined by the idealI theirreducible components on which the polynomialf is zero, is the following:Theprimary decomposition ofI:f{\displaystyle I:f^{\infty }}consists of the components of the primary decomposition of I that do not contain any power of f.

Computation of the saturation

[edit]

A Gröbner basis of the saturation byf of a polynomial ideal generated by a finite set of polynomialsF, may be obtained by eliminatingt inF{1tf},{\displaystyle F\cup \{1-tf\},} that is by keeping the polynomials independent oft in the Gröbner basis ofF{1tf}{\displaystyle F\cup \{1-tf\}} for an elimination ordering eliminatingt.

Instead of usingF, one may also start from a Gröbner basis ofF. Which method is most efficient depends on the problem. However, if the saturation does not remove any component, that is if the ideal is equal to its saturated ideal, computing first the Gröbner basis ofF is usually faster. On the other hand, if the saturation removes some components, the direct computation may be dramatically faster.

If one wants to saturate with respect to several polynomialsf1,,fk{\displaystyle f_{1},\ldots ,f_{k}} or with respect to a single polynomial which is a productf=f1fk,{\displaystyle f=f_{1}\cdots f_{k},} there are three ways to proceed which give the same result but may have very different computation times (it depends on the problem which is the most efficient).

Effective Nullstellensatz

[edit]

Hilbert's Nullstellensatz has two versions. The first one asserts that a set of polynomials has no common zeros over analgebraic closure of the field of the coefficients, if and only if 1 belongs to the generated ideal. This is easily tested with a Gröbner basis computation, because 1 belongs to an ideal if and only if 1 belongs to the Gröbner basis of the ideal, for any monomial ordering.

The second version asserts that the set of common zeros (in an algebraic closure of the field of the coefficients) of an ideal is contained in thehypersurface of the zeros of a polynomialf, if and only if a power off belongs to the ideal. This may be tested by saturating the ideal byf; in fact, a power off belongs to the ideal if and only if the saturation byf provides a Gröbner basis containing 1.

Implicitization in higher dimension

[edit]

By definition, an affinerational variety of dimensionk may be described byparametric equations of the form

x1=p1p0xn=pnp0,{\displaystyle {\begin{aligned}x_{1}&={\frac {p_{1}}{p_{0}}}\\&\;\;\vdots \\x_{n}&={\frac {p_{n}}{p_{0}}},\end{aligned}}}

wherep0,,pn{\displaystyle p_{0},\ldots ,p_{n}} aren+1 polynomials in thek variables (parameters of the parameterization)t1,,tk.{\displaystyle t_{1},\ldots ,t_{k}.} Thus the parameterst1,,tk{\displaystyle t_{1},\ldots ,t_{k}} and the coordinatesx1,,xn{\displaystyle x_{1},\ldots ,x_{n}} of the points of the variety are zeros of the ideal

I=p0x1p1,,p0xnpn.{\displaystyle I=\left\langle p_{0}x_{1}-p_{1},\ldots ,p_{0}x_{n}-p_{n}\right\rangle .}

One could guess that it suffices to eliminate the parameters to obtain the implicit equations of the variety, as it has been done in the case of curves. Unfortunately this is not always the case. If thepi{\displaystyle p_{i}} have a common zero (sometimes calledbase point), everyirreducible component of the non-empty algebraic set defined by thepi{\displaystyle p_{i}} is an irreducible component of the algebraic set defined byI. It follows that, in this case, the direct elimination of theti{\displaystyle t_{i}} provides an empty set of polynomials.

Therefore, ifk>1, two Gröbner basis computations are needed to implicitize:

  1. SaturateI{\displaystyle I} byp0{\displaystyle p_{0}} to get a Gröbner basisG{\displaystyle G}
  2. Eliminate theti{\displaystyle t_{i}} fromG{\displaystyle G} to get a Gröbner basis of the ideal (of the implicit equations) of the variety.

Algorithms and implementations

[edit]

Buchberger's algorithm is the oldest algorithm for computing Gröbner bases. It has been devised byBruno Buchberger together with the Gröbner basis theory. It is straightforward to implement, but it appeared soon that raw implementations can solve only trivial problems. The main issues are the following ones:

  1. Even when the resulting Gröbner basis is small, the intermediate polynomials can be huge. It results that most of the computing time may be spent inmemory management. So, specialized memory management algorithms may be a fundamental part of an efficient implementation.
  2. Theintegers occurring during a computation may be sufficiently large for makingfast multiplication algorithms andmultimodular arithmetic useful. For this reason, most optimized implementations use theGMPlibrary. Also,modular arithmetic,Chinese remainder theorem andHensel lifting are used in optimized implementations
  3. The choice of the S-polynomials to reduce and of the polynomials used for reducing them is devoted toheuristics. As in many computational problems, heuristics cannot detect most hidden simplifications, and if heuristic choices are avoided, one may get a dramatic improvement of the algorithm efficiency.
  4. In most cases most S-polynomials that are computed are reduced to zero; that is, most computing time is spent to compute zero.
  5. Themonomial ordering that is most often needed for the applications (pure lexicographic) is not the ordering that leads to the easiest computation, generally the orderingdegrevlex.

For solving 3. many improvements, variants andheuristics have been proposed before the introduction ofF4 and F5 algorithms byJean-Charles Faugère. As these algorithms are designed for integer coefficients or with coefficients in theintegers modulo a prime number, Buchberger's algorithm remains useful for more general coefficients.

Roughly speaking, F4 algorithm solves 3. by replacing many S-polynomial reductions by therow reduction of a single large matrix for which advanced methods oflinear algebra can be used. This solves partially issue 4., as reductions to zero in Buchberger's algorithm correspond to relations between rows of the matrix to be reduced, and the zero rows of the reduced matrix correspond to abasis of thevector space of these relations.

F5 algorithm improves F4 by introducing a criterion that allows reducing the size of the matrices to be reduced. This criterion is almost optimal, since the matrices to be reduced havefull rank in sufficiently regular cases (in particular, when the input polynomials form aregular sequence). Tuning F5 for a general use is difficult, since its performances depend on an order on the input polynomials and a balance between the incrementation of the working polynomial degree and of the number of the input polynomials that are considered. To date (2022), there is no distributed implementation that is significantly more efficient than F4, but, over modular integers F5 has been used successfully for severalcryptographic challenges; for example, for breakingHFE challenge.

Issue 5. has been solved by the discovery of basis conversion algorithms that start from the Gröbner basis for one monomial ordering for computing a Gröbner basis for another monomial ordering.FGLM algorithm is such a basis conversion algorithm that works only in the zero-dimensional case (where the polynomials have a finite number of complex common zeros) and has apolynomial complexity in the number of common zeros. A basis conversion algorithm that works is the general case is theGröbner walk algorithm.[6] In its original form, FGLM may be the critical step for solvingsystems of polynomial equations because FGML does not take into account thesparsity of involved matrices. This has been fixed by the introduction ofsparse FGLM algorithms.[7]

Most general-purposecomputer algebra systems have implementations of one or several algorithms for Gröbner bases, often also embedded in other functions, such as for solving systems of polynomial equations or for simplifying trigonometric functions; this is the case, for example, ofCoCoA,GAP,Macaulay 2,Magma,Maple,Mathematica,SINGULAR,SageMath andSymPy. When F4 is available, it is generally much more efficient than Buchberger's algorithm. The implementation techniques and algorithmic variants are not always documented, although they may have a dramatic effect on efficiency.

Implementations of F4 and (sparse)-FGLM are included in thelibraryMsolve.[8] Beside Gröbner algorithms, Msolve contains fast algorithms forreal-root isolation, and combines all these functions in an algorithm for the real solutions ofsystems of polynomial equations that outperforms dramatically the other software for this problem (Maple and Magma).[8] Msolve is available onGitHub, and is interfaced withJulia, Maple and SageMath; this means that Msolve can be used directly from within these software environments.

Complexity

[edit]

Thecomplexity of the Gröbner basis computations is commonly evaluated in term of the numbern of variables and the maximal degreed of the input polynomials.

In the worst case, the main parameter of the complexity is the maximal degree of the elements of the resulting reduced Gröbner basis. More precisely, if the Gröbner basis contains an element of a large degreeD, this element may containΩ(Dn){\displaystyle \Omega (D^{n})} nonzero terms whose computation requires a time ofΩ(Dn)>DΩ(n).{\displaystyle \Omega (D^{n})>D^{\Omega (n)}.} On the other hand, if all polynomials in the reduced Gröbner basis ahomogeneous ideal have a degree of at mostD, the Gröbner basis can be computed bylinear algebra on thevector space of polynomials of degree less than2D, which has a dimensionO(Dn).{\displaystyle O(D^{n}).}[1] So, the complexity of this computation isO(Dn)O(1)=DO(n).{\displaystyle O(D^{n})^{O(1)}=D^{O(n)}.}

The worst-case complexity of a Gröbner basis computation is doubly exponential inn. More precisely, the complexity is upper bounded by a polynomial ind2n.{\textstyle d^{2^{n}}.} Usinglittle o notation, it is therefore bounded byd2n+o(n).{\textstyle d^{2^{n+o(n)}}.} On the other hand, examples have been given of reduced Gröbner bases containing polynomials of degreed2Ω(n),{\textstyle d^{2^{\Omega (n)}},}or containingd2Ω(n){\textstyle d^{2^{\Omega (n)}}} elements. As every algorithm for computing a Gröbner basis must write its result, this provides a lower bound of the complexity.

Gröbner basis isEXPSPACE-complete.[9]

Generalizations

[edit]

The concept and algorithms of Gröbner bases have been generalized tosubmodules offree modules over a polynomial ring. In fact, ifL is a free module over a ringR, then one may consider thedirect sumRL{\displaystyle R\oplus L} as a ring by defining the product of two elements ofL to be0. This ring may be identified withR[e1,,el]/{eiej|1ijl}{\displaystyle R[e_{1},\ldots ,e_{l}]/\left\langle \{e_{i}e_{j}|1\leq i\leq j\leq l\}\right\rangle }, wheree1,,el{\displaystyle e_{1},\ldots ,e_{l}} is a basis ofL. This allows identifying a submodule ofL generated byg1,,gk{\displaystyle g_{1},\ldots ,g_{k}} with the ideal ofR[e1,,el]{\displaystyle R[e_{1},\ldots ,e_{l}]} generated byg1,,gk{\displaystyle g_{1},\ldots ,g_{k}} and the productseiej{\displaystyle e_{i}e_{j}},1ijl{\displaystyle 1\leq i\leq j\leq l}. IfR is a polynomial ring, this reduces the theory and the algorithms of Gröbner bases of modules to the theory and the algorithms of Gröbner bases of ideals.

The concept and algorithms of Gröbner bases have also been generalized to ideals over various rings, commutative or not, like polynomial rings over aprincipal ideal ring orWeyl algebras.

Areas of applications

[edit]

Error-Correcting Codes

[edit]

Gröbner bases have been applied in the theory of error-correcting codes for algebraic decoding. By using Gröbner basis computation on various forms of error-correcting equations, decoding methods were developed for correcting errors of cyclic codes,[10] affine variety codes,[11] algebraic-geometric codes and even general linear block codes.[12] Applying Gröbner basis in algebraic decoding is still a research area of channelcoding theory.

See also

[edit]

References

[edit]
  1. ^abcLazard, Daniel (1983). "Gröbner bases, Gaussian elimination and resolution of systems of algebraic equations".Computer Algebra. Lecture Notes in Computer Science. Vol. 162. pp. 146–156.doi:10.1007/3-540-12868-9_99.ISBN 978-3-540-12868-7.
  2. ^Renschuch, Bodo; Roloff, Hartmut; Rasputin, Georgij G.; Abramson, Michael (June 2003)."Contributions to constructive polynomial ideal theory XXIII: Forgotten works of Leningrad mathematician N. M. Gjunter on polynomial ideal theory"(PDF).ACM SIGSAM Bulletin.37 (2):35–48.doi:10.1145/944567.944569.S2CID 1819694.
  3. ^Cox, David A.; Little, John;O'Shea, Donal (1997).Ideals, Varieties, and Algorithms: An Introduction to Computational Algebraic Geometry and Commutative Algebra. Springer.ISBN 0-387-94680-2.
  4. ^Lazard, Daniel (2021)."Degree of a polynomial ideal and Bézout inequalities".
  5. ^Ene, Viviana; Herzog, Jürgen (2012).Gröbner Bases in Commutative Algebra. Graduate Studies in Mathematics. Vol. 130. Providence, RI: American Mathematical Society.ISBN 978-0-8218-7287-1.: Proposition 4.29
  6. ^Collart, Stéphane; Kalkbrener, Michael; Mall, Daniel (1997)."Converting bases with the Gröbner walk".Journal of Symbolic Computation.24 (3–4). Elsevier:465–469.doi:10.1006/jsco.1996.0145.
  7. ^Faugère, Jean-Charles; Chenqi, Mou (2017)."Sparse FGLM algorithms".Journal of Symbolic Computation.80. Elsevier:538–569.arXiv:1304.1238.doi:10.1016/j.jsc.2016.07.025.S2CID 149627.
  8. ^abBerthomieu \first1=Jérémy; Eder, Christian; Safey El Din, Mohab (2021).Msolve: a library for solving polynomial systems. 2021 International Symposium on Symbolic and Algebraic Computation. 46th International Symposium on Symbolic and Algebraic Computation. Saint Petersburg, Russia.arXiv:2104.03572.doi:10.1145/3452143.3465545.{{cite conference}}: CS1 maint: numeric names: authors list (link)
  9. ^Mayr, Ernst W. (September 1997), "Some Complexity Results for Polynomial Ideals",Journal of Complexity,13 (3):303–325,doi:10.1006/jcom.1997.0447
  10. ^Chen, X.; Reed, I.S.; Helleseth, T.; Truong, T.K. (1994)."Use of Gröbner bases to decode binary cyclic codes up to the true minimum distance".IEEE Transactions on Information Theory.40 (5):1654–61.doi:10.1109/18.333885.
  11. ^Fitzgerald, J.; Lax, R.F. (1998). "Decoding affine variety codes using Gröbner bases".Designs, Codes and Cryptography.13 (2):147–158.doi:10.1023/A:1008274212057.S2CID 2515114.
  12. ^Bulygin, S.; Pellikaan, R. (2009). "Decoding linear error-correcting codes up to half the minimum distance with Gröbner bases".Gröbner Bases, Coding, and Cryptography.Springer. pp. 361–5.ISBN 978-3-540-93805-7.

Further reading

[edit]

External links

[edit]
Bydegree
By properties
Tools and algorithms
International
National
Other
Retrieved from "https://en.wikipedia.org/w/index.php?title=Gröbner_basis&oldid=1296420479"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2025 Movatter.jp