OFFSET
1,1
COMMENTS
The associated prime(k) are in A136019.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1)=7 because 7 is smallest prime p such that (p+2)/3 is prime
a(2)=11 because 11 is smallest prime p such that (p+4)/5 is prime
a(3)=29 because 29 is smallest prime p such that (p+6)/7 is prime
MATHEMATICA
a = {}; Do[k = 1; While[ !PrimeQ[(Prime[k] + 2n)/(2n + 1)], k++ ]; AppendTo[a, Prime[k]], {n, 1, 100}]; a
Reap[Do[k = PrimePi[2n+1]; While[ !PrimeQ[((p=Prime[k]) + 2n)/(2n + 1)], k++ ]; Sow[p], {n, 100}]][[2, 1]] (* Zak Seidov, Mar 04 2013 *)
PROG
(PARI) a(n)=my(t); n+=n; forprime(p=2, default(primelimit), if(isprime(t=(n+1)*p-n), return(t))) \\ Charles R Greathouse IV, Feb 13 2011
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Artur Jasinski, Dec 10 2007
EXTENSIONS
Edited by R. J. Mathar, May 17 2009
STATUS
approved