2 To 100
2 To 100
Abstract. Sierpiński proved that there are infinitely many odd integers k
such that k · 2n +1 is composite for all n ≥ 0. These k are now called Sierpiński
numbers. We define a Sierpiński number base b to be an integer k > 1 for which
gcd(k+1, b−1) = 1, k is not a rational power of b, and k · bn +1 is composite
for all n > 0. We discuss ways that these can arise, offer conjectured least
Sierpiński number in each of the bases 2 < b ≤ 100 (34 are proven), and show
that all bases b admit Sierpiński numbers. We also show that under certain cir-
r
cumstances there are base b Sierpiński numbers k for which k, k2 , k3 , ..., k2 −1
are each base b Sierpiński numbers.
Key words and phrases. Sierpiński number, covering set, generalized Fermat number.
†Undergraduate student. The beginning of this work was partially supported by a University
of Tennessee at Martin College of Engineering and Natural Sciences undergraduate research grant.
1
2 BRUNNER, CALDWELL, KRYWARUCZENKO AND LOWNSDALE
gap by providing a definition and then extending the studied bases systematically
to include all of the bases up through 100.
We will prove that Sierpiński numbers exist for all bases b > 1, and offer conjec-
tured least Sierpiński numbers for the bases 2 < b ≤ 100. For 34 of these bases we
are able to prove that the conjectured values are indeed the least.
b k b k
6 6, 36, 216, 1296, 7776, 46656 46 46, 2116
8 2, 4, 8, 16, 32 48 48
10 10, 100, 1000 52 52, 2704
12 12, 144 58 58, 3364
16 16, 256, 4096, 65536 60 60, 3600
18 18, 324 64 4, 16
22 22, 484 66 66, 4356, 287496, 18974736
24 24, 576, 13824 70 70, 4900
26 26 72 72
28 28, 784 78 78, 6084
30 30 80 80
32 2, 4, 8 82 82, 6724
36 36, 1296 88 88
40 40, 1600, 64000 96 96, 9216
42 42, 1764 100 100
2 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048,
4096, 8192, 16384, 32768, 65536
4 4, 16, 64, 256, 1024, 4096, 16384, 65536
‡ Just those smaller than conjectured least base b Sierpiński and
with gcd(k + 1, b − 1) = 1.
numbers that are a power of 2. If the only Fermat numbers are the five known,
then 216 2n + 1 would be composite for n > 0, and therefore 216 = 65536 would be
the least Sierpiński number, not Selfridge’s 78557.
Since the existence of infinitely many Fermat primes is undecidable at this point
in time, it seems best to define generalized Sierpiński numbers in such a way as
to exclude the Fermat numbers and, for bases other than powers of 2, to exclude
n
the generalized Fermat numbers Fn (b) = b2 + 1 [12]. At the end of this section
(Theorem 2.3) we will show that this is equivalent to adding the requirement that
p
k is not a rational power of b (k 6= b q for integers p ≥ 0 and q > 0), and hence that
k > 1. Those values so omitted are listed in Table 2.
Combining these considerations we generalize Sierpiński numbers as follows.
Definition 2.2. Let b > 1 be an integer. A Sierpiński number base b (or b-
Sierpiński) is an integer k > 1 for which gcd(k+1, b−1) = 1, k is not a rational
power of b, and k · bn +1 is composite for all n > 0.
Notice that this definition extends the definition of Sierpiński numbers in base 2
as well as the larger integer bases—yet still 78557 likely remains the least possible
2-Sierpiński number.
We end this section by showing we have properly characterized those pairs k and
b which may generate infinitely many generalized Fermat numbers. Recall that the
order of b modulo a relatively prime integer p, denoted ordp (b), is the least positive
integer m for which p divides bm − 1. So in particular ordp (b) divides φ(p) (Euler’s
φ function of p).
4 BRUNNER, CALDWELL, KRYWARUCZENKO AND LOWNSDALE
Theorem 2.3. Let b > 1 and k > 0 be integers for which gcd(k + 1, b − 1) = 1.
There is an integer c > 1 for which k · bn + 1 = Fr (c) for infinitely many integer
values of r and n, if and only if k is a rational power of b.
Proof. Let b > 1 and k > 0 be fixed integers for which gcd(k + 1, b − 1) = 1.
Suppose there is an integer c for which k ·bn +1 (n > 1) is the generalized Fermat
number Fr (c) for infinitely many pairs of integers r and n. Choose two such pairs
(r, n) and (s, m) with n < m. Then
r s
k · bn + 1 = c2 + 1 and k · bm + 1 = c2 + 1.
s r 2s −2r m2r −n2s
Thus bm−n = c2 −2 , and it follows b = c m−n , k = c m−n , and therefore k is a
rational power of b (and both are rational powers of c).
Conversely, suppose k is a rational power of b, say k = be/f for relatively prime
integers e and f with e ≥ 0 and f > 0. Then because b is an integer, b = cf
and k = ce for some integer c. Write f = 2t f 0 where f 0 is an odd integer. Now
gcd(cf − 1, ce + 1) = 1, so by Lemma 2.1 c is even and the power of 2 which divides
e is at least as great as the power of 2 which divides f. So we may write e = 2t e0
for some integer (not necessarily odd) e0 . Note that if r is any positive multiple
of ordf 0 (2), then e0 ≡ e0 2r (mod f 0 ), so we may solve the following for a positive
integer n = n(r):
e0 + f 0 n = e0 2r .
So it follows
e + f n = 2t (e0 + f 0 n) = e0 2r+t ,
and there are infinitely many choices of r and n for which
0 r+t 0
k · bn + 1 = ce+f n + 1 = ce 2 + 1 = Fr+t (ce ).
Erdös apparently believed that all Sierpiński numbers arise from covers [16,
Section F13]. That is probably not the case [14]. In section 5 we will show that
not all b-Sierpiński numbers arise from covers. In practice, though, most small
examples do come from covers.
There are two basic ways of constructing covers: Sierpiński’s approach of using
the Fermat numbers (generalized Fermat numbers in our case) and Selfridge’s use
of factors of bn − 1. We begin with the latter.
Theorem 3.2. Every element of an N -cover S of k · bn +1 divides bN −1.
Proof. Choose p ∈ S. This p must divide k · bn +1 for some n ≤ N . It then also
divides k · bn+N +1, so divides their difference k · bn (bN −1). Since p does not divide
k · bn , this completes the proof.
For example, Selfridge’s Sierpiński 78557 arises from the cover {3, 5, 7, 13, 19, 37, 73}.
Each prime of this cover divides 236 −1, so to prove 78557 is a Sierpiński number base
2, it is sufficient to show that each of the first 36 terms in the sequence 78557 · 2n +1
(n > 0) are divisible by one of these seven primes.
The previous theorem also tells us that for every element p of an N -cover of
k · bn +1, ordp (b) divides N . It is easy to show N = lcmp∈S (ordp (b)).
Unless we say otherwise, in the rest of this article “N -cover” will mean non-
trivial N cover, that is N > 1. Note that if S is an N -cover for one k, b pair, then
it is also a cover for infinitely many other multipliers k and bases b.
Theorem 3.3. An N -cover S of k · bn + 1 is also a cover of K · B n + 1 for all
integers K ≡ k, B ≡ b (mod P ) where P is the product of the primes in the cover
S.
It follows by Dirichlet’s theorem that there are infinitely many prime multipliers,
and infinitely many prime bases covered by any given N -cover.
In what follows it will be helpful to recall the cyclotomic polynomials Φn (x).
These are defined by
Y Y
Φn (x) = (xd − 1)µ(n/d) and so xn − 1 = Φd (x).
d|n d|n
n
This makes Φn (x) the “primitive part” of x − 1 when factoring, and the φ(n) zeros
of Φn (x) are the primitive nth roots of unity. For integers n and b greater than one,
if a prime p divides Φn (b) but not n, then p ≡ 1 (mod n).
Theorem 3.2 can now be greatly sharpened for more specific values of N .
Theorem 3.4. Let p be a prime number. The sequence base b has a p-cover S if
and only if Φp (b) has at least p distinct prime divisors greater than p.
Proof. Suppose first S is a p-cover of k · bn +1. So there is an element of S which
divides each element of T = {k · b1 +1, k · b2 +1, . . . , k · bp +1}. If q is an element of
S, the ordq (b) must divide the period of S, which is p. This order can not be one,
or {q} would be a trivial cover, so ordq (b) = p and therefore p | q − 1. This means
q can only divide one element of T , hence there are at least p primes in S. Finally,
these primes do not divide b − 1, so they each divide (bp − 1)/(b − 1) = Φp (b).
On the other hand, if Φp (b) has the p distinct prime divisors: q1 , q2 , q3 , . . . , qp ,
each greater than p, then none divide b − 1 so we can use the Chinese Remainder
6 BRUNNER, CALDWELL, KRYWARUCZENKO AND LOWNSDALE
Theorem to show they form a p-cover by solving the system of linear equations
k · b1 + 1 ≡ 0 (mod q1 )
2
k·b +1 ≡ 0 (mod q2 )
..
.
k · bp + 1 ≡ 0 (mod qp )
For example, the base b has a 2-cover if and only if b + 1 has at least two distinct
odd prime divisors. Examples of this include b = 14, 20, 29, 32, 34, 38, 41, 44, 50,
54, 56, 59, 62, 64, 65, 68, 69, 74, 76, 77, 83, 84, 86, 89, 90, 92, 94, 98 . . . For all of
these listed bases except 68 and 86, we have proven the 2-cover generates the least
generalized Sierpiński number base b (see Table 5). Bowen [5] also used 2-covers to
address the bases b = 2s + 1 where s 6= 2m + 1 and s > 5.
Similarly, the base b has a 3-cover if and only if Φ3 (b) = b2 + b + 1 has at least
three distinct divisors greater than 3. The first such bases are b = 74, 81, 87, 100,
102, 107, 121, and 135. Of those bases b ≤ 100, only for 100 does the 3-cover yield
the least generalized Sierpiński number. Bases 74, 81, and 87 have 3-covers, but
these produce larger multipliers k than can be generated by other methods.
The minimal base for longer prime period covers grows quickly: (p, minimal base
b) = (2, 14), (3, 74), (5, 339), (7, 2601), (11, 32400), and (13, 212574).
The structure of composite period covers are more interesting. For example,
4-covers usually arise from an odd prime factor p for which the base b has order 2
(a divisor of Φ2 (b) = b + 1), and two primes q1 , q2 for which b has order 4 (divisors
of Φ4 (b) = b2 + 1). Then the terms of the sequence k · bn + 1 (n = 1, 2, . . .) are
divisible by the primes of the 4-cover in a pattern like
p, q , p, q2 , p, q1 , p, q2 , . . . .
| 1{z } | {z }
k · b1 + 1 ≡ 0 (mod p)
k · b2 + 1 ≡ 0 (mod q1 )
k · b4 + 1 ≡ 0 (mod q2 )
For 29 of the bases in Table 5, 4-covers provide the least known b-Sierpiński num-
bers.
Most 6-covers involve four primes. Often one prime in the cover, say p, has
period 2 (ordp (b) = 2), and there are three more of orders 3 or 6, say q1 , q2 , q3 ,
dividing the terms of k · bn + 1 in a pattern similar to
p, q , p, q , p, q3 , p, q1 , p, q2 , p, q3 , . . . .
| 1 {z 2 } | {z }
Most bases have a 12-cover. One way one of these can arise is if each of Φ2 (b),
Φ3 (b), Φ4 (b), Φ6 (b) and Φ12 (b) have a primitive divisor, call them p2 , p3 , p4 , p6 and
GENERALIZED SIERPIŃSKI NUMBERS BASE b 7
Theorem 3.8. If there is a non-trivial cover for k · bn +1, then k+1 has an odd
prime divisor.
The first of these was used by Stanton [32] in his analysis of possible covers for the
b = 2 case.
where the ith component is gcd(k · bi +1, bN −1) (1 ≤ i ≤ N ). From this it was a
simple hand calculation to find the actual covering set of primes.
This program was run on the 16 nodes of our Beowulf cluster for about 80-CPU
days to find the constants k and the associated covers in the first columns of Table 5
except for b = 3, 7 and 15. Some individual values (e.g., 71), required substantially
longer search times.
The program Sierpiński has several limitations. First, one must know some-
thing of N in advance because the program is set up to seek all covers with period
N dividing a specified constant. For Table 5 we usually sought periods dividing
7! = 5040. It is possible that we missed some covers for smaller k values.
Second, the program Sierpiński only seeks values of k belonging to covers.
Such k values are Sierpiński numbers base b, but there may be smaller b-Sierpiński
numbers that do not arise from covers. We will discuss this in the next section.
Finally, the program Sierpiński is too slow to find the least covers for bases
like 3. For those bases we may begin by factoring bN −1 for various small values of
N and construct covers as described in the previous section. For example, Brennen
[7] used this method to find 3574321403229074 (48-cover) for b = 3. (This improved
earlier results of Bowen [5] and Saouter [28].) With an improved algorithm Bosma
[6] reduced this to k = 125050976086 (144-cover).
are composite. These have been excluded by our definition, but we see no reason
that there could not be other examples of b-Sierpiński numbers that have neither
covers nor factorizations. (This same possibility is addressed for base 2 in [19, 14]).
2http://robert.gerbicz.googlepages.com/coveringsets
10 BRUNNER, CALDWELL, KRYWARUCZENKO AND LOWNSDALE
This cover has another interesting property: not only are the k so constructed
Sierpiński numbers, but so are k t for any odd integer t > 1. It turns out that by using
multiple different composite terms we may use essentially the same construction to
find k for which k t is also prime for t divisible by low powers of 2. This was done
by Filaseta et al. [14] for the regular Sierpiński numbers, and virtually the same
argument works here.
Theorem 7.1. Let b > 1 be an integer for which b + 1 is not a power of 2. If there
m
are at least r generalized Fermat numbers Fm (b) = b2 + 1 which are each divisible
by at least two distinct odd primes, then there are infinitely many integers k such
that k t is a b-Sierpiński number for all positive integers t not divisible by 2r .
0 0 0
Proof. Define the integer Fm (b) by Fm (b) = 2rm Fm (b) with rm ≥ 0 and Fm (b) odd.
0
In what follows we need the fact that Fm (b) has at least one prime factor. This
0 0
is the case unless Fm (b) = 1. If Fm (b) = 1, then b is odd. It follows that m = 0;
2m 0
otherwise b + 1 ≡ 2 (mod 8). So Fm (b) = 1 implies b + 1 = 2r0 . For the rest of
0
the proof we assume b + 1 is not a power of 2 (hence Fm (b) > 1).
Let m0 < m1 < · · · < mr−1 be non-negative integers for which the generalized
0
Fermat numbers Fmj (b), hence Fm j
(b), each have at least two distinct odd prime
0
factors, say pj and qj . Note these primes are all different as b, b − 1, and Fm (b)
(m ≥ 0) are pairwise relatively prime.
By the Chinese Remainder Theorem there are infinitely many solutions to the
following set of congruences.
0 (mod b − 1)
1 (mod b)
0
k≡ 1 (mod Fm (b)) for 0 ≤ m < mr−1 and m 6∈ {m0 , . . . , mr−1 }
1 (mod p j ) for 0≤j ≤r−1
2mj −j
b (mod qj ) for 0 ≤ j ≤ r − 1.
We further restrict k to those solutions which are greater than each of the moduli
above. The first of these modular restrictions guarantees gcd(k + 1, b − 1) = 1, and
the second guarantees that gcd(k, b) = 1, so k is not a rational power of b.
Given any positive integer t not divisible by 2r , say t = 2w t0 where t0 is odd and
0 ≤ w < r, we must show k t bn + 1 is composite for each positive integer n. Fix a
positive integer n and let n = 2i n0 where n0 is odd. We may complete this proof
by showing k t bn + 1 is divisible by
0
Fm (b) if i < mw and m 6∈ {m0 , . . . , mw }
d= pj if i = mj for some j with 0 ≤ j ≤ w
qw if i > mw .
Since d < k < k t bn + 1, this will show the latter term is composite.
i i 0
If i ≤ mw , then d divides b2 + 1 which divides b2 n + 1 = bn + 1. Because k ≡ 1
(mod d), it follows d divides k t bn + 1.
If instead i > mw , then
mw −w w 0 mw 0 0
k t ≡ (b2 )2 t
≡ (b2 )t ≡ (−1)t ≡ −1 (mod qw ).
2mw
Now d = qw divides b + 1 which divides
i i−1
b2 − 1 = (b − 1)(b + 1)(b2 + 1)(b4 + 1) · · · (b2 + 1).
GENERALIZED SIERPIŃSKI NUMBERS BASE b 11
k r
23140626796 1
3352282631064632411056 2
38454071854799507248067375352496 3
295612797233398523232282186442005794587542575896 4
1202250010386171287615458085\ 5
386724017477152933279927552922222324231610279296
4833\ 6
96281140918612511630787705875212985273405983905\
512852696056665671273849671134513427529509057456
18081740848967\ 7
53044039134711401288516658002520824319923798573\
210660688220428187289811356995735827761349820556
i
n0
Thus d divides b2 − 1 and it follows
k t · bn + 1 ≡ −(bn − 1) ≡ 0 (mod d).
This completes the proof of the theorem.
0
For example, when b = 5, Fm (b) is prime for m = 0, 1, and 2. It is composite
(with distinct prime divisors) for m0 = 3, m1 = 4, . . . , m10 = 13. If we let qj be
0
the smallest prime factor of Fm j
(b) and pj be the second smallest, then we get
the b-Sierpiński numbers in Table 3. Of course, as noted in the discussion before
the proof, rather than use prime factors, we may use any two relatively prime
(non-trivial) proper divisors. So it is sufficient to know any odd prime divisor and,
after checking that the given generalized Fermat is not a power of that prime,
use the cofactor as the second “prime.” Table 4 shows that there are 244 known
composite generalized Fermat numbers Fn (5), so there are 5-Sierpiński k for which
244
k, k 2 , k 3 , ..., k 2 −1 are all Sierpiński numbers (from [23]).
8. Conclusions
Of the many possible generalizations of the Sierpiński numbers, we have discussed
what seemed the most natural to us. It would be interesting, but difficult, to study
the generalized Fermat cases that we excluded in our definition. It seems likely that
bases b can be found so that the least Sierpiński number is arbitrarily large. One
can also ask the reverse question: given a value k, can we find a base b for which
12 BRUNNER, CALDWELL, KRYWARUCZENKO AND LOWNSDALE
k is a base b Sierpiński? A partial answer has been provided by one of the authors
[24].
Note that every cover of a sequence of the form k · bn + 1 (n > 0) is also a cover
of a sequence k 0 · bn − 1 (n > 0), and vice versa. Positive odd integers k for which
k · bn − 1 are composite for all n > 0 are called Riesel numbers after an article by
Riesel [27] in 1956 (so the Riesel numbers predate the Sierpiński numbers). Thus
another generalization to study would be the generalized Riesel numbers (k which
make k · bn − 1 composite for all n > 0 with suitable restrictions on k and b); as
well as the numbers that are both b-Riesels and b-Sierpińskis. Part of this work is
being done informally by Barnes and others [4], as are restrictive cases like seeking
the smallest b-Sierpiński numbers which are prime.
A. de Polignac conjectured (and quickly retracted) the guess that every positive
odd number can be written in the form 2n + p for a prime p and integer n > 0.
He did this even though Euler had previously shown that this was not the case for
127 or 929 [1]. Again every cover of k · bn + 1 (n > 0) is also a cover of a sequence
bn + k (n > 0), and vice versa.
References
1. L. Babai, C. Pomerance & P. Vértesi, The mathematics of Paul Erdös, Notices of the AMS,
45:1 (January 1998), 19–31.
2. R. Baillie, G. V. Cormack & H. C. Williams, The problem of Sierpiński concerning k·2n +1
Math. Comput., 37 (1981) 229–231; corrigendum, 39 (1982) 308.
3. A. S. Bang, Taltheoretiske Undersøgelser, Tidsskrift for Mat., 5(4) (1886), 70–80, 130–137.
4. G. Barnes, Sierpiński conjecture reservations, May 2008, http://gbarnes017.googlepages.
com/Sierp-conjecture-reserves.htm.
5. R. Bowen, The sequence kan + 1 composite for all n, Math. Monthly, 71:2 (1964), 175–176.
6. W. Bosma, Some computational experiments in number theory. In Discovering Mathematics
with Magma, volume 19 of Algorithms Comput. Math., pp. 1–30. Springer-Verlag, Berlin,
2006.
7. J. Brennen, PrimeForm e-mail discussion list, May 16, 2002, http://tech.groups.yahoo.com/
group/primenumbers/message/7147.
8. J. Brillhart, D.H. Lehmer & J.L. Selfridge, New primality criteria and factorizations of 2m ±1,
Math. Comp., 29 (1975) 620–647.
9. D. A. Buell & J. Young, Some large primes and the Sierpiński problem, SRC Technical Report
88-004, Supercomputing Research Center, Lanham, MD, May 1988.
10. Y. G. Chen, On integers of the forms kr −2n and kr 2n +1, J. Number Theory, 98:2 (2003),
310–319.
11. Y. G. Chen, On integers of the forms k ± 2n and k2n ± 1, J. Number Theory, 125:1 (2007),
14–25, MR 2333115.
12. H. Dubner & W. Keller, Factors of generalized Fermat numbers, Math. Comput., 64 (1995)
397-405, MR 1270618.
13. P. Erdös, On integers of the form 2k + p and some related problems, Summa Brasil. Math.,
2 (1950) 113-123, MR 0044558.
14. M. Filaseta, C. Finch & M. Kozek, On powers associated with Sierpiński numbers, Riesel
numbers and Polignac’s conjecture, J. Number Theory, 128:7 (2008)1916–1940, MR 2423742.
15. Y. Gallot, Proth.exe, primes.utm.edu/programs/gallot/, July 2005.
16. R. K. Guy, Unsolved Problems in Number Theory (3rd ed.), Problem Books in Mathematics,
Springer-Verlag, New York, 2004.
17. L. Helm, P. Moore, P. Samidoost & G. Woltman, Resolution of the mixed Sierpiński problem,
INTEGERS: Elec. J. Comb. Num. Th., to appear.
18. L. Helm & D. Norris, Seventeen or Bust–a distributed attack on the Sierpiński problem,
http://www.seventeenorbust.com/.
19. A. S. Izotov, A note on Sierpiński numbers, Fibonacci Quart., 33 (1995) 206–207; MR
96f:11020.
GENERALIZED SIERPIŃSKI NUMBERS BASE b 13
20. G. Jaeschke, On the smallest k such that all k · 2N + 1 are composite, Math. Comput., 40
(1983) 381–384; MR 84k:10006; corrigendum, 45 (1985) 637; MR 87b:11009.
21. L. Jones, Variations on a theme of Sierpiński, J. Integer Seq., 10 (2007).
22. W. Keller, Factors of Fermat numbers and large primes of the form k·2n +1, Math. Comput.,
41 (1983) 661–673; MR 85b:11119; II (incomplete draft, 92-02-19).
23. W. Keller, Factors of generalized Fermat numbers found after Björn & Riesel, http://www1.
uni-hamburg.de/RRZ/W.Keller/GFNfacs.html, Oct. 2008.
24. D. Krywaruczenko, A reverse Sierpiński number problem, Rose-Hulman Undergrad. Math. J.
(electronic) to appear.
25. C. Nash & J. Fougeron, OpenPFGW (Open source software), http://tech.groups.yahoo.
com/group/primeform/.
26. R. M. Robinson, A report on primes of the form k·2n +1 and on factors of Fermat numbers,
Proc. Amer. Math. Soc., 9 (1958) 673–681; MR 20 #3097.
27. H. Riesel, Några stora primtal (Swedish: Some large primes), Elementa, 39 (1956) 258–260.
28. Y. Saouter, A Fermat-like sequence and primes of the form 2h · 3n + 1, Research Report 2728,
Nov. 1995, citeseer.ist.psu.edu/saouter95fermatlike.html.
29. W. Sierpiński, Sur un problème concernant les nombres k·2n +1, Elem. Math., 15 (1960)
73–74; MR 22 #7983; corrigendum, 17 (1962) 85.
30. N. Sloane, The on-line encyclopedia of integer sequences, Conjectured smallest Sierpiński
numbers, www.research.att.com/~njas/sequences/A123159.
31. R. Smith, Sierpiński and Riesel bases 6 to 18, Conjectured smallest Sierpiński numbers, www.
mersenneforum.org/showthread.php?t=6895, August 2007.
32. R. G. Stanton, Further results on covering integers of the form 1 + k · 2N by primes, Com-
binatorial mathematics, VIII (Geelong, 1980), Springer Lecture Notes in Math., 884 (1981)
107–114; MR 84j:10009.
14 BRUNNER, CALDWELL, KRYWARUCZENKO AND LOWNSDALE