OFFSET
1,1
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 300 do if isprime( A001358(n)+2 ) then printf("%d, ", n) ; end if; end do: # R. J. Mathar, Apr 20 2010
MATHEMATICA
Position[Select[Range[1000], PrimeOmega[#]==2&], _?(PrimeQ[#+2]&)]//Flatten (* Harvey P. Dale, Dec 01 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Apr 14 2010
EXTENSIONS
Entries checked by R. J. Mathar, Apr 20 2010
STATUS
approved