OFFSET
1,1
COMMENTS
A259758(n) = (2*a(n) - 3) * (3*a(n) - 2).
Except for a(1)=3 this is the same sequence as primes p such that A288814(3*p) - A288814(2*p) = 5. - David James Sycamore, Jul 22 2017
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[Prime@ Range@ 270, Times @@ Boole@ Map[PrimeQ, {2 # - 3, 3 # - 2}] > 0 &] (* Michael De Vlieger, Jul 22 2017 *)
Select[Prime[Range[300]], AllTrue[{2#-3, 3#-2}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Mar 08 2020 *)
PROG
(Haskell)
import Data.List.Ordered (isect)
a259730 n = a259730_list !! (n-1)
a259730_list = a063908_list `isect` a088878_list
(PARI) lista(nn) = forprime(p=3, nn, if(isprime(2*p-3) && isprime(3*p-2), print1(p, ", "))); \\ Altug Alkan, Jul 22 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jul 05 2015
STATUS
approved