OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)
FORMULA
a(n) = A062326(n)^2-2. - Zak Seidov, Apr 29 2015
EXAMPLE
127^2 - 2 = 16127.
MATHEMATICA
f[n_]:=n^2-2; lst={}; Do[p=Prime[n]; If[PrimeQ[f[p]], AppendTo[lst, f[p]]], {n, 7!}]; lst (* Vladimir Joseph Stephan Orlovsky, Jul 16 2009 *)
Select[Prime[Range[150]]^2 - 2, PrimeQ] (* Vincenzo Librandi, Apr 29 2015 *)
PROG
(PARI) lista(nn) = forprime(p=1, nn, if (isprime(q=p^2-2), print1(q, ", "))); \\ Michel Marcus, Jan 08 2015
(Sage)
a = lambda p: p^2-2
[a(p) for p in primes(600) if is_prime(a(p))] # Bruno Berselli, Apr 29 2015
(Magma) [a: p in PrimesUpTo(1000) | IsPrime(a) where a is p^2-2 ]; // Vincenzo Librandi, Apr 29 2015
(Haskell)
a049002 n = a049002_list !! (n-1)
a049002_list = filter ((== 1) . a010051') a049001_list
-- Reinhard Zumkeller, Jul 30 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Herman H. Rosenfeld (herm3(AT)pacbell.net)
EXTENSIONS
More terms from James A. Sellers
STATUS
approved