OFFSET
1,1
COMMENTS
47-3=44=2^2*11, 47+3=50=2*5^2, ...
LINKS
John Cerkan, Table of n, a(n) for n = 1..10000
MATHEMATICA
<<NumberTheory`NumberTheoryFunctions` (* needed with older versions *)
lst={}; Do[p=Prime[n]; If[ !SquareFreeQ[p-3]&&!SquareFreeQ[p+3], AppendTo[lst, p]], {n, 6!}]; lst
Select[Prime[Range[500]], !SquareFreeQ[#-3]&&!SquareFreeQ[#+3]&] (* Harvey P. Dale, May 01 2012 *)
PROG
(PARI) lista(nn) = forprime(p=3, nn, if (! issquarefree(p-3) && ! issquarefree(p+3), print1(p, ", "))); \\ Michel Marcus, Jul 05 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Jan 21 2009
STATUS
approved