OFFSET
1,1
LINKS
John Cerkan, Table of n, a(n) for n = 1..10000
EXAMPLE
37 is here because neither 37+2 nor 37-2 is divisible by squares.
MAPLE
with(numtheory): A049282:=n->`if`(isprime(n) and issqrfree(n-2) and issqrfree(n+2), n, NULL): seq(A049282(n), n=1..10^3); # Wesley Ivan Hurt, Jun 25 2016
MATHEMATICA
lst={}; Do[p=Prime[n]; If[SquareFreeQ[p-2]&&SquareFreeQ[p+2], AppendTo[lst, p]], {n, 6!}]; lst (* Vladimir Joseph Stephan Orlovsky, Dec 20 2008 *)
PROG
(PARI) lista(nn) = forprime(p=2, nn, if (issquarefree(p-2) && issquarefree(p+2), print1(p, ", "))); \\ Michel Marcus, Jun 22 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved