OFFSET
1,2
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
MAPLE
isA001358 := proc(n) numtheory[bigomega](n) = 2 ; end proc:
A001358 := proc(n) option remember ; if n = 1 then return 4 ; else for a from procname(n-1)+1 do if isA001358(a) then return a; end if; end do; end if; end proc:
for n from 1 to 900 do if isprime(n+A001358(n)) then printf("%d, ", n) ; end if; end do: # R. J. Mathar, Apr 20 2010
MATHEMATICA
Module[{nn=2000, semi}, semi=Select[Range[nn], PrimeOmega[#]==2&]; Transpose[ Select[ Thread[{semi, Range[Length[semi]]}], PrimeQ[Total[#]]&]][[2]]] (* Harvey P. Dale, Jun 16 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Apr 14 2010
EXTENSIONS
Corrected (68 inserted, 156 inserted) by R. J. Mathar, Apr 20 2010
STATUS
approved