OFFSET
1,1
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
FORMULA
A154734(n+1) = k(a(n)) where k(m) = m*(m+1)/2+2*(-1)^m. - R. J. Mathar, Jul 16 2022
EXAMPLE
The initial values of k(m) are:
k(0) = 2
k(1) = 1-2 = -1
k(2) = 4+1 = 5
k(3) = 9-5 = 4 and both 3 and 5 are primes, so 3 is the first term in the sequence
k(4) = 16-4 = 12, and 11 & 13 are primes, so a(2) = 4
and so on - N. J. A. Sloane, Jul 14 2022
MATHEMATICA
k=2; lst={}; Do[k=n^2-k; If[PrimeQ[k-1]&&PrimeQ[k+1], AppendTo[lst, n]], {n, 8!}]; lst
(* Second program: *)
k = 2; Reap[Do[Set[k, m^2 - k]; If[AllTrue[k + {-1, 1}, PrimeQ], Sow[m]], {m, 4000}]][[-1, -1]] (* Michael De Vlieger, Jul 14 2022 *)
PROG
(PARI) a154736(upto, k0=2) = {my(k=k0); for(n=1, upto, my(kk=n^2-k); if(isprime(kk-1) && isprime(kk+1), print1(n, ", ")); k=kk)};
a154736(5000) \\ Hugo Pfoertner, Jul 14 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Jan 14 2009
EXTENSIONS
Better name from Pontus von Brömssen, Jul 14 2022
STATUS
approved