OFFSET
1,1
MAPLE
isA007510 := proc(n) if isprime(n) then not isprime(n-2) and not isprime(n+2) ; else false; end if ; end proc:
isA014574 := proc(n) isprime(n+1) and isprime(n-1) ; end proc:
isA167706 := proc(n) isA007510(n) or isA014574(n) ; end proc:
isA005117 := proc(n) n =1 or numtheory[issqrfree](n) ; end proc:
isA168270 := proc(n) isA005117(n) and isA167706(n) ; end proc:
for n from 1 to 600 do if isA168270(n) then printf("%d, ", n) ; fi; od: # R. J. Mathar, Dec 09 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Nov 22 2009
EXTENSIONS
Corrected (missing terms inserted) by R. J. Mathar, Dec 09 2009
STATUS
approved