Abstract 1 Introduction 2 Preliminaries 3 Proof of the Main Result 4 Concluding Remarks References

On Expansions of Monadic Second-Order Logic with Dynamical Predicates

Joris Nieuwveld ORCID Max Planck Institute for Software Systems, Saarland Informatics Campus, Saarbrücken, Germany Joël Ouaknine ORCID Max Planck Institute for Software Systems, Saarland Informatics Campus, Saarbrücken, Germany
Abstract

Expansions of the monadic second-order (MSO) theory of the structure ;< have been a fertile and active area of research ever since the publication of the seminal papers of Büchi and Elgot & Rabin on the subject in the 1960s. In the present paper, we establish decidability of the MSO theory of ;<,P, where P ranges over a large class of unary “dynamical” predicates, i.e., sets of non-negative values assumed by certain integer linear recurrence sequences. One of our key technical tools is the novel concept of (effective) prodisjunctivity, which we expect may also find independent applications further afield.

Keywords and phrases:
Monadic second-order logic, linear recurrence sequences, decidability, Baker’s theorem
Funding:
Joël Ouaknine: Also affiliated with Keble College, Oxford as emmy.network Fellow, and supported by ERC grant DynAMiCs (101167561) and DFG grant 389792660 as part of TRR 248.
Copyright and License:
[Uncaptioned image] © Joris Nieuwveld and Joël Ouaknine; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Logic and verification
; Theory of computation Formal languages and automata theory ; Mathematics of computing Discrete mathematics
Related Version:
Full Version: https://arxiv.org/abs/2507.16581
Editors:
Paweł Gawrychowski, Filip Mazowiecki, and Michał Skrzypczak

1 Introduction

The monadic second-order (MSO) theory of the structure ;< has been a foundational pillar of the field of automated verification, and more generally the area of logic in computer science, for many decades. Arguably the most important paper on the topic is due to Büchi [6], who in the early 1960s established decidability of this theory through a deep and fecund connection between logic and automata theory.

Shortly thereafter, in yet another seminal piece of work [11], Elgot and Rabin devised the contraction method to establish decidability of expansions of this base theory by various “arithmetic” unary predicates P.111In this paper, we adopt the convention that the set of natural numbers contains 0. We also use the adjective “positive” with the meaning of “non-negative”. In particular, they proved decidability of the MSO theory of ;<,P, where P could for example be taken to be the set 2 of powers of 2, or the set 𝖲𝗊 of perfect squares, or the set 𝖥𝖺𝖼 of factorial numbers, and so on.

Much progress followed in the ensuing decades, notably thanks to Semënov [24], who introduced the notion of effective almost periodicity, and Carton and Thomas [8], who substantially refined Elgot and Rabin’s contraction method into the notion of effective profinite ultimate periodicity. Other notable works in this area include articles by Rabinovich [22], Rabinovich and Thomas [23], and Berthé et al. [3].

In the present paper, we significantly extend this line of research by considering a large class of “dynamical” predicates derived from integer linear recurrence sequences (LRS).222The “dynamical” moniker was chosen to reflect the close kinship of these predicates with discrete-time linear dynamical systems; see, e.g., [2, 9, 16, 15, 18, 14, 1]. The complexity of an LRS can be measured in various ways; chief among them are the number of distinct dominant characteristic roots of the LRS, and whether the LRS is simple or not.333These notions are properly defined in Sec. 2.3. For P the set of positive values of an LRS having a single dominant root, such as the set of Fibonacci numbers, the decidability of the MSO theory of the structure ;<,P is readily established via Elgot and Rabin’s contraction method;444LRS having a single positive dominant root are particularly well behaved: they are either constant, or effectively ultimately monotonically increasing (or decreasing), effectively procyclic, and effectively sparse (see [3] for precise definitions of these notions and a proper account of these facts). In the case of a negative dominant root, essentially the same behaviour is observed by restricting to the positive terms. see also [3] in which expansions of ;< by adjoining multiple predicates obtained from such LRS are thoroughly investigated. Unfortunately, virtually nothing is known when considering LRS having more than a single dominant characteristic root.

Our main contribution is the following result:

Theorem 1.

Let P={un:n} for unn=0 a non-degenerate, simple, integer-valued linear recurrence sequence with two dominant roots. Then the MSO theory of the structure ;<,P is decidable.

An example of an LRS satisfying the hypotheses of Thm. 1 is the sequence unn=0 defined by

un+3=6un+213un+1+10un (1)

and u0=2,u1=4, and u2=7. Its subsequent values are: 10, 9, 6, 53, 150, 271, 206, 787, 4690, 15849, 41994, 92827, 169530, 230369, 106594, …. One can readily verify that unn=0 satisfies the formula

un=12(2+i)n+12(2i)n+2n, (2)

and that un is both infinitely often positive and infinitely often negative. Moreover, although limn|un|=+, |un| is not monotonically increasing: for all N there is an n such that |un+N|<|un|.

For our LRS unn=0, let P be as per Thm. 1. Viewing P as an ordered set, we have:

P={u0,u1,u2,u4,u3,u10,u11,u12,u13,u14,u17,u15,u16,u24,u25,u26,u27,u28,u30,}.

One immediately notices that there are two complications at play here. The first one is that we are “throwing away” all the negative values of our LRS (this is of course necessary since we are working over domain rather than ). This restriction is, however, entirely benign in view of the following result:

Corollary 2.

Let P¯={un:n} for unn=0 a non-degenerate, simple, integer-valued linear recurrence sequence with two dominant roots. Then the MSO theory of the structure ;0,<,P¯ is decidable.

This result is straightforwardly obtained from Thm. 1 via an application of Shelah’s celebrated “composition method” in model theory [25]. In the case at hand, one can directly invoke, for example, [26, Cor. 6], since the structure ;0,<,P¯ is isomorphic to the ordered sum ;<,P+;<,P, where the second summand is as per Thm. 1, and the first structure is obtained by setting P={un:n}(). Intuitively speaking, MSO sentences over ;<,P¯ can be faithfully decomposed into component subformulas dealing exclusively with either positive or negative values of our LRS, with the truth value of the original sentence obtained by appropriately piecing together truth values of each of the sub-sentences within their respective structures.

The second complication is that, as noted earlier, the ordering of the positive values of our LRS does not respect the index ordering of the LRS. This ostensibly precludes the direct application of classical techniques such as Elgot and Rabin’s contraction method (or more generally Carton and Thomas’s effective profinite ultimate periodicity criterion), or Semënov’s toolbox of effective almost periodicity.

In order to prove Thm. 1, we therefore rely instead on a new concept, that of (effective) prodisjunctivity, which we introduce shortly. Prodisjunctivity is itself predicated on the notion of disjunctivity, whose application to the decidability of logical theories (under the alternative terminology of “weak normality”) was pioneered by Berthé et al. [3]. In particular, Berthé et al. study the MSO theory of the structure ;<,2,𝖲𝗊, and establish decidability assuming that the binary expansion of 2 is disjunctive, i.e., contains every finite bit pattern as a factor infinitely often.

The hypothesis that 2 is disjunctive in base 2 is widely expected by number theorists to be true, but remains a major unsolved problem. Irrational algebraic numbers, along with most known transcendental numbers such as e and π, are in fact believed to satisfy a stronger property, that of being normal in every integer base: a real number α is normal in base b provided that every finite word w{0,,b1} appears with frequency b|w| in the base-b expansion of α. And although Borel [5] showed over a century ago that non-normal real numbers have null Lebesgue measure, establishing normality (or even disjunctivity) of everyday irrational numbers has remained fiercely elusive; see [13, 21, 17, 7] for more detailed accounts of results, research, and open problems in this area.

Given a finite alphabet Σ together with a subset SΣ, we say that an infinite sequence in Σω is disjunctive relative to S if every finite word over S appears infinitely often as a factor in the sequence. We are now in a position to define (effective) prodisjunctivity:

Definition 3.

Let pmm=0 be an infinite integer-valued sequence. For any integer M2, let SM be the set of residue classes modulo M that appear infinitely often in pmm=0:

SM={s{0,,M1}:m:pms(modM)}. (3)

We say that the sequence pmn=0 is prodisjunctive if, for all M2, the sequence of residues pmmodMm=0 is disjunctive relative to SM.

For effectiveness, we require in addition that the sequence pmmodMm=0 be computable and that, for each M, the set SM, together with the smallest index threshold NM beyond which all residue classes of the sequence lie in SM (i.e., for all mNM, (pmmodM)SM), also be computable.

Let us now sketch how prodisjunctivity relates to Thm. 1. Recall from (1) our LRS unn=0, along with its infinite set P of positive values (as per the notation of Thm. 1), and let pmm=0 be a strictly increasing enumeration of P; in other words, pm1 is the mth smallest element in P. We will establish the following instrumental result:

Theorem 4.

Let pmm=0 be as above, i.e., the strictly increasing sequence of positive terms of some non-degenerate, simple, integer-valued LRS having two dominant roots. Then pmm=0 is effectively prodisjunctive.

Let us return to our example and set M=5. Then one easily shows that

unmod5={0if n3(mod4)2if n=0, or if n2(mod4)4otherwise. (4)

It follows that S5={0,2,4}, N5=0, and

pmmod5m=0=2,4,2,4,0,2,0,4,4,2,4,0,4,4,4,2,0,4,2,4,2,0,4,4,4,2,0,0,4,4,2,

is in {0,2,4}ω.

Computing the first million terms of pmmod5m=0, one does not encounter the factor 0,0,0 once within that initial segment. Nevertheless, according to Thm. 4, it should appear infinitely often! Indeed, we prove this in Sec. 3.2, and moreover construct an index of roughly 2.181059 where this occurs. A more involved calculation (not shown here) can however establish that the very first occurrence of 0,0,0 within pmmod5m=0 is at index 8479226.

The above discussion suggests that, whilst the sequence pmm=0 is provably disjunctive, it is by all accounts not normal, i.e., a given factor w{0,2,4} does not necessarily appear with frequency 3|w|. This is however unsurprising in view of (4): the residue class 4 should statistically appear approximately twice as often as either of the other two residue classes, and indeed it is possible to prove that this is asymptotically the case.

Theorem 4 is the key technical device underpinning our main result, Thm. 1. One of the critical mathematical ingredients entering its proof is Baker’s theorem on linear forms in logarithms of algebraic numbers. We also make use of various automata-theoretic, topological, and combinatorial constructions and tools.

Let us now briefly comment on the motivation and scope of Thm. 1. As noted earlier, linear recurrence sequences and linear dynamical systems are intimately related; in particular, decision procedures for the former often underpin verification algorithms for the latter. Write P𝒖 to denote the set of positive values of a given integer LRS 𝒖. We already pointed out that Elgot and Rabin’s contraction method enables one to establish decidability of the MSO theory of ;<,P𝒖 provided that 𝒖 has a single dominant root, and that [3, Thm. 5.1] provides general conditions under which the MSO theory of ;<,P𝒖(1),,P𝒖(d) is decidable, once again under the assumption that each LRS 𝒖(i) possesses a single dominant root. On the other hand, the situation concerning predicates derived from LRS having more than one dominant root is entirely uncharted. The present paper can thus be seen as taking the first significant step (in some six decades!) in this direction.

Nevertheless, the hypotheses of Thm. 1 (the LRS must be non-degenerate, simple, and have exactly two dominant roots) may at first glance appear quite restrictive. It is worth noting, however, that if one were to pick LRS “at random”, then under most “reasonable” probability distributions, with probability tending to 1 such LRS would have either a single dominant root, or be non-degenerate, simple, and have exactly two dominant roots. In other words, the vast majority of LRS are of one of these two types, and thus fall under the scope of the Elgot-Rabin contraction method or our own Thm. 1. (This folklore result follows from the fact that a random polynomial with real coefficients will almost surely have either a single dominant real root, or two dominant complex conjugate roots; and that non-degeneracy and simplicity are null-measure conditions. For a more in-depth discussion of these matters, see the work of Dubickas and Sha [10] and citations therein.) In any event, we briefly return in Sec. 4 to the question of whether the hypotheses of Thm. 1 could be weakened in any way.

2 Preliminaries

2.1 Automata Theory

An alphabet Σ is a finite, non-empty set of letters. Σ and Σω denote respectively the sets of finite and infinite words over Σ. As words can be viewed as sequences, we freely identify the finite word w0w1wk with the finite sequence w0,w1,,wk, and the infinite word w0w1 with the infinite sequence w0,w1,, and conversely. A finite word w=w0w1wk is a factor of an (in)finite word w if there is an index n such that wn,,wn+k=w0,,wk. We let |w| denote the length a finite word w.

An infinite word w is recursive (or computable) if one can compute wj for every j0, and is disjunctive if each wΣ appears infinitely often as a factor of w. Other authors occasionally use the terminology weakly normal or rich instead of disjunctive. An infinite word w is periodic if w=vω for some vΣ+ and ultimately periodic if w=uvω for some u,vΣ. Here, for the smallest possible such u and v, |u| and |v| are referred to as the preperiod and period, respectively.

A deterministic Muller automaton 𝒜=(Σ,Q,qinit,δ,) consists of an alphabet Σ, a finite set of states Q, an initial state qinitQ, a transition function δ:Q×ΣQ, and an accepting family of sets 2Q. The run of wΣΣω on 𝒜 is the sequence of states visited while reading w starting in qinit and repeatedly updating the state using the transition function. We say that 𝒜 accepts wΣω if the set of states visited infinitely often upon reading w belongs to . The acceptance problem of a recursive word wΣω is the question of determining, given a deterministic Muller automaton 𝒜 with alphabet Σ, whether 𝒜 accepts w. We denote the acceptance problem by Accw.

Berthé et al. established the following [3, Thm. 4.16]:

Theorem 5.

If wΣω is recursive and disjunctive, then Accw is decidable.

A deterministic finite transducer =(Σin,Σout,Q,qinit,δ) is given by an input alphabet Σin, output alphabet Σout, set of states Q, initial state qinitQ, and transition function δ:Q×ΣinQ×Σout. The transducer starts in state qinit. It reads a word wΣinΣinω and upon reading letter a whilst in state q, it computes (q,w)=δ(q,a), moves to state q, and concatenates w to the output string. Write (w)ΣoutΣoutω to denote the output word thus computed upon reading wΣinΣinω.

We now recall a Turing-reducibility property between word acceptance problems [3, Lem. 4.5]:

Lemma 6.

Let wΣω and be a deterministic finite transducer. Then the problem Acc(w) reduces to Accw.

2.2 Monadic Second-Order Logic

A (unary or monadic) predicate P is a function P:{0,1}, which equivalently we can interpret as a subset P. For P an infinite predicate, let us write pmm=0 to denote the enumeration of P in increasing order, i.e., such that pm1 is the mth smallest element of P. The characteristic word w{0,1}ω of P is obtained by setting wn=P(n). We then have:

w=0p010p1p0110p2p1110p3p21. (5)

A predicate P is sparse if for every N0, there is M0 such that pm+1pmN for all mM. The predicate is effectively sparse if M can always be computed from N.

Monadic second-order logic (MSO) is an extension of first-order logic over signature {=,<,} that allows quantification over subsets of the universe . We also consider expansions of MSO by various fixed unary predicates P; in other words (abusing notation), the signature is expanded by a predicate symbol P, with interpretation the given subset P. We refer the reader to [4] for a thorough contemporary overview of MSO.

The deep connection between MSO and automata theory was brought to light in the seminal work of Büchi; see, for example, [27, Thms. 5.4 and 5.9].

Theorem 7.

The decidability of the MSO theory of the structure ;<,P is Turing equivalent to Accw, where w is the characteristic word of P.

As noted earlier, Elgot and Rabin devised the contraction method to establish decidability of the MSO theory of ;<,P, for various “arithmetic” predicates P. The following proposition is a variation on their method:

Proposition 8.

Let P be an infinite, recursive, and effectively sparse predicate with enumeration pmm=0. If, for each M2, AccpmmodMm=0 is decidable, the MSO theory of the structure ;<,P is decidable.

Proof.

By Thm. 7, it is sufficient to be able to decide whether a deterministic Muller automaton 𝒜=({0,1},Q,qinit,δ,) accepts the characteristic word (5). Let us restrict 𝒜 to a directed graph G with nodes Q and 0-transitions as arrows.

By construction, every node in G has outdegree 1 and so each state is in at most one cycle in G. We can therefore compute the least common multiple of the cycle lengths in G (call this number M) and the longest path to a cycle (call this number N). Then, for all states q and numbers nM+N and d1, reading 0n and 0n+dM leads to journeying through the exact same set of states and ending up in the same state.

Let K be such that pm+1pmM+N+1 for all mK (which can be computed as P is effectively sparse). We construct a deterministic finite transducer that hard-codes the initial segment winit:=0p010p1p0110pK+1pK11. For mK, after reading (pmmodM) and (pm+1modM), outputs 0km1, where M+N<km2M+N is congruent to pm+1pm1 modulo M. Then, by construction, a state q is visited infinitely often upon reading the characteristic word of P if and only if q is visited infinitely often when 𝒜 reads winit(pmmodMm=K+1).

Recall that AccpmmodMm=0 is assumed to be decidable. Then AccpmmodM)m=K+1 is also decidable (by hard-coding the initial segment) and thus Acc(pmmodMm=K+1) is decidable by Thm. 6. Therefore Accwinit(pmmodMm=K+1) is decidable (by again hard-coding the initial segment), and by construction, 𝒜 accepts the characteristic word of P if and only if 𝒜 accepts winit(pmmodMm=K+1). Hence Accw is decidable, as required.

2.3 Linear Recurrence Sequences

A number α is algebraic if F(α)=0 for some non-zero polynomial F[X]. The unique such polynomial (up to multiplication by a constant) of least degree is the minimal polynomial of α. We write ¯ to denote the field of algebraic numbers.

A linear recurrence sequence over a ring R (an R-LRS) is a sequence unn=0Rω such that there are numbers c1,,cdR, with cd0, having the property that, for all n,

un+d=c1un+d1++cdun.

The entire sequence is therefore completely specified by the 2d numbers c1,,cd and u0,,ud1. Although there may be several such recurrence relations, we shall assume that we are always dealing with the (unique) one for which d is minimal. In the remainder of the paper, whenever R is not specified, we are working over the ring of integers . The characteristic polynomial of the LRS is given by F(X)=Xdc1Xd1cdR[X]. The characteristic roots of the LRS are the roots of its characteristic polynomial, and the multiplicity of a characteristic root λ is its multiplicity as a root of F(X). Every ¯-LRS unn=0 with characteristic roots λ1,,λK admits a unique exponential-polynomial representation given by

un=k=1KQk(n)λkn,

where Qk¯[X] has degree the multiplicity of λk minus 1. A ¯-LRS is simple if all its characteristic roots have multiplicity 1, and is non-degenerate if no quotient of two distinct characteristic roots is a root of unity.

A characteristic root λ is dominant if |λ||μ| for all characteristic roots μ. We define the dominant part vnn=0 of unn=0 by writing

vn=λj dominantQj(n)λjn,

with the non-dominant part rnn=0 given by unvnn=0. One can compute positive real numbers r and R such that rRn>|rn| for all n and R<|λ| for λ a dominant characteristic root, see, e.g., [3, Lem. 2.5].

For every M2 and -LRS unn=0, the sequence unmodMn=0 is ultimately periodic and one can effectively compute its period and preperiod [3, Lem. 2.6].

Example 9.

Let unn=0 be the (-)LRS (1) from Sec. 1. The characteristic polynomial of unn=0 is

F(X)=X36X2+13X10=(X2)(X(2+i))(X(2i)).

The characteristic roots of unn=0 are 2, 2+i, and 2i. Using linear algebra, one easily recovers the exponential-polynomial representation (2). The LRS unn=0 is simple and non-degenerate. Its dominant part is the sequence defined by vn=12(2+i)n+12(2i)n, while its non-dominant part is given by rn=2n.

2.4 Number Theory

Our main number-theoretic tool is Baker’s theorem on linear forms in logarithms. Specifically, we make use of a flexible version due to Matveev along with some off-the-shelf applications due to Mignotte, Shorey, and Tijdeman.

Let α0 be an algebraic number of degree d with minimal polynomial F(X)=a0i=1d(Xαi). The logarithmic Weil height of α is defined as

h(α)=1d(log|a0|+i=1dlogmax{|αi|,1}).

Furthermore, set h(0)=0. For all algebraic numbers α1,,αk and n, we have the following properties:

h(α1++αk) logk+h(α1)++h(αk),
h(α1α2) h(α1)+h(α2),and
h(α1n) |n|h(α1).

A number field L is a field extension of such that the degree of L/ is finite. If α1,,αd¯, L=(α1,,αd) is a number field whose degree can be effectively computed.

Let L be a number field of degree D, M1, α1,,αML, and b1,,bM. Then set B=max{|b1|,,|bM|},

Λ=j=1Mαjbj1,andh(αj)=max{Dh(αj),|logαj|,0.16}.

Matveev [19] proved the following:

Theorem 10.

If Λ0, then

log|Λ|>330M+4(M+1)5.5D2(1+logD)(1+log(MB))j=1dh(αj).

In particular, there is a computable constant c such that when Λ0, |Λ|>Bc.

We will also need the following results from Mignotte, Shorey, and Tijdeman [28]:

Theorem 11.

Let unn=0 be a non-degenerate LRS with two dominant roots of magnitude |λ|. Then there are computable positive constants C1 and C2 such that

|un||λ|nnC1log(n)

whenever nC2.

Theorem 12.

Let unn=0 be a non-degenerate LRS with two dominant roots of magnitude |λ|. Then there are computable positive constants C3 and C4 such that

|un1un2||λ|n1n1C3log(n1)log(n2+2)

whenever n1>n2 and n1C4.

3 Proof of the Main Result

In this section, we prove Thm. 1: the MSO theory of the structure ;<,P is decidable whenever P is a predicate comprising the set of positive terms of some non-degenerate, simple, integer-valued LRS having two dominant characteristic roots.

We begin in Sec. 3.1 by untangling the definition of an LRS satisfying the above hypotheses and reduce Thm. 1 to Thm. 4. We then provide intuition underlying the proof of the latter through an extended example in Sec. 3.2, and consider a continuous version of our problem in Sec. 3.3. Finally, in the same section, we establish Thm. 4.

3.1 Reduction to Prodisjunctivity

Let unn=0 be an LRS satisfying the hypotheses of Thm. 1. We first record some elementary observations whose proof is in the full version.

Lemma 13.

Assume that unn=0 is an LRS satisfying the hypotheses of Thm. 1 and whose two dominant roots are λ1 and λ2. Then λ2=λ1¯, the argument of λ1 is not a rational multiple of π, and |λ1|>1.

When writing unn=0 in its exponential-polynomial form, we split it into its dominant part vnn=0 and non-dominant part rnn=0:

un=vn+rn=αλn+α¯λ¯n+rn. (6)

Here, α and λ are algebraic numbers with α0 (as otherwise λ and λ¯ would not be characteristic roots, i.e., roots of the polynomial corresponding to the minimal recurrence relation that unn=0 obeys), |λ|>1, and the argument of λ is not a rational multiple of π. We keep λ, α, and rnn=0 fixed for the remainder of the paper.

Recall that P={un:n} and pmm=0 is an enumeration of P in increasing order. To apply Prop. 8, we need the following lemma, whose proof relies heavily on the results of Mignotte, Shorey, and Tijdeman from Sec. 2.3 and can be found in the full version.

Lemma 14.

Let P be as per Thm. 1. Then P is infinite, recursive, and effectively sparse.

By Prop. 8, it now suffices to prove that for all M2, one can decide AccpmmodMm=0 (determine whether a given deterministic Muller automaton 𝒜 over alphabet {0,,M1} accepts pmmodMm=0). The next lemma shows how the effective prodisjunctivity of pmm=0 asserted by Thm. 4 enables us to do this. Its proof can be found in the full version.

Lemma 15.

Theorem 4 implies Thm. 1, i.e., if unn=0 is an LRS satisfying the conditions of Thm. 1, P:={un:n}, and the enumeration of P is prodisjunctive, then the MSO theory of the structure ;<,P is decidable.

Theorem 4 asserts that, for any M2, the sequence pmmodMm=NMSMω is disjunctive. Let us unpack this definition. We have that pmmodMm=NMSMω is disjunctive if for any 1, every pattern s1,,sSM appears infinitely often in pmmodMm=NM. That is, for all ,N with NNM and s1,,sSM, there are n1,,n such that

  1. 1.

    n1,,nN;

  2. 2.

    for all 1i, unisi(modM);

  3. 3.

    0un1<<un;

  4. 4.

    for all m0 such that un1umun, um{un1,,un}.

As the dominant part vnn=0 (where vn=αλn+α¯λ¯n, see (6)) only relies on two algebraic numbers, α and λ, it is easier to work with vnn=0 than unn=0. In the following lemma, we prove that we can actually effect this change.

Lemma 16.

Assume that for all natural numbers , N, T2, and t1,,t{0,,T1}, there are n1,,n such that

  1. 1.

    n1,,nN;

  2. 2.

    for all 1j, njtj(modT);

  3. 3.

    0<vn1<<vn;

  4. 4.

    for all m0 such that vn1<vm<vn, m{n2,,n1}.

Then the enumeration pmm=0 of P:={un:n} is effectively prodisjunctive.

Proof.

We claim that for some computable number N, um1<um2 if and only if vm1<vm2 whenever m1,m2N. Suppose not. Then, without loss of generality, m1>m2 and |vm1vm2|<|rm1|+|rm2|<2rRm1. But Thm. 12 implies that for m1C4,

log(2r)+m1logR >log|vm1vm2|
>m1log|λ|C3log(m1)2log(m2+2)
>m1log|λ|C3log(m1+1)3,

which cannot hold for m1N for some computable N as log|λ|>log|R|. Our claim therefore follows.

Assume that unmodMn=0 has period T (which can be effectively computed). If s1,,sS and 1j, there exists a tj{0,,T1} such that unT+tjsj(modM) whenever n is large enough. Therefore, if n1,,n are at least max{N,N} and satisfy the hypotheses of the lemma, it follows that 0un1<<un and for all m such that un1<um<un, m{n2,,n1}. In other words, if pm=un1, then for 2j, pm+j1=unj and pm+j1=unT+tjsj(modM) for some n.

As |α|0 and Lem. 16 is only concerned with inequalities vn1<vn2 and vn1>0 for natural numbers n1 and n2, we can scale vn by 1/|2α|. That is, we can assume that |α|=1/2. Write α=12eiϕ and λ=|λ|eiθ. Thus, vn=cos(θn+ϕ)|λ|n.

Let us now sketch the method we will use to establish the hypothesis of Lem. 16 whose proof relies heavily on the following notation.

Definition 17.

For integers d0 and real numbers 0<γ<δ, define 𝒥d(γ,δ)/(2π) as

𝒥d(γ,δ)={x/(2π):0<γcos(x)<cos(x+dθ)|λ|d<δcos(x)}.

Moreover, we define Jd(0,δ) to stand for the limit of Jd(γ,δ) as γ tends to 0.

In Lem. 18, we can quantify the size of these intervals. For a real number x, let |x|2π denote the distance of x to the nearest integer multiple of 2π. The proof of Lem. 18 can be found in the full version.

Lemma 18.

One can compute constants C5, C6, C7, and C8 such that for all 0<γ<δ<|λ| and all d1, 𝒥d(γ,δ) consists of a single interval,

C6δγ|λ|ddC7<|𝒥d(γ,δ)|<C5δγ|λ|d,

and |x(dθ±π/2)|2π<C8|λ|d for any x𝒥d(γ,δ).

We now continue with the most technical step in the proof of Thm. 1: establishing a continuous version of the hypothesis of Lem. 16. Items 3 and 4 of Lem. 16 involve terms of the form vn, and we can divide them by |λ|n1. Further, for 2j, we set bj:=njn1 and x=nθ+ϕ. Then, item 3 turns into

0<cos(x)<cos(x+b2θ)|λ|b2<<cos(x+bθ)|λ|b (7)

and item 4 turns into

d:cos(x)<cos(x+dθ)|λ|d<cos(bθ)|λ|bd{b2,,b1}. (8)

We want to find an interval /(2π) and numbers bj such that bjtjt1(modT) (for item 2), (7) holds for all x, and (8) hold for “most” x. This leads us to the following:

Lemma 19.

Let ,T2 and t1,,t{0,,T1}. Then there are an interval /(2π), natural numbers D,b2,,b1, and real numbers 1<δ2<<δ<|λ| such that

  1. (a)

    for all 2j, bjtjt1(modT);

  2. (b)

    for all x,

    0<cos(x)<cos(x+b2θ)|λ|b2<δ2cos(x)<cos(x+b3θ)|λ|b3<δ3cos(x)<cos(x+bθ)|λ|b<δcos(x); (9)
  3. (c)

    for all integers d<D not in the set {b2,,b1}, 𝒥d(1,δ)=;

  4. (d)

    d=D|𝒥d(1,δ)|<||.

Translating (b) in the notation of intervals 𝒥, we have that

(π/2,π/2)𝒥b2(1,δ2)𝒥b3(δ2,δ3)𝒥b3(δ1,δ).

The proof of Lem. 19 is in the full version. The remaining tools we need for this proof are listed in Sec. 3.3. Finally, we derive Thm. 1 from Lem. 19. But first, in the upcoming section we go through an example to recap the construction.

3.2 An Extended Example

Let unn=0 be the sequence (1) from Sec. 1 and assume M=5. Let us study the sequence un=12(2+i)n+12(2i)n+2n modulo 5.

From (4), we conclude that S5={0,2,4} and T=4 (that is, unmod5n=0 is ultimately periodic with period T=4, and {0,2,4} are exactly the congruence classes that appear infinitely often in unmod5n=0). Thus for all m0, pmmod5{0,2,4}, where pmm=0 is the enumeration of P:={un:n}. Then Thm. 4 asserts that every s1,,sS5 appears in pmmod5m=0 infinitely often as a factor. We will show that this indeed holds when =3 and s1,s2,s3=0,0,0.

Let us compute t1, t2, and t3 (the conjugacy classes modulo T=4 such that for large enough n, uTn+tisimod5). In view of (4), we are forced to take ti=3 for i=1,2,3 as si=0. Thus, to find 0,0,0 in pmm=0, we want to find n1,n2,n3 congruent to 3 modulo 4 such that 0<un1<un2<un3, and if un1<um<un3 for some m0, then m{n2} (i.e., m=n2).

By (2), the dominant part vnn=0 of unn=0 is given by vn=12(2+i)n+12(2i)n and the non-dominant part rnn=0 by rn=2n. As shown in Lem. 16, we can work with vn=αλn+α¯λ¯n=cos(nθ+ϕ)|λ|n instead of un for large enough n. Here, we have that λ=eiθ|λ|=2+i and α=12eiϕ=12. Thus, in this example, ϕ=0.

Hence, for a given N (for simplicity, let us take N=0) we wish to find n1,n2,n3N that are congruent to 3 modulo 4, satisfy

0<cos(n1θ)<cos(n2θ)|λ|n2n1<cos(n3θ)|λ|n3n1,

and such that m=n2 for any m satisfying the inequality cos(n1θ)<cos(mθ)|λ|mn1<cos(n3θ)|λ|n3n1. This is the statement of Lem. 16.

Next, we aim to find b2,b3 such that n1=n, b2=n2n, and b3=n3n for some nN congruent to 3 modulo 4 that satisfies our hypotheses. To ensure that for j{2,3}, njtj3(mod4), we have to require that bj0(mod4).

In order to solve this discrete problem (find natural numbers n1,n2,n3 meeting these constraints), we first want to solve a continuous variant of this problem: find b2 and b3 and an interval /(2π) such that for x, these properties hold “often”. We will construct an open, non-empty interval /(2π) such that for all x,

0<cos(x)<cos(x+b2θ)|λ|b2<cos(x+b3θ)|λ|b3. (10)

We cannot ensure that for all x and m, cos(x)<cos(x+mθ)|λ|m<cos(x+b3θ)|λ|b3 implies that m=b2. However, we can ensure that it happens for “many” x, including for infinitely many numbers of the form x=nθ, where n3(mod4).

We can translate the last part into the notation of the intervals 𝒥d(γ,δ). Then, arbitrarily setting δ2=1.95 and δ3=2, we strengthen (10) in the form of item (b) in Lem. 19: we want that for all x,

0<cos(x)<cos(x+b2θ)|λ|b2<1.95cos(x)<cos(x+b3θ)|λ|b3<2cos(x).

Thus, I(π/2,π/2) (as cos(x)>0 and x), I𝒥b2(1,1.95) (as 1cos(x)<cos(x+b2θ)|λ|b2<1.95cos(x)), and similarly, I𝒥b3(1.95,2).

Dealing with items (c) and (d) of Lem. 19 is harder. In Lem. 18, we have proven many useful results on the structure of the sets 𝒥d(γ,δ). For d and small enough δ, if 𝒥d(γ,δ) is non-empty, it is a single interval and |𝒥d(γ,δ))|=O((δγ)|λ|d). Hence, item (d) of Lem. 19 can easily be estimated using a geometric series:

d=D|𝒥d(1,δ)|O(δ|λ|D).

Moreover, every point in 𝒥d(γ,δ)) is at most O(δ|λ|d) away from dθ±π/2 in /(2π). We illustrate this in Fig. 1.

Figure 1: Let λ=1+2i. Then for d=1,2,3,4, 𝒥d(1,3) are drawn in cyan, green, magenta, and red, respectively, and |𝒥d(1,3)| are π/2, 0.464, 0.182, 0.073, respectively. As some of the intervals overlap, we have stacked them vertically for visual purposes. The points in black mark out dθ±π/2 for the corresponding value of d.

As we have guessed δ2=1.95 and δ3=2, we will construct inductively as

:=3:=𝒥b3(1.95,2)2:=𝒥b2(1,1.95)1.

At each step i=1,2,3, the items (a) and (b) hold for ji while items (c) and (d) also hold for the interval i.

First, for i=1, we take 1:=(1.1,1.1)(π/2,π/2). Then, for an integer d<0 and x1, we have that cos(x+θd)|λ|d<cos(x). So in item (c), we can take D=1. As δ3=2, we can prove that item (d) is satisfied:

|d=1𝒥d(1,δ3)|d=1|𝒥d(1,δ3)|<|1|. (11)

For i=2, recall that δ2=1.95 and that b2 has to satisfy 𝒥b2(1,1.95)1 and b20(mod4). We pick b2=4, being the smallest possible choice for b2. As shown in Fig. 2, we have that 2 is indeed in 1. For all integers d20 not equal to either 0 or 4, 2𝒥d(1,2) is empty. As d=21|𝒥d(1,2)|<|𝒥4(1,1.95)|, 2:=𝒥4(1,1.95) is not covered by the intervals 𝒥d(1,2) with d{0,4}. Hence, D=21 is valid choice in this step of our inductive procedure.

Figure 2: We drew the intersection of 1 and 𝒥d(1,1.95) for d=1,,20 in red when d0(mod4) and in green when d0(mod4). For ease of visibility, 𝒥d(1,2) is positioned higher for d=1,2,3, and for intervals 𝒥d(1,2) that are too small to draw, their position is marked out with a dot.

For b3, recall that b30(mod4). We find that 𝒥b(1.95,2)2 is non-empty for b=0,4,38,99,160,309,370,. The smallest such b that is not equal to 0 or 4 (which are already in use) and congruent to 0 modulo 4 is 160. Then,

d=1,d{4,160}|𝒥d(1,2)𝒥160(1.95,2)|<|𝒥160(1.95,2)|

allows us to take D=160 and :=𝒥160(1.95,2). This interval is tiny: it has length approximately 5.71058. However, it satisfies the constraints in Lem. 19.

Now we must show that there is some x=nθ such that n3(mod4) and item 4 of Lem. 16 holds, as the three other conditions are already satisfied.

Recall that for a real number x, |x|2π denotes the distance from x to the nearest integer multiple of 2π. If x=nθ and nN, then items 1, 2 and 3 of Lem. 16 are satisfied. Assume that for n as above item 4 in Lem. 16 is not satisfied. Then nθ𝒥d(1,2) for some d{0,4,160}. Hence, combining the information from Lem. 18 to the effect that 𝒥d(1,2) is close to dθ±π/2 modulo 2π, we have:

c1|λ|d>|nθ(dθ±π/2)|2π>|d+n|c2 (12)

for two constants c1 (derived from Lem. 18) and c2>0 (derived from Thm. 10).

As there are many nθ in that are fairly “evenly” distributed, we are able to prove that (12) cannot hold for all n. In other words, there must be some n for which

0<vn<vn+4<vn+160

and vn<vm<vn+160 implies that m=n+4. Then, using the strategy of Lemma 16, we translate this back to the sequence unn=0 to give the required result when n is large enough. In particular, we can calculate that when taking

n=218085867698737188268427463501308698889728969450963229999559,

we have that nθ, n3(mod4), and un<um<un+160 implies that m=4. Thus, for some value m, pm,pm+1,pm+2=un,un+4,un+160 where all three terms are divisible by 5. Thus, the pattern 0,0,0 does indeed appear in pmmod5m=0.

However, this construction is far from optimal. Although our choices of b2=4 and b3=160 were as small as possible with our choice of δ2 and δ3, δ2 and δ3 were not optimal choices. Still, that would give just a minor improvement as one has to be more careful with inequalities like (11). Then one can find that choosing b2=8, b3=28, and n=16958443, we indeed get that there is no m{0,8,28} such that um is between un, un+b2 and un+b3. Meanwhile, cos(nθ)0.404, cos((n+b2)θ)94.5 and cos((n+b3)θ)751, and so one should have had chosen much larger δ2 and δ3 for which the general bounds in Lem. 18 fail.

3.3 Proof of the Main Theorems

In this section, we establish Thms. 4 and 1. We first prove two technical lemmas. These two lemmas will both serve to prove these theorems but also to prove Lem. 19, the continuous version of disjunctivity.

For the first lemma, we want that to show that each interval in /(2π) contains intervals 𝒥d(0,|λ|) and numbers nθ+ϕ for relatively small numbers d and n in a certain congruent class. The proof of this lemma in in the full version.

Lemma 20.

Let T2. There is a number C9>0 such that for every t{0,,T1} and small enough interval (π/2,π/2)/(2π), there are ||C9n1,n22||C9 such that n1θ+ϕ, 𝒥n2(0,|λ|) and n1,n2t(modT).

For the fourth condition of Lem. 19, we would like D to be large. That is, the smallest d such that 𝒥d(1,δ) has to be quite large. This is hard to guarantee, and our solution is to show that we can find a subinterval for which this number D is large. We show this in the following lemma whose proof is in the full version.

Lemma 21.

Assume /(2π), b2,,b, 1<δ2<<δ<|λ| and D>0 satisfy the assertions of Lem. 19. Then for every small enough ε>0 there is a subinterval of of length ε for which the assertions of Lem. 19 hold for these b2,,b and δ1,,δ and some D>ε1/2.

With these last two lemmas, we have gathered all the tool needed to prove Lem. 19, which is also done in the full version. The idea of this proof is as described in Sec. 3.2: we inductively construct the numbers b2,,b and intervals

(π/2,π/2)12:=b2(1,δ2)3:=b3(δ2,δ3)b(δ1,δ)

such that on step i, we show that some bi exists for which the lemma holds for b2,,bi and i.

Having established the continuous version of our result in Lem. 19, we can solve the discrete version Thm. 4 and conclude Thm. 1: the MSO theory of the structure ;<,{un:n} is decidable. The proof shares some of the main ideas of the proof of Lem. 19.

Proof of Thms. 1 and 4.

As Thm. 4 implies Thm. 1 thanks to Lem. 15, it is sufficient to prove Thm. 4. In turn, Lem. 16 states that Thm. 4 is implied by the following statement: for all N, T2 and t1,,t{0,,T1}, we can find n1,,n such that

  1. 1.

    n1,,nN;

  2. 2.

    njtj(modT) for 1j;

  3. 3.

    0<vn1<<vn;

  4. 4.

    for all m such that vn1<vm<vn, m{n2,,n1}.

We will prove this statement for given N,T,, and t1,,t.

Lem. 19 gives numbers b2,,b, 1<δ<|λ| and an interval . We take n1=n and nj=n+bj for 2j.

When nN, nt1(modT) and nθ+ϕ, imply items 1, 2, and 3. Indeed, items 2 and 3 hold for n1=n, and for 2j, nj=n+bj>nN and

njn1+bjn1+(tjt1)tj(modT).

For item 3, we have that as nθ+ϕ,

0<cos(nθ+ϕ)<cos((n+b2)θ+ϕ)|λ|b2<<cos((n+b)θ+ϕ)|λ|b.

Multiplying these inequalities by |λ|n and noting that cos(mθ+ϕ)|λ|m=vm entails item 3.

Let 0<ε<|| be small enough. Lemma 21 implies that there is an interval ε of length ε such that ε𝒥d(1,δ)= for all d<ε1/2. By Lem. 20, there is a εC9<n<2εC9 such that nθ+ϕε and nt1(modT). Thus, nN for small enough ε and items 1, 2, and 3 are satisfied.

Now assume that m{n2,,n1} and vn1<vm<vn. Then, setting d=mn, we have that d{b2,,b1} and

0<cos(nθ+ϕ)|λ|n<cos((n+d)θ+ϕ)|λ|n+d<cos((n+b)θ+ϕ)|λ|n+b.

As nθ+ϕε, cos((n+b)θ+ϕ)|λ|b<δcos(nθ+ϕ). Thus, nθ+ϕ is in 𝒥d(1,δ) and so dε1/2 (which comes from Lem. 21).

By Lem. 18, nθ+ϕ𝒥d(1,δ) and dθ±π/2 are at most C8|λ|d apart. Thus, for a constant c1 derived from Thm. 10,

C8|λ|d |(nθ+ϕ)(dθ±π/2)|π/2
|ei(nθ+ϕ)(dθ±π/2)|
>π2|n+d|c1.

Taking logarithms,

log(C8)+c1log(n+d)>dlog|λ|.

As dε1/2 and nεC9, we have that n,d2 for small enough ε. Then, dnn+d and so

log(C8)+c1log(n)+c1log(d)>dlog|λ|.

Hence, either

2log(C8)+2c1log(d)>dlog|λ|or2c1log(n)>dlog|λ|.

The former is impossible for large enough d (and thus small enough ε) while for the latter, the upper and lower bounds for n and d give that

2c1log(2εC9)>2c1log(n)>dlog|λ|>log|λ|ε1/2,

which again is impossible for small enough ε. Hence item 4 also follows.

4 Concluding Remarks

Our main result, Thm. 1, significantly expands the decidability landscape of MSO theories of structures of the form ;<,P, by considering predicates P obtained as sets of positive terms of non-degenerate simple LRS having two dominant roots. A natural question is whether any of these constraints can be relaxed. It is conceivable that investigating simple LRS with three dominant roots might yield positive decidability results, although the present development would have to be significantly altered at various junctures. We also note that open instances of the Skolem Problem [12] can easily be reduced from in the presence of simple LRS having four or more dominant roots. Likewise, considering non-simple LRS having three dominant roots or more exposes one to Positivity-hardness [20]. Let us remark that non-degeneracy is essential, as lifting this restriction has the effect of voiding the constraint on the number of dominant roots.555More precisely, one can construct a degenerate LRS for which the two dominant roots “cancel” each other out at every odd index, and the sub-LRS that remains over odd indices (involving the non-dominant roots) can be arbitrarily complex. Finally, one might envisage expanding the present setup by adjoining further predicates. However, even in the simplest of cases, one expects quite formidable difficulties in attempting to follow that research direction; see [3].

References

  • [1] Rajab Aghamov, Christel Baier, Toghrul Karimov, Joël Ouaknine, and Jakob Piribauer. Linear dynamical systems with continuous weight functions. In HSCC, pages 22:1–22:11. ACM, 2024. doi:10.1145/3641513.3650173.
  • [2] Christel Baier, Florian Funke, Simon Jantsch, Toghrul Karimov, Engel Lefaucheux, Florian Luca, Joël Ouaknine, David Purser, Markus A. Whiteland, and James Worrell. The orbit problem for parametric linear dynamical systems. In CONCUR, volume 203 of LIPIcs, pages 28:1–28:17. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2021. doi:10.4230/LIPICS.CONCUR.2021.28.
  • [3] Valérie Berthé, Toghrul Karimov, Joris Nieuwveld, Joël Ouaknine, Mihir Vahanwala, and James Worrell. On the decidability of monadic second-order logic with arithmetic predicates. In LICS, pages 11:1–11:14. ACM, 2024. doi:10.1145/3661814.3662119.
  • [4] Achim Blumensath. Monadic Second-Order Model Theory. http://www.fi.muni.cz/˜blumens/, 2024. [Online; accessed on 24 January 2025].
  • [5] Émile M. Borel. Les probabilités dénombrables et leurs applications arithmétiques. Rendiconti del Circolo Matematico di Palermo (1884-1940), 27(1):247–271, 1909.
  • [6] J. Richard Büchi. Symposium on decision problems: On a decision method in restricted second order arithmetic. In Studies in Logic and the Foundations of Mathematics, volume 44, pages 1–11. Elsevier, 1966.
  • [7] Yann Bugeaud. Distribution modulo one and Diophantine approximation, volume 193. Cambridge University Press, 2012.
  • [8] Olivier Carton and Wolfgang Thomas. The monadic theory of morphic infinite words and generalizations. Information and Computation, 176(1):51–65, 2002. doi:10.1006/INCO.2001.3139.
  • [9] Julian D’Costa, Toghrul Karimov, Rupak Majumdar, Joël Ouaknine, Mahmoud Salamati, and James Worrell. The pseudo-reachability problem for diagonalisable linear dynamical systems. In MFCS, volume 241 of LIPIcs, pages 40:1–40:13. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2022. doi:10.4230/LIPICS.MFCS.2022.40.
  • [10] Artūras Dubickas and Min Sha. Positive density of integer polynomials with some prescribed properties. Journal of Number Theory, 159:27–44, 2016.
  • [11] Calvin C. Elgot and Michael O. Rabin. Decidability and undecidability of extensions of second (first) order theory of (generalized) successor. The Journal of Symbolic Logic, 31(2):169–181, 1966. doi:10.2307/2269808.
  • [12] Graham Everest, Alf van der Poorten, Igor Shparlinski, and Thomas Ward. Recurrence Sequences. Mathematical Surveys and Monographs. American Mathematical Society, 2003.
  • [13] Glyn Harman. One hundred years of normal numbers. In Surveys in Number Theory, pages 57–74. AK Peters/CRC Press, 2002.
  • [14] Toghrul Karimov, Edon Kelmendi, Joris Nieuwveld, Joël Ouaknine, and James Worrell. The power of positivity. In LICS, pages 1–11. IEEE, 2023. doi:10.1109/LICS56636.2023.10175758.
  • [15] Toghrul Karimov, Edon Kelmendi, Joël Ouaknine, and James Worrell. What’s decidable about discrete linear dynamical systems? In Principles of Systems Design, volume 13660 of Lecture Notes in Computer Science, pages 21–38. Springer, 2022. doi:10.1007/978-3-031-22337-2_2.
  • [16] Toghrul Karimov, Engel Lefaucheux, Joël Ouaknine, David Purser, Anton Varonka, Markus A. Whiteland, and James Worrell. What’s decidable about linear loops? Proc. ACM Program. Lang., 6(POPL):1–25, 2022. doi:10.1145/3498727.
  • [17] Lauwerens Kuipers and Harald Niederreiter. Uniform distribution of sequences. Courier Corporation, 2012.
  • [18] Florian Luca, Joël Ouaknine, and James Worrell. Algebraic model checking for discrete linear dynamical systems. In FORMATS, volume 13465 of Lecture Notes in Computer Science, pages 3–15. Springer, 2022. doi:10.1007/978-3-031-15839-1_1.
  • [19] Eugene M. Matveev. An explicit lower bound for a homogeneous rational linear form in the logarithms of algebraic numbers. II. Izvestiya: Mathematics, 64(6):1217, 2000.
  • [20] Joël Ouaknine and James Worrell. Positivity problems for low-order linear recurrence sequences. In SODA, pages 366–379. SIAM, 2014. doi:10.1137/1.9781611973402.27.
  • [21] Martine Queffélec. Old and new results on normality. Lecture Notes-Monograph Series, pages 225–236, 2006.
  • [22] Alexander Rabinovich. On decidability of monadic logic of order over the naturals extended by monadic predicates. Information and Computation, 205(6):870–889, 2007. doi:10.1016/J.IC.2006.12.004.
  • [23] Alexander Rabinovich and Wolfgang Thomas. Decidable theories of the ordering of natural numbers with unary predicates. In International Workshop on Computer Science Logic, pages 562–574. Springer, 2006.
  • [24] Aleksei Lvovich Semenov. Logical theories of one-place functions on the set of natural numbers. Mathematics of the USSR-Izvestiya, 22(3):587, 1984.
  • [25] S. Shelah. The monadic theory of order. Ann. Math., 102:379–419, 1975.
  • [26] Wolfgang Thomas. Ehrenfeucht games, the composition method, and the monadic theory of ordinal words. In Structures in Logic and Computer Science, volume 1261 of Lecture Notes in Computer Science, pages 118–143. Springer, 1997. doi:10.1007/3-540-63246-8_8.
  • [27] Wolfgang Thomas. Languages, automata, and logic. In Handbook of Formal Languages: Volume 3 Beyond Words, pages 389–455. Springer, 1997. doi:10.1007/978-3-642-59126-6_7.
  • [28] R. Tijdeman, M. Mignotte, and T. N. Shorey. The distance between terms of an algebraic recurrence sequence. J. Reine Angew. Math., 346:63–76, 1984.