OFFSET
5,1
COMMENTS
This sequence is well-defined per the work of Booker and Pomerance.
The number 4 in the congruence in the name could be replaced by any value; this number was chosen for being the first integer that is not squarefree.
LINKS
Rémy Sigrist, Table of n, a(n) for n = 5..10000
Andrew R. Booker, Carl Pomerance, Squarefree smooth numbers and Euclidean prime generators, arXiv:1607.01557 [math.NT], 2016-2017.
Andrew R. Booker and Carl Pomerance, Squarefree smooth numbers and Euclidean prime generators, Proceedings of the American Mathematical Society 145 (2017), 5035-5042.
Rémy Sigrist, Colored scatterplot of (n, a(n)) for n = 5..1000000 (where the color is function of (a(n)-4)/A000040(n)).
FORMULA
a(n) = A261144(n, k) for some k in 1..2^n.
EXAMPLE
For n = 7:
- the 7th prime is 17,
- the first squarefree 17-smooth integers s, alongside (s-4) mod 17, are:
s 1 2 3 5 6 7 10 11 13 14 15 17 21
------------ -- -- -- - - - -- -- -- -- -- -- --
(s-4) mod 17 14 15 16 1 2 3 6 7 9 10 11 13 0
- hence a(7) = 21.
MATHEMATICA
a[n_] := Module[{p = Prime[n], k = 4}, While[! SquareFreeQ[k] || FactorInteger[k][[-1, 1]] > p, k += p; Continue[]]; k]; Array[a, 100, 5] (* Amiram Eldar, Dec 08 2018 *)
PROG
(PARI) a(n) = my (p=prime(n)); forstep (v=4, oo, p, if (issquarefree(v), my (f=factor(v)); if (f[#f~, 1] <= p, return (v))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Rémy Sigrist, Dec 05 2018
STATUS
approved