OFFSET
1,1
MAPLE
A006530 := proc(n) numtheory[factorset](n) ; max(op(%)) ; end proc: isA166305 := proc(n) if type(n, 'even') then if numtheory[bigomega](n) = 2 then isprime(A006530(n)+8) ; else false; end if; else false; end if; end proc: for n from 4 to 3000 by 2 do if isA166305(n) then printf("%d, ", n) ; end if; end do: # R. J. Mathar, Jan 30 2010
MATHEMATICA
Select[2*Range[1500], PrimeOmega[#]==PrimeOmega[#+16]==2&] (* Harvey P. Dale, Dec 28 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Teofilatto, Oct 11 2009
EXTENSIONS
Extended by R. J. Mathar, Jan 30 2010
STATUS
approved