Algorithm Search⌘K

Proofs: F4 algorithm

These are the human proofs of F4 algorithm, at full length.

F4 terminates and returns the reduced Gröbner basis full

This proof shows: correctness, termination.

Let KK be a field, R=K[x1,…,xn]R = K[x_1, \ldots, x_n] with a monomial order >>, and I=⟨f1,…,fs⟩I = \langle f_1, \ldots, f_s \rangle. We say that ff reduces to zero modulo a finite set GG, written f→G0f \to_G 0, if f=∑kckukgkf = \sum_k c_k u_k g_k with ck∈Kc_k \in K, monomials uku_k, gk∈Gg_k \in G, and LM⁡(ukgk)≤LM⁡(f)\operatorname{LM}(u_k g_k) \le \operatorname{LM}(f) for all kk. This is a standard representation, and it remains one when GG grows. We use two facts.

  • Buchberger's criterion (CLO Chapter 2, Section 6, Theorem 6 and Chapter 2, Section 10). A finite G⊆IG \subseteq I that generates II is a Gröbner basis if the S-polynomial of every pair has a standard representation, or more generally an lcm-representation, with respect to GG.
  • Gebauer and Möller (1988). The pairs removed by UPDATE are covered: if all pairs that are never removed have S-polynomials with standard representations with respect to the final GG, then so do the removed pairs, through lcm-representations. This is the theorem behind the Gebauer-Möller criteria.

Lemma 1 (rows lie in the ideal). At all times G⊆IG \subseteq I, and GG generates II.

Proof. Line L1 puts the nonzero fif_i into GG. Every row of HH is a monomial multiple ugu g of some g∈Gg \in G: the halves at line L4 by construction, and the reducers at line L5 by construction. The row echelon form N~\tilde N consists of KK-linear combinations of rows of HH, so N~⊆I\tilde N \subseteq I. □\square

Lemma 2 (symbolic preprocessing). SYMBOLIC-PREPROCESSING terminates, and its output HH has this property: for every monomial mm that occurs in some row of HH and is divisible by LM⁡(g)\operatorname{LM}(g) for some g∈Gg \in G, there is a row of HH of the form ug′u g' with g′∈Gg' \in G and LM⁡(ug′)=m\operatorname{LM}(u g') = m.

Proof. Each step marks one monomial as handled and possibly adds a row ugu g whose leading monomial is that monomial mm, so every other monomial of the new row is smaller than mm. Consider the graph whose nodes are the handled monomials, with an edge from mm to each monomial of the row added for mm. Each node has finitely many successors, and every path is strictly decreasing for >>. A monomial order is a well-order, so there is no infinite path, and by König's lemma only finitely many monomials are reached from the finitely many monomials of LL. So the loop ends. At the end every monomial of every row is handled. For a monomial mm in LM⁡(L)\operatorname{LM}(L), the half itself is a row ugu g with leading monomial mm; for every other handled mm divisible by some LM⁡(g)\operatorname{LM}(g), the reducer was added. □\square

Lemma 3 (the span reduces to zero). Let GG be the basis before line L8, HH the rows of a round, N~\tilde N the row echelon form and N~+\tilde N^+ the new rows. Then every ff in the KK-span of HH satisfies f→G∪N~+0f \to_{G \cup \tilde N^+} 0.

Proof. By induction on LM⁡(f)\operatorname{LM}(f), which is well founded. If f=0f = 0 there is nothing to show. Otherwise N~\tilde N is a basis of the span with distinct leading monomials, so LM⁡(f)=LM⁡(h)\operatorname{LM}(f) = \operatorname{LM}(h) for some h∈N~h \in \tilde N, and m=LM⁡(f)m = \operatorname{LM}(f) is a monomial of some row of HH.

  • If m∈LM⁡(H)m \in \operatorname{LM}(H), then mm is the leading monomial of a row r=ugr = u g of HH, and every element of LM⁡(H)\operatorname{LM}(H) is divisible by the leading monomial of an element of GG (by Lemma 2 and the form of the halves). Then f′=f−(LC⁡(f)/LC⁡(r)) rf' = f - (\operatorname{LC}(f)/\operatorname{LC}(r))\, r is in the span of HH and has LM⁡(f′)<m\operatorname{LM}(f') < m.
  • If m∉LM⁡(H)m \notin \operatorname{LM}(H), then h∈N~+h \in \tilde N^+, and f′=f−(LC⁡(f)/LC⁡(h)) hf' = f - (\operatorname{LC}(f)/\operatorname{LC}(h))\, h is in the span and has LM⁡(f′)<m\operatorname{LM}(f') < m.

In both cases the subtracted term is c ugc\, u g with g∈G∪N~+g \in G \cup \tilde N^+ and LM⁡(ug)=m≤LM⁡(f)\operatorname{LM}(u g) = m \le \operatorname{LM}(f), and by induction f′f' has a standard representation with terms of leading monomial at most LM⁡(f′)<m\operatorname{LM}(f') < m. Adding the two gives one for ff. □\square

Lemma 4 (new leading monomials). For h∈N~+h \in \tilde N^+, LM⁡(h)\operatorname{LM}(h) is divisible by no LM⁡(g)\operatorname{LM}(g) with g∈Gg \in G.

Proof. LM⁡(h)\operatorname{LM}(h) is a monomial of some row of HH, since hh is a combination of rows. If LM⁡(g)∣LM⁡(h)\operatorname{LM}(g) \mid \operatorname{LM}(h) for some g∈Gg \in G, Lemma 2 gives a row of HH with leading monomial LM⁡(h)\operatorname{LM}(h), so LM⁡(h)∈LM⁡(H)\operatorname{LM}(h) \in \operatorname{LM}(H), which contradicts h∈N~+h \in \tilde N^+. □\square

Theorem. F4 terminates and returns the reduced Gröbner basis of II.

Proof. Termination. Let J=⟨LM⁡(G)⟩J = \langle \operatorname{LM}(G) \rangle at the start of a round. By Lemma 4, every h∈N~+h \in \tilde N^+ has LM⁡(h)∉J\operatorname{LM}(h) \notin J, so a round with N~+≠∅\tilde N^+ \ne \emptyset makes JJ strictly larger. By Dickson's lemma (equivalently, RR is Noetherian) a strictly ascending chain of monomial ideals is finite, so only finitely many rounds add elements, each adds finitely many, and GG stays finite. Pairs are created only by UPDATE, at most ∣G∣|G| for each new element, and each round removes at least one pair from PP at line L3. So the loop ends.

Correctness. Every pair (i,j)(i, j) that is selected at line L3 has its two halves in L⊆HL \subseteq H, so its S-polynomial is in the span of HH, and by Lemma 3 it reduces to zero modulo G∪N~+G \cup \tilde N^+, which is contained in the final GG. Every pair that is never selected was removed by UPDATE, and by the theorem of Gebauer and Möller its S-polynomial also has an lcm-representation with respect to the final GG. By Lemma 1 the final GG lies in II and generates it, so by Buchberger's criterion it is a Gröbner basis. Line L9 removes elements whose leading monomial is divisible by another one, makes the rest monic, and reduces their tails, which gives the reduced Gröbner basis; it is unique (CLO Chapter 2, Section 7). ■\blacksquare

Back to F4 algorithm