OFFSET
1,1
COMMENTS
All numbers of the form 2p, where p==2(mod 3) is prime, are in the sequence.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[3 Range@ 120 + 1, Mod[#, 3] == 2 &[FactorInteger[#][[-1, 1]]] &] (* Michael De Vlieger, Feb 07 2016 *)
PROG
(PARI) isok(n) = ((n % 3) == 1) && (n != 1) && ((vecmax(factor(n)[, 1]) % 3) == 2); \\ Michel Marcus, Feb 07 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Apr 13 2010
EXTENSIONS
Corrected and extended by Michel Marcus, Feb 07 2016
STATUS
approved