OFFSET
1,1
COMMENTS
The asymptotic density of this sequence is Sum_{p prime, p == 5 (mod 6)} ((p-1)/(Product_{q prime, q <= p} q)) = 0.1376937818... . - Amiram Eldar, Dec 04 2022
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) local p;
p:= 3;
do
p:= nextprime(p);
if n mod p <> 0 then return (p mod 6 = 5) fi
od
end proc:
select(filter, [seq(i, i=6..1000, 6)]); # Robert Israel, Dec 04 2023
MATHEMATICA
f[n_] := Module[{p = 2}, While[Divisible[n, p], p = NextPrime[p]]; p]; Select[Range[500], Mod[f[#], 6] == 5 &] (* Amiram Eldar, Dec 04 2022 *)
PROG
(PARI) isA358757(n) = A358755(n);
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Dec 03 2022
STATUS
approved