OFFSET
1,1
COMMENTS
Indices of 3-almost prime octagonal numbers.
Subset of A001358. - R. J. Mathar, Jun 07 2007
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
MAPLE
isA001358 := proc(n) if numtheory[bigomega](n) = 2 then true ; else false ; end ; end: isA129926 := proc(n) if isA001358(n) then isprime(3*n-2) ; else false ; fi ; end: for n from 1 to 1000 do if isA129926(n) then printf("%d, ", n) ; fi ; od ; # R. J. Mathar, Jun 07 2007
PROG
(PARI) list(lim)=my(v=List(), n); forprime(p=2, lim\2, forprime(q=2, min(lim\p, p), n=p*q; if(isprime(3*n-2), listput(v, n)))); Set(v) \\ Charles R Greathouse IV, Jan 31 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Teofilatto, Jun 06 2007
EXTENSIONS
More terms from R. J. Mathar, Jun 07 2007
STATUS
approved