Sylvester's sequence
![](http://upload.wikimedia.org/wikipedia/commons/thumb/a/a2/Sylvester-square.svg/300px-Sylvester-square.svg.png)
In number theory, Sylvester's sequence is an integer sequence in which each term is the product of the previous terms, plus one. Its first few terms are
- 2, 3, 7, 43, 1807, 3263443, 10650056950807, 113423713055421844361000443 (sequence A000058 in the OEIS).
Sylvester's sequence is named after James Joseph Sylvester, who first investigated it in 1880.[1] Its values grow doubly exponentially, and the sum of its reciprocals forms a series of unit fractions that converges to 1 more rapidly than any other series of unit fractions.[2] The recurrence by which it is defined allows the numbers in the sequence to be factored more easily than other numbers of the same magnitude,[3] but, due to the rapid growth of the sequence, complete prime factorizations are known only for a few of its terms.[4] Values derived from this sequence have also been used to construct finite Egyptian fraction representations of 1, Sasakian Einstein manifolds,[5] and hard instances for online algorithms.[6]
Formal definitions
[edit]Formally, Sylvester's sequence can be defined by the formula[7]
The product of the empty set is 1,[8] so this formula gives s0 = 2, without need of a separate base case.
Alternatively, one may define the sequence by the recurrence[3]
- with the base case s0 = 2.
It is straightforward to show by induction that this is equivalent to the other definition.[9]
Closed form formula and asymptotics
[edit]The Sylvester numbers grow doubly exponentially as a function of n. Specifically, it can be shown that
for a number E that is approximately 1.26408473530530...[10] (sequence A076393 in the OEIS). This formula has the effect of the following algorithm:
- s0 is the nearest integer to E 2; s1 is the nearest integer to E 4; s2 is the nearest integer to E 8; for sn, take E 2, square it n more times, and take the nearest integer.
This would only be a practical algorithm if we had a better way of calculating E to the requisite number of places than calculating sn and taking its repeated square root.[11]
The double-exponential growth of the Sylvester sequence is unsurprising if one compares it to the sequence of Fermat numbers Fn ; the Fermat numbers are usually defined by a doubly exponential formula, , but they can also be defined by a product formula very similar to that defining Sylvester's sequence:[12]
Connection with Egyptian fractions
[edit]The unit fractions formed by the reciprocals of the values in Sylvester's sequence generate an infinite series:[13]
The partial sums of this series have a simple form,
which is already in lowest terms.[14] This may be proved by induction, or more directly by noting that the recursion implies that
so the sum telescopes[14]
Since this sequence of partial sums (sj − 2)/(sj − 1) converges to one, the overall series forms an infinite Egyptian fraction representation of the number one:
One can find finite Egyptian fraction representations of one, of any length, by truncating this series and subtracting one from the last denominator:
The sum of the first k terms of the infinite series provides the closest possible underestimate of 1 by any k-term Egyptian fraction.[2] For example, the first four terms add to 1805/1806, and therefore any Egyptian fraction for a number in the open interval (1805/1806, 1) requires at least five terms.
It is possible to interpret the Sylvester sequence as the result of a greedy algorithm for Egyptian fractions, that at each step chooses the smallest possible denominator that makes the partial sum of the series be less than one.[15]
Uniqueness of quickly growing series with rational sums
[edit]As Sylvester himself observed, Sylvester's sequence seems to be unique in having such quickly growing values, while simultaneously having a series of reciprocals that converges to a rational number. This sequence provides an example showing that double-exponential growth is not enough to cause an integer sequence to be an irrationality sequence.[16]
To make this more precise, it follows from results of Badea (1993) that, if a sequence of integers grows quickly enough that
and if the series
converges to a rational number A, then, for all n after some point, this sequence must be defined by the same recurrence
that can be used to define Sylvester's sequence.[17]
Erdős & Graham (1980) conjectured that, in results of this type, the inequality bounding the growth of the sequence could be replaced by a weaker condition,[18]
Badea (1995) surveys progress related to this conjecture; see also Brown (1979).[19]
Divisibility and factorizations
[edit]If i < j, it follows from the definition that sj ≡ 1 (mod si ). Therefore, every two numbers in Sylvester's sequence are relatively prime. The sequence can be used to prove that there are infinitely many prime numbers, as any prime can divide at most one number in the sequence. More strongly, no prime factor of a number in the sequence can be congruent to 5 modulo 6, and the sequence can be used to prove that there are infinitely many primes congruent to 7 modulo 12.[20]
Much remains unknown about the factorization of the numbers in Sylvester's sequence. For instance, it is not known if all numbers in the sequence are squarefree, although all the known terms are.[21]
As Vardi (1991) describes, it is easy to determine which Sylvester number (if any) a given prime p divides: simply compute the recurrence defining the numbers modulo p until finding either a number that is congruent to zero (mod p) or finding a repeated modulus.[3] Using this technique he found that 1166 out of the first three million primes are divisors of Sylvester numbers,[22] and that none of these primes has a square that divides a Sylvester number. The set of primes that can occur as factors of Sylvester numbers is of density zero in the set of all primes:[23] indeed, the number of such primes less than x is .[24]
The following table shows known factorizations of these numbers (except the first four, which are all prime):[4]
n | Factors of sn |
---|---|
4 | 13 × 139 |
5 | 3263443, which is prime |
6 | 547 × 607 × 1033 × 31051 |
7 | 29881 × 67003 × 9119521 × 6212157481 |
8 | 5295435634831 × 31401519357481261 × 77366930214021991992277 |
9 | 181 × 1987 × 112374829138729 × 114152531605972711 × 35874380272246624152764569191134894955972560447869169859142453622851 |
10 | 2287 × 2271427 × 21430986826194127130578627950810640891005487 × P156 |
11 | 73 × C416 |
12 | 2589377038614498251653 × 2872413602289671035947763837 × C785 |
13 | 52387 × 5020387 × 5783021473 × 401472621488821859737 × 287001545675964617409598279 × C1600 |
14 | 13999 × 74203 × 9638659 × 57218683 × 10861631274478494529 × C3293 |
15 | 17881 × 97822786011310111 × 54062008753544850522999875710411 × C6618 |
16 | 128551 × C13335 |
17 | 635263 × 1286773 × 21269959 × C26661 |
18 | 50201023123 × 139263586549 × 60466397701555612333765567 × C53313 |
19 | 775608719589345260583891023073879169 × C106685 |
20 | 352867 × 6210298470888313 × C213419 |
21 | 387347773 × 1620516511 × C426863 |
22 | 91798039513 × C853750 |
As is customary, Pn and Cn denote prime numbers and unfactored composite numbers n digits long.
Applications
[edit]Boyer, Galicki & Kollár (2005) use the properties of Sylvester's sequence to define large numbers of Sasakian Einstein manifolds having the differential topology of odd-dimensional spheres or exotic spheres. They show that the number of distinct Sasakian Einstein metrics on a topological sphere of dimension 2n − 1 is at least proportional to sn and hence has double exponential growth with n.[5]
As Galambos & Woeginger (1995) describe, Brown (1979) and Liang (1980) used values derived from Sylvester's sequence to construct lower bound examples for online bin packing algorithms.[6] Seiden & Woeginger (2005) similarly use the sequence to lower bound the performance of a two-dimensional cutting stock algorithm.[25]
Znám's problem concerns sets of numbers such that each number in the set divides but is not equal to the product of all the other numbers, plus one. Without the inequality requirement, the values in Sylvester's sequence would solve the problem; with that requirement, it has other solutions derived from recurrences similar to the one defining Sylvester's sequence. Solutions to Znám's problem have applications to the classification of surface singularities (Brenton and Hill 1988) and to the theory of nondeterministic finite automata.[26]
Curtiss (1922) describes an application of the closest approximations to one by k-term sums of unit fractions, in lower-bounding the number of divisors of any perfect number, and Miller (1919) uses the same property to upper bound the size of certain groups.[27]
See also
[edit]Notes
[edit]- ^ Sylvester (1880).
- ^ a b This claim is commonly attributed to Curtiss (1922), but Miller (1919) appears to be making the same statement in an earlier paper. See also Rosenman & Underwood (1933), Salzer (1947), Soundararajan (2005), and Nathanson (2023).
- ^ a b c Vardi (1991).
- ^ a b All prime factors p of Sylvester numbers sn with p < 5×107 and n ≤ 200 are listed by Vardi. Ken Takusagawa lists the factorizations up to s9 and the factorization of s10. The remaining factorizations are from a list of factorizations of Sylvester's sequence maintained by Jens Kruse Andersen. Retrieved 2014-06-13.
- ^ a b Boyer, Galicki & Kollár (2005).
- ^ a b Galambos & Woeginger (1995); Brown (1979); Liang (1980).
- ^ Sloane, N. J. A. (ed.). "Sequence A000058 (Sylvester's sequence)". The On-Line Encyclopedia of Integer Sequences. OEIS Foundation.
- ^ Nešetřil & Matoušek (1998).
- ^ A proof by induction is given by Sylvester (1880), p. 333.
- ^ Graham, Knuth & Patashnik (1989), formula 4.17, p. 109, and exercise 4.37, p. 147; see also Golomb (1963).
- ^ Graham, Knuth & Patashnik (1989), p. 109.
- ^ Sloane, N. J. A. (ed.). "Sequence A000215 (Fermat numbers)". The On-Line Encyclopedia of Integer Sequences. OEIS Foundation.
- ^ This series is the starting point for Sylvester (1880)
- ^ a b Sylvester (1880), p. 334.
- ^ Nathanson (2023).
- ^ Guy (2004).
- ^ Badea (1993).
- ^ Erdős & Graham (1980).
- ^ Badea (1995); Brown (1979).
- ^ Guy & Nowakowski (1975).
- ^ Graham, Knuth & Patashnik (1989), research problem 4.65, p. 151; Vardi (1991); see also Chentouf (2020)
- ^ This appears to be a typo, as Andersen finds 1167 prime divisors in this range.
- ^ Jones (2006).
- ^ Odoni (1985).
- ^ In their work, Seiden and Woeginger refer to Sylvester's sequence as "Salzer's sequence" after the work of Salzer (1947) on closest approximation.
- ^ Domaratzki et al. (2005).
- ^ Curtiss (1922); Miller (1919).
References
[edit]- Badea, Catalin (1993). "A theorem on irrationality of infinite series and applications". Acta Arithmetica. 63 (4): 313–323. doi:10.4064/aa-63-4-313-323. MR 1218459.
- Badea, Catalin (1995). "On some criteria for irrationality for series of positive rationals: a survey" (PDF). Archived from the original (PDF) on 2008-09-11.
- Boyer, Charles P.; Galicki, Krzysztof; Kollár, János (2005). "Einstein metrics on spheres". Annals of Mathematics. 162 (1): 557–580. arXiv:math.DG/0309408. doi:10.4007/annals.2005.162.557. MR 2178969. S2CID 13945306.
- Brenton, Lawrence; Hill, Richard (1988). "On the Diophantine equation 1=Σ1/ni + 1/Πni and a class of homologically trivial complex surface singularities". Pacific Journal of Mathematics. 133 (1): 41–67. doi:10.2140/pjm.1988.133.41. MR 0936356.
- Brown, D. J. (1979). A lower bound for on-line one-dimensional bin packing algorithms. Tech. Rep. R-864. Coordinated Science Lab., Univ. of Illinois, Urbana-Champaign.
- Chentouf, A. Anas (2020). "On Sylvester's Sequence and some of its properties" (PDF). Parabola. 56 (2).
- Curtiss, D. R. (1922). "On Kellogg's diophantine problem". American Mathematical Monthly. 29 (10): 380–387. doi:10.2307/2299023. JSTOR 2299023.
- Domaratzki, Michael; Ellul, Keith; Shallit, Jeffrey; Wang, Ming-Wei (2005). "Non-uniqueness and radius of cyclic unary NFAs". International Journal of Foundations of Computer Science. 16 (5): 883–896. doi:10.1142/S0129054105003352. MR 2174328.
- Erdős, Paul; Graham, Ronald L. (1980). Old and new problems and results in combinatorial number theory. Monographies de L'Enseignement Mathématique, No. 28, Univ. de Genève. MR 0592420.
- Galambos, Gábor; Woeginger, Gerhard J. (1995). "On-line bin packing — A restricted survey". Mathematical Methods of Operations Research. 42 (1): 25. doi:10.1007/BF01415672. MR 1346486. S2CID 26692460.
- Golomb, Solomon W. (1963). "On certain nonlinear recurring sequences". American Mathematical Monthly. 70 (4): 403–405. doi:10.2307/2311857. JSTOR 2311857. MR 0148605.
- Graham, Ronald; Knuth, Donald E.; Patashnik, Oren (1989). Concrete Mathematics (2e ed.). Addison-Wesley. ISBN 978-0-201-55802-9.
- Guy, Richard K. (2004). "E24 Irrationality sequences". Unsolved Problems in Number Theory (3rd ed.). Springer-Verlag. p. 346. ISBN 0-387-20860-7. Zbl 1058.11001.
- Guy, Richard; Nowakowski, Richard (1975). "Discovering primes with Euclid". Delta (Waukesha). 5 (2): 49–63. MR 0384675.
- Jones, Rafe (2006). "The density of prime divisors in the arithmetic dynamics of quadratic polynomials". Journal of the London Mathematical Society. 78 (2): 523–544. arXiv:math.NT/0612415. Bibcode:2006math.....12415J. doi:10.1112/jlms/jdn034. S2CID 15310955.
- Liang, Frank M. (1980). "A lower bound for on-line bin packing". Information Processing Letters. 10 (2): 76–79. doi:10.1016/S0020-0190(80)90077-0. MR 0564503.
- Nešetřil, Jaroslav; Matoušek, Jiří (1998). Invitation to Discrete Mathematics. Oxford University Press. p. 12. ISBN 0-19-850207-9.
- Miller, G. A. (1919). "Groups possessing a small number of sets of conjugate operators". Transactions of the American Mathematical Society. 20 (3): 260–270. doi:10.2307/1988867. JSTOR 1988867.
- Nathanson, Melvyn B. (January 2023). "Underapproximation by Egyptian fractions". Journal of Number Theory. 242: 208–234. arXiv:2202.00191. doi:10.1016/j.jnt.2022.07.005.
- Odoni, R. W. K. (1985). "On the prime divisors of the sequence wn+1 =1+w1⋯wn". Journal of the London Mathematical Society. Series II. 32: 1–11. doi:10.1112/jlms/s2-32.1.1. Zbl 0574.10020.
- Rosenman, Martin; Underwood, F. (1933). "Problem 3536". American Mathematical Monthly. 40 (3): 180–181. doi:10.2307/2301036. JSTOR 2301036.
- Salzer, H. E. (1947). "The approximation of numbers as sums of reciprocals". American Mathematical Monthly. 54 (3): 135–142. doi:10.2307/2305906. JSTOR 2305906. MR 0020339.
- Seiden, Steven S.; Woeginger, Gerhard J. (2005). "The two-dimensional cutting stock problem revisited". Mathematical Programming. 102 (3): 519–530. doi:10.1007/s10107-004-0548-1. MR 2136225. S2CID 35815524.
- Soundararajan, K. (2005). "Approximating 1 from below using n Egyptian fractions". arXiv:math.CA/0502247.
- Sylvester, J. J. (1880). "On a point in the theory of vulgar fractions". American Journal of Mathematics. 3 (4): 332–335. doi:10.2307/2369261. JSTOR 2369261.
- Vardi, Ilan (1991). Computational Recreations in Mathematica. Addison-Wesley. pp. 82–89. ISBN 0-201-52989-0.
External links
[edit]- Irrationality of Quadratic Sums, from K. S. Brown's MathPages.
- Weisstein, Eric W. "Sylvester's Sequence". MathWorld.