OFFSET
1,2
COMMENTS
The asymptotic density of this sequence is 6 * Sum_{p prime, p == 5 (mod 6)} ((p-1)/(Product_{q prime, q <= p} q)) = 0.8261626908... . - Amiram Eldar, Dec 04 2022
FORMULA
{k | A053669(6*k) == 5 (mod 6)}.
EXAMPLE
35 is present as 6*35 = 210 = 2*3*5*7, and the first nondividing prime is 11, which is of the form 6m+5. This is the first multiple of 5 in this sequence.
385 is not present as 6*385 = 2310 = 2*3*5*7*11, and the first nondividing prime is 13, which is of the form 6m+1, not of 6m+5.
MATHEMATICA
f[n_] := Module[{p = 2}, While[Divisible[n, p], p = NextPrime[p]]; p]; Select[Range[100], Mod[f[6*#], 6] == 5 &] (* Amiram Eldar, Dec 04 2022 *)
PROG
(PARI) isA358848(n) = A358847(n);
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Dec 03 2022
STATUS
approved