OFFSET
1,2
MAPLE
A046315 := proc(n) option remember; if n = 1 then 9; else for a from procname(n-1)+2 by 2 do if numtheory[bigomega](a) = 2 then return a; end if; end do: end if; end proc:
A020639 := proc(n) numtheory[factorset](n) ; min(op(%)) ; end proc:
isA176550 := proc(n) os := A046315(n) ; p := A020639(os) ; q := os/p ; ( A046315(n+1) mod nextprime(p) ) = 0 or (A046315(n+1) mod nextprime(q) = 0 ) ; end proc:
for n from 1 to 300 do if isA176550(n) then printf("%d, ", n) ; end if; end do: # R. J. Mathar, May 30 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Apr 20 2010
EXTENSIONS
Corrected (39 removed, 43 inserted) and extended by R. J. Mathar, May 30 2010
STATUS
approved