Solving systems of equations in supernilpotent algebras
Erhard Aichinger
Institute for Algebra, Johannes Kepler University Linz, Linz, Austria
[email protected]
Abstract.
Recently, M. Kompatscher proved that for each finite supernilpotent algebra
A
in a congruence modular variety, there is a polynomial time algorithm
to solve polynomial equations over this algebra. Let μ be the maximal
arity of the fundamental operations of A, and let
[TABLE]
Applying a method that G. Károlyi and C. Szabó had used to solve
equations over finite nilpotent rings, we show that for A,
there is c∈N such that a solution of every system of s
equations in n variables
can be found
by testing
at most cnsd (instead of all ∣A∣n possible)
assignments to the variables. This also yields new information
on some circuit satisfiability problems.
Key words and phrases:
Supernilpotent algebras, polynomial equations, polynomial mappings,
circuit satisfiability
2010 Mathematics Subject Classification:
08A40 (68Q25)
Supported by the Austrian Science Fund (FWF):P29931.
1. Introduction
We study systems of polynomial equations over
a finite algebraic structure A. Such a system
is given by equations of the form p(x1,…,xn)≈q(x1,…,xn), where p,q
are polynomial terms of A; a polynomial term
of A is a term of the algebra A∗ which
is obtained by expanding A with one nullary
function symbol for each a∈A.
A solution to a system pi(x1,…,xn)≈qi(x1,…,xn) (i=1,…,s)
is an element \mathbfsla=(a1,…,an)∈An such that
piA(\mathbfsla)=qiA(\mathbfsla) for all
i∈{1,…,s}.
The problem to decide whether such a solution exists
has been called \textscPolSysSat(A), and \textscPolSat(A)
if the system consists of one single equation, and
the terms of the input are encoded as strings
over {x1,…,xn}∪A∪F, where
F is the set of function symbols of A.
A survey of results on the computational complexity
of this problem is given, e.g., in [IK18, Kom18].
In algebras such as groups, rings or Boolean algebras,
one can reduce an equation p(\mathbfslx)≈q(\mathbfslx) to an equation of the form f(\mathbfslx)≈y, where
y∈A. A system of equations of this form
then has the form fi(\mathbfslx)≈yi
(i=1,…,s).
For n∈N, let Poln(A) denote the n-ary polynomial
functions on A [MMT87, Definition 4.4].
For a finite nilpotent ring or group A,
[Hor11] establishes the existence of
a natural number dA such that
for every f∈Poln(A)
and for every \mathbfsla∈An, there exists
\mathbfslb such that fA(\mathbfsla)=fA(\mathbfslb) and \mathbfslb has at most
dA components that are different
from [math]. Hence the equation f(\mathbfslx)≈y
has a solution if and only if it has a solution
with at most dA nonzero entries.
Thus for the algebra A, testing only vectors with
at most dA nonzero entries is an algorithm,
which, given an equation f(\mathbfslx)≈y of length
n, takes at most c(\mathbfslA)⋅ndA+1 many steps
to find whether this equation is solvable:
there are at most
∑i=0dA(in)(∣A∣−1)i≤c1(A)⋅ndA
many evaluations to be done, each of them
taking at most c2(A)⋅n many steps.
The number dA in [Hor11]
is obtained from Ramsey’s Theorem and therefore
rather large.
In [Kom18], it is proved
that for every finite supernilpotent algebra in a congruence
modular variety, such a number dA exists, again
using Ramsey’s Theorem.
For rings, lower values
of dA have been obtained in [KS18] (cf. [KS15]).
In [Föl17, Föl18], A. Földvári provides polynomial
time algorithms for solving equations over finite nilpotent groups and rings
relying on the structure theory of these algebras.
In this paper, we extend the method developed in [KS18]
from finite nilpotent rings to arbitrary finite supernilpotent
algebras in congruence modular varieties. For such algebras,
we compute dA as ∣A∣log2(μ)+log2(∣A∣)+1 (Theorem 10).
The technique that allows to generalize Károlyi’s and Szabó’s method
is the coordinatization of nilpotent algebras of prime power order
by elementary abelian
groups from [Aic18, Theorem 4.2].
The method generalizes to systems of equations:
we show for a given finite supernilpotent
algebra A in a congruence modular variety, and a given s∈N0,
there is a polynomial time algorithm to test whether a system of at most s
polynomial equations over A has a solution.
If s is not fixed in advance, then [LZ06, Corollary 3.13]
implies that if A is not abelian,
\textscPolSysSat(A) is NP-complete.
Let us finally explain to which class of algebras our results applies:
A finite algebra A from a congruence modular variety
with finitely many fundamental operations is supernilpotent
if and only if it is a direct product of nilpotent algebras
of prime power order; modulo notational differences
explained, e.g., in [Aic18, Lemma 2.4],
this result has been proved in
[Kea99, Theorem 3.14].
Such an algebra is therefore always nilpotent,
has a Mal’cev term (cf.
[FM87, Theorem 6.2], [Kea99, Theorem 2.7]),
and hence generates a congruence permutable variety.
For a more detailed introduction to supernilpotency
and, for k∈N, to k-supernilpotency, we refer
to [AM10, AMO18, Aic18].
2. A theorem of Károlyi and Szabó
In this section, we state a special case of [KS18, Theorem 3.1].
Since their result is much more general than needed for our purpose, we
also include a self-contained proof, which is a reduction
Károlyi’s and Szabó’s proof to the case of elementary abelian groups.
For n∈N={1,2,3,…}, we denote the set {1,2,…,n} by n.
Let A be a set with an element 0∈A, and let J⊆n.
For \mathbfsla∈An, \mathbfsla(J) is defined by \mathbfsla(J)∈An, \mathbfsla(J)(j)=\mathbfsla(j) for
j∈J and \mathbfsla(J)(j)=0 for i∈n∖J.
Suppose that 1 is an element of A. Then by \mathbfsl1, we denote the vector
(1,1,…,1) in An, and for J⊆n, \mathbfsl1(J) is the
vector (v1,…,vn) with vj=1 if j∈J and vj=0 if j∈J.
For any sets C,D, we write C⊂D for (C⊆D and C=D).
We first need the following variation of [Bri11, Theorem 1] and
[KS18, Theorem 3.2], which is proved
using several arguments from the proof of [Alo99, Theorem 3.1] and
from [Bri11].
Lemma 1**.**
Let F be a finite field, let k,m,n∈N, let q:=∣F∣,
let p1,…,pm∈F[x1,…,xn] be polynomials
such that for each i∈m,
each monomial of pi contains at most k variables.
Then there exists J⊆n such that
∣J∣≤km(q−1) and pi(\mathbfsl1(J))=pi(\mathbfsl1)
for all i∈m.
Proof.
We proceed by induction on n. If n≤km(q−1), then
we take J:=n.
For the induction step, we assume that n>km(q−1).
We first produce a set J1⊂n such that
pi(\mathbfsl1(J1))=pi(\mathbfsl1) for all i∈m.
Seeking a contradiction, we suppose that no such J1
exists. Following an idea from the proof of [Alo99, Theorem 3.1],
we consider the polynomials
[TABLE]
We first show that for all \mathbfsla∈{0,1}n, q2(\mathbfsla)=0.
To this end, we first consider the case \mathbfsla=\mathbfsl1.
Then
q2(\mathbfsla)=1−∏i=1m1=0.
If \mathbfsla∈{0,1}n∖{\mathbfsl1}, then
by the assumptions, there is i∈m
such that pi(\mathbfsla)=pi(\mathbfsl1).
Then 1−(pi(\mathbfsla)−pi(\mathbfsl1))q−1=0.
Therefore q2(\mathbfsla)=0.
Hence the polynomial q2 vanishes at {0,1}n. By the Combinatorial
Nullstellensatz [Alo99, Theorem 1.1] applied to gj(xj):=xj2−xj,
q2
then lies in the ideal V of F[x1,…,xn] generated by
G={xj2−xj∣j∈n}. Hence
x1x2⋯xn−q1(x1,…,xn)∈V. Since
the leading monomials of the polynomials in G are coprime,
G is a Gröbner basis of V
(with respect to x1>x2>⋯>xn, lexicographic order, cf. [Eis95, p.337]).
Therefore,
reducing q1(x1,…,xn) modulo G, we must obtain x1x2⋯xn as
the remainder (as defined, e.g., in [Eis95, p.334]).
Because of the form of all polynomials in G
(all variables of gj occur in the leading term of gj),
none of the reduction steps
increases the number of variables in any monomial. Therefore,
q1(x1,…,xn) must contain a monomial that contains all n variables.
Computing the expansion of q1 by multiplying out all products from
its definition, we see
that each monomial in q1 contains at most km(q−1) variables.
Hence n≤km(q−1), which contradicts the assumption n>km(q−1).
This contradiction shows that there is set J1⊂n such that
pi(\mathbfsl1(J1))=pi(\mathbfsl1) for all i∈m.
Now we let n′:=∣J1∣, and we assume that
J1={j1,…,jn′} with j1<⋯<jn′.
For i∈m, we define pi′∈F[y1,…,yn′] by
[TABLE]
By the induction hypothesis, there exists J2⊆n′ with
∣J2∣≤km(q−1) such that
pi′(\mathbfsl1(J2))=pi′(\mathbfsl1) for all i∈m.
Now we define
J:={jt∣t∈J2}. We have J⊆J1, and therefore
\mathbfsl1(J)=\mathbfsl(\mathbfsl1(J))(J1).
Then pi(\mathbfsl1(J))=pi(\mathbfsl(\mathbfsl1(J))(J1))=pi′(\mathbfsl1(J)(j1),…,\mathbfsl1(J)(jn′))=pi′(\mathbfsl1(J2))=pi′(\mathbfsl1)=pi(\mathbfsl1(J1))=pi(\mathbfsl1), which completes the induction step. ∎
We will need the following special case of [KS18, Theorem 3.1].
Let Pk(n) denote the set
{I⊆n:∣I∣≤k} of subsets of n with at most k elements.
Theorem 2** (cf. [KS18, Theorem 3.1]).**
Let n∈N, let k∈N0, let p be a prime,
and let m∈N.
Let φ:Pk(n)→Zpm. Then there is
U⊆n with ∣U∣≤km(p−1) such that
[TABLE]
Proof.
We denote the vector φ(J) by ((φ(J))1,…,(φ(J))m), and
we define m polynomial functions f1,…,fm∈Zp[x1,…,xn]
by
[TABLE]
for i∈m.
By Lemma 1, there is a subset U of n
with ∣U∣≤km(p−1) such that
for all i∈m, we have
fi(\mathbfsl1)=fi(\mathbfsl1(U)). Hence
∑J∈Pk(n)(φ(J))i=fi(\mathbfsl1)=fi(\mathbfsl1(U))=∑J∈Pk(n),J⊆U(φ(J))i=∑J∈Pk(U)(φ(J))i. ∎
3. Absorbing components
Let A be a set, let 0A be an element of A, let B=(B,+,−,0)
be an abelian group, let n∈N, let f:An→B, and let I⊆n.
By Dep(f) we denote the set \{i\in\underline{n}\mid f\text{ depends on its i\,th argument}\}.
We say that f is absorbing in its j th argument if for all \mathbfsla=(\mathbfsla(1),…,\mathbfsla(n))∈An
with \mathbfsla(j)=0A we have f(\mathbfsla)=0. In the sequel, we will denote 0A simply by
[math].
We say that f is absorbing in I if Dep(f)⊆I and for every
i∈I, f is absorbing in its i th argument.
Lemma 3**.**
Let A be a set, let [math] be an element of A, let B=(B,+,−,0)
be an abelian group, let n∈N, and
let f:An→B. Then there is exactly one sequence (fI)I⊆n of functions from An to B such that
for each I⊆n, fI is absorbing in I and f=∑I⊆nfI.
Furthermore, each function fI lies in the subgroup F of BAn that is
generated by the functions \mathbfslx↦f(\mathbfslx(I)), where
I⊆n.
Proof.
We first prove the existence of such a sequence.
To this end, we define fI by recursion on ∣I∣.
We define f∅(\mathbfsla):=f(0,…,0) and for I=∅,
we let
[TABLE]
By induction on ∣I∣, we see that Dep(fI)⊆I and
that fI lies in the subgroup F.
We will
now show that each fI is absorbing in I, and we again proceed
by induction on ∣I∣. Let i∈I, and let \mathbfsla∈An be
such that \mathbfsla(i)=0. We have to show fI(\mathbfsla)=0. We compute
fI(\mathbfsla)=f(\mathbfsla(I))−∑J⊂IfJ(\mathbfsla).
By the induction hypothesis, we have fJ(\mathbfsla)=0 for those J with
i∈J. Hence f(\mathbfsla(I))−∑J⊂IfJ(\mathbfsla)=f(\mathbfsla(I))−∑J⊆I∖{i}fJ(\mathbfsla), and because
of \mathbfsla(I)=\mathbfsla(I∖{i}), this is equal to
f(\mathbfsla(I∖{i}))−∑J⊆I∖{i}fJ(\mathbfsla)=f(\mathbfsla(I∖{i}))−∑J⊂I∖{i}fJ(\mathbfsla)−fI∖{i}(\mathbfsla).
By the definition of fI∖{i}, the last expression is equal to
fI∖{i}(\mathbfsla)−fI∖{i}(\mathbfsla)=0. This completes the induction
proof; hence each fI is absorbing in I.
In order to show f=∑I⊆nfI, we choose \mathbfsla∈An and
compute ∑I⊆nfI(\mathbfsla)=fn(\mathbfsla)+∑I⊂nfI(\mathbfsla)=f(\mathbfsla(n))−∑J⊂nfJ(\mathbfsla)+∑I⊂nfI(\mathbfsla)=f(\mathbfsla).
This completes the proof of the existence of such a sequence.
For the uniqueness, assume that f=∑I⊆nfI=∑I⊆ngI and
that for all I, fI and gI are absorbing in I. We show by induction on ∣I∣ that
fI=gI. Let I:=∅.
First we notice that f(0,…,0)=∑J⊆nfJ(0,…,0)=∑J⊆ngJ(0,…,0).
Since fJ and gJ are absorbing, the summands with J=∅ are [math],
and thus f∅(0,…,0)=∑J⊆nfJ(0,…,0)=f(0,…,0)=∑J⊆ngJ(0,…,0)=g∅(0,…,0).
Since both f∅ and g∅ are constant functions, they are equal.
For the induction step, we assume ∣I∣≥1.
Let \mathbfsla∈An.
Then ∑J⊆nfJ(\mathbfsla(I))=∑J⊆ngJ(\mathbfsla(I)).
Only the summands with J⊆I can be nonzero, and therefore
∑J⊆IfJ(\mathbfsla(I))=∑J⊆IgJ(\mathbfsla(I)).
By the induction hypothesis, fJ=gJ for J⊂I. Therefore,
fI(\mathbfsla(I))=gI(\mathbfsla(I)). Since fI and gI depend only on the arguments
at positions in I, we obtain fI(\mathbfsla)=fI(\mathbfsla(I))=gI(\mathbfsla(I))=gI(\mathbfsla).
Thus fI=gI. ∎
Actually, the component fI can be computed by
fI(\mathbfsla)=∑J⊆I(−1)∣I∣+∣J∣f(\mathbfsla(J)).
Definition 4**.**
Let A be a set, let [math] be an element of A, let B=(B,+,−,0)
be an abelian group, let n∈N, let f:An→B, and let J⊆n.
Then we call the sequence (fI)I⊆n such that
for each I⊆n, fI is absorbing in I, and f=∑I⊆nfI the
absorbing decomposition of f, and fJ the J-absorbing component of f.
We define the absorbing degree of f by
adeg(f):=max({−1}∪{∣J∣:J⊆n and fJ=0}).
Theorem 5**.**
Let A be a set, let [math] be an element of A,
let p be a prime, let k∈N0, let n∈N, and let f1,…,fm:An→Zp.
We assume that each fi is of absorbing degree at most k.
Let \mathbfsla∈An. Then there is U with ∣U∣≤km(p−1) such that
for all i∈m, we have fi(\mathbfsla)=fi(\mathbfsla(U)).
Proof.
We define a function φ:Pk(n)→Zpm by
φ(J):=((f1)J(\mathbfsla),…,(fm)J(\mathbfsla)), where
for i∈m, \bigl{(}(f_{i})_{J}\bigr{)}_{J\subseteq\underline{n}}
is the absorbing
decomposition of fi.
Then Theorem 2 yields a subset U of n with
∣U∣≤km(p−1) such that
∑J∈Pk(n)φ(J)=∑J∈Pk(U)φ(J).
Since (fi)J=0 for all J with ∣J∣>k,
we have
∑J∈Pk(n)φ(J)=∑J∈Pk(n)((f1)J(\mathbfsla),…,(fm)J(\mathbfsla))=∑J⊆n((f1)J(\mathbfsla),…,(fm)J(\mathbfsla))=(f1(\mathbfsla),…,fm(\mathbfsla)) and
[TABLE]
∎
4. Polynomial mappings
In this section, we develop a property of polynomial mappings
of finite supernilpotent algebras in congruence modular varieties.
We call an algebra A=(A,+,−,0,(fi)i∈S) an
expanded group if its reduct A+=(A,+,−,0)
is a group, an expanded abelian group if
A+ is an abelian group, and an
expanded elementary abelian group if A+ is
elementary abelian, meaning that A+ is abelian and
all its nonzero elements have the same prime order.
For an algebra A and n,s∈N, we define the set Poln,s(A)
of
polynomial maps from An to As as the set
of all mappings \mathbfsla↦(f1(\mathbfsla),…,fs(\mathbfsla)) with
f1,…,fs∈Poln(A).
Lemma 6**.**
Let k,n∈N, let A be a k-supernilpotent
expanded abelian group, and let f∈Poln(A).
Then f is of
absorbing degree at most k.
Proof.
Let J⊆n with ∣J∣>k, and let
fJ be the J-absorbing component of f.
Let m:=∣J∣ and let J={i1,…,im}.
Using Lemma 3, we obtain that the
function g:Am→A defined by
g(ai1,…,aim):=fJ(\mathbfsla) for a∈An
is an absorbing function in Polm(A). Hence
[Aic18, Lemma 2.3] and the remark immediately
preceding that Lemma yield that g is the zero function.
Thus fJ=0. Hence the absorbing degree of f is at most k.
∎
We first consider polynomial mappings of supernilpotent
expanded elementary abelian groups of prime power order.
Theorem 7**.**
Let k,n,s,α∈N,
let p∈P,
and let
A be a k-supernilpotent
expanded elementary abelian group of
order pα.
Let F=(f1,…,fs)∈Poln,s(A), and let \mathbfsla∈An.
Then there is U⊆n with ∣U∣≤ksα(p−1) such that
F(\mathbfsla)=F(\mathbfsla(U)).
Proof.
We let π be a group isomorphism from (A,+,−,0) to
Zpα, and for a∈A, we denote
π(a) by (π1(a),…,πα(a)).
For each r∈s and each β∈α,
let fr,β:An→Zp be defined by
fr,β(\mathbfsla)=πβ(fr(\mathbfsla)); hence fr,β(\mathbfsla)
is the β th component of
fr(\mathbfsla). Since fr∈Poln(A) and
A is k-supernilpotent, Lemma 6 implies that
each of these fr,β
is of absorbing degree at most k.
Setting m:=sα, Theorem 5 yields
U with ∣U∣≤ksα(p−1) such that fr,β(\mathbfsla)=fr,β(\mathbfsla(U))
for all r∈s and β∈α.
Then clearly F(\mathbfsla)=F(\mathbfsla(U)). ∎
We apply this result to polynomial mappings of direct products
of finite supernilpotent
expanded elementary abelian groups.
For a vector \mathbfsla∈An, we call the number of its nonzero entries
the weight of a; formally,
wt(\mathbfsla):=∣{j∈n:\mathbfsla(j)=0}∣.
Theorem 8**.**
Let n,s,t,k1,…,kt∈N.
For
each i∈t,
let Bi a ki-supernilpotent
expanded elementary abelian group with ∣Bi∣=piαi,
where pi is a prime and αi∈N.
Let A:=∏i=1tBi, let
F∈Poln,s(A), and let \mathbfsla∈An.
Then there is \mathbfsly∈An with
wt(\mathbfsly)≤∑i=1tkisαi(pi−1) such that
F(\mathbfsla)=F(\mathbfsly).
Proof.
For i∈t, let
νi be the i th projection kernel.
Applying Theorem 7 to
A/νi, which is isomorphic to Bi, and
\mathbfslb:=\mathbfsla/νi, we obtain
Ui⊆n with ∣Ui∣≤kisαi(pi−1) such that
FA/νi(\mathbfslb(Ui))=FA/νi(\mathbfslb).
Lifting \mathbfslb(Ui) to A, we obtain (xi,1,…,xi,n)∈An such that
(xi,1,…,xi,n)/νi=\mathbfslb(Ui) and
xi,j=0 for j∈n∖Ui.
Now for every j∈n, we define yj∈A
by the equations
[TABLE]
For each i∈t, we have F(y1,…,yn)/νi=FA/νi(xi,1/νi,…,xi,n/νi)=FA/νi(\mathbfslb(Ui))=FA/νi(\mathbfslb)=FA/νi(\mathbfsla/νi)=F(\mathbfsla)/νi.
Hence F(\mathbfsly)=F(\mathbfsla).
For j∈n∖(U1∪⋯∪Ut),
and for all i∈t, we have xi,j=0, and therefore
yj=0. Hence the number of nonzero entries in \mathbfsly is
at most ∑i=1t∣Ui∣=∑i=1tkisαi(pi−1). ∎
Now we consider arbitrary finite supernilpotent algebras
in congruence modular varieties. In these algebras,
we can introduce group operations preserving nilpotency using [Aic18].
Lemma 9**.**
Let μ∈N,
let A=(A,(fi)i∈S) be a finite supernilpotent algebra in a congruence
modular variety all of whose fundamental operations have arity at most μ, and let z∈A.
Let t∈N0, let p1,…,pt be different primes, and let α1,…,αt∈N such that ∣A∣=∏i=1tpiαi.
For i∈t, let
ki:=(μ(piαi−1))αi−1.
Then there are operations + (binary), − (unary), [math] (nullary) on A such that
A′=(A,+,−,0,(fi)i∈S)
is isomorphic to
a direct product
∏i=1tBi′, where each Bi′ is a
ki-supernilpotent expanded elementary abelian group,
and 0A′=z.
Proof.
Since the result is true for ∣A∣=1, we henceforth assume ∣A∣≥2.
By [Kea99], A is isomorphic to a direct product
∏i=1tBi of nilpotent algebras of prime power order.
We let (π1(a),…,πt(a)) denote
the image of a of the underlying isomorphism.
As a finite supernilpotent algebra in a congruence
modular variety, A is nilpotent (cf. [Aic18, Lemma 2.4]) and
therefore has a Mal’cev term [FM87, Theorem 6.2].
We use
[Aic18, Theorem 4.2] to expand each Bi with operations
+i and −i such that
the expansion Bi′ is a nilpotent expanded elementary abelian group with
zero element πi(z). By [Aic18, Theorem 1.2],
Bi′ is ki-supernilpotent. ∎
We note that the supernilpotency degree of A′ may be strictly larger
than the supernilpotency degree of A.
Combining these results, we obtain the following result on polynomial mappings
on arbitrary finite supernilpotent algebras in congruence modular
varieties.
Theorem 10**.**
Let μ∈N,
let A be a finite supernilpotent algebra in a congruence
modular variety all of whose fundamental operations have arity at most μ.
Let p1,…,pt be distinct primes, and let α1,…,αt∈N such
that ∣A∣=∏i=1tpiαi.
Let F∈Poln,s(A) be a polynomial
map from An to As, and let z∈A.
Then for every \mathbfsla∈An there is \mathbfsly∈An
such that F(\mathbfsly)=F(\mathbfsla) and
∣{j∈n:\mathbfsly(j)=z}∣≤s∑i=1t(μ(piαi−1))αi−1αi(pi−1)≤sμ−1∣A∣log2(μ)+log2(∣A∣)log2(∣A∣)≤s∣A∣log2(μ)+log2(∣A∣)+1.
Proof.
Let A′=∏i=1tBi′ be the expansion of A produced by
Lemma 9, and for each i∈t, let pi∈P and αi∈N be such that
∣Bi∣=piαi.
Clearly, F is a also a polynomial map of
A′.
Let ki=(μ(piαi−1))αi−1.
Then Theorem 8 yields \mathbfsly∈An such that
∣{j∈n:\mathbfsly(j)=z}∣≤∑i=1tkisαi(pi−1).
Using the obvious estimate αi≤log2(∣A∣), we obtain
[TABLE]
∎
5. Systems of equations
We will now explain how these results give a polynomial time algorithm
for solving systems of a fixed number of equations over the finite
supernilpotent algebra A. The size m of a system of
polynomial equations is measured
as the length of the polynomial terms used to represent the system.
For measuring the “running time” of our algorithm, we count
the number of A-operations: each such A-operation, may,
for example, be done by looking up one value in the operation tables
defining A.
Theorem 11**.**
Let A be a finite supernilpotent algebra in a congruence modular
variety all of
whose fundamental operations are of arity at most μ,
and let s∈N. We consider the following algorithmic problem
s-\textscPolSysSat(A):
Given:* 2s polynomial terms f1,g1,…,fs,gs over A. *
Asked:* Does the system f1≈g1,…,fs≈gs have a solution
in A?*
Let m be the length of the input of this system, and
let
[TABLE]
Then we can decide s-\textscPolSysSat(A) using at most
O(me−1) evaluations of all terms occuring in the system.
Therefore, we have an algorithm
that determines whether a system of s polynomial equations over A
has a solution using O(me) many A-operations.
Proof.
Let n be the number of different variables that occur in the
given system. We may assume that these variables
are x1,…,xn, and that our system is
⋀i=1sfi(x1,…,xn)≈gi(x1,…,xn).
We choose an element z∈A, and we
will show: if this system has a solution in \mathbfsla∈An,
then it has a solution in
[TABLE]
For proving this claim, we first observe that
A is a finite nilpotent algebra in a congruence modular variety,
and it therefore has a Mal’cev term d.
We consider the polynomial map H=(h1,…,hs), where
hi(\mathbfslx):=d(fi(\mathbfslx),gi(\mathbfslx),z) for i∈s and \mathbfslx∈An.
Since \mathbfsla is a solution of the system, H(\mathbfsla)=(z,z,…,z).
By Theorem 10, there is \mathbfsly∈C such that
H(\mathbfsly)=H(\mathbfsla).
Then for every i∈s, we have
d(fi(\mathbfsly),gi(\mathbfsly),z)=z.
By [FM87, Corollary 7.4], the function x↦d(x,gi(\mathbfsly),z) is injective.
Since d(fi(\mathbfsly),gi(\mathbfsly),z)=z=d(gi(\mathbfsly),gi(\mathbfsly),z), this injectivity implies that fi(\mathbfsly)=gi(\mathbfsly).
Hence \mathbfsly is a solution that lies in C.
The algorithm for solving the system now simply evaluates the system at all places in
C; if a solution is found, the answer is “yes”. If we find no solution inside C , we answer “no”,
and by the argument above, we know that in this case, the system has no solution inside An at all.
We now estimate the complexity of this procedure:
There is a c∈N such that for all n∈N, ∣C∣≤c1ne−1,
hence we have to do O(ne−1) evaluations of all the terms fi,gi in the system.
Such an evaluation can be done using at most O(m) many A-operations.
Since the length of the input m is at least the number of variables n occuring in it, this solves s-\textscPolSysSat(A) using at most
O(me) many A-operations.
∎
6. Circuit satisfiability
With every finite algebra A,
[IK18] associates a number of computational problems that
involve circuits whose gates are taken from the fundamental operations of
A. One of these problems is \textscSCsat(A). It takes
as an input 2s circuits f1,g1,…,fs,gs over A with n input
variables, and asks whether there is a \mathbfsla∈An such that
the evaluations at \mathbfsla satisfy fi(\mathbfsla)=gi(\mathbfsla) for
all i∈s. For finite algebras in congruence modular varieties,
[LZ06, Corollary 3.13] implies that
\textscSCsat(A) is in P when A is abelian,
and NP-complete otherwise.
However, if we restrict the number s of circuits, we obtain a different
problem, which we call s-\textscSCsat(A) in the sequel.
Obviously, 1-\textscSCsat(A) is the circuit satisfiability
problem called \textscCsat(A) in [IK18].
The method used to prove Theorem 11 immediately
yields:
Theorem 12**.**
Let A be a finite supernilpotent algebra in a congruence
modular variety, and let s∈N. Then
s-\textscSCsat(A) is in P.
Hence a supernilpotent, but not abelian algebra A has
s-\textscSCsat(A) in P, whereas
SCsat is NP-complete.
In the converse direction, Theorem 9.1 from [IK18]
has the following corollary.
Corollary 13**.**
Let A be a finite algebra from a congruence modular variety.
If A has no homomorphic image A′ such that
2-\textscSCsat(A′) is NP-complete, then
A is nilpotent.
Proof.
Suppose that A has a homomorphic image A′ for
which \textscCsat(A′) is NP-complete.
Then also 2-\textscSCsat(A′) is NP-complete
because an algorithm solving 2-SCsat can be used
to solve an instance (∃\mathbfsla)(f(\mathbfsla)=g(\mathbfsla))
of \textscCsat(A′) by solving 2-SCsat on the input
(∃\mathbfsla)(f(\mathbfsla)=g(\mathbfsla)&f(\mathbfsla)=g(\mathbfsla)).
Thus the assumptions imply that for no homomorphic image A′ of A,
the problem \textscCsat(A′) is NP-complete.
Now by [IK18, Theorem 9.1], A is isomorphic to
N×D, where N is nilpotent and D
is a subdirect product of 2-element algebras each of which is
polynomially equivalent to a two element lattice.
If ∣D∣>1, then there is a homomorphic image
A2 of A such that A2 is polynomially equivalent to
a two element lattice. By [GK11], 2-\textscSCsat(A2)
is NP-complete, contradicting the assumptions.
Hence ∣D∣=1, and therefore A is nilpotent. ∎
Acknowledgements
The author thanks M. Kompatscher for dicussions on solving equations over
nilpotent algebras.
These discussions took place during a workshop organized by P. Aglianò at the
University of Siena in June 2018.
The author also thanks A. Földvári, C. Szabó, M. Kompatscher, and S. Kreinecker
for their comments on preliminary versions of the manuscript.